Ryokei Yoshimura

Preliminary results of multidisciplinary observations before, during and after the Kocaeli (Izmit) earthquake in the western part of the North Anatolian Fault Zone

On August 17, 1999, a destructive earthquake occurred in the western part of the North Anatolian Fault Zone, Turkey. The earthquake source region has been designated as a seismic gap and an M7-class earthquake has been supposed to occur someday in the future so as to fill this seismic gap. So far we have undertaken various kinds of observations in this area and we could obtain some valuable data before, during and after the mainshock. Here we report some of the preliminary results of our recent studies, which include field work started in late July this year and continued during and after the earthquake occurrence just in the earthquake source region and its vicinity, in addition to seismic observations carried out for several years before the mainshock. Much emphasis is put on magnetotelluric field data acquired during the mainshock; in fact, large variations caused by seismic waves were recorded. Such variations could be interpreted in terms of electromagnetic induction in the conducting crust caused by the velocity field interacting with the static magnetic field of the Earth. In particular, the first motion of seismic wave could be identified in the records and used for precise determination of the hypocenter of the mainshock.

Three‐dimensional resistivity structure and magma plumbing system of the Kirishima Volcanoes as inferred from broadband magnetotelluric data

Broadband magnetotelluric (MT) measurements were conducted in 2010 and 2011 in the vicinity of Shinmoe-dake Volcano in the Kirishima volcano group, Japan, where sub-Plinian eruptions took place 3 times during 26–27 January 2011. By combining the new observations with previous MT data, it is found that an anomalous phase in excess of 90° is commonly observed in the northern sector of the Kirishima volcano group. Because the anomalous phase is not explained by 1-D or 2-D structure with isotropic resistivity media, 3-D inversions were performed. By applying small errors to the anomalous phase, we successfully estimated a 3-D resistivity structure that explains not only the normal data but also the anomalous phase data. The final model shows a vertical conductor that is located between a deep-seated conductive body (at a depth greater than 10 km) and a shallow conductive layer. By applying the findings of geophysical and petrological studies of the 2011 sub-Plinian eruptions, we infer that the subvertical conductor represents a zone of hydrothermal aqueous fluids at temperatures over 400°C, in which a magma pathway (interconnected melt) is partially and occasionally formed before magmatic eruptions. To the north of the deep conductor, earthquake swarms occurred from 1968 to 1969, suggesting that these earthquakes were caused by volcanic fluids.

Magnetotelluric observations around the focal region of the 2007 Noto Hanto Earthquake (Mj 6.9), Central Japan

On 25 March 2007, a damaging earthquake (Mj 6.9) occurred near the west coast of the Noto Peninsula, Central Japan. A wideband magnetotelluric (MT) survey was carried out in the onshore area of the source region immediately after the mainshock, with the aim of imaging the heterogeneity of the crustal resistivity structure. The final observation network had consisted of 26 sites. As a preparatory step for imaging three-dimensional features of the resistivity around the focal region, we constructed two-dimensional resistivity models along five profiles using only the TM mode responses, in order to reduce three-dimensional effects. Four profiles are perpendicular to the fault strike, and a fifth profile is parallel to the strike through the mainshock epicenter. Significant characteristics of the resistivity models are: (1) beneath the mainshock hypocenter, there is a conductive body which spreads to the eastern edge of the active aftershock region; (2) a resistive zone is located in the gap of the aftershock distribution between the mainshock hypocenter and the largest eastern aftershock; (3) one of the largest aftershock occurred at the boundary of the resistive zone described above. These results suggest that the deep conductors represent fluid-filled zones and that the lateral heterogeneity could have controlled the slip distribution on the fault plane.

Three-dimensional electromagnetic imaging of upwelling fluids in the Kyushu subduction zone, Japan

A three-dimensional (3-D) lithospheric-scale electrical resistivity model, developed using network-magnetotelluric (network-MT) data, contains structures associated with arc magmatism beneath Kyushu Island in the Southwest Japan arc. Kyushu Island, where the Philippine Sea Plate (PSP) subducts beneath the Eurasian plate, can be divided into northern and southern volcanic regions separated by a nonvolcanic region. Many active Quaternary volcanoes occur along the volcanic front (VF) associated with the PSP in the two volcanic regions. Our 3-D electrical resistivity model shows three different shapes of upwelling fluid-like conductive anomalies, indicative of either slab-derived aqueous fluid and/or partial melt beneath the volcanic and nonvolcanic regions. A conductive anomaly in the northern volcanic region, located at some distance from the subducting PSP, extends from the surface to depths of 70 km depth, extends from the surface to depths of >100 km. In the nonvolcanic region, the upper region of a relatively conductive anomaly extends upward to a depth of ~50 km along the subducting plate. The degrees of magmatism and the relative contribution of slab-derived fluids to the magmatism vary spatially in the one nonvolcanic and two volcanic regions.

Modification of the Network-MT method and its first application in imaging the deep conductivity structure beneath the Kii Peninsula, southwestern Japan

The Network-Magnetotelluric (NMT) method is well-suited for investigating deep and large-scale conductivity structure; however, application of the method is strongly dependent on the availability of telecommunication facilities (specifically, metallic transmission cables). To overcome the problem posed by the progressive replacement of metallic transmission cables with fiber cables, we developed a modified NMT (modified NMT) method consisting of purpose-built electrodes, making use of local metallic telecommunication lines, without a transmission cable. We first applied this modified NMT method over the Kii Peninsula, southwestern Japan, undertaking two-dimensional conductivity modeling along a transect across the central part of the peninsula. The model is characterized by a large (∼20 km wide and depths of 10–60 km) and highly conductive (< 10 Ω m) zone in the central part of the peninsula between the Conrad discontinuity and the upper surface of the Philippine Sea slab. This zone contains the hypocenters of many deep low-frequency tremors but regular earthquakes are rare. The zone also corresponds to a high-Vp/Vs area. The presence of fluid in the zone plays a key role in the absence of regular earthquakes, occurrence of deep low-frequency tremors, and elevated Vp/Vs values, as well as enhancing conductivity.

Hydrothermal system in the Tatun Volcano Group, northern Taiwan, inferred from crustal resistivity structure by audio-magnetotellurics

The present study proposes an improved conceptual model for the hydrothermal system in the Tatun Volcano Group in northern Taiwan. In the study, audio-magnetotellurics (AMT) surveys were conducted to reveal the spatial distribution of resistivity, which is highly sensitive to fluids and hydrothermal alteration. By combining the obtained resistivity structure with other geophysical and geochemical evidence, the following hydrothermal system was inferred. Beneath Chishinshan, vapor-dominant hydrothermal fluids, supplied from a deeper part, are maintained in a low to relatively low resistivity region (5 to 20 Ω m) that is covered by a clay-rich cap, represented by an upper extremely low resistivity layer. Fluid ascent is suggested by a pressure source and clustered seismicity. Exsolved gases result in fumarolic areas, such as Siao-you-keng, while mixing of gases with shallow groundwater forms a shallow flow system of hydrothermal fluids in the Matsao area, represented by a region of less than 10 Ω m. The fumarole in the Da-you-keng area originates from vapor-dominant hydrothermal fluids that may be supplied from a deeper part beneath Cing-tian-gang, suggested by a pressure source and low to relatively low resistivity. Horizontally extended vapor-bearing regions also suggest the possibility of future phreatic eruptions. The proposed conceptual model may provide clues to detecting precursors of potential volcanic activity.

Magnetotelluric and temperature monitoring after the 2011 sub-Plinian eruptions of Shinmoe-dake volcano

Three sub-Plinian eruptions took place on 26–27 January 2011 at Shinmoe-dake volcano in the Kirishima volcanic group, Japan. During this event, GPS and tiltmeters detected syn-eruptive ground subsidence approximately 7 km to the WNW of the volcano. Starting in March 2011, we conducted broad-band magnetotelluric (MT) measurements at a site located 5 km NNW of the volcano, beneath which the Shinmoe-dake magma plumbing system may exist. In addition, temperature monitoring of fumaroles and hot-springs near the MT site was initiated in July 2011. Our MT data record changes in apparent resistivity of approximately ±5%, along with a ±1° phase change in the off-diagonal component of the impedance tensor (Zxy and Zyx ). Using 1-D inversion, we infer that these slight changes in resistivity took place at relatively shallow depths of only a few hundred meters, at the transition between a near-surface resistive layer and an underlying conductive layer. Resistivity changes observed since March 2012 are correlated with the observed temperature increases around the MT monitoring site. These observations suggest the existence beneath the MT site of pathways which enable volatile escape.

Resistivity characterisation of Hakone volcano, Central Japan, by three-dimensional magnetotelluric inversion

On 29 June 2015, a small phreatic eruption occurred at Hakone volcano, Central Japan, forming several vents in the Owakudani geothermal area on the northern slope of the central cones. Intense earthquake swarm activity and geodetic signals corresponding to the 2015 eruption were also observed within the Hakone caldera. To complement these observations and to characterise the shallow resistivity structure of Hakone caldera, we carried out a three-dimensional inversion of magnetotelluric measurement data acquired at 64 sites across the region. We utilised an unstructured tetrahedral mesh for the inversion code of the edge-based finite element method to account for the steep topography of the region during the inversion process. The main features of the best-fit three-dimensional model are a bell-shaped conductor, the bottom of which shows good agreement with the upper limit of seismicity, beneath the central cones and the Owakudani geothermal area, and several buried bowl-shaped conductive zones beneath the Gora and Kojiri areas. We infer that the main bell-shaped conductor represents a hydrothermally altered zone that acts as a cap or seal to resist the upwelling of volcanic fluids. Enhanced volcanic activity may cause volcanic fluids to pass through the resistive body surrounded by the altered zone and thus promote brittle failure within the resistive body. The overlapping locations of the bowl-shaped conductors, the buried caldera structures and the presence of sodium-chloride-rich hot springs indicate that the conductors represent porous media saturated by high-salinity hot spring waters. The linear clusters of earthquake swarms beneath the Kojiri area may indicate several weak zones that formed due to these structural contrasts.