Disequilibrium for Ca during growth of pelitic garnet

Compositional zoning in hundreds of almandine-rich garnets in amphibolite by facies micaceous quartzites from the Picuris Range, north-central New Mexico USA, indicates that although Mn, Mg and Fe achieve chemical equilibrium at hand-sample scale during garnet growth, Ca does not. Instead, Ca concentrations at the surface of growing garnets appear to depend strongly on kinetic factors that govern the local chemical environment, yielding disequilibrium for Ca at scales larger than the region immediately surrounding an individual porphyroblast. Detailed zoning profiles were obtained for 371 garnet crystals in a small volume of a single sample of garnetiferous quartzite, and core analyses were made of 97 additional crystals. Each analysis was made on a section that passed precisely through the morphological centre of the crystal, located by means of 3-D imagery from computed X-ray tomography. The data reveal strong correlations between crystal size and concentrations of Mn, Mg and Fe (but not Ca) in garnet cores; a relationship between crystal size and isolation; rigorous cross-correlations among concentrations of Mn, Mg and Fe (but not Ca); and systematic variations in Ca concentrations as a function of crystal size and core composition that are anomalous in comparison to the behaviour of the other divalent cations. We interpret these observations as the result of thermally accelerated diffusion-controlled garnet growth, in circumstances that promoted rapid intergranular diffusion and thus rock-wide equilibration of Mn, Mg and Fe, but that prevented equilibration at similar scale for Ca because of its more sluggish intergranular diffusion. The anomalous behaviour of Ca is made evident in these garnets by the presence of sharp spikes in Ca concentration, which are demonstrably not a consequence of any simultaneous rock-wide event, such as a change in pressure, temperature, or some other intensive parameter. Instead, Ca concentrations probably reflect the local extent of reaction in the immediate vicinity of each porphyroblast. To the degree that such kinetic factors introduce departures from chemical equilibrium for Ca, thermobarometric estimates that involve grossular contents of pelitic garnet will be in error.     

Fully normalized spherical cap harmonics: application to the analysis of sea-level data from TOPEX/POSEIDON and ERS-1

We introduce two sets of fully normalized harmonics for the spectral analysis of functions defined on a spherical cap. The harmonics are the products of Fourier functions and the fully normalized associated Legendre functions of non-integer degree. Using Sturm-Liouville theory for boundary-value problems, we present two convenient and stable formulae for computing the zeros of the associated Legendre functions that form two sets of orthogonal functions. Formulae for the stable numerical evaluation of the fully normalized associated Legendre functions of non-integer degree that avoid the gamma function are also derived. The result from the expansions of sea-level anomaly from altimetry into Set 2 fully normalized cap harmonics shows fast convergence of the series, and the degree variances decay rapidly without aliasing effects. The zero-degree coefficients (Set 2) of sea-level anomaly from TOPEX/ POSEIDON (T/P) and ERS-1 indicate an El Niño event during 1993 January-1993 July, and a La Niño event during 1993 November-1994 July, although the ERS-1 result is less obvious. Ocean circulations over the South China Sea and the Kuroshio area are clearly identified with the low-degree expansions of sea-surface topography (SST) from T/P and ERS-1. A cold-core eddy of 4° in diameter centred at 17.5°N, 118 E was detected with the expansion of SST from T/P cycle 47, and a property of the cap harmonics is used to compute this eddy's kinetic energy. The kinetic energy is at a low in winter and high in summer, and its variation seems to be periodic with an amplitude of 0.4 m² S^-2.     

A general correspondence principle for time-domain electromagnetic wave and diffusion fields

A general correspondence principle is presented that relates any time-domain electromagnetic diffusion field to an electromagnetic wavefield in a 'corresponding' configuration. The principle applies to arbitrarily inhomogeneous and anisotropic media and arbitrary transmitters and receivers. For the correspondence between the two types of electromagnetic fields to hold, the electric conductivity in the diffusive case and the permittivity in the wavefield case should have the same spatial variation, while the permeability distributions in space in the two cases are to be identical. Essential steps in the derivation of the correspondence principle are the use of the time Laplace transformation of causal signals, taken at real, positive values of the transform parameter, the Schouten-Van der Pol theorem in the theory of the Laplace transformation, and the reliance upon Lerch's theorem of the uniqueness of the interrelation between causal field quantities and their time-Laplace-transform representations at real, positive values of the transform parameter. Correspondence is then established between the tensorial Green's functions in the two cases, where the Green's functions are the point-receiver responses (either electric or magnetic field) to point-transmitter excitations (either electric- or magnetic-current source).
Through the correspondence principle, all transient electromagnetic wavefields (where losses are neglected) have as a counterpart a transient diffusive electromagnetic field (where the electric displacement current is neglected). The interrelation yields the tool to compare quantitatively the potentialities of the two types of fields in transient electromagnetic geophysical prospecting.
Finally, a general medium-parameter scaling law for time-domain electromagnetic wavefields is presented.     

The respiratory environment of the Namib Desert Golden Mole

We measured the partial pressure of oxygen (PO₂) in the interstitial gas surrounding the sand-swimming Namib moleEremitalpa granti namibensis. At a sand temperature of 26 °C, which produced a nearly maximal rate of oxygen consumption, thePO₂near the noses of the animals averaged only 0·9 kPa (6·7 Torr) below the level in the free atmosphere. High oxygen availability was a result of the notably low metabolic rate in the 20 g mammals and the dry, porous and metabolically inactive nature of dune sand. A mathematical model indicated that normal mammals weighing 200 g or more could comfortably exist completely encased in dune sand. We concluded that the moles' small size and low metabolic rate are not adaptations to hypoxia or hypercapnia underground but are probably related to low food availability and the energetic cost of foraging in their desert environment.     

地质速率计—矿物氧扩散的应用

概述了同位素封闭体系内的矿物氧扩散和同位素交换机制及其在地质速率计上的应用。火成岩从高温冷却或变质岩从高峰主质温度冷却过程中，由于冷却速度不同，扩散作用导致的矿物晶体内部及晶粒间氧同位素再平衡也有所不同。通过实测岩石中各组成矿物氧同位素比值，模式含量和颗粒半径，据矿物氧扩散和同位素交换模型，可以估算出岩石的冷却速率。     

GOCE Gravity Field Processing Strategy

A processing strategy and the corresponding software architecture for the processing of GOCE (Gravity field and steady-state Ocean Circulation Explorer) observables is presented and described, with the major objective to compute a high-accuracy, high-resolution spherical harmonic model of the Earth's gravity field. The combination of two numerical solution strategies, i.e. the rigorous solution of the corresponding large normal equation systems applying parallel processing (on a PC cluster) as the core solver, and the fast semianalytic approach as a quick-look gravity field analysis (QL-GFA) tool, is proposed. Such a method fusion benefits from the advantages of the individual components: the rigorous inversion of the system providing also the full variance-covariance information, and the quickness enabling the consecutive production of intermediate gravity field solutions, for the purpose to analyse partial and incomplete data sets and to derive a diagnosis of the performance of the GOCE measurement system. The functionality and operability of the individual components are demonstrated in the framework of a closed loop simulation, which is based on a realistic mission scenario both in terms of the orbit configuration and the coloured measuring noise. Special concern is given to the accuracy of the recovered coefficients, the numerical behaviour, the required computing time, and the particular role of the individual modules within the processing chain. In the case of the core solver, it is demonstrated that the assembling and rigorous solution of large normal equation systems can be handled by using Beowulf clusters within a reasonable computing time. The application of the quick-look tool to partial data sets with short-term data gaps is demonstrated on the basis of several case studies. Additionally, the spectral analysis of the residuals of the adjustment is presented as a valuable tool for the verification of the noise characteristics of the GOCE gradiometer.     

Multicomponent ionic diffusion in porewaters: Coulombic effects revisited

The diffusion of an ion in porewaters cannot occur independently of the other ions in solution as a result of Coulombic coupling, as well as from other effects not considered here. Unfortunately, a longstanding disagreement exists about the correct form and meaning of the equations that describe Coulombic coupling in porewaters, i.e., Ben-Yaakov [Am. J. Sci. 281 (1981) 974] vs. Lasaga [Am. J. Sci. 281 (1981) 981]. This paper re-examines this controversy by reformulating the problem starting from fundamental concepts of mass and charge conservation. We show that these antagonistic formulations are both valid and, in fact, equivalent, when the different interpretations of charge balance are resolved. Most of the disagreements between Ben-Yaakov and Lasaga are then shown to result from differing methods of solution, not fundamental disparities in their models. We note, however, that the explanation for the concept of “stationary” gradients of nonreacting ions as given Ben-Yaakov is inaccurate, and such gradients do lead to diffusive fluxes that are counterbalanced by electrochemical migrational fluxes to produce no net flux (excluding advective flux). We further find that the bicarbonate diffusive flux will not balance the diffusional charge flux of sulfate during its reduction if advection is present. This latter imbalance generates compensating fluxes in the other nonreacting ions. We have applied our theory to a simplified case of sulfate reduction in a marine sediment. The results show that nonreacting ions do diffuse and that with normally expected values of porewater advection, the ratio of the bicarbonate to the sulfate flux can be far different than the ideal value of −2.     

用混合和扩散云室对大气冰核浓度的观测分析

继1963年用混合云室对春季冰核浓度观测之后，1994～1996年，在同季节、同地点用混合和静力扩散云室又做了同步观测。对取得的资料进行整理分析，得出了凝华核浓度在不同天气背景下占总冰核浓度的百分比，两种云室测值的对应关系；分析了扩散云室的容积效应及冰核浓度随气压、能见度等气象因子的变化。     

时域瞬变场电磁场有限差分法

模拟电磁场在地下介质中的传播规律是理解地下介质响应的重要手段，直接应用有限差分法在时间域对二维、三维瞬变场进行数值分析是一种有效方法。应用反映电磁场基本规律的麦克斯韦方程组，导出时域场量的齐次扩散方程，对所研究的空间区域作差分离散，源作为初始条件加入，采用合适的边界条件，利用差分方程进行计算，从而得到场量的数值结果，展现瞬变场在地下随时间扩散的全过程，对电磁场的瞬态响应有更加直观的了解。