Numerical modelling of the three-dimensional stress field in southwestern Japan

It has been inferred from the focal mechanism of earthquakes and their hypocenter distribution (Shiono, 1977) that the stress field in southwestern Japan indicates complicated features; a NW-SE compression at shallow depths along the Nankai trough, an E-W or ESE-WNW compression in the inland crust, an extension parallel to the leading edge of the Philippine Sea plate at subcrustal depths in the region from the southern Chubu to northwestern Shikoku, and a down-dip tension beneath the Kyushu island.In order to investigate possible sources of these complex features of the stress state, a three-dimensional finite element method is employed to model the configuration of the subducting Philippine Sea plate, taking into consideration the following three possible types of forces:

1. (1) A negative buoyancy due to the density contrast between the subducting plate and the surrounding mantle.

2. (2) A northwestward compressive force generated by the movement of the Philippine Sea plate.

3. (3) A westward compressive force due to the movement of the Pacific plate.

For various combinations of different magnitudes of these forces, and of different elastic moduli, the stresses at a number of selected sites are calculated, and their directions are compared with those inferred from the focal mechanism and other geophysical information.It is found that the observed extensional stresses parallel to the leading edge of the subducting Philippine Sea plate may be caused mainly by the negative buoyancy. The northwestward compressive force may not play an important role in generating the complex stress field in southwestern Japan. The observed E-W compressional stress field prevailing in the inland region appears to result mainly from the westward-moving Pacific plate. The present results suggest that if a thin low-velocity transitional layer exists just above the subducting Philippine Sea plate, it could give appreciable effects on the mechanism of low-angled thrust faulting off the Kii peninsula and the Shikoku island.The magnitude of the shear stress in the continental crust and in the subducting plate is estimated to be of the order of several hundred bars, although the calculated shear stresses are considerably affected by the configuration of the subducting plate and also by the applied forces.It is interesting that the stress distribution appears to have some relations to seismicity of subcrustal earthquakes, and to the rupture process of large thrust earthquakes along the plate boundary.     

Rare earth elements in garnets and clinopyroxenes from garnet lherzolite nodules in kimberlites

Six pairs of coexisting garnets and clinopyroxenes were separated from the sheared and granular garnet lherzolite nodules in kimberlites and analyzed for rare earth elements (REE). The sheared and granular nodules can be distinguished in terms of REE pattern of both clinopyroxene and garnet. However, there are no significant differences in REE partitioning between clinopyroxene and garnet, indicating that the partitioning may be insensitive toP, T and composition. REE partition coefficients between garnet and liquid were estimated by using clinopyroxene-liquid partition coefficients found in the literature and clinopyroxene-garnet partitioning reported here. The estimated values agree with those reported by Philpotts et al. (1972). The estimated whole-rock REE pattern for the sheared nodules is similar to a chondritic pattern suggesting that the sheared nodules appear to be close to the primary mantle material. The REE data suggest that the granular nodules were originally garnet-free assemblages equilibrated with kimberlitic or nepheline-melilite basalt-like liquid, and later recrystallized as a garnet lherzolite assemblage.     

Rare earth element concentrations in a suite of basanitoids and alkali olivine basalts from Grenada,Lesser Antilles

A suite of basanitoids and alkali olivine basalts from Grenada, Lesser Antilles were analyzed for rare earth elements. The REE concentrations of these rocks are characterized by a small variation in the heavy REE (7 to 9 times chondrite) and a large variation in the light REE (17 to 93 times chondrite). Among the possible mechanisms to account for the REE variations, fractional crystallization processes at low and high pressures, and partial melting processes (both batch melting and fractional melting) were examined, using the partition relationships of REE among silicate minerals and melts. It is suggested that the observed REE variations are best explained by variable degrees of batch partial melting, in which garnet is present as one of the solid phases through 2 to 17% melting of a garnet lherzolite parent rock.     

Overview of the MITI Nankai Trough Wells: A Milestone in the Evaluation of Methane Hydrate Resources

Abstract. Bottom-simulating reflectors suggestive of the presence of methane hydrates are widely distributed below the ocean floor around Japan. In late 1999, drilling of the MITI Nankai Trough wells was conducted to explore this potential methane hydrate resource and a Tertiary conventional structure. The wells are located in the Northwest Pacific Ocean off Central Japan at a water depth of 945 m. A total of six wells were drilled, including the main well, two pilot wells, and three post survey wells at intervals of 10–100 m. All wells except the first confirmed the occurrence of hydrates based on logging-while-drilling, wire-line logging and/or coring using a pressure and temperature coring system in addition to conventional methods. Based on the various well profiles, four methane hydrate-bearing sand-rich intervals in turbidite fan deposits were recognized. Methane hydrates fill the pore spaces in these deposits, reaching saturation of up to 80 % in some layers. The methane hydrate-bearing turbiditic sand layers are less than 1 m thick, with a total thickness of 12–14 m. The bottom depth of high hydrate concentration correlates well with the depth of the bottom-simulating reflector. Based on these exploration results, the Japanese government inaugurated a 16-year methane hydrate exploitation program in 2001.     

Sensitivity of the central Asian climate to uplift of the Tibetan Plateau in the coupled climate model (MRI-CGCM1)

Abstract   The relationship between the altitude of the Tibetan Plateau and climate change in central Asia was investigated through a numeric experiment using the Meteorological Research Institute (MRI) coupled atmosphere–ocean general circulation model I (MRI-CGCM1). The results suggest that summer precipitation in central Asia decreased significantly as the Tibetan Plateau rose in height. Spring precipitation, however, increased during initial growth stages when the plateau height was up to 40% of its present-day height, and then decreased with further plateau growth. During the Tibetan Plateau uplift, the difference between precipitation and evaporation was minimal during spring. When the plateau attained a height exceeding 60% of its present height, relatively low precipitation but high evaporation in spring led to a lower amount of ground moisture. In the case of the high plateau, sensible heat flux during summer and fall largely exceeded latent heat flux. Change was particularly significant for cases when the plateau reached 40–60% of its present-day height. The duration of the predominant sensible heat flux became longer with the uplift of the Tibetan Plateau. The period in which latent heat exceeded sensible heat seems to have been restricted to winter and early spring. The numeric experiments suggest that a significant drying of central Asia corresponded to the period in which the Tibetan Plateau exceeded approximately half its present-day height.     

Sulfur Isotope Geochemistry of the Supergiant Xikuangshan Sb Deposit, Central Hunan, China: Constraints on Sources of Ore Constituents

Abstract. The supergiant Xikuangshan Sb deposit is located in the Middle to Upper Devonian limestone of central Hunan, China. Primary ores are composed of early-stage stibnite and calcite with rare pyrite, early main-stage stibnite and quartz, and late main-stage stibnite and calcite. New sulfur isotope data reveal the clustering of δ³⁴S values (+5 ∼ +8 %) for both early and late main-stage stibnite; a single early-stage stibnite exhibits δ³⁴S value (+7.5 %) identical to its main ore-stage counterparts and the coexisting calcite has almost unmodified carbon isotope composition (-4.4 %). The data suggest a probable common source of sulfur for stibnite that was deposited at different paragenetic stages. A much wider variation in δ³⁴S values for early main-stage stibnite (+3.5 to +16.3 %, av. +7.5 %) compared to that for late main-stage stibnite (+5.3 to +8.1 %, av. +6.2 %) can be interpreted to be due to local interaction of earlier ore fluid with Devonian host rocks. The previous studies show that the Precambrian basement contains elevated Sb concentrations, and two distinctive sulfur reservoirs with δ³⁴S_pyrite values at ca. +11 ∼ +24 % and -7.0 ∼-11 %. The homogenizing effect for sulfur hydrothermally leached from the two reservoirs might have provided ore constituents for the Xikuangshan fluids.     

The effective failure modes and stability of slopes in rock mass with two discontinuity sets

Summary The stability of rock slopes in discontinuous rock mass associated with the construction of power plants, highways and open pits is always of paramount importance during the lifetime of these structures. The likely forms of instabilities observed in the excavation of rock slopes and some mathematical methods for the stability analyses are well documented in literature. Since most of the mathematical approaches used are based on the limiting-equilibrium concept, there seems a need to check the validity of these approaches under some controlled conditions. In this paper, the authors describe methods for the stability of a blocky column and discontinuous rock slopes derived on the basis of dynamic equilibrium equations and compare the results calculated according to the developed method with those of experiments on model blocky columns and model slopes in the laboratory. Test results confirm that the limiting equilibrium approache is valid and an effective way of dealing with the stability problems in discontinuous rock mass as long as the likely forms of instability are properly treated in these approaches.     

Wind fluctuations near a cold vortex-tropopause funnel system observed by the MU radar

Vertical and temporal variations of three-dimensional wind velocity associated with an upper-tropospheric cold vortex-tropopause funnel system were observed by an MST radar in Japan (the MU radar). Marked changes of vertical velocity and horizontal wind direction were found between the inside and outside of the cold vortex. The vertical velocity activity outside the vortex was asymmetric; it was most active in a sector before the vortex. Unsaturated internal gravity waves in their generation stage contribute predominantly to the vertical velocity activity, suggesting that tropospheric occluded cyclones may be a possible source of middle-atmospheric gravity waves through the geostrophic adjustment process.     

Internal gravity wave selection in the upper troposphere and lower stratosphere observed by the MU radar: Preliminary results

Marked wavelike variations of the lower stratospheric wind observed on 7–10 May, 1985 by an MST radar in Japan (by the MU radar) are analyzed assuming that they are induced by monochromatic internal inertio-gravity waves. These variations are mainly composed of two modes (periods: 22 and 24 hours), both of which have zonal phase velocities (C _X ) slower than the mean westerly wind (). A statistical analysis of the zonal phase velocity shows thatC _X above andC _X below the tropopause jet stream, which is considered to be a vivid proof of wave selection due to the tropospheric mean flow and upward wave emission from the tropopause jet. A comparison between the MU radar results and routine meteorological observations leads to the conclusion that the marked waves appear when the jet stream takes a maximum wind speed.