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.     

Izu-Bonin arc subduction under the Honshu island, Japan: Evidence from geological and seismological aspect

The Izu-Bonin intra-oceanic arc with 20–35 km thick continental crust is being subducted under the Honshu, presumably since 17 Ma. Tomographic image clearly demonstrates that the whole Izu-Bonin arc is subducting under the Honshu arc. Geologic cross section and the thickness of continental crust do not support the accretion of thick crust in spite of the continued subduction over 17 Ma.     

Quasi-static slip on the plate boundary associated with the 2003 M8.0 Tokachi-oki and 2004 M7.1 off-Kushiro earthquakes, Japan

We analyzed small repeating earthquakes recorded over a 13-year period and GPS data recorded over an 8-month period to estimate interplate quasi-static slip associated with the 2003 Tokachi-oki earthquake (M8.0) and the 2004 off-Kushiro earthquake (M7.1). The repeating-earthquake analysis revealed that the slip rate near the source region of the Tokachi-oki earthquake was relatively low (< 5 cm/year) prior to the earthquake; however, in the last 3 years leading up to the event, a minor acceleration in slip occurred upon the deeper extension of the coseismic slip area of the earthquake. Repeating-earthquake and GPS data indicate that large amounts of afterslip occurred around the rupture area following the earthquake; the afterslip mainly propagated to the east of the coseismic slip area. We also infer that the occurrence of the 2004 off-Kushiro earthquake, located about 100 km northeast of the epicenter of the Tokachi-oki earthquake, was advanced by the afterslip associated with the Tokachi-oki earthquake.     

The dynamics of big mantle wedge, magma factory, and metamorphic–metasomatic factory in subduction zones

The East Asian continental margin is underlain by stagnant slabs resulting from subduction of the Pacific plate from the east and the Philippine Sea plate from the south. We classify the upper mantle in this region into three major domains: (a) metasomatic–metamorphic factory (MMF), subduction zone magma factory (SZMF), and the ‘big mantle wedge’ (BMW). Whereas the convection pattern is anticlockwise in the MMF domain, it is predominantly clockwise in the SZMF and BMW, along a cross section from the south. Here we define the MMF as a small wedge corner which is driven by the subducting Pacific plate and dominated by H₂O-rich fluids derived by dehydration reactions, and enriched in large ion lithophile elements (LILE) which cause the metasomatism. The SZMF is a zone intermediate between MMF and BMW domains and constitutes the main region of continental crust production by partial melting through wedge counter-corner flow. Large hydrous plume generated at about 200 km depth causes extensive reduction in viscosity and the smaller scale hydrous plumes between 60 km and 200 km also bring about an overall reduction in the viscosity of SZMF. More fertile and high temperature peridotites are supplied from the entrance to this domain. The domain extends obliquely to the volcanic front and then swings back to the deep mantle together with the subducting slab. The BMW occupies the major portion of upper mantle in the western Pacific and convects largely with a clockwise sense removing the eastern trench oceanward. Sporadic formation of hydrous plume at the depth of around 410 km and the curtain flow adjacent to the trench cause back arc spreading. We envisage that the heat source in BMW could be the accumulated TTG (tonalite–trondhjemite–granodiorite) crust on the bottom of the mantle transition zone. The ongoing process of transportation of granitic crust into the mantle transition zone is evident from the deep subduction of five intra-oceanic arcs on the subducting Philippine Sea plate from the south, in addition to the sediment trapped subduction by the Pacific plate and Philippine Sea plate. The dynamics of MMF, SZMF and BMW domains are controlled by the angle of subduction; a wide zone of MMF in SW Japan is caused by shallow angle subduction of the Philippine Sea plate and the markedly small MMF domain in the Mariana trench is due to the high angle subduction of Pacific plate. The domains in NE Japan and Kyushu region are intermediate between these two. During the Tertiary, a series of marginal basins were formed because of the nearly 2000 km northward shift of the subduction zone along the southern margin of Tethyan Asia, which may be related to the collision of India with Asia and the indentation. The volume of upper mantle under Asia was reduced extensively on the southern margin with a resultant oceanward trench retreat along the eastern margin of Asia, leading to the formation of a series of marginal basins. The western Pacific domain in general is characterized by double-sided subduction; from the east by the oldest Pacific plate and from the south by the oldest Indo-Australian plate. The old plates are hence hydrated extensively even in their central domains and therefore of low temperature. The cracks have allowed the transport of water into the deeper portions of the slab and these domains supply hydrous fluids even to the bottom of the upper mantle. Thus, a fluid dominated upper mantle in the western Pacific drives a number of microplates and promote the plate boundary processes.     

A model for the hydrogen coma of a comet

A model for the hydrogen coma of a comet on the basis of the Monte Carlo method is presented. In this model isotropic ejections of H atoms produced by photodissociation of H₂O and OH, thermalization of the H atoms due to collisions with ambient H₂O molecules, and the solar radiation pressure have been taken into account. A production spectrum of H atoms from OH is evaluated by using the predissociation rates and the level populations of OH, confirming that the spectrum has a sharp peak around 8.0 km sec^?1 with the standard deviation of 0.1 km sec^?1. Including the above effects, velocity distribution functions of the H atoms at various positions in the coma for the first time, as well as their density and outflow velocity profiles, have been calculated. It is pointed out that the collisional thermalization process in the inner coma is an important factor at small heliocentric distances in determining the density profiles and the velocity distributions. It is shown that thermalization leads to an increase in the H density in the inner coma larger than those expected from other models such as the vectorial model, in which collision is not taken into account. Lyman α isophotes and its line profiles in the optically thin region are computed by using the velocity distribution function.     

Rare gas isotopic compositions in natural gases of Japan

Isotopic and elemental compositions of rare gases in various types of gas samples collected in the Japanese Islands were investigated. Excess³He was found in most samples. Many samples showed a regionally uniform high³He/⁴He ratio of about 7 times the atmospheric ratio. The He concentrations varied from 0.6 to 1800 ppm, and they were low in CO₂-rich gases and high in N₂-rich gases. Ne isotopic deviations from the atmospheric Ne were detected in most volcanic gases. The deviations and the elemental abundance patterns in volcanic gases can be explained by a mixing between two components, one is mass fractionated rare gases and the other is isotopically atmospheric and is enriched in heavy rare gas elements. Ar was a mixture of mass fractionated Ar, atmospheric Ar and radiogenic Ar, and the contribution of radiogenic⁴⁰Ar was small in all samples. Except for He, elemental abundance patterns were progressively enriched in the heavier rare gases relative to the atmosphere. Several samples were highly enriched in Kr and Xe relative to the abundance pattern of dissolution equilibrium of atmospheric rare gases in water. The component which is highly enriched in heavy rare gases may be released from sedimentary materials in the crust.     

Partition of magnesium and iron between Olivine and spinel,and between pyroxene and spinel

Partition of Fe²⁺ and Mg between coexisting (Mg, Fe)₂SiO₄ spinel and (Mg, Fe)SiO₃ pyroxene was investigated at pressures 80 and 90 kbar and at temperatures 840 and 1050° C, using tetrahedral-anvil type of high pressure apparatus. Olivine-spinel solid solution equilibria in the system Mg₂SiO₄-Fe₂SiO₄ were discussed in the light of the partition reaction. Partition of Fe²⁺ and Mg in both olivine-spinel and pyroxene-spinel systems can not be regarded as that between ideal solid solutions. By applying the simple solution model for the partition of Fe²⁺ and Mg, sign of the heat of mixing was estimated to be positive for all olivine, spinel and pyroxene. Relative concentration of Fe²⁺ in spinel in the pyroxene-spinel system is likely to cause some change in the chemical composition of modified spinel () or spinel () in the transition zone of the mantle. A considerable change is also expected in the transition pressure of to ( + ) and ( + ) to .Presented at the symposium Recent Advances in the Studies of Rocks and Minerals at High Pressures and Temperatures held in Montreal, 1972. Jointly sponsored by the International Mineralogical Association and the Commission on Experimental Petrology.     

Step-like temperature fluctuations associated with inverted ramps in a stable surface layer

Ensemble averages of temperature before and after step-like temperature fluctuations reveal the presence of inverted ramps in a stable surface layer. Normalized frequency of upward steps increases with increasing stability, whereas normalized magnitude of the temperature step decreases with stability and becomes constant at about R _i = 1. These results suggest that the significance of temperature steps increases as stability increases. In moderate stability, the temperature pattern shows a gradual decrease after an upward step, which can be called a time-inverted ramp. Descending air and large downward heat flux are observed in a time-inverted ramp, suggesting a contribution from an ordered motion in wind. On the other hand, the temperature steps are related to gravity waves in strong stability.     

Imaging the heterogeneous source area of the 2003 M6.4 northern Miyagi earthquake, NE Japan, by double-difference tomography

A shallow M6.4 inland earthquake occurred on 26 July 2003 in the northern part of Miyagi Prefecture, northeastern Japan. This earthquake was a typical inland thrust earthquake, a type that is common in NE Japan. We obtained a detailed seismic velocity structure in the focal area of this earthquake by the double-difference tomography method. Arrival-time data came from temporary seismic stations deployed above the mainshock fault plane. Both the P-wave and S-wave velocities in the hanging wall were lower than those in the footwall. Aftershocks were aligned along a zone where the seismic velocity changes rapidly. This is consistent with the interpretation that the 2003 northern Miyagi earthquake occurred along a fault that acted as a normal fault in the Miocene and has been reactivated as a reverse fault under the present compressional stress regime. The large slip area by the main shock rupture (asperity) corresponds to an area with relatively high P- and S-wave velocities. A zone with low V_p/V_s was detected along the aftershock area. One of the possible causes of this low-V_p/V_s zone is the existence of high-aspect-ratio pores that contain water. Hypocenters of the main shock, largest foreshock, and largest aftershock are also located within the low-V_p/V_s zone.