Molecular dynamics simulations of the phase transition between low-temperature and high-temperature clinoenstatites

 Phase transition between low-temperature clinoenstatite (LT-CEn) and high-temperature clinoenstatite (HT-CEn) was studied by using molecular dynamics (MD) simulations, based on empirical potential parameters. Starting from LT-CEn, the MD calculations were carried out at atmospheric pressure and at elevated pressures (1–6 GPa). At elevated temperatures the transformation from the starting LT-CEn to HT-CEn occurred at any pressure. It was confirmed that the HT-CEn has the same space group C2/c as diopside but the M2 site is six-coordinated, unlike diopside. A significant difference in the MD-simulated cell volumes between LT-CEn and HT-CEn was also observed, showing a first-order transition. In addition, there were some temperature ranges where LT-CEN and HT-CEn would be coexistent and very small thermal hystereses between increasing and decreasing temperatures during the transition. These behaviors are consistent with the characteristic of a thermoelastic-martensitic transformation. The phase boundary between LT-CEn and HT-CEn was determined for the first time. Its positive dT/dP slope strongly shows that the high-pressure clinoenstatite is a significantly distinct phase from HT-CEn although the both phases have the same space group, C2/c. Received: 8 November 2000 / Accepted: 28 April 2001     

The effects of liquefaction-induced lateral spreading on pile foundations

Observations of pile foundation performance during previous earthquakes have shown that pile failure has been caused by lateral ground movements resulting from soil liquefaction. The recognition that lateral ground movements may play a critical role in pile performance during an earthquake has important implications for design and risk assessment, and requires that analytical models be devised to evaluate these potential problems.In this paper, parametric studies were conducted to estimate the maximum bending moments induced in piles subjected to lateral ground displacement. The results are summarized in charts using dimensionless parameters.The analyses reveal that the existence of a nonliquefiable layer at the ground surface can affect significantly the maximum bending moment of the pile. When a relatively thick nonliquefiable layer exists above a liquefiable layer, neither the material nonlinearity of the soil nor loss of soil stiffness within the liquefiable layer significantly affect the maximum bending moment. When the thickness of the liquefiable soils is greater than about three times that of an overlying intact layer, soil stiffness in the liquefiable layer must be chosen carefully when evaluating the maximum bending moment.     

Equivalent linear method considering frequency dependent characteristics of stiffness and damping

Equivalent linear dynamic response analysis of ground is based on complex moduli and Fourier series expansion; therefore, it is not an equivalent method but an approximate method. Two deficiencies in the conventional equivalent linear method represented by SHAKE are described first. The maximum shear strength is overestimated, resulting in overestimation of the peak acceleration under a strong ground motion, and the amplification is underestimated at high frequency. The latter sometimes results in underestimation of the peak acceleration under weak ground shaking, and gives an incident wave with unrealistic large accelerations or a divergence of analysis in deconvolution analysis under strong ground motion. Both deficiencies are shown to come from the same cause, i.e. computing the effective strain as a constant fraction of the maximum strain. Since this is a key concept of the equivalent linear analysis, one cannot overcome both deficiencies at the same time in the conventional method. An apparent frequency dependence in stiffness and damping is shown to appear in the dynamic response, although soil itself does not show frequency dependent characteristics. Following this observation, the effective strain is expressed in terms of frequency from the similarity concept of the strain–frequency relationship between time domain and frequency domain. This enables the reduction of both deficiencies at the same time, resulting in a marked improvement in the equivalent linear analysis. The accuracy of the proposed method is examined by the simulations of three vertical array records during large earthquakes. The proposed method always gives much better prediction than conventional equivalent linear methods for both convolution and deconvolution analyses, and it is confirmed to be applicable at more than 1% shear strain.     

Structural characteristics off Miyagi forearc region, the Japan Trench seismogenic zone, deduced from a wide-angle reflection and refraction study

The Japan Trench subduction zone, located east of NE Japan, has regional variation in seismicity. Many large earthquakes occurred in the northern part of Japan Trench, but few in the southern part. Off Miyagi region is in the middle of the Japan Trench, where the large earthquakes (M > 7) with thrust mechanisms have occurred at an interval of about 40 years in two parts: inner trench slope and near land. A seismic experiment using 36 ocean bottom seismographs (OBS) and a 12,000 cu. in. airgun array was conducted to determine a detailed, 2D velocity structure in the forearc region off Miyagi. The depth to the Moho is 21 km, at 115 km from the trench axis, and becomes progressively deeper landward. The P-wave velocity of the mantle wedge is 7.9–8.1 km/s, which is typical velocity for uppermost mantle without large serpentinization. The dip angle of oceanic crust is increased from 5–6° near the trench axis to 23° 150 km landward from the trench axis. The P-wave velocity of the oceanic uppermost mantle is as small as 7.7 km/s. This low-velocity oceanic mantle seems to be caused by not a lateral anisotropy but some subduction process. By comparison with the seismicity off Miyagi, the subduction zone can be divided into four parts: 1) Seaward of the trench axis, the seismicity is low and normal fault-type earthquakes occur associated with the destruction of oceanic lithosphere. 2) Beneath the deformed zone landward of the trench axis, the plate boundary is characterized as a stable sliding fault plain. In case of earthquakes, this zone may be tsunamigenic. 3) Below forearc crust where P-wave velocity is almost 6 km/s and larger: this zone is the seismogenic zone below inner trench slope, which is a plate boundary between the forearc and oceanic crusts. 4) Below mantle wedge: the rupture zones of thrust large earthquakes near land (e.g. 1978 off Miyagi earthquake) are located beneath the mantle wedge. The depth of the rupture zones is 30–50 km below sea level. From the comparison, the rupture zones of large earthquakes off Miyagi are limited in two parts: plate boundary between the forearc and oceanic crusts and below mantle wedge. This limitation is a rare case for subduction zone. Although the seismogenic process beneath the mantle wedge is not fully clarified, our observation suggests the two possibilities: earthquake generation at the plate boundary overridden by the mantle wedge without serpentinization or that in the subducting slab.     

Origin of three-dimensional shapes of chondrules: I. Hydrodynamics simulations of rotating droplet exposed to high-velocity rarefied gas flow

The origin of three-dimensional shapes of chondrules is an important information to identify their formation mechanism in the early solar nebula. The measurement of their shapes by using X-ray computed topography suggested that they are usually close to perfect spheres, however, some of them have rugby-ball-like (prolate) shapes [Tsuchiyama, A., Shigeyoshi, R., Kawabata, T., Nakano, T., Uesugi, K., Shirono, S., 2003. Lunar Planet. Sci. 34, 1271-1272]. We considered that the prolate shapes reflect the deformations of chondrule precursor dust particles when they are heated and melted in the high velocity gas flow. In order to reveal the origin of chondrule shapes, we carried out the three-dimensional hydrodynamics simulations of a rotating molten chondrule exposed to the gas flow in the framework of the shock-wave heating model for chondrule formation. We adopted the gas ram pressure acting on the chondrule surface of in a typical shock wave. Considering that the chondrule precursor dust particle has an irregular shape before melting, the ram pressure causes a net torque to rotate the particle. The estimated angular velocity is for the precursor radius of r₀=1 mm, though it has a different value depending on the irregularity of the shape. In addition, the rotation axis is likely to be perpendicular to the direction of the gas flow. Our calculations showed that the rotating molten chondrule elongates along the rotation axis, in contrast, shrinks perpendicularly to it. It is a prolate shape. The reason why the molten chondrule is deformed to a prolate shape was clearly discussed. Our study gives a complementary constraint for chondrule formation mechanisms, comparing with conventional chemical analyses and dynamic crystallization experiments that have mainly constrained the thermal evolutions of chondrules.     

A simplified 3D model for soil-structure interaction with radiation damping and free field input

This paper addresses the time history finite element analysis of rock-structure interaction. Modeled is not only the lateral energy dissipation, but also the interaction between the far field and the numerical model itself. This is accomplished by a preliminary analysis of the far field as a shear beam (for lateral excitation), and then velocities and displacements are transferred to the model as nodal forces through damping and stiffness matrices respectively. Details of the finite element implementation are given, along with an extensive series of simulations comparing this method, with the one of Lysmer for both 2D and 3D models. The model is derived from the principle of virtual work, and its implementation does not require any modification of existing finite element codes, only clever pre and postprocessing of results are needed.     

Cosmogenic nuclides in the Brenham pallasite

Abstract— Cosmic‐ray‐produced (cosmogenic) nuclides were studied in fragments of the Brenham pallasite, a large stony iron meteorite. The contents of light noble gases (He, Ne, and Ar) and long‐lived radionuclides (¹⁰Be, ²⁶Al, ³⁶Cl, and ⁵³Mn), produced by nuclear reactions with cosmic rays, were measured in the separated metal and olivine phases from numerous samples representing a wide range of shielding conditions in the meteoroid. The distribution of cosmogenic nuclide concentrations in the metal follows patterns similar to that observed in large iron meteorites. Shielding effects were estimated from the relative proportions of low‐ and high‐energy reaction products. The production rates varied, from surface to interior, by a factor of more than several hundred. The ³⁶Cl‐³⁶Ar cosmic‐ray exposure age of Brenham is 156 ± 8 Myr. This determination is based on a multiple nuclide approach that utilizes cosmogenic nuclide pairs. This approach not only yields a “shielding independent” exposure age but also demonstrates that the production of cosmogenic nuclides occurred in a single stage. The depth profiles of ¹⁰Be in the stone phase and ⁵³Mn in the metal phase are shown superimposed on corresponding profiles from the Apollo 15 long drill core. Surprisingly low abundances of lithophile elements, such as K, U, and Th, provided a unique opportunity to examine the production systematics of those nuclides whose inventories typically have significant contributions from non‐cosmogenic sources, particularly radiogenic contributions. The U and Th contents of the olivine samples are extremely low, allowing detection of cosmogenic ⁴He production from oxygen, magnesium, silicon, and iron.     

Chemical characteristics of the Cretaceous-Tertiary boundary layer at Gubbio, Italy

Platinum metals (Pd, Ir, Pt) and Au were determined by radiochemical neutron activation analysis (RNAA) on Cretaceous-Tertiary (K---T) boundary samples collected at Gubbio, Italy. In addition, about thirty elements were determined by instrumental NAA (INAA) for both the Cretaceous and the Tertiary layer samples, as well as the boundary clay samples. Iridium and other Pt metals were observed to be similarly enriched, relative to CI-chondrite, in the boundary layers, i.e., they occur in chondritic ratios. On the other hand, correlations among the metals are not very good, suggesting that Pt metals were not incorporated into clay materials in a single phase, but rather, they behaved separately according to their own solution chemistries. In the Gubbio section, the chemical composition does not change so much across the K---T boundary layer, except for Ir (and possibly other Pt group elements) and is very similar to that of the “North American shale composite (NASC)”. Arsenic and Sb were observed not to be anomalously enriched in the Gubbio K---T boundary layers. This implies that their abundances cannot be a clue in elucidating an event which must have happened at the K---T boundary period. NASC-normalized REE abundance patterns for the boundary samples show characteristic features due to some REE precipitation out of seawater. REE abundances are not variable among the Cretaceous, the Tertiary, and the boundary samples, suggesting that conditions for precipitation and sedimentation were similar across the K---T boundary at Gubbio.