Imaging of upper crustal structure beneath East Java-Bali, Indonesia with ambient noise tomography

The complex geological structures in East Java and Bali provide important opportunities for natural resource exploitation, but also harbor perils associated with natural disasters. Such a condition makes the East Java region an important area for exploration of the subsurface seismic wave velocity structure, especially in its upper crust. We employed the ambient noise tomography method to image the upper crustal structure under this study area. We used seismic data recorded at 24 seismographs of BMKG spread over East Java and Bali. In addition, we installed 28 portable seismographs in East Java from April 2013 to January 2014 for 2–8 weeks, and we installed an additional 28 seismographs simultaneously throughout East Java from August 2015 to April 2016. We constructed inter-station Rayleigh wave Green’s functions through cross-correlations of the vertical component of seismic noise recordings at 1500 pairs of stations. We used the Neighborhood Algorithm to construct depth profiles of shear wave velocity (Vs). The main result obtained from this study is the thickness of sediment cover. East Java’s southern mountain zone is dominated by higher Vs, the Kendeng basin in the center is dominated by very low Vs, and the Rembang zone (to the North of Kendeng zone) is associated with medium Vs. The existence of structures with oil and gas potential in the Kendeng and Rembang zones can be identified by low Vs.The data used in this study were taken from BMKG (Free) and portable seismometers installed independently by the help of funding from the Australian Department of Foreign Affairs and Trade Grant 71982

The Sumatra subduction zone: A case for a locked fault zone extending into the mantle

A current view is that the portion of the subduction interface that remains locked in the time interval between large interplate earthquakes, hereinafter referred to as the locked fault zone (LFZ), does not extend into the mantle because serpentinization of the mantle wedge would favor stable aseismic sliding. Here, we test this view in the case of the Sumatra subduction zone where the downdip end of the LFZ can be well constrained from the pattern and rate of uplift deduced from coral growth and from GPS measurements of horizontal deformation. These geodetic data are modeled from a creeping dislocation embedded in an elastic half-space and indicate that the LFZ extends 132 ± 10/7 km from the trench, to a depth between 35 and 57 km. By combining this information with the geometry of the plate interface as constrained from two-dimensional gravimetric modeling and seismicity, we show that the LFZ extends below the forearc Moho, which is estimated to lie at a depth of ~30 km, at a horizontal distance of 110 km from the trench. So, in this particular island arc setting, the LFZ most probably extends into the mantle, implying that either the mantle is not serpentinized, or that the presence of serpentine does not necessarily imply stable sliding. From thermal modeling, the temperature at the downdip end of the LFZ is estimated to be 260 ± 100°C. This temperature seems too low for thermally activated ductile flow, so that aseismic slip is most probably due to pressure and/or temperature induced steady state brittle sliding, possibly favored by fluids released from the subducting slab

Shallow structure beneath the Central Volcanic Complex of Tenerife from new gravity data: implications for its evolution and recent reactivation

We present a new local Bouguer anomaly map of the Central Volcanic Complex (CVC) of Tenerife, Spain, constructed from the amalgamation of 323 new high precision gravity measurements with existing gravity data from 361 observations. The new anomaly map images the high-density core of the CVC and the pronounced gravity low centred in the Las Cañadas caldera in greater detail than previously available. Mathematical construction of a sub-surface model from the local anomaly data, employing a 3D inversion based on 'growing' the sub-surface density distribution via the aggregation of cells, enables mapping of the shallow structure beneath the complex, giving unprecedented insights into the sub-surface architecture. We find the resultant density distribution in agreement with geological and other geophysical data. The modelled sub-surface structure supports a vertical collapse origin of the caldera, and maps the headwall of the ca. 180 ka Icod landslide, which appears to lie buried beneath the Pico Viejo–Pico Teide stratovolcanic complex. The results allow us to put into context the recorded ground deformation and gravity changes at the CVC during its reactivation in spring 2004 in relation to its dominant structural building blocks. For example, the areas undergoing the most significant changes at depth in recent years are underlain by low-density material and are aligned along long-standing structural entities, which have shaped this volcanic ocean island over the past few million years

Post-supereruption recovery at Toba Caldera

Large calderas, or supervolcanoes, are sites of the most catastrophic and hazardous events on Earth, yet the temporal details of post-supereruption activity, or resurgence, remain largely unknown, limiting our ability to understand how supervolcanoes work and address their hazards. Toba Caldera, Indonesia, caused the greatest volcanic catastrophe of the last 100 kyr, climactically erupting ~74 ka. Since the supereruption, Toba has been in a state of resurgence but its magmatic and uplift history has remained unclear. Here we reveal that new 14 C, zircon U-Th crystallization and (U-Th)/He ages show resurgence commenced at 69.7±4.5 ka and continued until at least ~2.7 ka, progressing westward across the caldera, as reflected by post-caldera effusive lava eruptions and uplifted lake sediment. The major stratovolcano north of Toba, Sinabung, shows strong geochemical kinship with Toba, and zircons from recent eruption products suggest Toba's climactic magma reservoir extends beneath Sinabung and is being tapped during eruptions

Receiver function study in northern Sumatra and the Malaysian peninsula

International audienceIn this receiver function study, we investigate the structure of the crust beneath six seismic broadband stations close to the Sunda Arc formed by subduction of the Indo-Australian under the Sunda plate. We apply three different methods to analyse receiver functions at single stations. A recently developed algorithm determines absolute shear-wave velocities from observed frequency-dependent apparent incidence angles of P waves. Using waveform inversion of receiver functions and a modified Zhu and Kanamori algorithm, properties of discontinuities such as depth, velocity contrast, and sharpness are determined. The combination of the methods leads to robust results. The approach is validated by synthetic tests. Stations located on Malaysia show high-shear-wave velocities () near the surface in the range of 3.4-3.6 km s attributed to crystalline rocks and 3.6-4.0 km s in the lower crust. Upper and lower crust are clearly separated, the Moho is found at normal depths of 30-34 km where it forms a sharp discontinuity at station KUM or a gradient at stations IPM and KOM. For stations close to the subduction zone (BSI, GSI and PSI) complexity within the crust is high. Near the surface low of 2.6-2.9 km s indicate sediment layers. High of 4.2 km s are found at depth greater than 6 and 2 km at BSI and PSI, respectively. There, the Moho is located at 37 and 40 km depth. At station GSI, situated closest to the trench, the subducting slab is imaged as a north-east dipping structure separated from the sediment layer by a 10 km wide gradient in between 10 and 20 km depth. Within the subducting slab ≈ 4.7 km s. At station BSI, the subducting slab is found at depth between 90 and 110 km dipping 20° ± 8° in approximately N 60° E. A velocity increase in similar depth is indicated at station PSI, however no evidence for a dipping layer is found

Magma plumbing systems: a geophysical perspective

Over the last few decades, significant advances in using geophysical techniques to image the structure of magma plumbing systems have enabled the identification of zones of melt accumulation, crystal mush development, and magma migration. Combining advanced geophysical observations with petrological and geochemical data has arguably revolutionised our understanding of, and afforded exciting new insights into, the development of entire magma plumbing systems. However, divisions between the scales and physical settings over which these geophysical, petrological, and geochemical methods are applied still remain. To characterise some of these differences and promote the benefits of further integration between these methodologies, we provide a review of geophysical techniques and discuss how they can be utilised to provide a structural context for and place physical limits on the chemical evolution of magma plumbing systems. For example, we examine how Interferometric Synthetic Aperture Radar (InSAR), coupled with Global Positioning System (GPS) and Global Navigation Satellite System (GNSS) data, and seismicity may be used to track magma migration in near real-time. We also discuss how seismic imaging, gravimetry and electromagnetic data can identify contemporary melt zones, magma reservoirs and/or crystal mushes. These techniques complement seismic reflection data and rock magnetic analyses that delimit the structure and emplacement of ancient magma plumbing systems. For each of these techniques, with the addition of full-waveform inversion (FWI), the use of Unmanned Aerial Vehicles (UAVs) and the integration of geophysics with numerical modelling, we discuss potential future directions. We show that approaching problems concerning magma plumbing systems from an integrated petrological, geochemical, and geophysical perspective will undoubtedly yield important scientific advances, providing exciting future opportunities for the volcanological community