Experimental study of a low-alkali tholeiite at 1–5 kbar: optimal condition for the crystallization of high-An plagioclase in hydrous arc tholeiite

We conducted melting experiments on a low-alkali tholeiite (SiO₂ ~52 wt%, MgO ~6.5 wt%, CaO/Na₂O~4.4, Al₂O₃/SiO₂ ~0.33) under both H₂O-undersaturated and H₂O-saturated conditions to investigate the effect of H₂O on the Ca–Na partitioning between plagioclase and melt. Experiments were performed in the temperature and pressure ranges of 1,000–1,300°C and 1–5 kbar, respectively, with varying H₂O contents of 0–12wt%. Redox condition was 0–2 log unit above NNO (nickel–nickel oxide) buffer. Temperature-bulk H₂O diagrams for the low-alkali tholeiite are constructed at 1, 2, and 5 kbar, and compositions of near-liquidus plagioclase and coexisting melt are determined. To exclude the effect of melt composition (CaO/Na₂O and Al₂O₃/SiO₂ ratios) on plagioclase composition and to reveal the effect of H₂O on An (=100×Ca/(Ca+Na)) content and (=(Ca/Na)_pl/(Ca/Na)_melt), we focused on the composition of near-liquidus plagioclases which crystallized from melts with nearly constant CaO/Na₂O and Al₂O₃/SiO₂ ratios. Our experimental results show that, at each experimental pressure, An content of the near-liquidus plagioclase and the K_D^Ca-Na almost linearly increases as H₂O content in melt increases. Each of the An content and the variations in a low-alkali tholeiitic system (CaO/Na₂O~4.0–4.5, Al₂O₃/SiO₂ ~0.27–0.33) can be described by one equation using temperature, pressure, and melt H₂O content as parameters. An content and of liquidus plagioclase increases with increasing melt H₂O and with decreasing pressure, elucidating that nearly H₂O-saturated conditions of 2–3 kbar is optimal for the crystallization of the most An-rich plagioclase (>An₈₈). We suggest this pressure condition of 2–3 kbar, corresponding to depth of 7–11 km, plays an important role for the origin of An-rich plagioclase in H₂O-rich low-alkali tholeiite. At pressures more than ca. 4 kbar, crystallization of liquidus Ca-rich clinopyroxene decreases the CaO/Na₂O ratio of liquid, thus prohibiting the crystallization of high-An plagioclase from hydrous tholeiite.     

Collision of comets with meteoroids

Statistical analysis of the quantity of dust in the cometary atmosphere in relation to the direction of motion of the comet about the Sun suggests an excess of dust for the retrograde comets. This excess is analyzed in the light of Harwit's theory of the cloud of “boulders” and of Öpik's impact theory. A comparison is also made between these excesses and other cometary phenomena such as splittings and outbursts.     

Polarization changes with time during solar microwave impulsive bursts

The evolution with time of circular polarization (t) from solar bursts at 7 GHz presents, in the majority of cases, a polarization degree peak before the maximum flux time. The subsequent evolution of (t) is continuous and usually increasing. The changes could be caused by superimposed polarization effects, due to the fast emissive electrons (dominant in the first phase), and to the propagation effects caused by the coronal condensation where the event occurred (dominant in the second phase). In an approximate approach, (t) is connected to the movement of the source in the second phase, being qualitatively sound, but limited to the lack of knowledge on acceleration processes and on magnetic field topology in the active region where the flares take place.     

Continuing observation of disturbed reproductive cycle and ovarian spermatogenesis in the giant abalone, Haliotis madaka from an organotin-contaminated site of Japan

Histological examination of gonads as well as chemical analysis of organotin compounds in tissues of the giant abalone, Haliotis madaka, was conducted to evaluate continuing endocrine disruption in abalone populations in Japan. Abalone specimens were collected from two different areas, Tsushima as a reference site and Jogashima as a site representative of declining abalone populations where serious organotin contamination had been observed, each month from January 1998 to March 1999. Scores were given to the development stages of reproductive cells in the ovary and testis. The degree of sexual maturation was evaluated by calculating the mean value of a histogram of these scores for the reproductive cells of each abalone. The temporal variations in degree of sexual maturation showed that female and male abalone from Tsushima matured synchronously, while those from Jogashima did not, which were similar to results of the previous study during September 1995–November 1996. Approximately 19% of the female abalone from Jogashima were masculinized with an ovo-testis, which was also similar to the result of the previous study. The masculinization of female abalone is reported to be induced by tributyltin (TBT) and triphenyltin (TPhT) from antifouling paints. Concentrations of the sum of butyltins (TBT, dibutyltin (DBT) and monobutyltin (MBT): ΣBTs) and the sum of phenyltins (TPhT, diphenyltin (DPhT) and monophenyltin (MPhT): ΣPhTs) in the muscle of abalone from Jogashima (n = 73) of 7.8 ± 9.0 ng/g wet wt. and 4.5 ± 6.8 ng/g wet wt., respectively, were significantly higher than those from Tsushima (n = 87) of 4.7 ± 4.9 ng/g wet wt. and 0.8 ± 1.7 ng/g wet wt., respectively (p < 0.05 for ΣBTs; p < 0.001 for ΣPhTs), although concentrations of TBT and TPhT in the muscle of abalone from Jogashima (n = 73) of 2.2 ± 2.5 ng/g wet wt. and 5.8 ± 5.1 ng/g wet wt., respectively, were insignificantly and significantly higher than those from Tsushima (n = 87) of 0.4 ± 0.6 ng/g wet wt. and 0.5 ± 0.9 ng/g wet wt., respectively, (p > 0.05 for TBT; p < 0.001 for TPhT). Thus, endocrine disruption as well as contamination by organotins in the giant abalone from Jogashima is still persisting.     

The light curve of the semiregular variable L2 Puppis – I. A recent dimming event from dust

The nearby Mira-like variable L₂ Pup is shown to be undergoing an unprecedented dimming episode. The stability of the period rules out intrinsic changes to the star, leaving dust formation along the line of sight as the most likely explanation. Episodic dust obscuration events are fairly common in carbon stars but have not been seen in oxygen-rich stars. We also present a 10-μm spectrum, taken with the Japanese Infrared Telescope in Space satellite, showing strong silicate emission that can be fitted with a detached, thin dust shell, containing silicates and corundum.     

The impulsiveness of microwave bursts and its association with sunspot types

A statistical analysis on 220 impulsive microwave bursts showed a marked preference to occur in sunspots of type D, E and F (after Zürich classification). Assuming a shock mechanism for accelerating electrons, this tendency seems to be related to the size of the magnetic loops.     

On the use of dimensioned measures of error to evaluate the performance of spatial interpolators

Spatial cross‐validation and average‐error statistics are examined with respect to their abilities to evaluate alternate spatial interpolation methods. A simple cross‐validation methodology is described, and the relative abilities of three, dimensioned error statistics—the root‐mean‐square error (RMSE), the mean absolute error (MAE), and the mean bias error (MBE)—to describe average interpolator performance are examined. To illustrate our points, climatologically averaged weather‐station temperatures were obtained from the Global Historical Climatology Network (GHCN), Version 2, and then alternately interpolated spatially (gridded) using two spatial‐interpolation procedures. Substantial differences in the performance of our two spatial interpolators are evident in maps of the cross‐validation error fields, in the average‐error statistics, as well as in estimated land‐surface‐average air temperatures that differ by more than 2°C. The RMSE and its square, the mean‐square error (MSE), are of particular interest, because they are the most widely reported average‐error measures, and they tend to be misleading. It (RMSE) is an inappropriate measure of average error because it is a function of three characteristics of a set of errors, rather than of one (the average error). Our findings indicate that MAE and MBE are natural measures of average error and that (unlike RMSE) they are unambiguous.     

Precise estimation of pressure–temperature paths from zoned minerals using Markov random field modeling: theory and synthetic inversion

The chemical zoning profile in metamorphic minerals is often used to deduce the pressure–temperature (P–T) history of rock. However, it remains difficult to restore detailed paths from zoned minerals because thermobarometric evaluation of metamorphic conditions involves several uncertainties, including measurement errors and geological noise. We propose a new stochastic framework for estimating precise P–T paths from a chemical zoning structure using the Markov random field (MRF) model, which is a type of Bayesian stochastic method that is often applied to image analysis. The continuity of pressure and temperature during mineral growth is incorporated by Gaussian Markov chains as prior probabilities in order to apply the MRF model to the P–T path inversion. The most probable P–T path can be obtained by maximizing the posterior probability of the sequential set of P and T given the observed compositions of zoned minerals. Synthetic P–T inversion tests were conducted in order to investigate the effectiveness and validity of the proposed model from zoned Mg–Fe–Ca garnet in the divariant KNCFMASH system. In the present study, the steepest descent method was implemented in order to maximize the posterior probability using the Markov chain Monte Carlo algorithm. The proposed method successfully reproduced the detailed shape of the synthetic P–T path by eliminating appropriately the statistical compositional noises without operator’s subjectivity and prior knowledge. It was also used to simultaneously evaluate the uncertainty of pressure, temperature, and mineral compositions for all measurement points. The MRF method may have potential to deal with several geological uncertainties, which cause cumbersome systematic errors, by its Bayesian approach and flexible formalism, so that it comprises potentially powerful tools for various inverse problems in petrology.     

Quasi-Static Strain and Stress Fields due to a Moment Tensor in Elastic–Viscoelastic Layered Half-Space

We derived explicit expressions in the time domain for 3-D quasi-static strain and stress fields, due to a point moment tensor source in an elastic surface layer overlying viscoelastic half-space under gravity. The expressions of strain in the elastic surface layer were directly obtained from the expressions of displacement in our previous paper. The conversion of strain into stress is easy, because the stress–strain relation of elastic material is linear. In the viscoelastic substratum, the expressions of strain were obtained by applying the correspondence principle of linear viscoelasticity to the associated elastic solution. The strain–stress conversion is not straightforward, as the stress–strain relation of viscoelastic material is usually given in a differential form. To convert strain into stress, we used an integral form of the stress–strain relation instead of the usual differential form. The expressions give the responses of elastic half-space at \( t = 0 \) , and the responses of an elastic plate floating on non-viscous liquid at \( t = \infty \) . The moment tensor is rationally decomposed into the three independent force systems, corresponding to isotropic expansion, shear faulting and crack opening, and so the expressions include the strain and stress fields for these force systems as special cases. As the first numerical example, we computed the temporal changes in strain and stress fields after the sudden opening of an infinitely long vertical crack cutting the elastic surface layer. Here, we observe that the stress changes caused by the sudden crack opening gradually decay with time and vanish at \( t = \infty \) everywhere. After the completion of stress relaxation, a characteristic pattern of shear strain remains in the viscoelastic substratum. Since the strain and stress fields at \( t = \infty \) can be read as the strain- and stress-rate fields caused by steady crack opening, respectively, this numerical example demonstrates the realization of a steady stress state supported by steady viscous flow in the asthenosphere, associated with steady seafloor spreading at mid-ocean ridges. For the second numerical example, we computed the temporal changes in strain and stress fields after the 2011 Tohoku-oki mega-thrust earthquake, which occurred at the North American-Pacific plate interface. In this numerical example, the stress changes caused by coseismic fault slip vanish at \( t = \infty \) in the viscoelastic substratum, but remain in the elastic surface layer. The coseismic stress changes (and also strain changes) in the elastic surface layer diffuse away from the source region with time, due to gradual stress relaxation in the viscoelastic substratum.