SUBDUCTION ZONE SEGMENTATION ALONG THE SUNDA MARGIN, INDONESIA Dissertation zur Erlangung des Doktorgrades der Mathematisch-Naturwissenschaftlichen Fakultät der Christian-Albrechts-Universität zu Kiel vorgelegt von Alexey Shulgin Kiel, 2012 Referent: Prof. Dr. H. Kopp Koreferent: Prof. Dr. E.R. Flüh Date of the defense: 21 June 2012 2 Contents Contents ........................................................................................................................ 3 Abstract ......................................................................................................................... 4 Zusammenfassung ......................................................................................................... 6 1. Motivation ............................................................................................................ 8 2. Factors controlling subduction dynamics: the Sunda margin ............................... 8 2.1 Rate and direction of convergence ...................................................................... 9 2.2 Sediment thickness and the seafloor age ............................................................. 9 2.3 Anomalous relief (geometry, topography) of subducting plate ........................ 10 2.4 Forearc high ...................................................................................................... 11 2.5 Structure of the overriding plate ....................................................................... 13 2.6 Seismicity along the Sunda margin ................................................................... 13 2.7 Summary ........................................................................................................... 15 3. The objectives of the study ................................................................................. 20 Part 1: the Sindbad Project offshore eastern Java ................................................... 20 Part 2: the SeaCause project offshore northern Sumatra ........................................ 21 4. Data and methodology ........................................................................................ 21 5. Major Results and Novelty ................................................................................. 23 6. Prospects in subduction system research ............................................................ 27 7. Acknowledgements ............................................................................................ 29 8. References .......................................................................................................... 30 9. Peer-reviewed papers produced during the study ............................................... 32 3 Abstract The Sunda margin marks the eastern termination of the Indian Ocean along the Indonesian archipelago, where the subduction of the Indo-Australian plate below Eurasia is taking place. The convergence direction changes from orthogonal in the eastern part of the margin to a more oblique convergence towards the north-west offshore Sumatra. Along the margin, a wide range of subduction scenarios is observed, which makes this margin a perfect natural laboratory to study the role of various physical and geological parameters on subduction dynamics. This work builds on geophysical and tectonic interpretations of data acquired along eight seismic/gravity and several MCS marine profiles, distributed along the Sunda margin. The data include wide-angle seismic records, multi-channel seismic data, and gravity data. Data interpretation and modeling include a joint reflection/refraction seismic tomography used to constrain geophysical models across the subduction complex at several locations, further constrained by geological information from the multichannel seismic data. Additionally, other geophysical information including gravity modeling and bathymetry data is employed to refine geodynamic models based on seismic interpretations. Excellent data quality allows detailed structural models of the subduction complex to be constrained in several locations along the Sunda margin. The geophysical models obtained in this study show significant variations of the crustal and upper mantle structure of the subduction complex along-strike and across-strike of the margin. The study revealed an increased thickness of the crystalline crust in the Savu Sea, attributed to the approach of the 4 Australian shelf to the trench. Offshore Lombok Island, the oceanic crust thickness is found to be 7 km, and heavily fractured by normal faults. The crustal structure of the Roo Rise oceanic plateau was modeled for the first time, revealing the crustal thickness of 15 km. Its subduction is causing inhomogeneous deformation of the forearc and complex evolution of the entire subduction processes. In addition, large variability in the amount of sediment on the incoming plate is found in adjacent sectors of the margin, ranging from ~3 km in the Savu Sea to almost absence in the Lombok and East Java areas. Offshore north Sumatra, comparison of the crustal-scale profiles located on different sides of the segment boundary between 2004 and 2005 earthquakes rupture areas, showed variability of the structure of the oceanic plate and the amount of sediment present at the trench, resulting in the principally different sizes of the accretionary prisms and the width of the seismogenic zone. A comparison of the subduction complex along the adjacent profiles shows that changes in the structure of the incoming oceanic plate result in drastic changes in the structure of the subduction complex and its evolution. New crustal-scale geophysical models enable analysis of the geodynamic effects caused by the transition from oceanic subduction to continental shelf collision, the presence of anomalous relief on the oceanic plate, variations in the morphology of the incoming oceanic plate, and the presence of segment boundaries, and thus provide a solid basis for geodynamic interpretations of subduction processes and the associated geohazards. 5 Zusammenfassung Der Sunda Kontinentalrand bildet die östliche Begrenzung des Indischen Ozeans entlang des indonesischen Archipels, wo die Subduktion der Indo-Australischen Platte unter Eurasien statt findet. Die Konvergenzrichtung variiert von orthogonal im östlichen Bereich des Kontinentalrandes zu einer schrägen Konvergenzzone im Bereich vor Sumatra. Die Variation von Subduktionsszenarien entlang dieses Kontinentalrandes ermöglicht das Studieren von unterschiedlichen Subduktions-Parametern und deren Auswirkung auf die Dynamik der Subduktionszone. Diese Dissertation integriert geophysikalische und tektonische Interpretationen von acht Seismik/Gravimetrie- und etlichen marinen Mehrkanalseismik - Profilen, die den Sunda Kontinentalrand erfassen. Die Daten umfassen Weitwinkelseismik, Mehrkanalseismik, und Gravimetriedaten. Die Modellierung einer refraktionsseismischen Tomographie inkludiert geologische Informationen aus Mehrkanalseismiken und ermöglicht die Erstellung von geophysikalischen Modellen entlang mehreren Lokalitäten des Subduktionskomplexes. Zusätzliche geophysikalische Daten aus Gravimetriemodellen und bathymetrische Datensätze erlauben eine Verbesserung des ursprünglichen geodynamischen Modells. Durch die exzellente Datenqualität konnten hochauflösende strukturelle Modelle des Subduktionskomplexes im Streichen und quer zum Sunda Kontinentalrand erstellt werden. Diese Modelle zeigen signifikante Unterschiede im strukturellen Aufbau der Kruste und des oberen Erdmantels innerhalb des Subduktionskomplexes. Die Studie zeigt eine Zunahme der kristallinen Krustenmächtigkeit in der Savu See, aufgrund der Annäherung des Australischen 6 Schelfs in Richtung Tiefseerinne. Die ozeanische Kruste vor der Insel Lombok zeigt eine Mächtigkeit von 7 km, die durch starke tektonische Abschiebungsbrüche charakterisiert ist. Die Krustenstruktur des ozeanischen Roo Rise Plateaus wurde in dieser Arbeit erstmalig modelliert und weist eine Mächtigkeit von 15 km auf. Die Subduktion des Plateaus führt zur inhomogenen Deformation des Forearcs sowie zu einer komplexen Entwicklung des gesamten Subduktionsprozesses. Weiterhin wurden große Unterschiede in der Sedimentmächtigkeit auf der einfahrenden Platte beobachtet. Diese reichen von ~3 km in der Savu See bis nur geringfügiger Sedimentbedeckung in den Bereichen um Lombok und Ost-Java. Ein Vergleich verschiedener Krustenprofile vor Nord-Sumatra, entlang der Segmentgrenze der Erdbeben von 2004 und 2005, zeigt Unterschiede in der Struktur der ozeanischen Platte sowie Unterschiede in der Sedimentmächtigkeit in der Tiefseerinne. Dies resultiert in den unterschiedlich großen Akkretionskeilen und der Breite der seismogenen Zone. Ein Vergleich der unterschiedlichen Profile des Subduktionskomplexes zeigt, dass strukturelle Variationen der subduzierten ozeanischen Platte dramatische Veränderungen in der Struktur und Entwicklung des Subduktionskomplexes hervorbringen. Unsere krustenmaßstäblichen geophysikalischen Modelle zeigen die geodynamischen Effekte, die durch den Übergang von ozeanischer Subduktion zu kontinentaler Schelfkollision, das anomale Relief der ozeanischen Platte, Veränderungen in der Morphologie der subduzierten ozeanischen Kruste, sowie die Präsenz von Segmentgrenzen hervorgerufen werden. Diese Modelle bilden eine solide Basis für geodynamische Interpretationen von Subduktionsprozessen und den damit verbundenen Naturgefahren und Georisiken. 7 1. Motivation Subduction zones, as an important end-member in the dynamics of