EM Monitoring of Crustal Processes Including the Use of the Network-MT Observations

There are several kinds of coupling mechanisms which can convert mechanical, chemical or thermal energies due to seismic or volcanic activities into electromagnetic energies. As a result of concentrated efforts in laboratory and theoretical research, the basic relationship between the intensity of electromagnetic sources and changes in mechanical, chemical and thermal state is becoming established. Also with the progress of the electromagnetic simulation techniques, it has been possible to evaluate in situ sensitivity. Based on this progress and also due to extensive improvement in measuring techniques, many field experiments have been performed to elucidate subsurface geophysical processes underlying the preparation stage, onset, and subsequent healing stage of earthquakes and volcanic eruptions. In volcanic studies, many studies have reported the measurement of electromagnetic signals which were successfully interpreted in terms of various driving mechanisms. Although there have been numerous reports about the existence of precursory electromagnetic signals in seismic studies, only a few of them could be successfully explained by the proposed mechanisms, whereas coseismic phenomena are often consistent with those mechanisms including the absence of detectable signals. In many cases, one or two orders of higher sensitivity were required, especially for precursory signals. Generally, electromagnetic methods are more sensitive to near-surface phenomena. It will be necessary to discriminate electromagnetic signals due to these near-surface sources, which often possess no relationship with the crustal activities. Further efforts to enhance in situ sensitivity through improvements in observation techniques and in data processing techniques are recommended. At the same time, multi-disciplinary confirmation against the validity of electromagnetic phenomena will inevitably be necessary. A Network-MT observation technique has been developed to determine large-scale deep electrical conductivity structure. In the method, a telephone line network or purpose-built long baseline cables are utilized to measure voltage differences with long electrode separations. Because of the averaging effect of the electric fields, static shift problems due to small-scale, near-surface lateral heterogeneities can be alleviated. Several field experiments revealed regional scale deep electrical conductivity structures related to slab subduction or its stagnation, which enable us to elucidate underlying physical processes caused by the slab motion. The technique can also be applied to monitor the electric potential field related to crustal activities. The annual variation of the potential field and electrical conductivity in the French Alps were interpreted to be caused by the annual variation of lake water level. The method was also used to monitor the regional scale spatio-temporal variation of the SP field and electrical conductivity before and at the onset of earthquakes and volcanic eruptions.     

Finite‐element analysis of laminated rubber bearing of building frame under seismic excitation

Finite element analysis is carried out for a building frame supported by laminated rubber bearings to simultaneously investigate global displacement and local stress responses under seismic excitation. The frame members and the rubber bearings are discretized into hexahedral solid elements with more than 3 million degrees of freedom. The material property of rubber is represented by the Ogden model, and the frame is assumed to remain in elastic range. It is shown that the time histories of non‐uniform stress distribution and rocking behavior of the rubber bearings under a frame subjected to seismic excitation can be successfully evaluated, and detailed responses of base and frame can be evaluated through large‐scale finite element analysis. Copyright © 2015 John Wiley & Sons, Ltd.     

Dynamic FE simulation of four‐story steel frame modeled by solid elements and its validation using results of full‐scale shake‐table test

Dynamic finite element analyses of a four‐story steel building frame modeled as a fine mesh of solid elements are performed using E‐Simulator, which is a parallel finite element analysis software package for precisely simulating collapse behaviors of civil and building structures. E‐Simulator is under development at the National Research Institute for Earth Science and Disaster Prevention (NIED), Japan. A full‐scale shake‐table test for a four‐story frame was conducted using E‐Defense at NIED, which is the largest shaking table in the world. A mesh of the entire structure of a four‐story frame with approximately 19 million degrees of freedom is constructed using solid elements. The density of the mesh is determined by referring to the results of elastic–plastic buckling analyses of a column of the frame using meshes of different densities. Therefore, the analysis model of the frame is well verified. Seismic response analyses under 60, 100, and 115% excitations of the JR Takatori record of the 1995 Hyogoken‐Nanbu earthquake are performed. Note that the simulation does not reproduce the collapse under the 100% excitation of the Takatori record in the E‐Defense test. Therefore, simulations for the 115% case are also performed. The results obtained by E‐Simulator are compared with those obtained by the E‐Defense full‐scale test in order to validate the results obtained by E‐Simulator. The shear forces and interstory drift angles of the first story obtained by the simulation and the test are in good agreement. Both the response of the entire frame and the local deformation as a result of elastic–plastic buckling are simulated simultaneously using E‐Simulator. Copyright © 2014 John Wiley & Sons, Ltd.     

Tidal residual circulation produced by a tidal vortex

“TIDAL VORTEX” is a term for a kind of starting vortex formed as a pair of vortices at the head of a tidal jet emanating from a narrow entrance into a bay. In this study, its formation and movement have been investigated by means of a hydraulic experiment and an analytical model. A tidal vortex is formed as a result of flow separation at the abrupt widening of a channel entrance followed by rolling up of the discontinuity surface around its free end. The vortex shows three types of life-history (type I, II and III), which are characterized by the Strouhal number and the aspect ratio of the horizontal shape of the entrance channel. In the case of type-I, the tidal vortex proceeds toward the inner region of the bay and there amalgamates with successive vortex cores into a core of tidal residual circulation. In the case of type-II, the tidal vortex core flows out into the entrance channel on the ebb but returns back into the bay on the subsequent flood. And, in the case of type-III, the tidal vortex core which was formed on the bay-side opening of the entrance channel flows out to the open sea and never comes back, whereas the core which was formed on the open-sea side of the entrance flows into the bay and never flows out. The circulation of a tidal vortex core is proportional to the reciprocal of the Strouhal number. The movement of the core near the bay entrance is determined by interaction between the cores and transportation due to the irrotational component of the tidal current. There are three types of tidal residual circulation, corresponding to three life-history types of tidal vortices. In the case of type-I, a strong tidal residual circulation is formed, but in type-II a small and weak circulation is formed. While, in type-III, the circulation having an inverse sense of rotation to that of type-I is formed.     

SAR- and gravity change-based characterization of the distribution pattern of pyroclastic flow deposits at Mt. Merapi during the past 10 years

Mt. Merapi, Indonesia, is one of the most active and dangerous volcanoes in the Torrid Zone. This volcano has erupted frequently and has produced pyroclastic flows following the collapse of the summit lava dome. We used Synthetic Aperture Radar (SAR) data acquired by JERS-1 and RADARSAT-1 satellites from April 1996 to July 2006 to clarify the distribution patterns of the pyroclastic flow deposits. The extent of the deposits, termed P-zones, was accurately extracted by ratio operation and low-level feature extraction from SAR intensity images. These images highlighted temporal changes of the distribution area, perimeter, flow distance, included angle, and collapse direction. To validate the image-processing results, reflectance spectra of the rock samples collected after the eruption in June 2006 were measured in a laboratory. The reflectance spectra of all samples showed similar characteristics to the reference spectra, which were derived from atmospheric correction of Hyperion sensor image data covering the lava dome at the summit. Therefore, P-zones were confirmed to be the pyroclastic flow deposits originating from destruction of the lava dome at the summit. The image-processing results clarified that the extent of the distribution areas, perimeter, flow distances, and included angle of the P-zones were variable among the eruptions, while the collapse direction had a constant pattern. The collapse pattern followed a clockwise change from the south toward the west. By comparing the ratio maps of Bouguer gravity anomaly data in two periods, the change was interpreted to originate from the inclination of the conduit and the formation of shallow and deep magma reservoirs.     

Spatial and temporal patterns of seagrass habitat use by fishes at the Ryukyu Islands,Japan

To investigate whether or not regional–temporal patterns of seagrass habitat use by fishes existed at the Ryukyu Islands (southern Japan), visual surveys were conducted in seagrass beds and adjacent coral reefs in northern, central, and southern Ryukyu Islands, in November 2004, and May, August, and November 2005, the northern region having less extensive seagrass beds compared with the central and southern regions. During the study period, the seagrass beds were utilized primarily by 31 species, the densities of some of the latter differing significantly among regions. With the exception of Apogonidae and Holocentridae, all species were diurnal and could be divided into 6 groups based on seagrass habitat use patterns; (1) permanent residents A (10 species, e.g. Stethojulis strigiventer), juveniles and adults living in seagrass beds as well as other habitats; (2) permanent residents B (5 species, e.g. Calotomus spinidens), juveniles and adults living only or mainly in seagrass beds; (3) seasonal residents A (4 species, e.g. Cheilodipterus quinquelineatus), juveniles living in seagrass beds as well as other habitats; (4) seasonal residents B (6 species, e.g. Lethrinus atkinsoni), juveniles living only or mainly in seagrass beds; (5) transients (5 species, e.g. Parupeneus indicus), occurring in seagrass beds in the course of foraging over a variety of habitats; and (6) casual species (1 species, Acanthurus blochii), occurring only occasionally in seagrass beds. Regarding temporal differences, juvenile densities in each group were high in May and August compared with November in each region, whereas adult densities did not differ drastically in each month. For regional differences, juvenile and adult densities of permanent residents A and B were higher in the southern and central regions than in the northern region. Moreover, some seasonal residents showed possible ontogenetic habitat shift from seagrass beds to coral reefs in each region. These results indicated that seagrass habitat use patterns by fishes changed temporally and regionally and there may be habitat connectivity between seagrass beds and coral reefs via ontogenetic migration in the Ryukyu Islands.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of time-space distributions of submarine ground water discharge

Submarine ground water discharge (SGD) rates were measured continuously by automated seepage meters to evaluate the process of ground water discharge to the ocean in the coastal zone of Suruga Bay, Japan. The ratio of terrestrial fresh SGD to total SGD was estimated to be at most 9% by continuous measurements of electrical conductivity of SGD. Semidiurnal changes of SGD due to tidal effects and an inverse relation between SGD and barometric pressure were observed. Power spectrum density analyses of SGD, sea level, and ground water level show that SGD near shore correlated to ground water level changes and SGD offshore correlated to sea level changes. SGD rates near the mouth of the Abe River are smaller than those elsewhere, possibly showing the effect of the river on SGD. The ratio of terrestrial ground water discharge to the total discharge to the ocean was estimated to be 14.7% using a water balance method.     

Petrogenesis of the Kaikomagatake granitoid pluton in the Izu Collision Zone, central Japan: implications for transformation of juvenile oceanic arc into mature continental crust

The Miocene Kaikomagatake pluton is one of the Neogene granitoid plutons exposed in the Izu Collision Zone, which is where the juvenile Izu-Bonin oceanic arc is colliding against the mature Honshu arc. The pluton intrudes into the Cretaceous to Paleogene Shimanto accretionary complex of the Honshu arc along the Itoigawa-Shizuoka Tectonic Line, which is the collisional boundary between the two arcs. The pluton consists of hornblende–biotite granodiorite and biotite monzogranite, and has SiO₂ contents of 68–75 wt%. It has high-K series compositions, and its incompatible element abundances are comparable to the average upper continental crust. Major and trace element compositions of the pluton show well-defined chemical trends. The trends can be interpreted with a crystal fractionation model involving the removal of plagioclase, biotite, hornblende, quartz, apatite, and zircon from a potential parent magma with a composition of ~68 wt% SiO₂. The Sr isotopic compositions, together with the partial melting modeling results, suggest that the parent magma is derived by ~53% melting of a hybrid lower crustal source comprising ~30% Shimanto metasedimentary rocks of the Honshu arc and ~70% K-enriched basaltic rocks of the Izu-Bonin rear-arc region. Together with previous studies on the Izu Collision Zone granitoid plutons, the results of this study suggest that the chemical diversity within the parental magmas of the granitoid plutons reflects the chemical variation of basaltic sources (i.e., across-arc chemical variation in the Izu-Bonin arc), as well as a variable contribution of the metasedimentary component in the lower crustal source regions. In addition, the petrogenetic models of the Izu Collision Zone granitoid plutons collectively suggest that the contribution of the metasedimentary component is required to produce granitoid magma with compositions comparable to the average upper continental crust. The Izu Collision Zone plutons provide an exceptional example of the transformation of a juvenile oceanic arc into mature continental crust.