Seismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan

Transcription

1 Seismic structure along transitions from flat to normal subduction: central Mexico, southern Peru, and southwest Japan Thesis by Sara L. Dougherty In Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy California Institute of Technology Pasadena, California 2014 (Defended May 23, 2014)

2 ii c 2014 Sara L. Dougherty All Rights Reserved

3 iii Acknowledgments First and foremost, I would like to sincerely thank my advisors, Rob Clayton and Don Helmberger, without whom this thesis would not have been possible. Don taught me to appreciate every wiggle of a seismogram and introduced me to complex waveforms and the fine-scale structures they can resolve. Rob taught me to appreciate the bigger picture and to explore the possible tectonic implications of observed fine-scale structures. Their support and guidance throughout my time at Caltech has been immensely valuable. I greatly appreciate Rob s patience, encouragement, and accessibility, especially as this thesis was coming together. I would like to thank my academic advisor, Jennifer Jackson, for her assistance with selecting courses and applying for candidacy and for the thought-provoking questions that she poses in my thesis committee meetings. I am also grateful to Joann Stock for serving on my committee and always providing insightful comments. The former grad students and postdocs who helped me in the initial stages of my research and along the way deserve special thanks: Alex Song, Risheng Chu, Min Chen, YoungHee Kim, and Ting Chen. Shengji Wei has also been a great help. I would like to thank Xyoli Pérez-Campos and Steve Skinner for insightful comments and valuable discussions. I am grateful for the support staff in the Seismo Lab and Tectonics Observatory, especially Rosemary Miller, Donna Mireles, Viola Carter, Evelina Cui, and Lisa Christiansen. Thanks also to Mike Black for all the IT support. Special thanks to Jason Saleeby and Jean-Philippe Avouac for leading the Pahoehoe field trip to Hawaii and the division enrichment trip to Greece, respectively, both of which were great experiences. Finally, my deep thanks to my husband, Walt Stelz, for his continuous support and understand-

4 iv ing, especially during the busy times. His ability to always make me laugh has been invaluable, and the occasional coding advice didn t hurt either.

5 v Abstract The fine-scale seismic structure of the central Mexico, southern Peru, and southwest Japan subduction zones is studied using intraslab earthquakes recorded by temporary and permanent regional seismic arrays. The morphology of the transition from flat to normal subduction is explored in central Mexico and southern Peru, while in southwest Japan the spatial coincidence of a thin ultra-slow velocity layer (USL) atop the flat slab with locations of slow slip events (SSEs) is explored. This USL is also observed in central Mexico and southern Peru, where its lateral extent is used as one constraint on the nature of the flat-to-normal transitions. In western central Mexico, I find an edge to this USL which is coincident with the western boundary of the projected Orozco Fracture Zone (OFZ) region. Forward modeling of the 2D structure of the subducted Cocos plate using a finite-difference algorithm provides constraints on the velocity and geometry of the slab s seismic structure in this region and confirms the location of the USL edge. I propose that the Cocos slab is currently fragmenting into a North Cocos plate and a South Cocos plate along the projection of the OFZ, by a process analogous to that which occurred when the Rivera plate separated from the proto-cocos plate 10 Ma. In eastern central Mexico, observations of a sharp transition in slab dip near the abrupt end of the Trans Mexican Volcanic Belt (TMVB) suggest a possible slab tear located within the subducted South Cocos plate. The eastern lateral extent of the USL is found to be coincident with these features and with the western boundary of a zone of decreased seismicity, indicating a change in structure which I interpret as evidence of a possible tear. Analysis of intraslab seismicity patterns and focal mechanism orientations and faulting types provides further support for a possible tear in the South Cocos slab. This potential tear, together with the tear along the projection of the

6 vi OFZ to the northwest, indicates a slab rollback mechanism in which separate slab segments move independently, allowing for mantle flow between the segments. In southern Peru, observations of a gradual increase in slab dip coupled with a lack of any gaps or vertical offsets in the intraslab seismicity suggest a smooth contortion of the slab. Concentrations of focal mechanisms at orientations which are indicative of slab bending are also observed along the change in slab geometry. The lateral extent of the USL atop the horizontal Nazca slab is found to be coincident with the margin of the projected linear continuation of the subducting Nazca Ridge, implying a causal relationship, but not a slab tear. Waveform modeling of the 2D structure in southern Peru provides constraints on the velocity and geometry of the slab s seismic structure and confirms the absence of any tears in the slab. In southwest Japan, I estimate the location of a possible USL along the Philippine Sea slab surface and find this region of low velocity to be coincident with locations of SSEs that have occurred in this region. I interpret the source of the possible USL in this region as fluids dehydrated from the subducting plate, forming a high pore-fluid pressure layer, which would be expected to decrease the coupling on the plate interface and promote SSEs.

7 vii Contents Acknowledgments iii Abstract v List of Figures x List of Tables xiii 1 Introduction 1 2 Seismic structure in central Mexico: Implications for fragmentation of the subducted Cocos plate Abstract Introduction Data Analysis Data D Velocity Modeling Ultra-slow Velocity Layer D Velocity Modeling Seismicity and Slab Dip Across USL Edge Discussion Conclusions Supplemental Figures References

8 viii 3 Seismicity and structure in central Mexico: Evidence for a possible slab tear in the South Cocos plate Abstract Introduction Tectonic Setting Data Analysis Data Slab Dip Ultra-slow Velocity Layer Seismicity D Velocity Modeling D Velocity Modeling Discussion Conclusions Supplemental Figures References Appendix Seismic structure in southern Peru: Evidence for a smooth contortion between flat and normal subduction of the Nazca plate Abstract Introduction Data Analysis Data Seismicity Slab Dip Slab Transition Ultra-slow Velocity Layer

9 ix D Velocity Modeling Discussion Conclusions Supplemental Figures References Further evidence for P-wave complexity in a region with slow slip events Abstract Introduction Data and Method Results Discussion and Conclusions Future Work References Conclusions 173

10 x List of Figures 1.1 Overview map of study locations Map of western central Mexico with event and station locations D P- and S-wave velocity models tested for western central Mexico Comparison of 1D modeling results for the five models tested in western central Mexico D modeling results for the ncm velocity model Profile across the MARS array Raypath schematic and examples of complex P waveforms in western central Mexico Maps of local S-to-P conversion points in western central Mexico D velocity model across the MARS array D modeling results Comparison of 2D modeling results for five different models Seismicity and slab dip across the western USL edge D schematic of two-tear plate model Supplemental: 1D modeling results for the SoCal velocity model Supplemental: 1D modeling results for the SoCalx velocity model Supplemental: 1D modeling results for the Kim et al. (2010) velocity model Supplemental: 1D modeling results for the Song et al. (2009) velocity model Map of eastern central Mexico with event and station locations Topographic-bathymetric map of eastern central Mexico Lateral variation in slab dip across eastern central Mexico

11 xi 3.4 Raypath schematic and examples of complex P waveforms in eastern central Mexico Eastern lateral extent of ultra-slow velocity layer in central Mexico Intraslab seismicity in eastern central Mexico Focal mechanisms for intraslab earthquakes in eastern central Mexico Cross-sections of CAP focal mechanisms in eastern central Mexico D velocity models tested in eastern central Mexico Comparison of 1D modeling results in eastern central Mexico D modeling results for ncmc velocity model D velocity model cross-section locations and results for eastern central Mexico Example 2D velocity model for eastern central Mexico Comparison of 2D modeling results for example event D schematic of slab rollback process and plate tearing through time Supplemental: Map of MASE, VEOX, SSN, and OXNET station locations and names Supplemental: Enlarged map of the eastern lateral extent of the USL in central Mexico with event labels Supplemental: Additional 1D velocity models tested in eastern central Mexico Supplemental: 1D modeling results for ncmc velocity model without station overlap Supplemental: Additional example 2D velocity model tested for eastern central Mexico Supplemental: Comparison of 2D modeling results for additional example event Appendix: Summary map of USL extent in central Mexico Map of southern Peru with event and station locations Topographic-bathymetric map of southern Peru Seismicity pattern in southern Peru Focal mechanisms for intermediate depth earthquakes in southern Peru Lateral variation in slab dip across southern Peru Map of along-strike seismicity in southern Peru Along-strike variation in seismicity in southern Peru

12 xii 4.8 3D view of the flat-to-normal transition in the Nazca slab Raypath schematic and examples of complex P waveforms in southern Peru Lateral extent of ultra-slow velocity layer in southern Peru S-to-P conversion point locations from the top of the Nazca slab D velocity models tested in southern Peru D velocity model cross-section locations D velocity models for event Comparison of 2D modeling results for event D velocity models for event Comparison of 2D modeling results for event Comparison of slab strain in central Mexico and southern Peru Supplemental: Cross-section of seismicity illustrating shallow cluster Supplemental: Complete set of trench-parallel seismicity cross-sections Map of southwest Japan with event and station locations Raypath schematic and examples of complex P waveforms in southwest Japan Raypath schematic and examples of teleseismic underside reflections Map of teleseismic short-period array and broadband station locations Map of local S-to-P conversion points and underside reflection bounce points

13 xiii List of Tables 2.1 Events used in western central Mexico and their source parameters Events used in eastern central Mexico and their source parameters Events used in southern Peru and their source parameters Model parameters and results for each 2D profile in southern Peru