Seismicity changes related to a water injection experiment in the Nojima Fault Zone

Abstract A water injection experiment was carried out by the scientific drilling program named the 'Nojima Fault Zone Probe' during the two periods 9–13 February and 16–25 March 1997. The pumping pressure at the surface was approximately 4 MPa. The total amount of injected water was 258 m³. The injection was made between depths of 1480 m and 1670 m in the Disaster Prevention Research Institute, Kyoto University (DPRI) 1800 m borehole drilled into the Nojima Fault zone. A seismic observation network was deployed to monitor seismic activity related to the water injections. Seismicity suddenly increased in the region not far from the injection hole 4 or 5 days after the beginning of each water injection. These earthquakes were likely to be induced by the water injections. Most of the earthquakes had magnitudes ranging from −2 to +1. Numerous earthquakes occurred during the first injection, but only one could be reliably located and it was approximately 2 km north of the injection site. Between the two injection periods, earthquakes concentrated in the region approximately 1 km northwest of the injection site. During and after the second injection experiment, earthquakes were located approximately 1.5 km west of the injection site. Those earthquakes were located approximately 3 km or 4 km from the injection point and between 2 km and 4 km in depth. Values of intrinsic permeability of 10⁻¹⁴–10⁻¹⁵ m² were estimated from the time lapse of the induced seismic activity. The coefficient of friction in the area where the induced earthquakes occurred was estimated to be less than 0.3.     

Oxygen and carbon isotope zonations of wall rocks around the Kamioka Pb---Zn skarn deposits, central Japan: application to prospecting

Mapping of the oxygen and carbon isotopic composition of hydrothermally altered wall rocks was conducted during blind ore prospecting for Pb---Zn skarn deposits in the Kamioka mining district, central Japan. The wall rocks consist of heterogeneous rock units. Oxygen and carbon isotope ratios were determined for 35 limestones and 33 silicate rocks from the area around the Mozumi deposit (3 km × 3 km) in the Kamioka mining district. The results (δ¹⁸O_spsmow of − 1.1 to + 17.3% and δ¹³C_sppdb of − 5.0 to +4.8% for limestones, and δ¹⁸O_spsmow of −0.8 to + 12.5% for silicate rocks) show isotope zonations of the wall rocks, with lighter isotopic compositions toward the center of mineralization. The isotope zonations likely formed by interaction of thermal waters with the wall rocks during skarn mineralization. The isotopically light zone indicates a higher paleotemperature or higher water-to-rock ratios, and occurs in the footwall of the 7Gohi fault. This structure is spatially related to the distribution of orebodies, indicating that the fault was the main conduit of the ore-forming fluids. The oxygen and carbon isotopic compositions of limestones vary regularly across limestone bodies hosted by the silicate wall rocks, suggesting that the thermal waters were pervasive throughout the wall rocks at the time of mineralization.An isotopically light zone was also found in the southeastern corner of the study area, where significant mineralization had yet to be identified. This suggested an extension of the extinct hydrothermal system to this area, and the possibility of hidden orebodies underneath. Recent drilling in this area has intercepted a zone 45 m thick with a grade of 13.4% Zn, 0.03% Pb and 8 g per metric ton at about 380 m depth.     

Oxygen isotope study of the Kamioka Zn-Pb skarn deposits,Central Japan

The Kamioka Zn-Pb deposits consist of clinopyroxene-rich skarns replacing limestone lenses in Hida gneissose rocks. Paleozoic metabasite, early Jurassic Funatsu granitic rocks and late Cretaceous porphyritic dikes and stocks, are the known igneous rocks in the mining area. Oxygen isotopic compositions of eight separates each of clinopyroxene and quartz from the deposits (Table 1), and some minerals and rocks from the major rock units around the deposits (Table 2), have been analyzed. The estimated oxygen isotopic ratios of the fluids responsible for the formation of the deposits, range as low as -4 to +3 (SMOW). These values are considerably lower than expected from magmatic waters, and demonstrate that none of the igneous rocks in the area could be a direct source of the skarn- and ore-forming fluids. The data instead show that the deposits are formed by a huge convective circulation of solutions of meteoric water origin, promoted by a hidden batholithic intrusion at the time of the late Cretaceous.     

Sulfur Isotope Study on Hg and Sb Deposits in Japan

Abstract. Sulfur isotope ratios of cinnabar from Hg deposits and stibnite, jamesonite and berthierite from Sb deposits in Japan are examined in order to understand metallogeneses of Hg and Sb deposits in Japanese island arcs. The studied Hg and Sb deposits include the Hg deposit at Yamato‐suigin (Honshu) and the Sb deposit at Ichinokawa (Shikoku) in the Southwest Japan arc. In addition, Hg deposits including Itomuka and Ryushoden in central Hokkaido and Hg and Sb mineralizations in Northeast Japan arc are examined. The δ³⁴S values of cinnabar from the Hidaka‐Kitami district, central Hokkaido, including the Itomuka and Ryushoden deposits range widely, from ‐10 to +16 %o, the highest values encountered at the Samani deposit. The δ³⁴S values of cinnabar from other areas in Japan range from ‐12 to +5 %o, having δ³⁴S values higher than +2 %o from southwestern Hokkaido (Meiji deposit), Shikoku (Suii deposit) and Kyushu (Hasami and Yamagano deposits). On the other hand, the δ³⁴S values of stibnite from all areas in Japan range from ‐14 to +5 %o, having positive δ³⁴S values higher than +2 %o up to +5 %o from southwestern Hokkaido (Yakumo, Toyotomi and Teine deposits) and eastern‐central Honshu (Hachiman and Daikoku deposits). The variation in δ³⁴S values of Hg and Sb deposits may reflect the variation in δ³⁴S values of country rocks or variation in mixing ratio of sulfur extracted from the country rocks, sulfur derived from seawater sulfate, and sulfur derived from magmatic emanations. The relatively high δ³⁴S values of cinnabar and stibnite higher than +2 %o from southwestern Hokkaido, eastern‐central Honshu and Kyushu are probably caused by contribution of volcanic emanation from arc magmas having positive σδ³⁴S values, whereas the positive δ³⁴S values of cinnabar higher than +2 %o from Suii deposit in Shikoku may be attributed to structurally substituted sulfate in limestone country rocks and/or sulfur derived from seawater sulfate. However, the wide range of the δ³⁴S values of cinnabar from the Hidaka‐Kitami district, central Hokkaido, is difficult to explain at this moment. Other relatively low, negative δ³⁴S values of cinnabar and stibnite, berthierite from other areas in Japan may be attributed to 1) incorporation of isotopically light sedimentary sulfur or sulfur derived from ilmenite‐series silicic magma, or 2) less contribution of volcanic emanation from arc magmas having positive σδ³⁴S values.     

Intra-plate seismicity gap along the Median Tectonic Line and oblique plate convergence in southwest Japan

The intra-plate seismicity map for southwest Japan, based on fairly complete historical data for the past four hundred years, reveals an inverse correlation between the seismic activity along the island arc and the slip-rate along the Median Tectonic Line during the Late Quaternary. In the eastern part, the tectonic line is geologically inactive but regional historic seismicity has been high. The intra-plate seismic activity is probably related with the well developed mosaic-like conjugate system of strike-slip faults there. Conversely, the historic seismicity has been low in the western part, especially low in an area along the most geologically active segment of the Median Tectonic Line. Since no creep movement has been found there, energy greater than that of the Mino-Owari earthquake of 1891 (M = 8.0) seems to be stored in this seismicity gap. The difference in seismic released energy between the two regions for the last four hundred years would be balanced by the strain energy accumulated in the seismic gap. The fairly uniform strain release is conformable to the idea, proposed on the basis of the trend of maximum compression axes, that the Philippine Sea plate is dragging southwest Japan southwestward along the Nankai trough.     

Evaluation of various satellite sensors for waterline extraction in a coral reef environment: Majuro Atoll, Marshall Islands

The ability of five satellite sensor bands (IKONOS band 4, Terra ASTER bands 3 and 4, and Landsat ETM+ bands 4 and 5) was examined to extract the waterline at coral reef coasts (Majuro Atoll, Marshall Islands) using different wavelength regions (near infrared [NIR] and shortwave infrared [SWIR]) and different spatial resolutions (4, 15, and 30 m). After performing georeferencing and normalization of the images, density slicing was used to extract the waterline. Comparisons of extracted waterline positions with ground-level data for eight transects and global positioning system (GPS) tracks of the island shorelines showed that NIR bands were superior to SWIR bands because of the characteristics of the coral reef coasts, including a lack of foam and suspended sediments (which can affect the NIR wavelength region, if present) and the presence of remnant water on reef flats during low tide (which can affect the SWIR wavelength region). A linear relationship was found between the estimation errors of waterline positions and the spatial resolutions of the NIR sensors. Analysis on estimation errors and image costs showed that Terra ASTER band 3 was the most cost-effective sensor for extracting waterlines with reasonable accuracy. The results serve as general guidelines for using satellite-derived data to estimate intertidal topography and detect and monitor shorelines in coral reef environments.     

Three-dimensional P-wave velocity structure beneath the Indonesian region

We present the P-wave seismic tomography image of the mantle to a depth of 1200 km beneath the Indonesian region. The inversion method is applied to a dataset of 118,203 P-wave travel times of local and teleseismic events taken from ISC bulletins. Although the resolution is sufficient for detailed discussion in only a limited part of the study region, the results clarify the general tectonic framework in this region and indicate a possible remnant seismic slab in the lower mantle.

Structures beneath the Philippine Islands and the Molucca Sea region are well resolved and high-velocity zones corresponding to the slabs of the Molucca Sea and Philippine Sea plates are well delineated. Seismic zones beneath the Manila, Negros and Cotabato trenches are characterized by high-velocity anomalies, although shallow structures were not resolved. The Molucca Sea collision zone and volcanic zones of the Sangihe and Philippine arcs are dominated by low-velocity anomalies. The Philippine Sea slab subducts beneath the Philippine Islands at least to a depth of 200 km and may reach depths of 450 km. The southern end of the slab extends at least to about 6°N near southern Mindanao. In the south, the two opposing subducting slabs of the Molucca Sea plate are clearly defined by the two opposing high-velocity zones. The eastward dipping slab can be traced about 400 km beneath the Halmahera arc and may extend as far north as about 5°N. Unfortunately, resolution is not sufficient to reveal detailed structures at the boundary region between the Halmahera and Philippine Sea slabs. The westward dipping slab may subduct to the lower mantle although its extent at depth is not well resolved. This slab trends N-S from about 10°N in the Philippine Islands to northern Sulawesi. A NE-SW-trending high-velocity zone is found in the lower mantle beneath the Molucca Sea region. This high-velocity zone may represent a remnant of the former subduction zone which formed the Sulawesi arc during the Miocene.

The blocks along the Sunda and Banda arcs are less well resolved than those in the Philippine Islands and the Molucca Sea region. Nevertheless, overall structures can be inferred. The bowl-shaped distribution of the seismicity of the Banda arc is clearly defined by a horseshoe-shaped high-velocity zone. The tomographic image shows that the Indian oceanic slab subducts to a depth deeper than 300 km i.e., deeper than its seismicity, beneath Andaman Islands and Sumatra and may be discontinuous in northern Sumatra. Along southern Sumatra, Java and the islands to the east, the slab appears to be continuous and can be traced down to at least a depth of the deepest seismicity, where it appears to penetrate into the lower mantle.     

Pre-seismic crustal deformation caused by an underthrusting oceanic plate,in eastern Hokkaido,Japan

Eastern Hokkaido is being considered as one of the prime candidates for a future great earthquake in Japan. Geodetic work and seismological data suggest that the continental plate is being compressed and dragged down into the asthenosphere by the underthrusting Pacific plate.A quantitative examination of this idea was carried out by an application of the finite-element method to a static-field problem of crustal deformation. The effect of heterogeneous medium with Poisson's ratio differing from 0.25 is significant. An inverse problem was solved by applying a quadratic programming method to the observed vertical displacement.A model, which assumes a nearly uniform tangential displacement along the interface of the oceanic and continental lithospheres down to a depth of 100 km, explains a gravity-change as well as the general features of crustal deformation observed in eastern Hokkaido. The displacement rate is estimated to be 2.7 cm/year. It indicates that a rebound of the continental lithosphere now brings about a thrust-type earthquake whose reverse slip is comparable to that estimated for great shallow earthquakes along the Kuril trench. A possibility of aseismic rebound is suggested at the lower part of the interface.     

State of stress within a thin elastic wedge: a model of internal deformation of the continental plate at arc—arc junctions

Summary. Anomalous stress distributions near arc—arc junctions have been found in several parts of the world. As a simple two-dimensional model, a thin elastic wedge is introduced to model the continental side of the arc—arc junction, and the stress distribution within it is calculated to interpret those anomalies. Uniform normal and/or shear stress boundary conditions are given separately to each edge to approximate the effect of the two neighbouring converging plate boundaries; Stress components are obtained by numerical integration. Among the various results, the following features are predominant and seem to be actually taking place at arc—arc junctions.
(1) The anomaly in stress pattern within the wedge is basically due to the difference in stress conditions between the two boundaries.
(2) A large stress concentration takes place when symmetrical shear stresses are applied to the two edges.