Tectonic evolution along the northernmost border of the Philippine Sea plate since about 1 Ma

Recent structural, tephrochronologic and magnetostratigraphic studies conducted along the northernmost border of the Philippine Sea (PHS) plate enable us to reconstruct the precise tectonic evolution along the convergent boundary between the PHS plate and the Northeast Japan (NEJ) plate or the North American (NAM) plate since about 1 Ma. The authors of the present study split the tectonic evolution into five stages and present the characteristics of each stage. A plate tectonic interpretation is proposed, based upon the tectonic evolution, with special reference to the mode of convergent motion. In brief, our interpretations are as follows: the relative motion between the PHS and the NEJ plates was not recognized geologically within the area studied from about 1.0 to 0.9 Ma (Stage 1), suggesting either none or small influence from the coupling between the two plates during that period of time. Convergence between the PHS and the NEJ plates was possibly in N-S direction from 0.9 to 0.5 Ma (Stage II), and probably north-northwestward since 0.5 Ma (Stages III to V). The mode of the convergent motion was that of buoyant subduction in Stages II and III. The mode changed gradually from buoyant subduction during Stage IV to collision in Stage V (0.07 Ma to the present).     

Geophysical studies of crustal structure of the Ongul Islands and the Northern Mizuho Plateau, East Antarctica

Explosion seismic experiments, gravity measurements and aeromagnetic surveys were made in the northern Mizuho Plateau including the Ongul Islands, East Antarctica, from 1979 to 1982 by the Japanese Antarctic Research Expeditions. The objective of these field operations was to determine the crustal structure along the 300 km-long oversnow traverse route between Syowa and Mizuho Stations. Three big shots were fired; at sea near Syowa Station, in an ice hole near Mizuho Station and in an ice hole between both stations. Twenty-seven temporal seismic stations were set up along the route. Gravity measurements were carried out at 30 points along this route. Aeromagnetic surveys over the area were made four times.In the seismic experiments, clear refracted waves from the Conrad (estimated depth 30 km) and the Moho (estimated depth 40 km) discontinuities were recorded. No layer with a velocity of less than 6 km/s was found in the Ongul Islands nor beneath the ice sheet in the surveyed area. The P-wave velocity in the upper layer varies with depth from 6.0 km/s on the surface to 6.4 km/s at a depth of 13 km. Comparing the observed record section with synthetic seismograms, it was derived that the Conrad was not associated with a sharp velocity discontinuity, but a linear velocity increase of 0.55 km/s in a transition zone of 2.4 km thick. Velocities of P^* and P_n were determined as 6.95 km/s and 7.93 km/s assuming a flat layered structure.Bouguer gravity anomalies could not be calculated along the whole profile because of a lack of data on bedrock topography, so reduced gravity anomalies were calculated. These anomalies indicate no abrupt changes of the bedrock topography.     

Paleoposition of Southwest Japan at 16 Ma: Implication from paleomagnetism of the Miocene Ichishi Group

New paleomagnetic data from shallow-marine sediments of the Ichishi Group suggest a clockwise tectonic rotation of Southwest Japan in the Middle Miocene. Samples have been collected from mud or tuff layers at 17 sites. Stability of remanent magnetization has been examined by using alternating field and thermal demagnetization. The polarity sequence, composed of four normal and seven reversed polarity sites, is correlated to Polarity Epoch 16 (15.2–17.6 Ma), based on micropaleontological assignment of the upper Ichishi Group to Blow's Zone N8. The mean paleomagnetic direction of the 11 sites shows an anomalous declination toward the northeast. This result suggests that Southwest Japan was subjected to a clockwise rotation through 45° since 16 Ma. The clockwise rotation can be explained by the drift of Southwest Japan associated with the spreading of the Japan Sea during the Middle Miocene.     

Role of Stress-controlled Flow Pathways in HDR Geothermal Reservoirs

— We addressed effects of in situ stress on the formation of flow pathways in fractured rocks in geothermal reservoirs, especially for HDR projects. Here we focused on fractures which are critically-stressed, causing shear slip in a current stress field. The sliding is likely to break sealing in the fractures and, as a result, to increase their permeability. Such a mechanism is possibly significant under high-temperature conditions at geothermal fields because of temperature enhancement on chemical reactions for the sealing. We present a procedure to estimate the orientation of the critically-stressed fractures relative to axes of in situ principal stress with the aid of the Mohr diagram. The procedure allows us to evaluate intuitively how the orientation changes with factors such as magnitude of in situ principal stresses and pore pressure. We applied the procedure to estimate possible orientations of the critically-stressed fractures in major HDR test sites. Results show that overall alignments of microseismicity during hydraulic stimulation are within predicted ranges for possible orientations of the critically-stressed fractures. Furthermore, it was found that if the state of in situ stress is not favorable to cause sliding of natural fractures, it tends to lead a high wellhead pressure at hydraulic stimulation and a high recovery rate at circulation. On the other hand, if the state of in situ stress is favorable for sliding, it tends vice versa to lead a low wellhead pressure at hydraulic stimulation and a low recovery rate at circulation.     

Insights into operation of the subduction factory from the oxygen isotopic values of the southern Izu–Bonin–Mariana Arc

Abstract Oxygen is the most abundant element in the earth, and isotopic analysis of this element in island arc lavas potentially provides sensitive constraints on the proportion of oxygen recycled from subducted material, relative to that extracted from the mantle. Here we report on 225 new oxygen isotopic analyses of whole‐rock and glass samples, and clinopyroxene separates, from lavas collected from the southernmost 1500 km of the Izu–Bonin–Mariana (IBM) convergent margin. Whole‐rock samples clustered around a mean of 6.11 ± 0.47‰, whereas Mariana Trough glasses and mafic melts, calculated to be in equilibrium with mafic phenocrysts, clustered narrowly around a mean of 5.7‰. These data demonstrate that unequivocal identification of magmatic oxygen requires analysis of fresh glass or mafic minerals, and that the source of southern IBM Arc melts is entirely, or almost entirely, in equilibrium with normal mantle oxygen. If the elemental enrichments characteristic of the subduction component originate in subducted materials, these oxygen isotopic data are most consistent with the interaction of a small amount of sediment melt (<4%; mostly less than 1%) with mantle peridotite to yield the hybrid mantle that melts to form IBM Arc magmas.     

新疆北部富蕴县沙尔布拉克玻安岩的地球化学特征及构造意义

详细的地球化学研究表明，新疆北部富蕴县沙尔布拉克中泥盆统北塔山组的辉石安山岩具有中等的SiO2含量（52.13%-60.45%)；高MgO(3.91%-9.31%)和相容元素（包括Ni、Cr、Co等）含量，以及高Mg^#和Al2O3／TiO2比值；低TiO2含量（0.21%-0.35%)；同时具有U型稀土元素分布模式；并强烈亏损高场强元素。这些特征均与玻安岩的地球化学特征一致，应属于玻安岩，形成于前弧环境。沙尔布拉克玻安岩的存在表明，古亚洲洋在早－中泥盆世向南发生了洋壳俯冲作用。     

Multi-resolution clustering analysis and 3-D visualization of multitudinous synthetic earthquakes

We use modern and novel techniques to study the problems associated with detection and analysis of multitudinous seismic events, which form the background for isolated great earthquakes. This new approach involves multivariate analysis of low and large magnitude events recorded in space over a couple of centuries in time. We propose here the deployment of the clustering scheme both for extracting small local structures and large-scale trends in synthetic data obtained from four numerically simulated models with: uniform properties (U), a Parkfield-type asperity (A), fractal brittle properties (F), and multi-size-heterogeneity fault zone (M). The mutual nearest neighbor (mnn) clustering scheme allows for extraction of multi-resolutional seismic anomalies in both the spatio-temporal and multi-dimensional feature space. We demonstrate that the large earthquakes are correlated with a certain pathway of smaller events. Visualization of the anomalies by using a recently introduced visualization package Amira reveals clearly the spatio-temporal relationships between clusters of small, medium and large earthquakes, indicating significant stress relaxation across the entire fault region. We demonstrate that this mnn scheme can extract distinct clusters of the smallest events, which precede and follow a singularly large magnitude earthquake. These clusters form larger spatio-temporal structures comprising a series of large earthquakes. The link between the large and medium magnitude events is not so clearly understood. Short-ranged correlations are dominated by strong spatio-temporal anomalies, thus reflecting the global seismic properties of the entire fault zone.Electronic Supplementary Material Supplementary material is available for this article if you access the article at . A link in the frame on the left on that page takes you directly to the supplementary material.     

Garnet-ilmenite-perovskite transitions in the system Mg4Si4O12-Mg3Al2Si3O12 at high pressures and high temperatures: phase equilibria, calorimetry and implications for mantle structure

Phase relations in the system Mg₄Si₄O₁₂-Mg₃Al₂Si₃O₁₂ were examined at pressures of 19-27 GPa and relatively low temperatures of 800-1000 °C using a multianvil apparatus to clarify phase transitions of pyroxene-garnet assemblages in the mantle. Both of glass and crystalline starting materials were used for the experiments. At 1000 °C, garnet solid solution (s.s.) transforms to aluminous ilmenite s.s. at 20-26 GPa which is stable in the whole compositional range in the system. In Mg₄Si₄O₁₂-rich composition, ilmenite s.s. transforms to a single-phase aluminous perovskite s.s., while Mg₃Al₂Si₃O₁₂-rich ilmenite s.s. dissociates into perovskite s.s. and corundum s.s. These newly determined phase relations at 1000 °C supersede preliminary phase relations determined at about 900 °C in the previous study. The phase relations at 1000 °C are quite different from those reported previously at 1600 °C where garnet s.s. transforms directly to perovskite s.s. and ilmenite is stable only very close to Mg₄Si₄O₁₂. The stability field of Mg₃Al₂Si₃O₁₂ ilmenite was determined at 800-1000 °C and 25-27 GPa by reversed phase boundaries. In ilmenite s.s., the a-axis slightly increases but the c-axis and molar volume decrease substantially with increasing Al₂O₃ content. Enthalpies of ilmenite s.s. were measured by differential drop-solution calorimetry method using a high-temperature calorimeter. The excess enthalpy of mixing of ilmenite s.s. was almost zero within the errors. The measured enthalpies of garnet-ilmenite and ilmenite-perovskite transitions at 298 K were 105.2±10.4 and 168.6±8.2 kJ/mol, respectively, for Mg₄Si₄O₁₂, and 150.2±15.9 and 98.7±27.3 kJ/mol, respectively, for Mg₃Al₂Si₃O₁₂. Thermodynamic calculations using these data give rise to phase relations in the system Mg₄Si₄O₁₂-Mg₃Al₂Si₃O₁₂ at 1000 and 1600 °C that are generally consistent with those determined experimentally, and confirm that the single-phase field of ilmenite expands from Mg₄Si₄O₁₂ to Mg₃Al₂Si₃O₁₂ with decreasing temperature. The earlier mentioned phase relations in the simplified system as well as those in the Mg₂SiO₄-Fe₂SiO₄ system are applied to estimate mineral proportions in pyrolite as a function of depth along two different geotherms: one is a horizontally-averaged temperature distribution in a normal mantle, and the other being 600 °C lower than the former as a possible representative geotherm in subducting slabs. Based on the previously described estimated mineral proportions versus depth along the two geotherms, density and compressional and shear wave velocities are calculated as functions of depth, using available mineral physics data. Along a normal mantle geotherm, jumps of density and velocities at about 660 km corresponding to the post-spinel transition are followed by steep gradients due to the garnet-perovskite transition between 660 and 710 km. In contrast, along a low-temperature geotherm, the first steep gradients of density and velocities are due to the garnet-ilmenite transition between 610 and 690 km. This is followed by abrupt jumps at about 690 km for the post-spinel transition, and steep gradients between 700 and 740 km that correspond to the ilmenite-perovskite transition. In the latter profile along the low-temperature geotherm, density and velocity increases for garnet-ilmenite and ilmenite-perovskite transitions are similar in magnitude to those for the post-spinel transition. The likely presence of ilmenite in cooler regions of subducting slabs is suggested by the fact that the calculated velocity profiles along the low-temperature geotherm are compatible with recent seismic observations indicating three discontinuities or steep velocity gradients at around 600-750 km depth in the regions of subducting slabs.     

Initial rupture process of microearthquakes recorded by high sampling borehole seismographs at the Nojima fault, central Japan

We analyze high sampling waveforms of the initial part of P-wave recorded at the 1800-m-deep borehole seismographs at the Nojima fault from December 1999 to May 2000 to clarify the initial rupture process of microearthquakes. We select 12 events with high S/N, whose magnitudes range from −0.3 to 2.2 and hypocentral distances from 1 to 11 km. We adopt the two different source models by Sato and Hirasawa (1973) and by Sato and Kanamori (1999). The former (model by Sato and Hirasawa (SH model)) generates only a ramp-like onset of velocity pulse. The later (model by Sato and Kanamori (SK model)) is able to generate a weak initial phase that is controlled by a trigger factor and the length of pre-existing crack. We perform the waveform inversion to estimate the optimum source parameters of each model. Waveforms of 5 of the 12 events are clearly reproduced by both SH model and SK model with a large trigger factor and a small length of pre-existing crack. The others are explained by not SH model but only SK model with a small trigger factor and a large length of the pre-existing crack, indicating that the weak initial phase is a nucleation phase and reflects the source process. These seven events satisfy roughly a relation that a large event has a large length of the pre-existing crack; the final crack length is proportional to the length of the pre-existing crack.