Magma plumbing beneath Anak Krakatau volcano, Indonesia: evidence for multiple magma storage regions

Understanding magma plumbing is essential for predicting the behaviour of explosive volcanoes. We investigate magma plumbing at the highly active Anak Krakatau volcano (Indonesia), situated on the rim of the 1883 Krakatau caldera by employing a suite of thermobarometric models. These include clinopyroxene-melt thermobarometry, plagioclase-melt thermobarometry, clinopyroxene composition barometry and olivine-melt thermometry. Petrological studies have previously identified shallow magma storage in the region of 2–8 km beneath Krakatau, while existing seismic evidence points towards mid- to deep-crustal storage zone(s), at 9 and 22 km, respectively. Our results show that clinopyroxene in Anak Krakatau lavas crystallized at a depth of 7–12 km, while plagioclase records both shallow crustal (3–7 km) and sub-Moho (23–28 km) levels of crystallization. These magma storage regions coincide with well-constrained major lithological boundaries in the crust, implying that magma ascent and storage at Anak Krakatau is strongly controlled by crustal properties. A tandem seismic tomography survey independently identified a separate upper crustal (<7 km) and a lower to mid-crustal magma storage region (>7 km). Both petrological and seismic methods are sensitive in detecting magma bodies in the crust, but suffer from various limitations. Combined geophysical and petrological surveys, in turn, offer increased potential for a comprehensive characterization of magma plumbing at active volcanic complexes.     

The Krakatau islands: The geotectonic setting

The Sunda Strait is located in a transitional zone between two different modes of subduction, the Java frontal and Sumatra oblique subductions. Western Java and Sumatra are, however, geologically continuous.The Krakatau complex lies at the intersection of two graben zones and a north-south active, shallow seismic belt, which coincides with a fracture zone along this seismic belt with fissure extrusion of alkali basaltic rocks commencing at Sukadana and continuing southward as far as the Panaitan island through Rajabasa, Sebuku and Krakatau.Paleomagnetic studies suggest that the island of Sumatra has been rotating clockwise relative to Java from at least 2.0 Ma to the present at a rate of 5–10h/Ma, and therefore the opening of the Sunda Strait might have started before 2 Ma (Nishimura et al. 1986).From geomorphological and seismological studies, it is estimated that the western part of Sumatra has been moving northward along the Semangko fault and the southern part of Sunda Strait has been pulled apart.Assuming that the perpendicular component (58 mm/yr; Fitch 1972) of the oblique subduction has not changed, we can estimate that the subduction started at 7–10 Ma. Huchon and LePichon (1984) also estimated that the subduction started at 13 Ma.Recent crustal earthquakes in the Sunda Strait area are clustered into three groups: (1) beneath the Krakatau complex where they are typically of tectonic origin, (2) inside a graben in the western part of the strait, and (3) in a more diffuse zone south of Sumatra. The individual and composite focal mechanisms of the events inside the strait show an extensional regime. A stress tensor, deduced from the individual focal mechanisms of the Krakatau group shows that the tensional axis has a N 130°E orientation (Harjono et al. 1988).These studies confirm that the Sunda Strait is under a tensional tectonic regime as a result of clockwise rotation along the continental margin and northward movement of the Sumatra sliver plate along the Semangko fault zone.     

Complex slab structure and arc magmatism beneath the Japanese Islands

A dense nationwide seismic network recently constructed in Japan has resulted in the production of a large amount of high-quality data that have enabled the high-resolution imaging of deep seismic structures in the Japanese subduction zone. Seismic tomography, precise locations of earthquakes, and focal mechanism research have allowed the identification of the complex structure of subducting slabs beneath Japan, revealing that the subducting Philippine Sea slab underneath southwestern Japan has an undulatory configuration down to a depth of 60–200 km, and is continuous from Kanto to Kyushu without disruption or splitting, even within areas north of the Izu Peninsula. Analysis of the geometry of the Pacific and Philippine Sea slabs identified a broad contact zone beneath the Kanto Plain that causes anomalously deep interplate and intraslab earthquake activity. Seismic tomographic inversions using both teleseismic and local events provide a clear image of the deep aseismic portion of the Philippine Sea slab beneath the Japan Sea north of Chugoku and Kyushu, and beneath the East China Sea west of Kyushu down to a depth of ∼450 km. Seismic tomography also allowed the identification of an inclined sheet-like seismic low-velocity zone in the mantle wedge beneath Tohoku. A recent seismic tomography work further revealed clear images of similar inclined low-velocity zones in the mantle wedge for almost all other areas of Japan. The presence of the inclined low-velocity zones in the mantle wedge across the entirety of Japan suggests that it is a common feature to all subduction zones. These low-velocity zones may correspond to the upwelling flow portion of subduction-induced convection systems. These upwelling flows reach the Moho directly beneath active volcanic areas, suggesting a link between volcanism and upwelling.     

西南日本火山区三维地震成像研究

利用日本气象厅(JMA)以及日本国立大学联合地震观测台网(JUNEC)记录到的3218个地震事件的231918条P波到时资料,反演求得西南日本160km深度范围内的三维P波速度结构。研究表明,在九州地区,俯冲的菲律宾海板块以高速为主要特征,该海洋板块在30～60km深度处的脱水使得弧前地幔楔顶端的橄榄石蛇纹岩化,在120km深度处的脱水使得地幔楔中的岩石局部熔融,融体上升引起该区的火山活动。在本州西部地区大山火山之下,低速异常显著,并伴随低频地震活动,说明该火山可能是个潜在的活火山,将来有喷发的可能性。     

Structural characteristics off Miyagi forearc region, the Japan Trench seismogenic zone, deduced from a wide-angle reflection and refraction study

The Japan Trench subduction zone, located east of NE Japan, has regional variation in seismicity. Many large earthquakes occurred in the northern part of Japan Trench, but few in the southern part. Off Miyagi region is in the middle of the Japan Trench, where the large earthquakes (M > 7) with thrust mechanisms have occurred at an interval of about 40 years in two parts: inner trench slope and near land. A seismic experiment using 36 ocean bottom seismographs (OBS) and a 12,000 cu. in. airgun array was conducted to determine a detailed, 2D velocity structure in the forearc region off Miyagi. The depth to the Moho is 21 km, at 115 km from the trench axis, and becomes progressively deeper landward. The P-wave velocity of the mantle wedge is 7.9–8.1 km/s, which is typical velocity for uppermost mantle without large serpentinization. The dip angle of oceanic crust is increased from 5–6° near the trench axis to 23° 150 km landward from the trench axis. The P-wave velocity of the oceanic uppermost mantle is as small as 7.7 km/s. This low-velocity oceanic mantle seems to be caused by not a lateral anisotropy but some subduction process. By comparison with the seismicity off Miyagi, the subduction zone can be divided into four parts: 1) Seaward of the trench axis, the seismicity is low and normal fault-type earthquakes occur associated with the destruction of oceanic lithosphere. 2) Beneath the deformed zone landward of the trench axis, the plate boundary is characterized as a stable sliding fault plain. In case of earthquakes, this zone may be tsunamigenic. 3) Below forearc crust where P-wave velocity is almost 6 km/s and larger: this zone is the seismogenic zone below inner trench slope, which is a plate boundary between the forearc and oceanic crusts. 4) Below mantle wedge: the rupture zones of thrust large earthquakes near land (e.g. 1978 off Miyagi earthquake) are located beneath the mantle wedge. The depth of the rupture zones is 30–50 km below sea level. From the comparison, the rupture zones of large earthquakes off Miyagi are limited in two parts: plate boundary between the forearc and oceanic crusts and below mantle wedge. This limitation is a rare case for subduction zone. Although the seismogenic process beneath the mantle wedge is not fully clarified, our observation suggests the two possibilities: earthquake generation at the plate boundary overridden by the mantle wedge without serpentinization or that in the subducting slab.     

Three-dimensional P- and S-wave velocity structures beneath the Japan Islands obtained by high-density seismic stations by seismic tomography

We construct fine-scale 3D P- and S-wave velocity structures of the crust and upper mantle beneath the whole Japan Islands with a unified resolution, where the Pacific (PAC) and Philippine Sea (PHS) plates subduct beneath the Eurasian (EUR) plate. We can detect the low-velocity (low-V) oceanic crust of the PAC and PHS plates at their uppermost part beneath almost all the Japan Islands. The depth limit of the imaged oceanic crust varies with the regions. High-V_P/V_S zones are widely distributed in the lower crust especially beneath the volcanic front, and the high strain rate zones are located at the edge of the extremely high-V_P/V_S zone; however, V_P/V_S at the top of the mantle wedge is not so high. Beneath northern Japan, we can image the high-V subducting PAC plate using the tomographic method without any assumption of velocity discontinuities. We also imaged the heterogeneous structure in the PAC plate, such as the low-V zone considered as the old seamount or the highly seismic zone within the double seismic zone where the seismic fault ruptured by the earthquake connects the upper and lower layer of the double seismic zone. Beneath central Japan, thrust-type small repeating earthquakes occur at the boundary between the EUR and PHS plates and are located at the upper part of the low-V layer that is considered to be the oceanic crust of the PHS plate. In addition to the low-V oceanic crust, the subducting high-V PAC plate is clearly imaged to depths of approximately 250 km and the subducting high-V PHS zone to depths of approximately 180 km is considered to be the PHS plate. Beneath southwestern Japan, the iso-depth lines of the Moho discontinuity in the PHS plate derived by the receiver function method divide the upper low-V layer and lower high-V layer of our model at depths of 30–50 km. Beneath Kyushu, the steeply subducting PHS plate is clearly imaged to depths of approximately 250 km with high velocities. The high-V_P/V_S zone is considered as the lower crust of the EUR plate or the oceanic crust of the PHS plate at depths of 25–35 km and the partially serpentinized mantle wedge of the EUR plate at depths of 30–45 km beneath southwestern Japan. The deep low-frequency nonvolcanic tremors occur at all parts of the high-V_P/V_S zone—within the zone, the seaward side, and the landward side where the PHS plate encounters the mantle wedge of the EUR plate. We prove that we can objectively obtain the fine-scale 3D structure with simple constraints such as only 1D initial velocity model with no velocity discontinuity.     

Deep slab subduction and dehydration and their geodynamic consequences: Evidence from seismology and mineral physics

We present new pieces of evidence from seismology and mineral physics for the existence of low-velocity zones in the deep part of the upper mantle wedge and the mantle transition zone that are caused by fluids from the deep subduction and deep dehydration of the Pacific and Philippine Sea slabs under western Pacific and East Asia. The Pacific slab is subducting beneath the Japan Islands and Japan Sea with intermediate-depth and deep earthquakes down to 600 km depth under the East Asia margin, and the slab becomes stagnant in the mantle transition zone under East China. The western edge of the stagnant Pacific slab is roughly coincident with the NE–SW Daxing'Anling-Taihangshan gravity lineament located west of Beijing, approximately 2000 km away from the Japan Trench. The upper mantle above the stagnant slab under East Asia forms a big mantle wedge (BMW). Corner flow in the BMW and deep slab dehydration may have caused asthenospheric upwelling, lithospheric thinning, continental rift systems, and intraplate volcanism in Northeast Asia. The Philippine Sea slab has subducted down to the mantle transition zone depth under Western Japan and Ryukyu back-arc, though the seismicity within the slab occurs only down to 200–300 km depths. Combining with the corner flow in the mantle wedge, deep dehydration of the subducting Pacific slab has affected the morphology of the subducting Philippine Sea slab and its seismicity under Southwest Japan. Slow anomalies are also found in the mantle under the subducting Pacific slab, which may represent small mantle plumes, or hot upwelling associated with the deep slab subduction. Slab dehydration may also take place after a continental plate subducts into the mantle.     

Morphology of the andean wadati-benioff zone, andesitic volcanism, and tectonic features of the nazca plate

The morphology of the Andean Wadati-Benioff zone south of 10° S was established on the basis of the distribution of earthquake foci. The existence of an intermediate aseismic gap, closely connected with the Andean andesitic volcanism, was confirmed. The gap, interpreted as a partially melted zone, is supposed to be the source of primary magma for active andesitic volcanoes. A clear correlation between the depth range of the Andean Wadati-Benioff zone and the major structural units of the Nazca plate was found. It implicates a non-uniform rate of subduction along the Peru—Chile trench due to the hampering effect related to the main tectonic features of the subducting oceanic plate. The Andean deep earthquakes with focal depths greater than 500 km are interpreted as a remnant of the foregoing cycle of subduction.     

East Asia: Seismotectonics,magmatism and mantle dynamics

In this article, we review the significant recent results of geophysical studies and discuss their implications on seismotectonics, magmatism, and mantle dynamics in East Asia. High-resolution geophysical imaging revealed structural heterogeneities in the source areas of large crustal earthquakes, which may reflect magma and fluids that affected the rupture nucleation of large earthquakes. In subduction zone regions, the crustal fluids originate from the dehydration of the subducting slab. Magmatism in arc and back-arc areas is caused by the corner flow in the mantle wedge and dehydration of the subducting slab. The intraplate magmatism has different origins. The continental volcanoes in Northeast Asia (such as Changbai and Wudalianchi) seem to be caused by the corner flow in the big mantle wedge (BMW) above the stagnant slab in the mantle transition zone and the deep dehydration of the stagnant slab as well. The Tengchong volcano in Southwest China is possibly caused by a similar process in BMW above the subducting Burma microplate (or Indian plate). The Hainan volcano in southernmost China seems to be a hotspot fed by a lower-mantle plume associated with the Pacific and Philippine Sea slabs’ deep subduction in the east and the Indian slab’s deep subduction in the west down to the lower mantle. The occurrence of deep earthquakes under the Japan Sea and the East Asia margin may be related to a metastable olivine wedge in the subducting Pacific slab. The stagnant slab finally collapses down to the bottom of the mantle, which may trigger upwelling of hot mantle materials from the lower mantle to the shallow mantle beneath the subducting slabs and cause the slab–plume interactions. Some of these issues, such as the origin of intraplate magmatism, are still controversial, and so further detailed studies are needed from now.     

What triggers Koyna region earthquakes? Preliminary results from seismic tomography digital array

The cause for prolific seismicity in the Koyna region is a geological enigma. Attempts have been made to link occurrence of these earthquakes with tectonic strain as well as the nearby reservoirs. With a view to providing reliable seismological database for studying the earth structure and the earthquake process in the Koyna region, a state of the art digital seismic network was deployed for twenty months during 1996–97. We present preliminary results from this experiment covering an area of 60 × 80 km² with twenty seismic stations. Hypocentral locations of more than 400 earthquakes confined to 11×25 km² reveal fragmentation in the seismicity pattern — a NE — SW segment has a dip towards NW at approximately 45°, whilst the other two segments show a near vertical trend. These seismic segments have a close linkage with the Western Ghat escarpment and the Warna fault. Ninety per cent of the seismicity is confined within the depth range of 3–10 km. The depth distribution of earthquakes delimits the seismogenic zone with its base at 10 km indicating a transition from an unstable to stable frictional sliding regime. The lack of shallow seismicity between 0 and 3 km indicates a mature fault system with well-developed gouge zones, which inhibit shallow earthquake nucleation. Local earthquake travel time inversion for P- and S-waves show ≈ 2% higher velocity in the seismogenic crust (0–10 km) beneath the epicentral tract relative to a lower velocity (2–3%) in the adjoining region. The high P- and S-wave velocity in the seismogenic crust argues against the presence of high pressure fluid zones and suggests its possible linkage with denser lithology. The zone of high velocity has been traced to deeper depths (≈ 70 km) through teleseismic tomography. The results reveal segmented and matured seismogenic fault systems in the Koyna region where seismicity is possibly controlled by strain build up due to competent lithology in the seismic zone with a deep crustal root.     

Tomographic imaging of hydrated crust and mantle in the subducting Pacific slab beneath Hokkaido,Japan: Evidence for dehydration embrittlement as a cause of intraslab earthquakes

We estimate detailed three-dimensional seismic velocity structures in the subducting Pacific slab beneath Hokkaido, Japan, using a large number of arrival-time data from 6902 local earthquakes. A remarkable low-velocity layer with a thickness of ~ 10 km is imaged at the uppermost part of the slab and is interpreted as hydrated oceanic crust. The layer gradually disappears at depths of 70–80 km, suggesting the breakdown of hydrous minerals there. We find prominent low-velocity anomalies along the lower plane of the double seismic zone and above the aftershock area of the 1993 Kushiro-oki earthquake (M7.8). Since seismic velocities of unmetamorphosed peridotite are much higher than the observations, hydrous minerals are expected to exist in the lower plane as well as the hypocentral area of the 1993 earthquake. On the other hand, regions between the upper and lower planes, where seismic activity is not so high compared to the both planes, show relatively high velocities comparable to those of unmetamorphosed peridotite. Our observations suggest that intermediate-depth earthquakes occur mainly in regions with hydrous minerals, which support dehydration embrittlement hypothesis as a cause of earthquake in the subducting slab.     

Earthquake forerunner as probable precursor — an example from north burma subduction zone

The Burmese Arc seismic activity is not uniform for its ∼ 1100 km length; only the Northern Burmese Arc (NBA) is intensely active. Six large earthquakes in the magnitude range 6.1–7.4 have originated from the NBA Benioff zone between 1954–2011, within an area of 200 × 300 km² where the Indian plate subducts eastward to depths beyond 200 km below the Burma plate. An analysis on seismogenesis of this interplate region suggests that while the subducting lithosphere is characterized by profuse seismicity, seismicity in the overriding plate is rather few. Large earthquakes occurring in the overriding plate are associated with the backarc Shan-Sagaing Fault (SSF) further east. The forecasting performance of the Benioff zone earthquakes in NBA as forerunner is analysed here by: (i) spatial earthquake clustering, (ii) seismic cycles and their temporal quiescence and (iii) the characteristic temporal b-value changes. Three such clusters (C1–C3) are identified from NBA Benioff Zones I & II that are capable of generating earthquakes in the magnitude ranges of 7.38 to 7.93. Seismic cycles evidenced for the Zone I displayed distinct quiescence (Q₁, Q₂ and Q₃) prior to the 6^th August 1988 (M 6.6) earthquake. Similar cycles were used to forecast an earthquake (Dasgupta et al. 2010) to come from the Zone I (cluster C1); which, actually struck on 4 February 2011 (M 6.3). The preparatory activity for an event has already been set in the Zone II and we speculate its occurrence as a large event (M > 6.0) possibly within the year 2012, somewhere close to cluster C3. Temporal analysis of b-value indicates a rise before an ensuing large earthquake.     

Plate subduction,and generation of earthquakes and magmas in Japan as inferred from seismic observations: An overview

A dense nationwide seismic network recently constructed in Japan has been yielding large volumes of high-quality data that have made it possible to investigate the seismic structure in the Japanese subduction zone with unprecedented resolution. In this article, recent studies on the subduction of the Philippine Sea and Pacific plates beneath the Japanese Islands and the mechanism of earthquake and magma generation associated with plate subduction are reviewed. Seismic tomographic studies have shown that the Philippine Sea plate subducting beneath southwest Japan is continuous throughout the entire region, from Kanto to Kyushu, without disruption or splitting even beneath the Izu Peninsula as suggested in the past. The contact of the Philippine Sea plate with the Pacific plate subducting below has been found to cause anomalously deep interplate and intraslab earthquake activity in Kanto. Detailed waveform inversion studies have revealed that the asperity model is applicable to interplate earthquakes. Analyses of dense seismic and GPS network data have confirmed the existence of episodic slow slip accompanied in many instances by low-frequency tremors/earthquakes on the plate interface, which are inferred to play an important role in stress loading at asperities. High-resolution studies of the spatial variation of intraslab seismicity and the seismic velocity structure of the slab crust strongly support the dehydration embrittlement hypothesis for the generation of intraslab earthquakes. Seismic tomography studies have shown that water released by dehydration of the slab and secondary convection in the mantle wedge, mechanically induced by slab subduction, are responsible for magma generation in the Japanese islands. Water of slab origin is also inferred to be responsible for large anelastic local deformation of the arc crust leading to inland crustal earthquakes that return the arc crust to a state of spatially uniform deformation.     

Tomographic imaging of hydrated crust and mantle in the subducting Pacific slab beneath Hokkaido, Japan: Evidence for dehydration embrittlement as a cause of intraslab earthquakes

We estimate detailed three-dimensional seismic velocity structures in the subducting Pacific slab beneath Hokkaido, Japan, using a large number of arrival-time data from 6902 local earthquakes. A remarkable low-velocity layer with a thickness of ~ 10 km is imaged at the uppermost part of the slab and is interpreted as hydrated oceanic crust. The layer gradually disappears at depths of 70–80 km, suggesting the breakdown of hydrous minerals there. We find prominent low-velocity anomalies along the lower plane of the double seismic zone and above the aftershock area of the 1993 Kushiro-oki earthquake (M7.8). Since seismic velocities of unmetamorphosed peridotite are much higher than the observations, hydrous minerals are expected to exist in the lower plane as well as the hypocentral area of the 1993 earthquake. On the other hand, regions between the upper and lower planes, where seismic activity is not so high compared to the both planes, show relatively high velocities comparable to those of unmetamorphosed peridotite. Our observations suggest that intermediate-depth earthquakes occur mainly in regions with hydrous minerals, which support dehydration embrittlement hypothesis as a cause of earthquake in the subducting slab.     

An interpretation of the 1883 cataclysmic eruption of Krakatau from geochemical studies on the partial melting of granite

Pumice flow from the 1883 Krakatau eruption significantly differs in both mineral and chemical compositions from any other volcanic rocks or ejecta of the Krakatau group, which belong to the tholeiitic series. Lithic fragments of granitic Rock, discovered in the pumice flow, are similar to West Malayan granitic rocks. No other granitic rock occurs throughout the Krakatau group, therefore, we consider that the granitic fragments came from the underlying complex at depths, where they were captured as foreign materials by the magma.It is possible that sialic crustal materials plunged into depths along a peculiar tectonic structure located at the Sunda Strait, which appears to be a sheared portion caused by deformation of the Sunda arc due to differential movement between the Indo-Australian oceanic plate and the Eurasian continental crust. The crustal materials were partially melted and produced a magma of granitic composition. The magma was mixed with or assimilated by an ascending basaltic magma originating probably from the upper mantle. This resulted in a dacitic magma distinctly dominant in silica, alkalis and volatile components, and the 1883 Krakatau eruption, characterized by the pumice flow of dacitic composition, took place.     

Seismic slip rates,potential subsurface rupture areas and seismic potential of the Vienna Basin Transfer Fault

The Vienna Basin Transfer Fault (VBTF) is a slow active fault with moderate seismicity (I _max~8–9, M _max~5.7) passing through the most vulnerable regions of Austria and Slovakia. We use different data to constrain the seismic potential of the VBTF including slip values computed from the seismic energy release during the 20th century, geological data on fault segmentation and a depth-extrapolated 3-D model of a generalized fault surface, which is used to define potential rupture zones. The seismic slip of the VBTF as a whole is in the range of 0.22–0.31 mm/year for a seismogenic fault thickness of 8 km. Seismic slip rates for individual segments vary from 0.00 to 0.77 mm/year. Comparing these data to geologically and GPS-derived slip velocities (>1 mm/year) proofs that the fault yields a significant seismic slip deficit. Segments of the fault with high seismic slip contrast from segments with no slip representing locked segments. Fault surfaces of segments within the seismogenic zone (4–14 km depth) vary from 55 to 400 km². Empirical scaling relations show that these segments are sufficiently large to explain both, earthquakes observed in the last centuries, and the 4th century Carnuntum earthquake, for which archeo-seismological data suggest a magnitude of M ≥ 6. Based on the combination of all data (incomplete earthquake catalog, seismic slip deficits, locked segments, potential rupture areas, indications of strong pre-catalog earthquakes) we argue, that the maximum credible earthquake for the VBTF is in the range M _max = 6.0–6.8, significantly larger than the magnitude of the strongest recorded events (M = 5.7).     

P-wave tomography and origin of the Changbai intraplate volcano in Northeast Asia

We present the first seismic image of the upper mantle beneath the active intraplate Changbai volcano in Northeast Asia determined by teleseismic travel time tomography. The data are measured at a new seismic network consisting of 19 portable stations and 3 permanent stations. Our results show a columnar low-velocity anomaly extending to 400-km depth with a P-wave velocity reduction of up to 3%. High velocity anomalies are visible in the mantle transition zone, and deep-focus earthquakes occur at depths of 500–600 km under the region, suggesting that the subducting Pacific slab is stagnant in the transition zone, as imaged clearly by global tomography. These results suggest that the intraplate Changbai volcano is not a hotspot like Hawaii but a kind of back-arc volcano related to the deep subduction and stagnancy of the Pacific slab under Northeast Asia.     

Deep slab subduction and dehydration and their geodynamic consequences: Evidence from seismology and mineral physics

We present new pieces of evidence from seismology and mineral physics for the existence of low-velocity zones in the deep part of the upper mantle wedge and the mantle transition zone that are caused by fluids from the deep subduction and deep dehydration of the Pacific and Philippine Sea slabs under western Pacific and East Asia. The Pacific slab is subducting beneath the Japan Islands and Japan Sea with intermediate-depth and deep earthquakes down to 600 km depth under the East Asia margin, and the slab becomes stagnant in the mantle transition zone under East China. The western edge of the stagnant Pacific slab is roughly coincident with the NE–SW Daxing'Anling-Taihangshan gravity lineament located west of Beijing, approximately 2000 km away from the Japan Trench. The upper mantle above the stagnant slab under East Asia forms a big mantle wedge (BMW). Corner flow in the BMW and deep slab dehydration may have caused asthenospheric upwelling, lithospheric thinning, continental rift systems, and intraplate volcanism in Northeast Asia. The Philippine Sea slab has subducted down to the mantle transition zone depth under Western Japan and Ryukyu back-arc, though the seismicity within the slab occurs only down to 200–300 km depths. Combining with the corner flow in the mantle wedge, deep dehydration of the subducting Pacific slab has affected the morphology of the subducting Philippine Sea slab and its seismicity under Southwest Japan. Slow anomalies are also found in the mantle under the subducting Pacific slab, which may represent small mantle plumes, or hot upwelling associated with the deep slab subduction. Slab dehydration may also take place after a continental plate subducts into the mantle.     

Role of crustal heterogeneity beneath Andaman–Nicobar Islands and its implications for coastal hazard

The Andaman–Nicobar (A–N) Islands region has attracted many geo-scientists because of its unique location and complex geotectonic settings. The recent occurrence of tsunamis due to the megathrust tsunamigenic north Sumatra earthquake (Mw 9.3) with a series of aftershocks in the A–N region caused severe damage to the coastal regions of India and Indonesia. Several pieces of evidence suggest that the occurrence of earthquakes in the A–N region is related to its complex geodynamical processes. In this study, it has been inferred that deep-seated structural heterogeneities related to dehydration of the subducting Indian plate beneath the Island could have induced the process of brittle failure through crustal weakening to contribute immensely to the coastal hazard in the region. The present study based on 3-D P-wave tomography of the entire rupture zone of the A–N region using the aftershocks of the 2004 Sumatra–Andaman earthquake (Mw 9.3) clearly demonstrates the role of crustal heterogeneity in seismogenesis and in causing the strong shakings and tsunamis. The nature and extent of the imaged crustal heterogeneity beneath the A–N region may have facilitated the degree of damage and extent of coastal hazards in the region. The 3-D velocity heterogeneities reflect asperities that manifest what type of seismogenic layers exist beneath the region to dictate the size of earthquakes and thereby they help to assess the extent of earthquake vulnerability in the coastal regions. The inference of this study may be used as one of the potential inputs for assessment of seismic vulnerability to the region, which may be considered for evolving earthquake hazard mitigation model for the coastal areas of the Andaman–Nicobar Islands region.     

Genesis of a new slab tear fault in the Indo-Australian plate,offshore northern Sumatra,Indian ocean

Following the December 2004 and March 2005 major shallow foci inter-plate earthquakes in the north Sumatra region, a slab-tear fault located within the subducting Indian plate ruptured across the West Sunda Trench (WST) within the marginal intra-plate region. Trend, length and movement pattern of this New Tear Fault (NTF) segment is almost identical to another such slab-tear fault mapped previously by Hamilton (1979), located around 160 km south of NTF. Seismic activity along the NTF remained quasi-stable till the end of the year 2011, when an earthquake of magnitude 7.2 occurred on 10.01.2012 just at the tip of NTF, only around ～100 km within the intra-plate domain west of WST. The NTF rupture propagated further towards SSW with the generation of two more large earthquakes on 11.04.2012. The foreshock (10.01.12; M7.2) — mainshock (11.04.12; M 8.6) — aftershock (11.04.12; M 8.2) sequence along with numerous smaller magnitude aftershocks unmistakably define the extension of NTF, a slab-tear fault that results tectonic segmentation of the convergent plate margin. Within the intra-plate domain most earthquakes display consistent left-lateral strike slip mechanism along NNE trending fault plane.