Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
PUBLICATIONS Geophysical Research Letters RESEARCH LETTER 10.1002/2016GL068590 Key Points: • Subducting Indian ocean plate has been imaged beneath the Andaman-Nicobar trench • Evidence of tear in subducting Indian ocean plate using Ps and Sp conversions • The lithospheric thickness varies from ~40 to 70 km all along the Island from north to south Supporting Information: • Supporting Information S1 Correspondence to: P. Kumar, [email protected] Citation: Kumar, P., G. Srijayanthi, and M. Ravi Kumar (2016), Seismic evidence for tearing in the subducting Indian slab beneath the Andaman arc, Geophys. Res. Lett., 43, 4899–4904, doi:10.1002/ 2016GL068590. Received 8 MAR 2016 Accepted 28 APR 2016 Accepted article online 2 MAY 2016 Published online 18 MAY 2016 Seismic evidence for tearing in the subducting Indian slab beneath the Andaman arc Prakash Kumar1, G. Srijayanthi1, and M. Ravi Kumar1 1 CSIR-National Geophysical Research Institute, Hyderabad, India Abstract Segmentation of a subduction zone through tearing is envisaged as an inevitable consequence of the differential rate of slab rollback along the strike of convergent plate boundaries. It is a key feature that controls plate tectonics and seismogenesis in a subduction setting. Globally, lithospheric tears are mostly recognized by seismic tomography and seismicity trends. However, such an intriguing feature has never been imaged with high resolution. Here we present seismological evidence for tearing of the Indian oceanic plate at shallow depths along the Andaman arc. Our image of the subducted plate using the shear-wave receiver function technique reveals three distinct plate segments. The middle lithospheric chunk has an abrupt offset of ~20 km relative to the northern and southern segments along the entire stretch of Andaman-Nicobar Islands. We interpret that this abrupt offset in the base of the lithosphere is caused by the tearing of the subducted oceanic plate. For the plate age of ~80 to 60 Myr, the lithospheric thickness varies from ~40 to 70 km. 1. Introduction Oceanic plates are created at mid-oceanic ridges and consumed along the trenches by subduction. Along the trenches, the oceanic plates are in continuous motion. In such convergent tectonic settings, the subducting slabs are migrated backward [Jarrard, 1986; Schellart et al., 2007] with a variable rate along the strike. As a consequence of such rollback, the subducted slabs experience deformation in terms of bending and/or tearing. This creates a discontinuity in the geometry of the subducting slab along the trench [Spence, 1977] and usually occurs in the case of a retreating lithospheric segment [Rosenbaum and Piana Agostinetti, 2015]. The important factors, which control slab rollback, are the buoyancy of the lithosphere and the contrasting physical properties with respect to the ambient asthenosphere [Garfunkel et al., 1986]. Other factors that control the segmentation of the slab are the topographic anomalies within the subducting plate, such as ridges, fracture zones, and seamount chains [Kodaira et al., 2000], the major structures abutting the over-riding plate [Ryan and Scholl, 1993] and large-scale variations in the buoyancy of the subducting plate related to its thermal age [Yáñez and Cembrano, 2004]. Wilson [1965] opined that tearing is inevitable in a tectonic setting where subduction terminates, thus creating an offset between the adjacent segments of the oceanic lithospheric plate. Slab-tear is an important entity that controls plate tectonics, in particular the tectonic activities in convergent plate boundaries. It possibly provides pathways for ascending magmas [Gvirtzman and Nur, 1999; Rosenbaum et al., 2008; Gasparon et al., 2009]. Tearing of oceanic plates has been documented in many subduction zones of the world [Isacks et al., 1969; Forsyth, 1975; Chatelain et al., 1992], based on the observations or interpretations of the abrupt break in seismicity [Reyners et al., 1991] or the sudden change in the velocity anomalies evidenced in tomographic models [Lucente et al., 1999; Benoit et al., 2011; Giacomuzzi et al., 2012]. Examples include the South Sandwich trench [Barker, 2001], Tonga trench [Millen and Hamburger, 1998], Lesser Antilles [Clark et al., 2008], Gibraltar arc [Meighan et al., 2013], and north Apennines [Rosenbaum and Piana Agostinetti, 2015]. ©2016. American Geophysical Union. All Rights Reserved. KUMAR ET AL. Other than the observation of seismicity distribution and seismic images from low-resolution tomographic studies, such an intriguing structural feature of the Earth has never been imaged with high resolution, using body waves. We designed a seismic experiment to image the nature and geometry of the subducting Indian plate beneath the Andaman-Nicobar Island arc (as shown in Figures 1 and S1 in the supporting information) by installing broadband seismic instruments. In the aftermath of the mega 2004-Sumatra earthquake (Mw 9.3) followed by the deadly Tsunami, this region has attracted the attention of seismologists world over. Further, in the years 2005, 2009, and 2010, other mega earthquakes having Mw greater than 7.5 rocked this region. In order to understand the genesis of the mega thrust earthquake of 26 December 2004, Shapiro et al. [2008] produced a shear velocity model and imaged the highly oblique subducting slab. However, this study TEARING IN THE INDIAN OCEANIC PLATE 4899 Geophysical Research Letters 10.1002/2016GL068590 did not focus on tearing of the slab. North of the Andaman Islands, slab tear between the Burma and Andaman arc regions has been documented based on the vertical distribution of seismicity [Richards et al., 2007]. Pesicek et al. [2010] reported the lack of a clear subduction in the same region and further suggested that the vertical offset of lithosphere might be responsible for the northward termination of the rupture zone of the 2004 mega earthquake. The Andaman region is characterized by the existence of an active spreading center in the back arc Figure 1. Topographic map with the station disposition and seismic profiles. beneath the Andaman Sea [e.g., The inverted triangles (white) are the locations of broadband seismic stations. Khan and Chakraborty, 2005] and AB, CD, and EF are the three profiles along which the migrated receiver the interaction of the linear 90E function images are constructed in Figures 2 and 3. The seismicity (gray dots) and major earthquakes (thick black circles) are also shown. The pink dots are the ridge with the trench (Figure 1). The 2005 swarm earthquakes observed after the 26 December 2004 earthquake complexity of the region imposes a (Mw, 9.0–9.3). The plate boundary and the contours (blue lines) of oceanic major challenge in mapping this zone. crustal age (gray contours) (http://www.ngdc.noaa.gov/mgg/ocean_age/ocean_age_2008.html) of the Indian ocean region are also shown. The red triangles Seismologically, this region is one of the least studied in the world. In order are the locations of the volcano [taken from Curray, 2005]. to image the plate geometry beneath the Andaman-Nicobar Islands, we applied the shear-to-compressional (Sp) converted wave technique to a large number of teleseismic earthquake waveforms. We also present the compressional-to-shear (Ps) images as an independent support to our findings (supporting information). 2. Data and Method The seismological data from the Andaman-Nicobar Islands are from an 11-station seismic network being operated by the Council for Scientific and Industrial Research-National Geophysical Research Institute and the India Meteorological Department. The operation periods of our experiment are different for different stations (Table-S1 in the supporting information), and the stations are located only on the islands covering the entire Andaman-Nicobar region. The mean elevation of the region is close to or below the mean sea level (Figure 1). In order to generate images of the plate architecture using S-receiver functions (SRF) [Farra and Vinnik, 2000], we use teleseismic waveforms registered at these broadband seismic stations. The SRFs are generated from the vertical component of ground motion recorded at the seismic station. The advantage of SRFs is that they are devoid of S wave multiple reflections within the time window of arrival of converted phases from the shallower discontinuities. This provides a clean image of the subsurface. In recent times, the SRF technique has become an effective tool to map the uppermost mantle discontinuities such as the lithosphere-asthenosphere boundary (LAB) and has been applied in diverse tectonic settings. The distance ranges of the teleseismic earthquakes whose waveforms are selected for the analysis are 60°–85° for S waves and 85°–120° for SKS waves:). To ensure good signal-to-noise ratios, the earthquakes having magnitude greater than 5.7 Mb from all the available back-azimuths are chosen (Figures S2 and S3 in supporting information). Further, the S phases on the radial component of the waveforms are visually inspected and only those which show good SV arrivals are retained. The waveforms are then filtered with a low-pass filter with a corner frequency of 0.25 Hz for S wave data and rotated into the radial and transverse components using the back-azimuth information. To isolate the weak Sp phase from the parent S phase, we further rotate the Z-R-T components into L-Q-T components using an incidence angle derived from the minimum SV energy on the P component at zero time, in order to have an optimal isolation of P-SV-SH [Kumar et al., 2006; Kumar and Kawakatsu, 2011]. In the presence of large noise and strong lateral heterogeneity, the incidence KUMAR ET AL. TEARING IN THE INDIAN OCEANIC PLATE 4900 Geophysical Research Letters Figure 2. S-receiver function data along two horizontal profiles (a) AB and (b) CD (Figure 1). The bathymetry with station locations (black inverted triangles) is also shown at the top of each profile (see wiggle plots in Figure S3). The earthquakes (Mb > 5 ) are well within the lithosphere. The observation shows two separate slab segments from the Indian ocean side and from the Andaman Sea side dipping in opposite directions. 10.1002/2016GL068590 angles determined in this way may be erroneous, and therefore, we discard the events with incidence angles >52°. This coordinate rotation ensures a clean decomposition of Sto-p conversions. Once the traces are rotated into P and SV components, we employ source normalization using deconvolution adopting a time domain deconvolution technique [Berkhout, 1977] to generate the S-receiver functions. In order to achieve a time to depth migration, all the receiver function amplitudes are placed at their proper conversion points by tracing the ray from the source to the receiver using the global IASP91 Earth model [Kennett and Engdahl, 1991]. 3. Seismic Images Figures 2 and 3 show the migrated images below the Andaman-Nicobar Islands down to the uppermost mantle. While generating these images, we have taken care of the artifacts arising due to the lateral heterogeneity in structure, since the S-to-p conversion points sample a wide spatial extent (see Figure S4 for piercing points for S-to-p conversions). Figure 2, containing near east-west sections along the profiles AB and CD (see Figure 1), reveals that there are at least two prominent phases corresponding to conversions from the Moho (red-positive polarity, crust-mantle boundary) and the LAB (blue-negative polarity, marked in the figure) representing a velocity increase across the Moho and decrease across the LAB. The wiggle plots of these figures are presented in Figure S5 to demonstrate the robustness of the observation. A close examination of the section down to 20–30 km depth reveals that there are two oppositely dipping features, i.e., one from the Indian ocean side and the other from the Andaman back arc side. The thickness of the 80–90 Ma old Indian ocean plate subducting below the Andaman-Nicobar Island region is found to be ~50 km (Figure 2), which is somewhat lower than that predicted by the thermal model. The primary reason for this discrepancy in the thickness might be due the deviation from a normal oceanic mantle due to the presence of sea-mounts and ridges reported here [e.g., Kamesh Raju et al., 2004]. The presence of Nintyeast Ridge might affect the plate structurally. The nature and size of the velocity jump across the oceanic LAB have been discussed by Kawakatsu et al. [2009] in the context of Pacific and Philippine sea plates. In the context of converting conversion times to the depths to the interfaces, issues such as (i) thickness of the warmer younger plate tending to be overestimated by using a 1-D reference velocity model, (ii) presence of azimuthal anisotropy, and (iii) regional variational in Vp/Vs ratios [Kumar and Kawakatsu, 2011] cannot be ruled out. Further, the Indian oceanic plate in this region may not be a representative of the normal oceanic plate. The migrated images (Figure 2) show a complex lithospheric architecture beneath the Andaman islands. In both the images (Figures 2a and 2b), two separate lithospheric fragments are seen and the seismicity (Mb > 5) along the respective profiles is well within the lithosphere. In order to image the nature of the Indian plate subduction beneath the Andaman trench, we constructed north-south images along a curved profile (Figure 1). The imaging (Figures 3) brings out the complex tectonic scenario, especially at the crustal level. At a deeper level, the base of the lithospheric plate is evident at a depth of about 50 km, except between 7 and 8° latitude, where it deepens down to ~70 km (see Figure S5 for wiggle plots and Figure S6 for comparison with PRF images). In Figure S5 binning and stacking of the traces are performed based on the piercing points of Sp conversions at depths of 20 and 60 km with a bin step of 0.05° and 0.1° with overlap of 0.1°, respectively, and bins which contain more than 5 traces are only stacked. The north-south section reveals that there are three distinct segments of the subducting Indian plate along the trench that is (i) to the north of the intersection of the 90°E ridge with the Andaman trench, (ii) KUMAR ET AL. TEARING IN THE INDIAN OCEANIC PLATE 4901 Geophysical Research Letters 10.1002/2016GL068590 Figure 3. Receiver function image along the north-south profile EF parallel to the trench, with the station locations (red inverted triangles). The average thickness of the lithosphere (blue color) is ~50 km from north to south along the arc, except in the vicinity of ~8 deg latitude, where it deepens to ~80 km. Interestingly, the positive phase (red) which is at ~20 km all along the profile suddenly deepens to ~40 km, where the base of the lithosphere is deepest. The seismicity (Mb > 5, green dots) with mega earthquakes of 2005 and 2010 is well within the lithosphere. This abrupt offset of the middle part of the lithosphere has been interpreted as caused by the tearing in oceanic slab (interpretive model in bottomright inset). The numerals 1, 2, and 3 denote the three distinct lithospheres observed in our seismic sections corresponding to those in the inset cartoon. between the Andaman and Nicobar islands, and (iii) on the Nicobar side. An abrupt offset of ~20 km is observed in the middle part of the lithosphere with respect to the adjacent segments. In order to demonstrate the robustness of the observation, we further show the stacked receiver functions from three representative stations, which are in the north, middle, and south segments of the arc. The station stacks displayed in Figure 4 (also see Figure S7 for Ps stacks) show that the phase, which we identified as a Sp conversion from the LAB, has a substantial amplitude above the ±2 SE. 4. Discussion The seismic images reveal clear evidence for an offset in the subducted lithosphere relative to the adjacent downgoing slabs along the arc. We interpret this offset as being caused by tearing in the oceanic lithosphere. Seismicity and tomography studies have already indicated tearing in other places along the Burma-AndamanSumatra trench [e.g., Dasgupta et al., 2003], which plays a significant role in shaping the tectonics in this region. The seismicity trends in our study region also vindicate the possibility of a slab tear (Figure S8). The hypocentral distribution in the tear region shows somewhat deeper seismicity, i.e., tending toward the trench. Such tears could have served as Figure 4. S-receiver function stacks from three stations CMBY, KAMO, and vents for the outpouring of the asthenoHUTB plotted at their geographical locations. The stacks are generated spheric material. This is substantiated by after applying moveout correction with a reference slowness of 6.4 s/° the presence of numerous volcanoes in using the IASP91 global Earth model. The short vertical bars are the ±2 SE these regions [Curray, 2005]. Also, lava in time. The thin lines on either side of the wiggles indicate the error in amplitudes. A similar stack section has also been shown for P-to-s conversions samples from the Barren islands reveal their affinity toward the mantle or deep in Figure S7. KUMAR ET AL. TEARING IN THE INDIAN OCEANIC PLATE 4902 Geophysical Research Letters 10.1002/2016GL068590 mantle origin [e.g., Alam et al., 2004]. The cartoon shown as an inset in Figure 3 summarizes our observations. The thickness of the oceanic plate as mapped by the Sp conversions (also by Ps conversions) varies between ~40 and 70 km along the arc. The two slabs, one on the northern side and the other on the southern side of 8° latitude, have gentler dips, whereas the middle segment has a steeper subduction angle with an offset of ~20 km relative to the other two adjacent segments. This is typically governed by the surface features of the oceanic slab like the age-variability of the Indian plate along the trench (Figure 1), the existence of active spreading in the back arc beneath the Andaman Sea, and an anomalously strong strain partitioning [e.g., Socquet et al., 2006], in which the oblique Sumatra-Andaman subduction is accommodated. Since the downgoing slabs are in constant motion, a Subduction Transform Edge Propagator [Govers and Wortel, 2005] develops in the downgoing slab as a result of slab roll back. It is possible to observe such features with structural variability in the gradients of the gravity anomaly map (see Figure S9), where, along the arc, the anomaly suddenly diminishes around the place where we found tearing in the slab. Such gaps in slabs have been observed at deeper depths [Hafkenscheid et al., 2001] and mostly attributed to slab detachment. The downgoing rate of the slab is controlled by the presence of 90°E ridge, which meets the trench in the latitude range of about 8–10°N. The differential dips in the Mariana slabs have been attributed to the tearing in the slab [van der Hilst and Seno, 1993]. However, this is the first ever demonstration of tear in the subducting slab along the Andaman-Nicobar Island trench. Such a feature not only affects the local seismicity pattern in the region but also modulates the stress regimes. Results from GPS studies affirm this notion by demonstrating that the largest horizontal displacement (Figure S10) [Chlieh et al., 2007; Gahalaut et al., 2006; Subarya et al., 2006] due to the rupture generated by the 26 December 2004 Sumatra-Andaman earthquake (Mw 9.0–9.3) occurred exactly in the area where we have imaged the tearing of the lithosphere. It seems that such close correlation between the slip and the plate configuration may indicate the irregularity of the surface that defines the top of the subducting Indian plate. The interesting point is the correlation between the two independent observations; however, at present their relationship cannot be ascertained. It can only be authenticated by more complex modeling incorporating realistic structure of the region. Further, such tearing of the plate may provide clues about the anatomy of the Indian oceanic plate. The plate reconstruction models suggest eastward motion of the plate with varying degree of obliquity and rate of subduction [Sieh and Natawidjaja, 2000]. The subduction initiated due to the convergence between the Australia/Indian plate and the Eurasian/Sunda plates orthogonally. However, the plate progressively moved obliquely westward in order to accommodate the arc [Curray et al., 1979; McCaffrey, 1992]. Such oblique subduction pattern implies the southward migration of the bottom tip of the northern slab diminishing the slab gap over time in the deeper part of the mantle. We further propose that such slab discontinuities could be the locales for the intense seismicity. For example, the 2005-swarm activity, due to the rupturing generated by the 26 December 2004 earthquake (Mw 9.0–9.3), occurred exactly in the region of tearing. Acknowledgments We thank INCOIS (Ministry of Earth Sciences) for funding the project. Seismic data analysis was performed in Seismic Handler (K. Stammler, http://www.seismic-handler.org/). Plots are generated using Generic Mapping Tool [Wessel and Smith, 1965]. We thank. V.K. Gahalaut for fruitful discussions. The manuscript has been benefited immensely from the constructive comments by two anonymous reviewers. Michael Wysession has handled the manuscript excellently. KUMAR ET AL. References Alam, M. A., D. Chandrasekharam, O. Vaselli, B. Capaccioni, P. Manetti, and A. P. Santo (2004), Petrology of the prehistoric flows and dyke of the Barren Island, Andaman Sea, Indian Ocean, in Proceedings of the Indian Academy of Sciencs (Earth and Planetary Sciences), Special Issue: Magmatism in India Through Time, vol. 113, pp. 715–721. Barker, P. F. (2001), Scotia Sea regional tectonic evolution: Implications for mantle flow and palaeocirculation, Earth Sci. Rev., 55(1–2), 1–39. Benoit, M. H., M. Torpey, K. Liszewski, V. Levin, and J. Park (2011), P and S wave upper mantle seismic velocity structure beneath the northern Apennines: New evidence for the end of subduction, Geochem. Geophys. Geosyst., 12, Q06004, doi:10.1029/2010GC003428. Berkhout, A. J. (1977), Least square inverse filtering and wavelet decomposition, Geophysics, 42, 1369–1383. Chatelain, J.-L., P. Molnar, R. Prevot, and B. Isacks (1992), Detachment of part of the downgoing slab and uplift of the New Hebrides (Vanuatu) islands, Geophys. J. Int., 30(19), 1507–1510. Clark, S. A., M. Sobiesiak, C. A. Zelt, M. B. Magnani, M. S. Miller, M. J. Bezada, and A. Levander (2008), Identification and tectonic implications of a tear in the South American plate at the southern end of the Lesser Antilles, Geochem. Geophys. Geosyst., 9, Q11004, doi:10.1029/2008GC002084. Chlieh, M., et al. (2007), Coseismic slip and afterslip of the great Mw 9.15 Sumatra-Andaman earthquake of 2004, Bull. Seismol. Soc. Am., 97(1A), S152–S173, doi:10.1785/0120050631. Curray, J. R. (2005), Tectonics and history of the Andaman Sea region, J. Asian Earth Sci., 25, 187–232. Curray, J. R., D. Moore, L. Lawyer, F. Emmel, R. Raitt, M. Henry, and R. Kieckhefer (1979), Tectonics of the Andaman sea and Burma, in Geological and Geophysical Investigations of Continental Margins, edited by J. S. Watkins, L. Montadert, and P. Dickerson, Am. Assoc. Petrol. Geol. Mem., 29, 189–198, AAPG, Tulsa, Okla. Dasgupta, S., M. Mukhopadhyay, A. Bhattacharya, and T. K. Jana (2003), The geometry of the Bur-mese-Andaman subducting lithosphere, J. Seismolog., 7, 155–174. Farra, V., and L. Vinnik (2000), Upper mantle stratification by P and S receiver functions, Geophys. J. Int., 141, 699–712. Forsyth, D. W. (1975), Fault plane solutions and tectonics of the South Atlantic and Scotia Sea, J. Geophys. Res., 80, 1429–1443, doi:10.1029/ JB080i011p01429. TEARING IN THE INDIAN OCEANIC PLATE 4903 Geophysical Research Letters 10.1002/2016GL068590 Gahalaut, V. K., B. Nagarajan, J. K. Catherine, and S. Kumar (2006), Constraints on 2004 Sumatra-Andaman earthquake rupture from GPS measurements in Andaman-Nicobar Islands, Earth Planet. Sci. Lett., 242, 365–374. Garfunkel, Z., C. A. Anderson, and G. Schubert (1986), Mantle circulation and the lateral migration of subducted slabs, J. Geophys. Res., 91, 7205–7223, doi:10.1029/JB091iB07p07205. Gasparon, M., G. Rosenbaum, J. Wijbrans, and P. Manetti (2009), The transition from subduction arc to slab tearing: Evidence from Capraia Island, northern Tyrrhenian Sea, J. Geodyn., 47, 30–38. Giacomuzzi, G., M. P. De Gori Civalleri, and C. Chiarabba (2012), A 3D Vs model of the upper mantle beneath Italy: Insight on the geodynamics of central Mediterranean, Earth Planet. Sci. Lett., 335–336, 105–120. Govers, R., and M. J. R. Wortel (2005), Lithosphere tearing at STEP-faults: Response to edges of subduction zones, Earth and Planet. Sci. Lett., 236, 505–523. Gvirtzman, Z., and A. Nur (1999), The formation of Mount Etna as the consequence of slab rollback, Nature, 401, 782–785. Hafkenscheid, E., S. J. H. Buiter, M. J. R. Wortel, W. Spakman, and H. Bijwaard (2001), Modelling the seismic velocity structure beneath Indonesia: A comparison with tomography, Tectonophysics, 333, 35–46. Isacks, B., L. Sykes, and J. Oliver (1969), Focal mechanisms of deep and shallow earthquakes in the Tonga-Kermadec region and the tectonics of island arcs, Geol. Soc. of Am. Bull., 80, 1443–1470. Jarrard, R. D. (1986), Relations among subduction parameters, Rev. Geophys., 24, 217–284, doi:10.1029/RG024i002p00217. Kamesh Raju, K. A., T. Ramprasad, P. S. Rao, B. R. Rao, and J. Varghese (2004), New insights into the tectonic evolution of the Andaman Basin, northeast Indian Ocean, Earth Planet. Sci. Lett., 221, 145–162. Kawakatsu, H., P. Kumar, Y. Takei, M. Shinohara, T. Kanazawa, E. Araki, and K. Suyehiro (2009), Seismic evidence for sharp lithosphereasthenosphere boundaries of oceanic plates, Science, 324, 499–502, doi:10.1126/science.1169499. Kennett, B. L. N., and E. R. Engdahl (1991), Traveltimes for global earthquake location and phase identification, Geophys. J. Int., 105, 429–465. Khan, P. K., and P. P. Chakraborty (2005), Two-phase opening of Andaman Sea: A new seismotectonic insight, Earth Planet. Sci. Lett., 229, 259–271. Kodaira, S., N. Takahashi, A. Nakanishi, S. Miura, and Y. Kaneda (2000), Subducted seamount imaged in the rupture zone of the 1946 Nankaido earthquake, Science, 289(5476), 104–106, doi:10.1126/science.289.5476.104. Kumar, P., and H. Kawakatsu (2011), Imaging the seismic lithosphere-asthenosphere boundary of the oceanic plate, Geochem. Geophys. Geosyst., 12, Q01006, doi:10.1029/2010GC003358. Kumar, P., X. Yuan, R. Kind, and J. Ni (2006), Imaging the colliding Indian and Asian lithospheric plates beneath Tibet, J. Geophys. Res., 111, B06308, doi:10.1029/2005JB003930. Lucente, F. P., C. Chiarabba, G. B. Cimini, and D. Giardini (1999), Tomographic constraints on the geodynamic evolution of the Italian region, J. Geophys. Res., 104, 20,307–20,327, doi:10.1029/1999JB900147. McCaffrey, R. (1992), Oblique plate convergence, slip vectors, and forearc deformation, J. Geophys. Res., 97, 8905–8915, doi:10.1029/ 92JB00483. Meighan, H. E., U. ten Brink, and J. Pulliam (2013), Slab tears and intermediate-depth seismicity, Geophys. Res. Lett., 40, 1–5, doi:10.1002/ grl.50830. Millen, D. W., and M. W. Hamburger (1998), Seismological evidence for tearing of the Pacific plate at the northerm termination of the Tonga subduction zone, Geology, 26(7), 659–662. Pesicek, J. D., C. H. Thurber, H. Zhang, H. R. DeShon, E. R. Engdahl, and S. Widiyantoro (2010), Teleseismic double-difference relocation of earthquakes along the Sumatra-Andaman subduction zone using a 3-D model, J. Geophys. Res., 115, B10303, doi:10.1029/2010JB007443. Reyners, M., K. Gledhill, and D. Waters (1991), Tearing of the subducted Australian plate during the Te Anau, New Zealand, earthquake of 1988 June 3, Geophys. J. Int., 104, 105–115. Richards, S., G. Lister, and B. Kennett (2007), A slab in depth: Three-dimensional geometry and evolution of the Indo-Australian plate, Geochem. Geophys. Geosyst., 8, Q12003, doi:10.1029/2007GC001657. Rosenbaum, G., and N. Piana Agostinetti (2015), Crustal and upper mantle responses to lithospheric segmentation in the northern Apennines, Tectonics, 34, 648–661, doi:10.1002/2013TC0034981. Rosenbaum, G., M. Gasparon, F. P. Lucente, A. Peccerillo, and M. S. Miller (2008), Kinematics of slab tear faults during subduction segmentation and implications for Italian magmatism, Tectonics, 27, TC2008, doi:10.1029/2007TC002143. Ryan, H. F., and D. W. Scholl (1993), Geologic implications of great interplate earthquakes along the Aleutian arc, J. Geophys. Res., 98, 22,135–22,146, doi:10.1029/93JB02451. Schellart, W. P., J. Freeman, D. R. Stegman, L. Moresi, and D. May (2007), Evolution and diversity of subduction zones controlled by slab width, Nature, 446, 308–311, doi:10.1038/nature05615. Shapiro, N. M., M. H. Ritzwoller, and E. R. Engdahl (2008), Structural context of the great Sumatra-Andaman Islands earthquake, Geophys. Res. Lett., 35, L05301, doi:10.1029/2008GL033381. Sieh, K., and D. Natawidjaja (2000), Neotectonics of the Sumatran fault, Indonesia, J. Geophys. Res., 105, 295–326, doi:10.1029/2000JB900120. Socquet, A., C. Vigny, N. Chamot-Rooke, W. Simons, C. Rangin, and B. Ambrosius (2006), India and Sunda plates motion and deformation along their boundary in Myanmar determined by GPS, J. Geophys. Res., 111, B05406, doi:10.1029/2005JB003877. Spence, W. (1977), The Aleutian are: tectonic blocks, episodic subduction, strain diffusion, and magma generation, J. Geophys. Res., 82, 213–230, doi:10.1029/JB082i002p00213. Subarya, C., M. Chlieh, L. Prawirodirdjo, J.-P. Avouac, Y. Bock, K. Sieh, A. J. Meltzner, Natawidjaja, H. Danny, and R. McCaffrey (2006), Plateboundary deformation associated with the great Sumatra-Andaman earthquake, Nature, 440, 46–51. van der Hilst, R. D., and T. Seno (1993), Effects of relative plate motion on the deep structure and penetration depth of slabs below the Izu-Bonin and Mariana island arcs, Earth Planet. Sci. Lett., 120, 395–407. Wessel, P., and W. H. F. Smith (1998), New improved version of generic mapping tools released, Eos Trans. AGU, 79(47), 579, doi:10.1029/ 98EO00426. Wilson, J. T. (1965), A new class of faults and their bearing on continental drift, Nature, 207, 343–347. Yáñez, G., and J. Cembrano (2004), Role of viscous plate coupling in the late Tertiary Andean tectonics, J. Geophys. Res., 109, B02407, doi:10.1029/2003JB002494. KUMAR ET AL. TEARING IN THE INDIAN OCEANIC PLATE 4904