Parameter identification method for determination of elastic modulus of rock based on adjoint equation and blasting wave measurements

The purpose of this paper is to present a parameter identification method to determine the force of a blast and the elastic modulus of the ground using the measurements of a dynamic elastic wave, the adjoint equation method of optimal control theory, and the finite element method. Before the excavation of rocky ground, it is important to estimate the ground properties. In this paper, the elastic modulus is determined as the performance function is minimized using a technique based on the first‐order adjoint method. The performance function is a square sum of the discrepancies between the computed and the observed values of the velocities. After the determination of the magnitude of the blasting force, we can determine the elastic modulus of the rock. As the basic equation to calculate the velocities of dynamic elastic body, elastic equilibrium equations with linear viscosity are employed. The adjoint equation method has been utilized in order to calculate the gradient of the performance function with respect to the parameters. The gradient of the performance function is calculated using the first‐order adjoint equation. The weighted gradient method is applied for minimization. In order to solve the state equations in space and time, the finite element method and the Newmark $\frac{1}{4}$ method are used. In this paper, we tested the practical application of our proposed method for determination of the elastic modulus of rock at the Ikawa tunnel located in the Tokushima prefecture, Japan. Copyright © 2009 John Wiley & Sons, Ltd.     

Alteration and mass transfer inferred from the Hirabayashi GSJ drill penetrating the Nojima Fault, Japan

Abstract Mineralogical and geochemical studies on the fault rocks from the Nojima–Hirabayashi borehole, south-west Japan, are performed to clarify the alteration and mass transfer in the Nojima Fault Zone at shallow depths. A complete sequence from the hornblende–biotite granodiorite protolith to the fault core can be observed without serious disorganization by surface weathering. The parts deeper than 426.2 m are in the fault zone where rocks have suffered fault-related deformation and alteration. Characteristic alteration minerals in the fault zone are smectite, zeolites (laumontite, stilbite), and carbonate minerals (calcite and siderite). It is inferred that laumontite veins formed at temperatures higher than approximately 100°C during the fault activity. A reverse component in the movement of the Nojima Fault influences the distribution of zeolites. Zeolite is the main sealing mineral in relatively deep parts, whereas carbonate is the main sealing mineral at shallower depths. Several shear zones are recognized in the fault zone. Intense alteration is localized in the gouge zones. Rock chemistry changes in a different manner between different shear zones in the fault zone. The main shear zone (MSZ), which corresponds to the core of the Nojima Fault, shows increased concentration of most elements except Si, Al, Na, and K. However, a lower shear zone (LSZ-2), which is characterized by intense alteration rather than cataclastic deformation, shows a decreased concentration of most elements including Ti and Zr. A simple volume change analysis based on Ti and Zr immobility, commonly used to examine the changes in fault rock chemistry, cannot account fully for the different behaviors of Ti and Zr among the two gouge zones.     

Evaluation of sheet pile-ring countermeasure against liquefaction for oil tank site

The evaluation of a countermeasure against liquefaction which uses a sheet pilering for oil tank sites is presented. The simulation of earthquake responses observed at tank sites with and without sheet pile-ring is first performed to validate the three-dimensional finite element numerical model. Using the numerical model, liquefaction analysis is performed and the excess pore water pressure generated in the soil and the settlement of tank are investigated. The comparison of two- and three-dimensional models is also conducted to assess the applicability of two-dimensional analysis. The results show that the numerical model could simulate the observed earthquake responses of tank-ring-soil system, and that the excess pore water pressure and the settlement of the tank could be significantly reduced using a sheet pile-ring. The two-dimensional analysis proves to be capable of representing the main features of the dynamic response of the three-dimensional tank-ring-soil system.     

Helium Isotopes in Pacific Waters from Adjacent Region of Honshu, Japan

We have measured helium isotopic ratios of thirty-seven Pacific water samples from various depths collected in adjacent regions of Honshu, Japan. The ³He/⁴He ratios vary significantly from 0.989 R _atm to 1.208 R _atm where R _atm is the atmospheric ratio of 1.39 × 10⁻⁶. The mid-depth (750–1500 m) profile of ³He/⁴He ratios at ST-1 located Northwestern Pacific Ocean east of Japan (Off Joban; 37°00′ N, 142°40′ E) is significantly different from that at ST-2 of the Northern Philippine Sea south of Japan (Nankai Trough; 33°07′ N, 139°59′ E), suggesting that these waters were separated by a topographic barrier, the Izu-Ogasawara Ridge. Taking ³He/⁴He data of the Geosecs expeditions in the western North Pacific, an extensive plume of 15% excess ³He relative to the air may be traced at ST-1 over 12,000 kilometers to the northwest of the East Pacific Rise where the mantle helium may originate. The 20% excess found at ST-2 may be attributable to the additional source of the subduction-type mantle helium in the Okinawa Trough. A 15% excess of ³He has also been discovered at a depth of about 1000∼1500 m at ST-3 adjacent to Miyakejima Island (33°57′ N, 139°22′ E) and ST-4 of Sagami Bay (35°00′ N, 139°22′ E). It is confirmed that mid-depth all over the western North Pacific water is affected by the mantle helium with a high ³He/⁴He ratio. This revised version was published online in July 2006 with corrections to the Cover Date.     

Ion microprobe U–Pb dating of zircon with a 15 micrometer spatial resolution using NanoSIMS

We developed a ²³⁸U–²⁰⁶Pb and ²⁰⁷Pb–²⁰⁶Pb zircon dating method using a Cameca NanoSIMS NS50 ion microprobe. A 7-to 9-nA O⁻ primary beam was used to sputter a 15-μm crater, and secondary positive ions were extracted for mass analysis using the Mattauch–Herzog geometry. The multicollector system was modified to detect ⁹⁰Zr⁺, ²⁰⁴Pb⁺, ²⁰⁶Pb⁺, ²³⁸U¹⁶O⁺, and ²³⁸U¹⁶O₂⁺ ions simultaneously. A mass resolution of about 4000 at 10% peak height and with a flat peak top was attained, and the sensitivity of Pb was about 4 cps·nA^− 1·ppm^− 1. A multicrystal zircon standard (QGNG) from South Australia with a U–Pb age of 1842 Ma was used as a reference for Pb⁺/UO⁺–UO₂⁺/UO⁺ calibration, and on the basis of the positive correlation between these ratios, we determined the sample ²⁰⁶Pb/²³⁸U ratios. ²⁰⁷Pb/²⁰⁶Pb ratios were measured by magnetic scanning in single-collector mode. The standard zircons 91500, from Canada, and SL13, from Sri Lanka, were analyzed against QGNG. Observed ²³⁸U–²⁰⁶Pb and ²⁰⁷Pb–²⁰⁶Pb ages agreed well with published ages within experimental error. Then, 16 zircon grains in a metamorphic rock from Nagasaki, Japan, were analyzed. Observed ages were compatible with SHRIMP ages, suggesting that the NanoSIMS with a 15-μm probe diameter is suitable for ion microprobe U–Pb zircon dating.     

Geographical distribution of helium isotope ratios in northeastern Japan

In order to study the precise geographical distribution of helium isotope ratios in northeastern Japan and compare it with geophysical data, we collected 43 gas and water samples from hot and mineral springs in the region where the ratio had never been reported, and measured the ³He/⁴He and ⁴He/²⁰Ne ratios of these samples. It was found that the ³He/⁴He ratios show clear contrasts between the forearc and the back-arc regions in the Tohoku district in northeastern Japan. In the forearc region, the ratios are smaller than 1 R_A (1 R_A = 1.4 × 10⁻⁶; R_A means the ³He/⁴He ratio of the atmosphere). On the other hand, those along the volcanic front and in the back-arc region are apparently higher. Moreover, we found a variation in the ³He/⁴He ratios along the volcanic front. In Miyagi Prefecture (38–39°N), the ratios range from 2 to 5 R_A. On the other hand, the ratios are less than 1 R_A in and around the southern border between Iwate and Akita Prefectures (39–39.5°N). Comparing the distribution of helium isotope ratios to results of recent geophysical studies, we found that the features in geographical distribution of helium isotope ratios are similar to those of seismic low-velocity zone distributions and high Qp⁻¹ distributions in the uppermost mantle. These observations strongly suggest that the helium isotope ratios reflect the distribution of melts in the uppermost mantle and are a useful tool for investigating the origin, behavior, and distribution of deep fluids and melts.