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.     

Collision structure in the upper crust beneath the southwestern foot of the Hidaka Mountains, Hokkaido, Japan as derived from explosion seismic observations

The P-wave velocity structure of the upper crust beneath a profile ranging from Niikappu to Samani in the southwestern foot of the Hidaka Mountains, Hokkaido, Japan was obtained through analysis of refraction and wide-angle reflection data. The mountains are characterized by high seismicity and a large gravity anomaly. The present profile crosses the source region of the 1982 Urakawa-oki earthquake (M_s 6.8). The length of the profile is 66 km striking northwest and southeast. Along the profile, 64 vertical geophones were set up and 5 shot points were chosen. For each shot, a 400–600 kg charge of dynamite was detonated. The studied area is composed of four major geological belts: Neogene sedimentary rocks, the Kamuikotan belt, the Yezo Group, and the Hidaka belt. The measurement line crosses these geological trend at an oblique angle. The structure obtained is characterized by remarkable velocity variations in the lateral direction and reflects the surface geological characteristics. A thin, high-velocity layer (HVL) was found between low-velocity materials in the central part of the profile, beneath the Kamuikotan Metamorphic Belt, at a depth ranging from 0.5 to 6 km, overthrusting toward the west on the low-velocity materials consisting of Neogene sedimentary rocks, and forming gentle folds. Outlines of the velocity structure of the Hidaka Mountains yielded by other studies have shown a large-scale overthrust structure associated with the collision of the Outer Kurile and the Outer Northern Honshu Arcs. The shallow velocity structure inferred by the present study showed a similar (although small scale) overthrust structure. The obtained structure shows that the composite tectonic force, comprising westward movement of the Outer Kurile Arc and northward movement of the Outer Northern Honshu Arc, plays an important role in the evolution of the tectonic features of the crust and upper mantle in a wide depth range beneath the Hidaka Mountains.     

Phase F,a new hydrous magnesium silicate synthesized at 1000 °C and 17 GPa: Crystal structure and estimated bulk modulus

The structure of phase F, a non-stoichiometric hydrous silicate synthesized in a uniaxial, split-sphere, multi-anvil apparatus at conditions of 17 GPa and 1000° C, has been solved and refined in space group P6₃cm, using synchrotron X-ray data for a single crystal of a size 18 × 24 × 30 μm. The composition and unit cell for phase F are Mg_3.35Si_5.51H_7.26O₁₈, a=5.073(3) Å, c=14.013(9) Å, 7= 312.3(5) ų. The structure contains layers with many similarities to superhydrous phase B. The layers of oxygen atoms are stacked in the ABCBAC-type double cubic closest packing arrangement. The bulk modulus of phase F was estimated from the structural and compositional relationship to superhydrous phase B and periclase.     

Crustal movement along a collision boundary of plates (case of eastern Taiwan)

Recent crustal movements in eastern Taiwan were investigated based on geodetic and seismological data as a case study of the tectonics along a colliding plate boundary. Secular horizontal strain along the Longitudinal Valley, a probable arc-continent collision boundary between the Philippine Sea and Eurasion plates, is compressive with a large rate of 2 microstrains/year, and the maximum compression axis trends in NW-SE direction. This means that nearly two-thirds of the convergence rate between the two plates is consumed in the narrow zone along the valley and confirms the view that the valley is the plate boundary between the two. The amount of uplift in the Central Range, deduced from trigonometric leveling, reaches more than one meter during the past 60 years. Such active orogeny is caused by the high strain rate in this area. Horizontal displacement of the Coastal Range, a leading edge of the Philippine Sea plate, has a left-lateral component, suggesting that the boundary is a region of oblique collision and the range moves toward the north, producing subduction of the Philippine Sea plate off the Coast of northeastern Taiwan.An extensional strain pattern is found northeastern Taiwan near Hualien, in northern Taiwan near Ilan, and in southeastern Taiwan south Taitung. In northeastern Taiwan, the direction of the maximum extension axis is nearly parallel to the isodipth contours of deep-focus earthquakes that occurred under the southern part of the Ryukyu Trench. From a model calculation using a finite-element method, this pattern of strain is interpreted as an edge effect of the collision-subduction junction in northeastern Taiwan. The direction of the maximum extension in northern Taiwan forms high angles with the direction of the tectonic line near Ilan stretching from the Okinawa Trough. This fact is harmonious with the view that the Okinawa Trough is spreading. Extensional strain observed in southeastern Taiwan may be related to the eastward subduction since the area is situated on a transition zone, from collision to subduction.     

Restoring forces on vertical circular cylinders forced by earthquakes

An eigenfunction solution is presented for the dynamic response of vertical circular cylinders to earthquake excitation in a compressible fluid of finite depth. This single eigenseries expansion eliminates the need for a double summation over both the eigenfunctions and the trial functions as required by Rayleigh-Ritz methods. Revised definitions for the added mass and hydrodynamic radiation damping coefficients per unit length are derived from the hydrodynamic fluid pressures. Based on comparisons between these newly defined coefficients, the compressibility of the fluid is found to be relatively more important at dimensionless frequencies greater than unity (ω > 1.0) when analysing both rigid and flexible cylinders having relatively large diameter to water depth ratios,r₀/h > 0.25 (squatty type). This conclusion regarding the relative importance of the fluid compressibility is derived from a comparison between the relative magnitudes and the vertical distributions over depth of both the added mass and radiation damping coefficients per unit length for both rigid and flexible squatty cylinders. From additional comparisons with Rayleigh-Ritz solutions that require trial functions, the results for totally immersed flexible slender cylinders (r₀/h< 0.10) are shown to be equivalent; but the results for totally immersed flexible squatty cylinders (r₀/h > 0.25) are not. The reason for this difference appears to be in the truncation of the trial function series in the Rayleigh-Ritz methods, which excludes the higher mode shapes, and in the definitions of the added mass coefficients. Comparisons with laboratory data for both rigid and flexible cylinders confirm the accuracy of the solutions obtained by the eigenseries in the limited frequency interval above the highest frequency for surface gravity waves (f > 1.0 Hz) and below the first dimensionless cut-off frequency for acoustic waves (ω< 1.0).     

Computational model of the structural and elastic properties of the ilmenite and perovskite phases of MgSiO3

The structural and elastic properties of the ilmenite and perovskite phases of MgSiO₃ are investigated with a computational model based on energy minimization. The potential energies of these two crystals are approximated by the sum of Coulomb, van der Waals, and repulsion terms between atoms. Required energy parameters are derived by fitting the parameters to the observed crystal structures of these two phases as well as to the measured elastic constants of the ilmenite phase. The resulting potential model is applied to predicting the elastic constants of the perovskite phase. The calculated bulk modulus of the perovskite phase compares favorably with the data obtained from volume-compression experiments as well as the values estimated from empirical elasticity systematics of perovskite type compounds. The predicted shear modulus of the perovskite phase is also in reasonable agreement with the values proposed from similar empirical elasticity systematics. Subsequently, the model is used to simulate the high pressure behaviors of the crystal structures and elastic constants of these two phases.     

Vortex structures around a hemispheric hump

The vortex structure formed in front of and behind a hemispheric surface hump has been studied in a water recirculating flume. In the recirculation zone behind the hemisphere, arched vortex tubes were formed. Up to a certain critical Reynolds number, each tube was shed individually from the recirculation zone. However, when the Reynolds number exceeded a critical value, several of the tubes coalesced in the recirculation zone before shedding. In front of the hemisphere, some parabolic vortex tubes formed. Their number increased with Reynolds number below its critical value and decreased above this value. The dimensionless height of the parabolic vortex tube decreased with Reynolds number.