Petrochronology of oxidized granulites from southern Peru

Sapphirine–quartz granulites from the Cocachacra region of the Arequipa Massif in southern Peru record early Neoproterozoic ultrahigh‐temperature metamorphism. Phase equilibrium modelling and zircon petrochronology are used to quantify timing and pressure–temperature (P–T) conditions of metamorphism. Modelling of three magnetite‐bearing sapphirine–quartz samples indicates peak temperatures of >950°C at ~0.7 GPa and a clockwise P–T evolution. Elevated concentrations of Al in orthopyroxene are also consistent with ultrahigh‐temperature conditions. Neoblastic zircon records ages of c. 1.0–0.9 Ga that are interpreted to record protracted ultrahigh‐temperature metamorphism. Th/U ratios of zircon of up to 100 reflect U‐depleted whole‐rock compositions. Concentrations of heavy rare earth elements in zircon do not show systematic trends with U–Pb age but do correlate with variable whole‐rock compositions. Very large positive Ce anomalies in zircon from two samples probably relate to strongly oxidizing conditions during neoblastic zircon crystallization. Low concentrations of Ti‐in‐zircon (<10 ppm) are interpreted to result from reduced titania activities due to the strongly oxidized nature of the granulites and the sequestration of titanium‐rich minerals away from the reaction volume. Whole‐rock compositions and oxidation state have a strong influence on the trace element composition of metamorphic zircon, which has implications for interpreting the geological significance of ages retrieved from zircon in oxidized metamorphic rocks.     

相对论天体力学和后牛顿运动方程

叙述了与Astrod工程有关的相对论天体力学基础内容。包括相对论天体力学、广义相对论基本原理、PPN方法体系、PPN多体问题、PPN二体问题。高阶PN二体问题等     

Asteroid families: evidence of ageing of the proper elements

In recent times it has been emphasized that the present kinematical structures of asteroid families should be evolved with respect to the original post-impact situations, according to numerical simulations performed taking into account also the previously neglected Yarkovsky effect. In this paper we show that also a “classical” approach based on an analysis of the current kinematical properties of families leads to conclude that the distributions of proper eccentricities and semimajor axes of family members exhibit evidence of an evolution. The importance of this approach is that it yields a fully independent and quantitative estimate of an evolutionary spreading of the proper elements. In particular, we find that the original post-impact families had to be on the average about twice more compact than the families we observe now, when considering family members down to about 5 km in size. This result can be used in future analyses to derive estimates of the ages of different families, and to better constrain the typical values of the ejection velocities of the fragments in family-forming events.     

A continuum‐discrete hydromechanical analysis of granular deposit liquefaction

A coupled continuum‐discrete hydromechanical model was employed to analyse the liquefaction of a saturated loose deposit of cohesionless particles when subjected to a dynamic base excitation. The pore fluid flow was idealized using averaged Navier–Stokes equations and the discrete element method was employed to model the solid phase particles. A well established semi‐empirical relationship was utilized to quantify the fluid–particle interactions. The conducted simulations revealed a number of salient micro‐mechanical mechanisms and response patterns associated with the deposit liquefaction. Space and time variation of porosity was a major factor which affected the coupled response of the solid and fluid phases. Pore fluid flow was within Darcy's regime. The predicted response exhibited macroscopic patterns consistent with experimental results and case histories of the liquefaction of granular soil deposits. Copyright © 2004 John Wiley & Sons, Ltd.     

The problem of clock synchronization: A relativistic approach

The problem of synchronization of the Earth-based clocks has been discussed in the framework of General Relativity Theory. The synchronization is considered as the transformation of the observers' proper time scales to the coordinate time scale of local inertial geocentric reference system, which is single for all the observers. The formulas for the relativistic corrections occurring in some methods of Earth-based clock synchronization (transported clock, duplex communication via geostationary satellite and meteor-burst link, LASSO experiments) have been derived enabling one to attain the accuracy of 0.1 ns.     

Minimal Energy and Gravitation

Quantum theory in Robertson – Walker spacetime suggests the existence of a minimal energy ε of the order of 10⁻⁴⁵ erg. Reasonable forms for ε give the expansion factor R=R(t)(t= the cosmic time) with no need of gravitational field equations.Einstein's theory should be modified in gravitational fields of strength less than ε c/ħ ∼ 10⁻⁸ cm/s² where c is the speed of light and ħ is the reduced Planck constant. The cosmological term λ is expected to decrease as the universe expands.In the Appendix, ε is derived from a big bang – big crunch Newtonian cosmology.     

An HLLC Riemann solver for relativistic flows – II. Magnetohydrodynamics

An approximate Riemann solver for the equations of relativistic magnetohydrodynamics (RMHD) is derived. The Harten–Lax–van Leer contact wave (HLLC) solver, originally developed by Toro, Spruce and Spears, generalizes the algorithm described in a previous paper to the case where magnetic fields are present. The solution to the Riemann problem is approximated by two constant states bounded by two fast shocks and separated by a tangential wave. The scheme is Jacobian-free, in the sense that it avoids the expensive characteristic decomposition of the RMHD equations and it improves over the HLL scheme by restoring the missing contact wave.
Multidimensional integration proceeds via the single step, corner transport upwind (CTU) method of Colella, combined with the constrained transport (CT) algorithm to preserve divergence-free magnetic fields. The resulting numerical scheme is simple to implement, efficient and suitable for a general equation of state. The robustness of the new algorithm is validated against one- and two-dimensional numerical test problems.