H-adaptivity applied to liquefiable soil in nonlinear analysis of soil–pile interaction

This paper is concerned with application of the h-adaptive finite element method to dynamic analysis of a pile in liquefiable soil considering large deformation. In finite element analysis of pile behavior in liquefiable soil during an earthquake, especially considering large deformation of liquefied ground, error due to discretization in the zone near the pile becomes very large. Our purpose was to refine the approximation of the finite element method. The updated Lagrangian formulation and a cyclic elasto-plastic model based on the kinematic hardening rule were adopted to deal with the nonlinearity of the soil. The mixed finite element and finite difference methods together with the u-p formulation and Biot's two-phase mixture theory were used. To improve the accuracy and increase the efficiency of finite element analysis, an h-adaptive scheme that included a posteriori error estimation and h-version mesh refinement was applied to the analysis. The calculated results of effective stress were smoothed locally by the extrapolation method and smoothed stress was used to calculate the L₂ norm of the effective stress error in the last step of the calculation of each time increment. The mesh was refined by a fission procedure based on the indication of the error estimate As a numerical example, a soil–pile interaction system loaded cyclically was analyzed by our method.     

Fault system of the 2004 Mid Niigata Prefecture Earthquake and its aftershocks

Data recorded by a seismic network deployed the day after the 2004 Mid Niigata Prefecture Earthquake (M6.8) in central Japan are used to determine the major source faults responsible for the mainshock and major aftershocks. Using this high-resolution seismic data, three major source faults are identified: two parallel faults dipping steeply to the west located 5 km apart, and the other dipping eastward and oriented perpendicular to the west-dipping faults. The analysis also reveals that the lateral variation in seismic velocity observed at the surface extends to a depth of 15 km, encompassing the source area of the mainshock. This strong heterogeneity of the crust, related to the complex geological and tectonic evolution of the area, is considered to be responsible for the prominent aftershock activity following the 2004 Niigata event.     

Tidal triggering of earthquakes in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, Japan

We found a characteristic space–time pattern of the tidal triggering effect on earthquake occurrence in the subducting Philippine Sea plate beneath the locked zone of the plate interface in the Tokai region, central Japan, where a large interplate earthquake may be impending. We measured the correlation between the Earth tide and earthquake occurrence using microearthquakes that took place in the Philippine Sea plate for about two decades. For each event, we assigned the tidal phase angle at the origin time by theoretically calculating the tidal shear stress on the fault plane. Based on the distribution of the tidal phase angles, we statistically tested whether they concentrate near some particular angle or not by using Schuster's test. In this test, the result is evaluated by p-value, which represents the significance level to reject the null hypothesis that earthquakes occur randomly irrespective of the tidal phase angle. As a result of analysis, no correlation was found for the data set including all the earthquakes. However, we found a systematic pattern in the temporal variation of the tidal effect; the p-value significantly decreased preceding the occurrence of M ≥ 4.5 earthquakes, and it recovered a high level afterwards. We note that those M ≥ 4.5 earthquakes were considerably larger than the normal background seismicity in the study area. The frequency distribution of tidal phase angles in the pre-event period exhibited a peak at the phase angle where the tidal shear stress is at its maximum to accelerate the fault slip. This indicates that the observed small p-value is a physical consequence of the tidal effect. We also found a distinctive feature in the spatial distribution of p-values. The small p-values appeared just beneath the strongly coupled portion of the plate interface, as inferred from the seismicity rate change in the past few years.     

230Th/234U and 14C dating of a lowstand coral reef beneath the insular shelf off Irabu Island, Ryukyus, southwestern Japan

Abstract   High-resolution seismic reflection profiles delineated the distribution of mound-shaped reflections, which were interpreted as reefs, beneath the insular shelf western off Irabu Island, Ryukyus, southwestern Japan. A sediment core through one of the mounded structures was recovered from the sea floor at a depth of −118.2 m by offshore drilling and was dated by radiometric methods. The lithology and coral fauna of the core indicate that the mounded structure was composed of coral–algal boundstone suggesting a small-scaled coral reef. High-precision α-spectrometric ²³⁰Th/²³⁴U dating coupled with calibrated accelerator mass spectrometric ¹⁴C ages of corals obtained reliable ages of this reef ranging from 22.18 ± 0.63 to 30.47 ± 0.98 ka. This proves that such a submerged reef was formed during the lowstand stage of marine oxygen isotope stages 3–2. The existence of low-Mg calcite in the aragonitic coral skeleton of 22.18 ± 0.63 ka provides evidence that the reef had once been exposed by lowering of the relative sealevel to at least −126 m during the last glacial maximum in the study area. There is no room for doubt that a coral reef grew during the last glacial period on the shelf off Irabu Island of Ryukyus in the subtropical region of western Pacific.     

Fracture Toughness Evaluation Based on Tension-softening Model and its Application to Hydraulic Fracturing

This paper discusses the applicability of the tension-softening model in the determination of the fracture toughness of rocks, where the fracture toughness evaluated based on the tension-softening model is compared with the crack growth resistance deduced from laboratory-scale hydraulic fracturing tests. It is generally accepted that the fracture process is dominated by the growth of a fracture process zone for most types of rocks. In this study, the J-integral based technique is employed to determine the fracture toughness of Iidate granite on the basis of the tension-softening model, where compact tension specimens of different dimensions were tested in order to examine the specimen size effect on the measured fracture toughness. It was shown that the tension-softening relation deduced from the J-integral based technique allowed us to determine the specimen size independent fracture toughness K_c of Iidate granite. Laboratory-scale hydraulic fracturing tests were performed on cubic specimens (up to a 10 m sized specimen), where cyclic pressurization was conducted using a rubber-made straddle packer to observe the extent of the hydraulically induced crack. The experimental results of pressure and crack length were then used to construct the crack growth resistance curve based on the stress intensity factor K. The crack growth resistance obtained from the hydraulic fracturing tests was observed to initially increase and then level off, giving a constant K value for a long crack extension stage. The plateau K value in the crack growth resistance curve was found to be in reasonable agreement with the fracture toughness K_c deduced from the tension-softening relation. It was demonstrated that the tension-softening model provides a useful tool to determine the appropriate fracture toughness of rocks, which may be applicable for the analysis of the process of large-scale crack extension in rock masses.     

Analysis of a guided wave along a conducting structure in a borehole

A guided wave along a borehole is often observed in borehole radar measurements. These guided waves deform the antenna pattern and can cause artefacts in radar measurements. A water-filled borehole or a conducting logging cable can function as a waveguide for electromagnetic waves under some conditions. We describe the theoretical characteristics of such a guided wave in a borehole and compare them with our experiments. The measured signal discussed was obtained with a directional borehole radar. This radar uses a cylindrical conformal array antenna as receiver and is a model of a conducting structure in a borehole. The induced field around the borehole was compared with the theory. The most fundamental symmetrical and asymmetrical modes were TM₀₁ and HE₁₁, and they were identified in the measured signals using time–frequency distribution analysis and by observation of the azimuthal field distribution of the magnetic field.     

Three-dimensional Mapping of Magma Source and Transport Regions from Seismic Data: The Mantle Wedge beneath Northeastern Japan

—We investigate the distribution of partial melt in island arc using the seismic velocity structure of the mantle wedge beneath northeastern Japan. The comparison of the seismic tomography with laboratory velocity data on a partially-molten mantle rock yields estimates of melting zones in three dimensions. We employ experimental data on the degree of partial melt in hydrous peridotite to give constraints on the melt fraction and temperature. Melting and magma-rich zones derived from the velocity structure coincide with observed low Q zones. The results of the three-dimensional mapping indicate that the source of magma in island arc is diapir-like melting patches localized within the low velocity zones of the mantle wedge. Extensive volcanic activity along the volcanic front is due to the presence of vast magma-rich zones just beneath the Moho. Those melting zones in the uppermost mantle may, in turn, cause melting of lower crustal materials and produce felsic magma. Melt appears to stay at and beneath the Moho, where crystallization fractionation may proceed. Melt exists at greater depths in the back-arc region, which may correlate with across-arc variations of chemical compositions of the volcanic rocks observed in northeastern Japan. We suggest that magma migration in the ductile lower crust may cause low-frequency microearthquakes, and magma penetration into the brittle upper crust may produce mid-crustal S-wave reflectors.     

The Holocene development of Kushu Lake on Rebun Island in Hokkaido, Japan

Diatom assemblages and sulfur content in sediments were analyzed to clarify changes in the sedimentary environment of Kushu Lake, a coastal lake on Rebun Island in Hokkaido, Japan. Salinity variations were assessed by means of a diatom-based index of paleosalinity and the sedimentary sulfur content. This paper discusses the Holocene development of the lake, in relation to Holocene relative sea-level change. For paleoenvironmental interpretation of the lake development, the rationale of the threshold method (Anundsen et al., 1994) was applied.At ca. 8000 yr BP, a coastal embayment (paleo-Kushu Bay) resulted from marine ingression. The threshold elevation at the mouth of the paleo-Kushu Bay kept pace with the rising sea-level, resulting in its enclosure at the culmination of Holocene marine transgression (ca. 6500–5000 yr BP). From predicted relative sea-level at ca. 6000 yr BP for Rebun Island (Nakada et al., 1991), the threshold may have been at least above –3 to –5 m altitude. A freshwater lake environment with strongly anoxic bottom conditions may have occurred from ca. 5500 to 5100 yr BP. After an important episode of marine ingression, the lake was isolated completely from the open sea at ca. 4900 yr BP. The diatom record suggests that the maximum lacustrine extent occurred at ca. 4900–3100 yr BP. Thereafter, water depth decreased at the lake margins.In Kushu Lake, the threshold elevation, due to a build-up of a coastal barrier, prevents us from determining the amplitude of sea-level changes, even though the age of isolation contacts corresponds to periods of regression and climatic deterioration. In spite of isostatic subsidence, the effective protection provided by the well-developed barrier did not allow registration of any relative sea-level fluctuations since its isolation.