由三元系相图计算活度的方法Ⅱ.有二元同分熔点化合物的体系

利用熔化自由能 ,由具有同分熔点二元化合物的三元系相图计算活度。计算的范围遍及三元相图的液相面、液相线、共熔线和三元低共熔点     

Anelastic structure of the upper mantle beneath the northern Philippine Sea

We investigated the upper mantle anelastic structure beneath the northern Philippine Sea region, including the Izu-Bonin subduction zone and the Shikoku Basin. We used regional waveform data from 69 events in the Pacific and the Philippine Sea slabs, recorded on F-net and J-array network broadband stations in western Japan. Using the S–P phase pair method, we obtained differential attenuation factors, δt^*, which represent the relative whole path Q. We conducted a tomographic inversion using 978 δt^* values to invert for a fine-scale (50–100 km) three-dimensional anelastic structure.

The results shows two high-Q regions (Q_P>1000) which are consistent with the locations of the Pacific and the Philippine Sea slabs. Also there is a low-Q (Q_P110) area extending to the deeper parts (350–400 km) of the model just beneath the old spreading center and the Kinan Seamount Chain in the Shikoku Basin. A small depth dependence of the laterally averaged Q_P was found, with values of 266 (0–250 km), 301 (250–400 km), and 413 (400–500 km).     

Practical VTI approximations: a systematic anatomy

Transverse isotropy (TI) with a vertical symmetry axis (VTI) often provides an appropriate earth model for prestack imaging of steep-dip reflection seismic data. Exact P-wave and SV-wave phase velocities in VTI media are described by complicated equations requiring four independent parameters. Estimating appropriate multiparameter earth models can be difficult and time-consuming, so it is often useful to replace the exact VTI equations with simpler approximations requiring fewer parameters. The accuracy limits of different previously published VTI approximations are not always clear, nor is it always obvious how these different approximations relate to each other. Here I present a systematic framework for deriving a variety of useful VTI approximations. I develop first a sequence of well-defined approximations to the exact P-wave and SV-wave phase velocities. In doing so, I show how the useful but physically questionable heuristic of setting shear velocities identically to zero can be replaced with a more precise and quantifiable approximation. The key here to deriving accurate approximations is to replace the stiffness a₁₃ with an appropriate factorization in terms of velocity parameters. Two different specific parameter choices lead to the P-wave approximations of Alkhalifah (Geophysics 63 (1998) 623) and Schoenberg and de Hoop (Geophysics 65 (2000) 919), but there are actually an infinite number of reasonable parametrizations possible. Further approximations then lead to a variety of other useful phase velocity expressions, including those of Thomsen (Geophysics 51 (1986) 1954), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Harlan (Stanford Exploration Project Report 89 (1995) 145), and Stopin (Stopin, A., 2001. Comparison of v(θ) equations in TI medium. 9th International Workshop on Seismic Anisotropy). Each P-wave phase velocity approximation derived this way can be paired naturally with a corresponding SV-wave approximation. Each P-wave or SV-wave phase velocity approximation can then be converted into an equivalent dispersion relation in terms of horizontal and vertical slownesses. A simple heuristic substitution also allows each phase velocity approximation to be converted into an explicit group velocity approximation. From these, in turn, travel time or moveout approximations can also be derived. The group velocity and travel time approximations derived this way include ones previously used by Byun et al. (Geophysics 54 (1989) 1564), Dellinger et al. (Journal of Seismic Exploration 2 (1993) 23), Tsvankin and Thomsen (Geophysics 59 (1994) 1290), Harlan (89 (1995) 145), and Zhang and Uren (Zhang, F. and Uren, N., 2001. Approximate explicit ray velocity functions and travel times for P-waves in TI media. 71st Annual International Meeting, Society of Exploration Geophysicists, Expanded Abstracts, 106–109).     

Computational study of tetrahedral Al–Si and octahedral Al–Mg ordering in phengite

 As part of a wider study of the nature and origins of cation order–disorder in micas, a variety of computational techniques have been used to investigate the nature of tetrahedral and octahedral ordering in phengite, K₂ ^[6](Al₃Mg)^[4](Si₇Al)O₂₀(OH)₄. Values of the atomic exchange interaction parameters J _n used to model the energies of order–disorder were calculated. Both tetrahedral Al–Si and octahedral Al–Mg ordering were studied and hence three types of interaction parameter were necessary: for T–T, O–O and T–O interactions (where T denotes tetrahedral sites and O denotes octahedral sites). Values for the T–T and O–O interactions were taken from results on other systems, whilst we calculated new values for the T–O interactions. We have demonstrated that modelling the octahedral and tetrahedral sheets alone and independently produces different results from modelling a whole T–O–T layer, hence justifying the inclusion of the T–O interactions. Simulations of a whole T–O–T layer of phengite indicated the presence of short-range order, but no long-range order was observed. Received: 8 August 2002 / Accepted: 14 February 2003 Acknowledgements The authors are grateful to EPSRC (EJP) and the Royal Society (CIS) for financial support. Monte Carlo simulations were performed on the Mineral Physics Group's Beowulf cluster and the University of Cambridge's High Performance Computing Facility.     

Global thermal and dynamical perturbations due to Cretaceous mantle avalanche

The numerical models of mantle convection agree to depict avalanches behaviour according to the level of endothermicity of the spinel → perovskite phase change. Their potential effects on the global thermal and dynamical states of the mantle have been computed thanks to a numerical code, which takes into account both the 400-km exothermic and the 660-km endothermic phase changes. The cycle followed by the avalanches is: local layering, destabilization of the 660-km thermal layer, travelling and spreading on the core, and reappearing of the local layering. Therefore, mantle convection is characterized by quiet periods of partial layering embedded in catastrophic events. During the avalanche, the amplitude of the surface velocity is multiplied by two, which would imply an enhanced plate tectonic and ridge activities. The global thermal effects of the avalanche are compatible with a high mantle temperature and an acceleration of Earth's rotation during the Cretaceous. They also offer a coherent explanation to locate the origin of mantle plumes both within the CMB and just below the transition zone.     

确定精密重力仪相位滞后方法的研究

从系统辨识的原理出发,给出了确定精密重力仪相位滞后的一种方法,即记录精密重力仪在阶跃信号作用下的输出数据,利用这些数据对仪器的传递函数进行辨识,从而确定重力仪对各潮汐波产生的相位滞后角度.用这种方法对ET-21地潮重力仪进行了相位滞后的确定,并给出这台仪器对15个主要潮汐波的相位滞后角度.     

Observation of gravity waves in a boreal forest

In this paper we report the results of the analysis of two 60-min wave events that occurred in a boreal aspen forest during the 1994 BOREAS (Boreal Ecosystems-Atmosphere Study) field experiment. High frequency wind and temperature data were provided by three 3-D sonic anemometer/thermometers and fourteen fine-wire thermocouples positioned within and above the forest. Wave phase speeds, estimated from information revealed by spectral analysis and linear plane wave equations, are 2.2 and 1.3 m s^-1 for the two events. The wavelengths are 130 m and 65 m respectively and are much larger than the vertical wave displacements. There is strong evidence from the present analysis and from the literature supporting our postulate that these waves are generated by shear instability. We propose that wind shear near the top of the stand is often large enough to reduce the gradient Richardson number below the critical value of 0.25 and thus is able to trigger the instability. When external conditions are favorable, the instability will grow into waves.     

滑坡监测的新方法——PMP测量方法

文章从国内外滑坡灾害、滑坡灾害的危害及现有滑坡监测方法等方面，论述了滑坡灾害带来的危害性，从而验证了滑坡监测的重要性。但目前现行的滑坡监测手段各有不足之处。针对这种情况，提出了新的测量方法改进其不足之处。位相测量剖面术的测量原理是将LCD投影仪安装于滑坡体上，在滑坡体外稳定地带放置监测系统，CCD光束中心垂直于滑坡体位移方向，监测系统由望远镜头和CCD组成，CCD输出的视频信号进入计算机。通过软件对图像进行处理，绘制出滑坡整体的三维面貌，并计算出滑动的距离，从而有效的观测滑坡体的运动过程。     

High-T phase transition of synthetic ANaB(LiMg)CMg5Si8O22(OH)2 amphibole: an X-ray synchrotron powder diffraction and FTIR spectroscopic study

The synthetic amphibole Na_0.95(Li_0.95Mg_1.05)Mg₅Si₈O₂₂(OH)₂ was studied in situ at high-T, using IR OH-stretching spectroscopy and synchrotron X-ray powder diffraction. At room-T the sample has P2₁ /m symmetry, as shown by the FTIR spectrum. It shows in the OH region two well-defined and intense absorptions at 3,748 and 3,712 cm⁻¹, respectively, and two minor bands at 3,667 and 3,687 cm⁻¹. The main bands are assigned to the two independent O–H groups in the primitive structure. The two minor bands evidencing the presence of small amount of vacant A-site (^A□_0.05). With increasing T, these bands shift continuously and merge into a unique absorption at high temperature. A change as a function of increasing T is revealed by the evolution of the refined unit-cell parameters, whose trend shows a transition to C2/m at about 320–330°C. The spontaneous scalar strain, fitted with a tricritical 2–6 Landau potential, gives a T _c of 325(10)°C (β parameter = 0.27). Comparison with the second-order P2₁ /m ⇔ C2/m phase transition at 255°C for synthetic amphibole ^ANa_0.8^B(Na_0.8Mg_1.2)^CMg₅Si₈O₂₂(OH)₂ indicates that the substitution of Na with Li at the B-sites strongly affects the thermodynamic character and the T _c of the phase transition. The comparison of LNMSH amphiboles with cummingtonitic ones shows that the high-T thermodynamic behaviour is affected by A-site occupancy.