Pressure dependence on partition coefficients for trace elements between olivine and the coexisting melts

Partition coefficients between olivine and melt at upper mantle conditions, 3 to 14 GPa, have been determined for 27 trace elements (Li, Be, B, Na, Mg, Al, Si, P, K, Ca, Sc, Ti, V, Cr, Mn, Fe, Co, Ni, Cu, Rb, Sr, Y, Zr, Cs, Ba, La and Ce) using secondary-ion mass-spectrometry (SIMS) and electron-probe microanalysis (EPMA). The general pattern of olivine/melt partitioning on Onuma diagrams resembles those reported previously for natural systems. This agreement strongly supports the argument that partitioning is under structural control of olivine even at high pressure. The partition coefficients for mono- and tri-valent cations show significant pressure dependence, both becoming larger with pressure, and are strongly correlated with coupled substitution into cation sites in the olivine structure. The dominant type of trace element substitution for mono- and tri-valent cations into olivine changes gradually from (Si, Mg)↔(Al, Cr) at low pressure to (Si, Mg)↔(Al, Al) and (Mg, Mg)↔(Na, Al) at high pressure. The change in substitution type results in an increase in partition coefficients of Al and Na with pressure. An inverse correlation between the partition coefficients for divalent cations and pressure has been observed, especially for Ni, Co and Fe. The order of decreasing rate of partition coefficient with pressure correlates to strength of crystal field effect of the cation. The pressure dependence of olivine/melt partitioning can be attributed to the compression of cation polyhedra induced by pressure and the compensation of electrostatic valence by cation substitution. Received: March 6, 1997 / Revised, accepted: March 12, 1998     

Response analysis of the Higashi–Kobe Bridge and surrounding soil in the 1995 Hyogoken–Nanbu Earthquake

This paper presents results of observation and analysis of the response of one of the longest cable-stayed bridges in the world to the Hyogoken–Nanbu (Kobe) Earthquake of 17 January 1995. It is determined that interaction of the foundations of the bridge towers with the supporting soil plays a decisive role in the overall structural behaviour. The key factor governing the changes of the soil properties at this site is pore-water pressure buildup, which results in liquefaction of the saturated surface soil layers under large dynamic loads. Models of the soil and structure are created and initially validated by accurately simulating the system response to a small earthquake. Soil parameters reflecting the pore-water pressure buildup in the strong earthquake are determined by an advanced non-linear effective stress analysis, combining the Ramberg–Osgood model of stress–strain dependence with a pore pressure model based on shear work concept. They are utilized to investigate and simulate the interaction of the foundation and the supporting soil using the program SASSI with the flexible volume substructuring approach. The results show a good agreement with the observations and have useful implications to the scientific and engineering practice. © 1998 John Wiley & Sons, Ltd.     

Seasonal variations of spectra of wind speed and air temperature in the mesoscale frequency range

Seasonal variations of the spectra of wind speed and air temperature in the mesoscale frequency range from 1.3 × 10^-4 to 1.5 × 10^-3 Hz (10 min to 2 h periods) have been studied through observations over land for one year. Spectrographs [time series of isopleths of spectral densities, f · S(f) vs f] of wind speed and air temperature contain occasional peaks that are attributed to short-lived mesoscale atmospheric activity with narrow frequency bands. Significant spectral peaks of wind speed were found in 19% of the total observations in winter, and in 15–16% in the other seasons; for air temperature, they occured in 12% of observations in autumn, and in 16–19% in the other seasons. The peaks most often occurred in the period range from 30 min to 1 h; most had durations less than 24 h. Mesoscale fluctuations of wind speed and air temperature were highly correlated, and in most cases, phase differences were 90–180 ° with air temperature leading wind speed. Significant spectral peaks of wind speed often occurred during northerly seasonal cold winds in winter, and accompanied tropical and/or mid-latitude cyclones in the other seasons. When the peaks occurred, wind speed was usually relatively high and the atmospheric surface layer was unstable.     

Simulations of seasonal variations of sulfur compounds in the remote marine atmosphere

A photochemical box model is used to simulate seasonal variations in concentrations of sulfur compounds at latitude 40° S. It is assumed that the hydroxyl radical (OH) addition reaction to sulfur in the dimethyl sulfide (DMS) molecule is the predominant pathway for methanesulfonic acid (MSA) production, and that the rate constant increases as the air temperature decreases. Concentration of the nitrate radical (NO₃) is a function of the DMS flux, because the reaction of DMS with NO₃ is the most important loss mechanism of NO₃. While the diurnally averaged concentration of OH in winter is a factor of about 8 smaller than in summer, due to the weak photolysis process, the diurnally averaged concentration of NO₃ in winter is a factor of about 4–5 larger than in summer, due to the decrease of DMS flux. Therefore, at middle and high latitudes in winter, atmospheric DMS is mainly oxidized by the reaction with NO₃. The calculated ratio of the MSA to SO₂ production rates is smaller in winter than in summer, and the MSA to non-sea-salt sulfate (nssSO₄ ^2-) molar ratio varies seasonally. This result agrees with data on the seasonal variation of the MSA/nssSO₄ ^2- molar ratio obtained at middle and high latitudes. The calculations indicate that during winter the reaction of DMS with NO₃ is likely to be a more important sink of NO_x (NO+NO₂) than the reaction of NO₂ with OH, and to serve as a significant pathway of the HNO₃ production. If dimethyl sulfoxide (DMSO) is produced through the OH addition reaction and is heterogeneously oxidized in aqueous solutions, half of the nssSO₄ ^2- produced in summer may be through the oxidation process of DMSO. It is necessary to further investigate the oxidation products by the reaction of DMS with OH, and the possibility of the reaction of DMS with NO₃ during winter.     

From an initial transient-state to a steady-state in metamorphic reactions: An experimental approach in the system dolomite-quartz-H2O

Abstract We carried out hydrothermal experiments in the system dolomite‐quartz‐H₂O to track the temporal change in reaction rates of simultaneous reactions during the development of reaction zones. Two types of configurations for the starting materials were prepared: dolomite single crystals + quartz powder + water and quartz single crystals + dolomite powder + water, both sealed separately in gold capsules. Runs at 0.1GPa and 600°C with cold seal pressure vessels gave the following results. (i) In short duration (45–71 h) runs metastable layer sequences involving wollastonite and talc occur in the reaction zone, whereas they disappear in longer duration (168–336 h) runs. (ii) The layer sequence of the reaction zones in short duration runs differs from place to place on the dolomite crystal even in the same run. (iii) The diversity of layer sequences in the short duration runs merges into a unique layer sequence in the longer duration runs. (iv) The reaction zone develops locally on the dolomite crystal, but no reaction zone was observed on the quartz crystal in any of the runs. The lines of evidence (i)–(iii) show that the system evolves from an initial transient‐ to a steady‐state and that the kinetic effect is important in the development of reaction zones. A steady diffusion model for the unique layer sequence Qtz/Di/Fo + Cal/Dol + Cal/Dol shows that the Dol + Cal layer cannot be formed by diffusion‐controlled process and that the stability of the layer sequence Qtz/Di/Fo + Cal/Dol depends not only on L‐ratios (a = /L_CaOCaO and b = /L_MgOMgO) but also on the relative rate P = (−2ξ₁ − ξ₂)/(–ξ₁ − 2ξ₂) of competing reactions: Dol + 2Qtz = Di + 2CO₂ (ξ₁) and 2Dol + Qtz = Fo + 2Cal + 2CO₂ (ξ₂). For smaller P the stability field will shift to higher values of a and b. The steady diffusion model also shows that the apparent‐non‐reactivity on the quartz surface can be attributed to void formation in a large volume fraction in the diopside layer.     

Wave energy in surface layers for energy-based damage evaluation

Seismic wave energy in surface layers is calculated based on vertical array records at four sites during the 1995 Hyogo-ken Nambu earthquake by assuming vertical propagation of SH waves. The upward energy generally tends to decrease as it goes up from the base layer to the ground surface particularly in soft soil sites. Theoretical study on 1D multi-layers model to investigate the basic energy flow mechanism indicates that the energy at the ground surface can be smaller on softer soils due to high soil damping during strong shaking even if resonance effect is considered. A simple calculation for a shear-vibrating structure resting on foundation ground shows that induced strain in the structure is directly related to the energy or the energy flux of surface layers. Hence, a general perception that soft soil sites tend to suffer heavier damage than stiff sites should be explained not by greater incident energy but by other reasons such as degree of resonance. Furthermore, it is recommended that not only acceleration or velocity but also S-wave velocity should be specified at a layer where a design seismic motion is given, so that the seismic wave energy can clearly be quantified in seismic design practice.     

Symmetrical fft technique and its applications to earthquake engineering

The objective of this paper is to propose an FFT technique symmetrical for time and frequency, in which the symmetrization is done by introducing generating functions. By the symmetrical treatment, it becomes possible to derive several complementary relations in the time and frequency domains. Based on the idea of the symmetrical FFT technique, the phase problem of seismic waves is discussed not only in the frequency domain but also in the time domain. Furthermore, the deconvolution technique to evaluate the minimum-phase-shift and all-pass functions of seismic waves, which is not discussed sufficiently in earthquake engineering, is developed in the two domains. Numerical examples to present these basic ideas are also illustrated in the symmetrical framework.     

Finite element analysis and experimental verification of the scrap tire rubber pad isolator

A new base isolation system using scrap tire rubber pads (STRP) has been introduced for seismic mitigation of ordinary residential buildings. The rubber and the steel reinforcing cords used in manufacturing the tire are the alternative materials of the proposed base isolation system. The steel reinforcing cords represent the steel plates used in conventional laminated rubber bearings. These steel reinforcing cords shall prevent the lateral bulging of the rubber bearing. The proposed base isolation system has no bonding between the superstructure and the foundation beam which allows for rollover deformation. In the first part of the study, the STRP layers were just stacked one on top of another without applying the adhesive. This paper presents loading test as well as finite element analysis (FE analysis) of strip STRP isolators that are subjected to any given combination of static vertical and lateral loads. The results of the static vertical and horizontal loading test conducted on STRP isolators were used to calculate the mechanical properties of the isolators, including stiffness and damping values. The load–displacement relationship of STRP isolators were compared between experimental and FE analysis results and the results were found to be in close agreement. The stress state within the STRP isolators was also analyzed in order to estimate the maximum stress demand in the rubber and steel reinforcing cords. These STRP isolators have several advantages over conventional laminated rubber bearings including superior damping properties, lower incurred cost, light weight and easily available material. This study suggests that using the STRP as low cost base isolation device for ordinary residential buildings is feasible.     

Groundwater-level Changes Due to Pressure Gradient Induced by Nearby Earthquakes off Izu Peninsula, 1997

Anomalous water level changes were observed at two wells associated with seismic swarm activity off Izu Peninsula on March, 1997. These are coseismic water level drops followed by gradual postseismic water level rise at the time of large earthquakes during the swarm activity. The post-seismic water level rises, which can be fitted by an exponential function with a time constant of about six hours, are explained in terms of the horizontal pressure diffusion due to the pressure gradient in the aquifer induced by the coseismic static strain.     

Resonance capacity method for earthquake response analysis of hysteretic structures

The Resonance Capacity Method is proposed for the earthquake response analysis of hysteretic structures. Resonance Capacity is a physical quantity of structures which is related to the hysteretic energy absorbed by structures in one cycle and is equated to the acceleration, velocity and displacement amplitudes α₀, d₀ and d₀ of earthquake ground motions at resonance.¹ According to the idealized trapezoidal approximation of earthquake ground motions in the logarithmic period–velocity plane as proposed by Veletsos and Newmark,⁸ the Resonance Capacity property applies in each period range, short, medium and long, where α₀, v₀ and d₀ respectively are approximately constant. In the medium range of periods, the energy dissipated in hysteretic loops and the deformation amplitudes of a single-degree system with elasto–plastic force–deformation relationships are calculated for the case of El Centro 1940, 18 May earthquake, by this Resonance Capacity Method. The result is compared with results from conventional numerical response analyses obtained by Berg and Thomaides,¹⁴ Kato and Akiyama¹² and Veletsos and Newmark,⁸ and the general agreement is seen to be good. Therefore, it may be possible to apply this Resonance Capacity Method over the entire range of periods. By means of this method the earthquake response analysis of hysteretic systems can be performed easily, and the hysteretic energy and fatigue characteristics of structures may be taken into account directly, up to the point of fracture.