Seismic evidence of pop-up tectonics beneath the Shillong Plateau area of Northeast India

Our detailed 3-D seismic tomographic assimilation using high-quality phase arrival time data recorded by the local seismographic network demonstrated that heterogeneities in the crustal faults have contributed significantly to the pop-up tectonics beneath the Shillong Plateau, characterized by high-V and low-σ. The major seismogenic faults, namely, the north-dipping Dapsi thrust in association with Dauki fault in the south and south dipping Brahmaputra fault in the north, located either side of the Shillong Plateau that acted as the causative factors for the pop-up, which attributed to the lithostatic (high-V, low-σ) and sedimentary (low-V, high-σ) load, respectively. Seismicity is found confined to a depth ≤ 60 km. Uneven distribution of structural heterogeneities in the upper crust is responsible for earthquake genesis of varying strengths. It is intriguing to note that high-velocity anomalies and low-ϭ in the uppermost crust, interpreted as the Shillong Plateau that acted as a geometric asperity and the juxtaposition of high-V and low-V became the source zone of the 1897 Shillong earthquake (Ms 8.7) as a novel observation for the region. Structural heterogeneities are distinctly distributed between low-V, high-σ and high-V, low-σ in the lower crust plays a major role for future intense seismogenesis due to differential strain accumulation.

Uncertainty of the 2D Resistivity Survey on the Subsurface Cavities

We examined the uncertainty of the two-dimensional (2D) resistivity method using conceptual cavity models. The experimental cavity study was conducted to validate numerical model results. Spatial resolution and sensitivity to resistivity perturbations were also assessed using checkerboard tests. Conceptual models were simulated to generate synthetic resistivity data for dipole-dipole (DD), pole-dipole (PD), Wenner–Schlumberger (WS), and pole-pole (PP) arrays. The synthetically measured resistivity data were inverted to obtain the geoelectric models. The highest anomaly effect (1.46) and variance (24,400 Ω·m) in resistivity data were recovered by the DD array, whereas the PP array obtained the lowest anomaly effect (0.60) and variance (2401 Ω·m) for the shallowest target cavity set at 2.2 m depth. The anomaly effect and variance showed direct dependency on the quality of the inverted models. The DD array provided the highest model resolution that shows relatively distinct anomaly geometries. In contrast, the PD and WS arrays recovered good resolutions, but it is challenging to determine the correct anomaly geometries with them. The PP array reproduced the lowest resolution with less precise anomaly geometries. Moreover, all the tested arrays showed high sensitivity to the resistivity contrasts at shallow depth. The DD and WS arrays displayed the higher sensitivity to the resistivity perturbations compared to the PD and PP arrays. The inverted models showed a reduction in sensitivity, model resolution, and accuracy at deeper depths, creating ambiguity in resistivity model interpretations. Despite these uncertainties, our modeling specified that two-dimensional resistivity imaging is a potential technique to study subsurface cavities. We inferred that the DD array is the most appropriate for cavity surveys. The PD and WS arrays are adequate, while the PP array is the least suitable for cavity studies.

Mantle wedge diapirs detected by a dense seismic array in Northern Taiwan

It is conventionally believed that magma generation beneath the volcanic arc is triggered by the infiltration of fluids or melts derived from the subducted slab. However, recently geochemical analyses argue the arc magma may be formed by mélange diapirs that are physically mixed by sediment, altered oceanic crust, fluids, and mantle above the subducted slab. Further numerical modeling predicts that the mantle wedge diapirs have significant seismic velocity anomalies, even though these have not been observed yet. Here we show that unambiguously later P-waves scattered from some obstacles in the mantle wedge are well recorded at a dense seismic array (Formosa Array) in northern Taiwan. It is the first detection of seismic scattering obstacles in the mantle wedge. Although the exact shape and size of the scattered obstacles are not well constrained by the arrival-times of the later P-waves, the first order approximation of several spheres with radius of ~ 1 km provides a plausible interpretation. Since these obstacles were located just beneath the magma reservoirs around depths between 60 and 95 km, we conclude they may be mantle wedge diapirs that are likely associated with magma generation beneath active volcanoes.

3D anisotropic structure of the Japan subduction zone

How mantle materials flow and how intraslab fabrics align in subduction zones are two essential issues for clarifying material recycling between Earth's interior and surface. Investigating seismic anisotropy is one of a few viable technologies that can directly answer these questions. However, the detailed anisotropic structure of subduction zones is still unclear. Under a general hexagonal symmetry anisotropy assumption, we develop a tomographic method to determine a high-resolution three-dimensional (3D) P wave anisotropic model of the Japan subduction zone by inverting 1,184,018 travel time data of local and teleseismic events. As a result, the 3D anisotropic structure in and around the dipping Pacific slab is firstly revealed. Our results show that slab deformation plays an important role in both mantle flow and intraslab fabric, and the widely observed trench-parallel anisotropy in the forearc is related to the intraslab deformation during the outer-rise yielding of the subducting plate.

New Insights of Geomorphologic and Lithologic Features on Wudalianchi Volcanoes in the Northeastern China from the ASTER Multispectral Data

Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imaging system onboard NASA’s (National Aeronautics and Space Administration’s) Terra satellite is capable of measuring multispectral reflectance of the earth’s surface targets in visible and infrared (VNIR) to shortwave infrared (SWIR) (until 2006) as well as multispectral thermal infrared (TIR) regions. ASTER VNIR stereo imaging technique can provide high-resolution digital elevation models (DEMs) data. The DEMs data, three-dimensional (3D) perspective, and ratio images produced from the ASTER multispectral data are employed to analyze the geomorphologic and lithologic features of Wudalianchi volcanoes in the northeastern China. Our results indicate that the 14 major conical volcanic craters of Wudalianchi volcanoes are arranged as three sub-parallel zones, extending in a NE (Northeast) direction, which is similar to the direction of regional fault system based on the ASTER DEMs data. Among the 14 volcanic craters in Wudalianchi, the Laoheishan, and Huoshaoshan lavas flows, after the historic eruptions, pouring down from the crater, partially blocked the Baihe River, which forms the Five Large Connected Pools, known as the Wudalianchi Lake. Lithologic mapping shows that ASTER multispectral ratio imagery, particularly, the Lava Flow Index (LFI) (LFI = B10/B12) imagery, can clearly distinguish different lava flow units, and at least four stages of volcanic eruptions are revealed in the Wudalianchi Quaternary volcano cluster. Thus, ASTER multispectral TIR data can be used to determine relative dating of Quaternary volcanoes in the semi-arid region. Moreover, ASTER 3D perspective image can present an excellent view for tracking the flow directions of different lavas of Wudalianchi Holocene volcanoes.

Seismic imaging of mantle wedge corner flow and arc magmatism

I reviewed studies on the inhomogeneous seismic structure of the mantle wedge in subduction zones, in relation to corner flow and its implications for arc magmatism. Seismic studies in Tohoku clearly imaged the descending flow portion of the corner flow as a thin seismic low-velocity layer right above the slab. Slab-derived H₂O is fixed to the layer as hydrous minerals, which are brought down by the slab and eventually decompose. The released H₂O rises and encounters the ascending flow, formed to fill the gap caused by the descending flow. The combination of H₂O addition and adiabatic decompression causes partial melting within the ascending flow. For many subduction zones, seismic tomography has distinctly imaged the ascending flow of the corner flow as a seismic low-velocity and/or high-attenuation layer in the mantle wedge inclined nearly parallel to the slab. These observations indicate that the volcanic front in subduction zones is formed both by the ascending flow and the addition of slab-derived H₂O.

Volcanology. The Yellowstone magmatic system from the mantle plume to the upper crust

The Yellowstone supervolcano is one of the largest active continental silicic volcanic fields in the world. An understanding of its properties is key to enhancing our knowledge of volcanic mechanisms and corresponding risk. Using a joint local and teleseismic earthquake P-wave seismic inversion, we revealed a basaltic lower-crustal magma body that provides a magmatic link between the Yellowstone mantle plume and the previously imaged upper-crustal magma reservoir. This lower-crustal magma body has a volume of 46,000 cubic kilometers, ~4.5 times that of the upper-crustal magma reservoir, and contains a melt fraction of ~2%. These estimates are critical to understanding the evolution of bimodal basaltic-rhyolitic volcanism, explaining the magnitude of CO2 discharge, and constraining dynamic models of the magmatic system for volcanic hazard assessment.

Tomography of the Source Area of the 1995 Kobe Earthquake: Evidence for Fluids at the Hypocenter?

Seismic tomography revealed a low seismic velocity (-5%) and high Poisson's ratio (+6%) anomaly covering about 300 square kilometers at the hypocenter of the 17 January 1995, magnitude 7.2, Kobe earthquake in Japan. This anomaly may be due to an overpressurized, fluid-filled, fractured rock matrix that contributed to the initiation of the Kobe earthquake.

Lithospheric Contributions to Arc Magmatism: Isotope Variations Along Strike in Volcanoes of Honshu, Japan

Major chemical exchange between the crust and mantle occurs in subduction zone environments, profoundly affecting the chemical evolution of Earth. The relative contributions of the subducting slab, mantle wedge, and arc lithosphere to the generation of island arc magmas, and ultimately new continental crust, are controversial. Isotopic data for lavas from a transect of volcanoes in a single arc segment of northern Honshu, Japan, have distinct variations coincident with changes in crustal lithology. These data imply that the relatively thin crustal lithosphere is an active geochemical filter for all traversing magmas and is responsible for significant modification of primary mantle melts.