Reconstruction of Subsurface Structures from Ambient Seismic Noise in Jakarta Area

Similar documents
Identification of engineering bedrock in Jakarta by using array observations of microtremors

Supplementary Figure 1. Distribution of seismic event locations determined using the final 3-D velocity model. We separate the crust-related

Tomographic imaging of P wave velocity structure beneath the region around Beijing

Preliminary Estimation of Engineering Bedrock Depths from Microtremor Array Measurements in Solo, Central Java, Indonesia

Seismic Velocity Structures beneath the Guntur Volcano Complex, West Java, Derived from Simultaneous Tomographic Inversion and Hypocenter Relocation

Ambient Noise Tomography of Merapi Complex, Central Java, Indonesia: A Preliminary Result

Magnitude 7.5 PALU, INDONESIA

Estimation of S-wave Velocity Structures by Using Microtremor Array Measurements for Subsurface Modeling in Jakarta

A magnitude 7.4 earthquake struck 255 km (158 miles) southwest of Tonga, according to the US Geological Survey, but there were no reports of damage.

We A10 12 Common Reflection Angle Migration Revealing the Complex Deformation Structure beneath Forearc Basin in the Nankai Trough

INVESTIGATION OF SUBSURFACE CHARACTERISTICS BY USING A V s30 PARAMETER AND A COMBINATION OF THE HVSR AND SPAC METHODS FOR MICROTREMOR ARRAYS

V. Sajama, Bolivia. Kevin M. Ward, George Zandt, Susan L. Beck, Ryan C. Porter, Lara S. Wagner, Estela Minaya, Hernando Tavera

Seismic and Tsunami Observations of Indonesia

Supporting Online Material for

Locating Events using Borehole Microseismic Monitoring by Inclusion of Particle Motion Analysis

Geophysical Journal International

STUDY ON TSUNAMIGENIC EARTHQUAKE CRITERIA FOR THE INDONESIAN TSUNAMI EARLY WARNING SYSTEM

M. Koch and T.H. Münch. Department of Geohydraulics and Engineering Hydrology University of Kassel Kurt-Wolters-Strasse 3 D Kassel

Earthquake patterns in the Flinders Ranges - Temporary network , preliminary results

NUMERICAL SIMULATIONS FOR TSUNAMI FORECASTING AT PADANG CITY USING OFFSHORE TSUNAMI SENSORS

SHORT PERIOD SURFACE WAVE DISPERSION FROM AMBIENT NOISE TOMOGRAPHY IN WESTERN CHINA. Sponsored by National Nuclear Security Administration 1,2

Seismic imaging and petrology explain highly explosive eruptions of Merapi Volcano, Indonesia

Ambient Noise Tomography in the Western US using Data from the EarthScope/USArray Transportable Array

2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

Magnitude 7.1 PHILIPPINES

Geological structure analysis in Central Java using travel time tomography technique of S waves

Crustal Velocity Structure from Surface Wave Dispersion Tomography in the Indian Himalaya

Supporting Information for An automatically updated S-wave model of the upper mantle and the depth extent of azimuthal anisotropy

Imaging Moho topography beneath the Alps by multdisciplinary seismic tomography

Tomography of the 2011 Iwaki earthquake (M 7.0) and Fukushima

Contents of this file

USArray the first five years

The Japanese University Joint Seismic Observations at the Niigaka-Kobe Tectonic Zone

Effects of Surface Geology on Seismic Motion

Dynamic Crust Practice

SURFACE WAVE GROUP VELOCITY MEASUREMENTS ACROSS EURASIA

SOEE3250/5675/5115 Inverse Theory Lecture 10; notes by G. Houseman

Regional 3D velocity structure

5. Gravity. 5.1 Geoid Variations. The Australian Continent: A Geophysical Synthesis Gravity

Dynamic Triggering Semi-Volcanic Tremor in Japanese Volcanic Region by The 2016 Mw 7.0 Kumamoto Earthquake

9th Workshop on Three-Dimensional Modelling of Seismic Waves Generation, Propagation and their Inversion

Th N Seismic Imaging in Gas Obscured Areas - Q Anomaly Detection and Q Migration Applied to Broadband Data

The stratification of seismic azimuthal anisotropy in the western US

The Structure of the Earth and Plate Tectonics

7.1 FIJI 1, :57:22 UTC

Coupling of Wave and Hydrodynamic Models for Predicting Coastal Inundation: A case study in Jakarta and Semarang

Determination of the seismicity and peak ground acceleration for Lombok Island: An evaluation on tectonic setting

Tu N Fault Shadow Removal over Timor Trough Using Broadband Seismic, FWI and Fault Constrained Tomography

ASSESSMENT OF A POTENTIAL HYDROCARBON RESERVOIR IN A YOGYAKARTA BASIN USING BROADBAND-NOISE ANALYSIS

GEO-DEEP9300 Lithosphere and Asthenosphere: Composition and Evolution

Lithospheric structure of West Greenland

DR

Magnitude 7.0 PAPUA, INDONESIA

Data Repository Item

DOUBLE-PLANED STRUCTURE OF INTERMEDIATE- DEPTH SEISMIC ZONE AND THERMAL STRESS IN THE DESCENDING PLATE. (Received December 20, 1983)

Magnitude 7.3 NEPAL. Tuesday, May 12, 2015 at 07:05:19 UTC

SOCIAL STUDIES GENERAL

Magnetic Case Study: Raglan Mine Laura Davis May 24, 2006

Regional Resource Area Mapping in Nevada Using the USArray Seismic Network

Seismic Activity near the Sunda and Andaman Trenches in the Sumatra Subduction Zone

Lecture 4: Earthquakes and Seismic Waves

Magnitude 8.2 NORTHWEST OF IQUIQUE, CHILE

Global geophysics and wave propagation

Long-period Ground Motion Characteristics of the Osaka Sedimentary Basin during the 2011 Great Tohoku Earthquake

Magnitude 7.0 PERU. This region of the Andes is a sparsely populated area, there were no immediate reports of injuries or damage.

Development of Nationwide Vs30 Map and Calibrated Conversion Table for Indonesia using Automated Topographical Classification

2008 Monitoring Research Review: Ground-Based Nuclear Explosion Monitoring Technologies

Determination of prominent crustal discontinuities from waveforms of local earthquakes Pavla Hrubcová Václav Vavryčuk Alena Boušková Josef Horálek

The Coso Geothermal Area: A Laboratory for Advanced MEQ Studies for Geothermal Monitoring

Seismic tomography: Art or science?

Seismic tomography: Art or science? Frederik J Simons Princeton University

Rayleigh Wave Group Velocity Dispersion Across Northern Africa, Southern Europe and the Middle East

Seismic Recording Station TA_109C Summary Report

The Structure of the Earth and Plate Tectonics

Geophysical Journal International

The Subsurface Soil Effects Study Using the Short and Long Predominant Periods From H/V Spectrum In Yogyakarta City

Borehole Strong Motion Observation along the Itoigawa-Shizuoka Tectonic Line

Auxiliary Material. Subduction of oceanic asthenosphere: evidence from sub-slab seismic anisotropy. Teh-Ru Alex Song, Hitoshi Kawakatsu

On the Limitation of Receiver Functions Method: Beyond Conventional Assumptions & Advanced Inversion Techniques

G 3. AN ELECTRONIC JOURNAL OF THE EARTH SCIENCES Published by AGU and the Geochemical Society

Investigation of long period amplifications in the Greater Bangkok basin by microtremor observations

Dynamic Subsidence and Uplift of the Colorado Plateau. Supplementary Material

Seismic Recording Station AZ_PFO Summary Report

NUMERICAL SIMULATION AS GUIDANCE IN MAKING TSUNAMI HAZARD MAP FOR LABUAN ISLAND

Dynamic Earth A B1. Which type of plate boundary is located at the Jordan Fault? (1) divergent (3) convergent (2) subduction (4) transform

Double-difference relocations of the 2004 off the Kii peninsula earthquakes

ESTIMATION OF SEDIMENT THICKNESS BY USING MICROTREMOR OBSERVATIONS AT PALU CITY, INDONESIA. Pyi Soe Thein. 11 November 2013

Seismic scatterers within subducting slab revealed from ambient noise autocorrelation

An Earthquake is a rapid vibration or shaking of the Earth s crust created by a release in energy from sudden movement of a part of a plate along a

Magnitude 7.9 SE of KODIAK, ALASKA

1. University of Ottawa 2. Dublin Institute for Advanced Studies 3. University of Texas at Austin

Seismic Reflection Views of the 1994 Northridge Earthquake Hypocentral Region Using Aftershock Data and Imaging Techniques

THREE-DIMENSIONAL FINITE DIFFERENCE SIMULATION OF LONG-PERIOD GROUND MOTION IN THE KANTO PLAIN, JAPAN

Magnitude 7.3 OFFSHORE EL SALVADOR

Chapter 3. Geology & Tectonics

Magnitude 7.0 VANUATU

Banda Aceh December 26th Earthquake monitored by GPS

FOCAL MECHANISMS OF SUBDUCTION ZONE EARTHQUAKES ALONG THE JAVA TRENCH: PRELIMINARY STUDY FOR THE PSHA FOR YOGYAKARTA REGION, INDONESIA

4-D Geodynamic Modeling With Data Assimilation: Subduction and Continental Evolution

Magnitude 7.0 NEW CALEDONIA

Transcription:

Reconstruction of Subsurface Structures from Ambient Seismic Noise in Jakarta Area Sri Widiyantoro 1 Andri Dian Nugraha 1 Zulfakriza 2 Bayu Pranata 3 1 Global Geophysics Research Group, Faculty of Mining and Petroleum Engineering, ITB 2 Ph. D. student, Faculty of Earth Science and Technology, ITB 3 Meteorological, Climatological and Geophysical Agency (BMKG), Jakarta

Outline Introduction Motivation Data and Method Presentation of Tomograms Discussion Conclusions Further Work

Introduction Indonesia is one of regions with high seismic activity in the world. Understanding of shallow subsurface structure is necessary to reduce risk of seismic hazard. Seismic Ambient Noise Tomography (ANT) represents an efficient tool to image the shallow subsurface structure, especially for Jakarta (as no active seismic sources needed). Dense seismographic stations are available in the Jakarta area.

Motivation Capital city of the Republic of Indonesia Most densely populated territory Figure 1. Tectonic complexity of Jakarta and its vicinity (courtesy of Dr. D. H. Natawidjaja, 2013). Thin red lines depict faults in the study region. Thick red line in the ocean depicts the Sunda trench, where the Indian oceanic lithosphere is being subducted beneath the Sunda arc.

Data Figure 2. Seismographic station network in Jakarta and its surrounding areas that includes 4 permanent stations (black squares) and 8 temporary stations (colored triangles) that were moved around. The difference in the color of triangles indicates the difference in the period of instrument installation.

Waveform Data Figure 3. Example of three-component seismographic recordings at station JN01 for 24 hours; N-S component (top), E-W component (middle), and vertical component (bottom).

Instrument Figure 4. Example of seismometer used in the observation.

Method of Cross Correlation

A Waveform Cross Correlation B Figure 5. A) Example of waveform cross correlation between station JN02 and other stations. B) Example of dispersion curve of stations JN02 and JN26.

Resolution Test A B Figure 6. A) Input pattern for checkerboard test using 60 x 60 grid points (cell size = 1.85 km x 1.85 km). B) Recovery of the checkerboard pattern for period = 1.0 sec; ray path coverage (for period = 1.0 sec) is shown in the inset.

Presentation of Tomograms Figure 7. Tomographic images of ambient seismic noise cross correlation beneath Jakarta and its vicinity produced using conventional method for period = 1 sec, 2.0 sec, 3.0 sec and 4.0 sec.

Discussion A B Figure 8. A) Tomographic images of ambient seismic noise cross correlation beneath Jakarta and its vicinity produced using Bayesian trans-dimensional inversion for periods = 1,0 sec, 2.0 sec, 3.0 sec and 4.0 sec. B) Geological map for Jakarta and surrounding area (Modified from Turkandi et al., 1992 and Effendi et al., 1998).

Discussion (cont.) (a) (b) Legend Layer Depth (m) Vs > 750 m/sec 350-375 375-400 400-425 425-450 450-475 475-500 500-525 525-550 550-575 575-600 600-625 625-650 650-675 675-700 700-725 (c). Legend Layer Depth (m) Vs = 500 m/sec 50-60 60-70 70-80 80-90 90-100. 100-110 110-120 120-130 130-140 140-150 (d) Legend Layer Depth (m) Vs = 700 m/sec 150-175 Legend Layer Depth (m) Vs = 900 m/sec 175-200 200-225 225-250 250-275 275-300 300-325 325-350 350-375. 375-400 400-425 425-450 450-475 475-500 500-525 525-550 350-375 375-400 400-425 425-450 450-475 475-500 500-525 525-550. 550-575 575-600 600-625 625-650 650-675 675-700 700-725 Figure 9. Tomographic images of ambient seismic noise cross correlation vs results of microtremor studies (Ridwan and co-workers, on going project); maps of layer depths for (a) Vs > 750 m/sec, (b) Vs = 500 m/sec, (c) Vs = 700 m/sec, and (d) Vs = 900 m/sec.

Conclusions Resulting tomographic images depict that the study region (Jakarta and its vicinity) is mostly underlain by pronounced low velocity anomalies. The observed low velocity anomalies are related to thick sedimentary layers, in particular in the northern part of Jakarta.

Futher Work Add more seismographic stations in the Jakarta area in order to achieve higher resolution images. Conduct more detailed analysis by applying azimuthal anisotropy using Rayleigh wave data.

Aknowledgments We thank the Asahi Glass Foundation for supporting this research. Thanks also go to P. R. Cummins and E. Saygin for their kind help with the data acquisition and fruitful discussions; M. Ridwan for providing his results from microtremor studies.

Publications (in preparation) Upper crustal structure of Jakarta and its surrounding areas from trans-dimensional seismic ambient noise tomography Azimuthal anisotropic structure beneath Jakarta and its vicinity from ambient noise tomography

Thank you for your kind attention

Supplementary Materials Tomograms for periods ranging from 0.1 sec to 10.0 sec and their associated ray path coverage Checkerboard tests using different cell size, i.e., 1.85 km x 1.85 km, and 2.22 km x 2.22 km

Tomograms and ray path coverage (using 50 x 50 grid points)

Tomograms and ray path coverage (using 70 x 70 grid points)

Checkerboard Test Cell size = 1.85 km x 1.85 km

Checkerboard Test Cell size = 2.22 km x 2.22 km