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).     

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.     

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.     

Optimal placement of sensors on buildings subjected to intermediate‐storey excitation according to QN control method

An algorithm to calculate direct velocity feedback gain with limited number of sensors is developed in a simple way such that a certain performance index is minimized according to QN control method. If a limited number of velocity outputs can be measured, full velocity responses of the whole structure can be interpolated based on the mode shapes. By defining the performance index function as a combination of the structure's velocity responses and control forces only, feedback gain can be determined according to QN control method with the external excitation being taken into account throughout the entire algorithm. Control forces are then regulated by the time‐invariant feedback gain matrix. The effective location of the active control devices for a building structure subjected to intermediate‐storey excitation has been determined to be in the three floors adjacent to the vibration source. Hence for the purpose of this paper, only the optimal placement of sensors is verified. It is shown in this paper that if the dynamic behaviour of the structure is well described by a mathematical model, sufficient response reduction effect can be achieved according to the new DVFC algorithm, and the degradation of control performance due to time delay can also be verified. Copyright © 2001 John Wiley & Sons, Ltd.     

Detrital zircon multi‐chronology,provenance, and low‐grade metamorphism of the Cretaceous Shimanto accretionary complex,eastern Shikoku,Southwest Japan: Tectonic evolution in response to igneous activity within a subduction zone

Detrital zircon multi‐chronology combined with provenance and low‐grade metamorphism analyses enables the reinterpretation of the tectonic evolution of the Cretaceous Shimanto accretionary complex in Southwest Japan. Detrital zircon U–Pb ages and provenance analysis defines the depositional age of trench‐fill turbidites associated with igneous activity in provenance. Periods of low igneous activity are recorded by youngest single grain zircon U–Pb ages (YSG) that approximate or are older than the depositional ages obtained from radiolarian fossil‐bearing mudstone. Periods of intensive igneous activity recorded by youngest cluster U–Pb ages (YC1σ) that correspond to the younger limits of radiolarian ages. The YC1σ U–Pb ages obtained from sandstones within mélange units provide more accurate younger depositional ages than radiolarian ages derived from mudstone. Determining true depositional ages requires a combination of fossil data, detrital zircon ages, and provenance information. Fission‐track ages using zircons estimated YC1σ U–Pb ages are useful for assessing depositional and annealing ages for the low‐grade metamorphosed accretionary complex. These new dating presented here indicates the following tectonic history of the accretionary wedge. Evolution of the Shimanto accretionary complex from the Albian to the Turonian was caused by the subduction of the Izanagi plate, a process that supplied sediments via the erosion of Permian and Triassic to Early Jurassic granitic rocks and the eruption of minor amounts of Early Cretaceous intermediate volcanic rocks. The complex subsequently underwent intensive igneous activity from the Coniacian to the early Paleocene as a result of the subduction of a hot and young oceanic slab, such as the Kula–Pacific plate. Finally, the major out‐of‐sequence thrusts of the Fukase Fault and the Aki Tectonic Line formed after the middle Eocene, and this reactivation of the Shimanto accretionary complex as a result of the subduction of the Pacific plate.     

Geological and Hydrodynamical Examination of the Bathyal Tsunamigenic Origin of Miocene Conglomerates in Chita Peninsula, Central Japan

A conglomerate appears on a rocky coast called ??Tsubutega-ura Coast??, located on the southwestern coast near the southern tip of the Chita Peninsula, Aichi Prefecture, central Japan. The conglomerate belongs to Miocene sedimentary rocks termed the Morozaki Group. The conglomerate includes meter-scale boulders, indicating that it was formed by an extraordinary event. In the geological investigation, we observed that the conglomerate shows alternate changes of paleocurrent directions between seaward and landward. This feature is supposed to be formed by tsunami currents. In the hydrodynamical investigation, we obtained following results: (1) the lowest limit of a current velocity to move a boulder of about 3?m in diameter would be about 2?C3?m/s, (2) the speed of tsunami currents reproduced by tsunami simulation exceeds 3?m/s at 300?m in depth when the tsunami is generated by a gigantic earthquake with magnitude 9.0 or more, (3) the transport distance of the boulder would be several hundred?meters to several kilometers by one tsunami event caused by a gigantic earthquake. We conclude that tsunamis best explain the formation of the conglomerate deposited in upper bathyal environments about 200?C400?m depth, both from geological and hydrodynamical viewpoints.     

Hydrogen isotopic substitution of solid methylamine through atomic surface reactions at low temperatures: A potential contribution to the D/H ratio of methylamine in molecular clouds

We experimentally studied hydrogen (H)–deuterium (D) substitution reactions of solid methylamine (CH₃NH₂) under astrophysically relevant conditions. We also calculated the potential energy surface for the H–D substitution reactions of methylamine isotopologues using quantum chemical methods. Despite the relatively large energy barrier of more than 18 kJ mol^?1, CH₃NH₂ reacted with D atoms to yield deuterated methylamines at 10 K, suggesting that the H–D substitution reaction proceeds through quantum tunneling. Deuterated methylamines reacted with H atoms as well. On the basis of present results, we propose that methylamine has potential for D enrichment through atomic surface reactions on interstellar grains at very low temperatures in molecular clouds. D enrichment would occur in particular in the methyl group of methylamine.     

Low-temperature single crystal reflection spectra of forsterite

The infrared reflectivities of crystalline forsterite (Mg₂SiO₄) were measured for the temperature range 295–50 K for each crystal axis, between wavenumber 5000 and 100 cm⁻¹. The reflection spectra show clear dependence of temperature; most of the bands become more intense, sharper and their peak positions shift to higher wavenumber with decreasing temperature. Reflection spectra were fitted with dispersion formula of damped oscillator model of the dielectric constants and the oscillator parameters in the model were derived. The absorption spectra of forsterite particle are calculated with the derived dielectric constants to show that the forsterite features are good thermal indicator for cold temperature range below 295 K.