Origin of the hydrocarbon component of carbonaceous chondrites: the star-meteorite connection

We report the results of an experiment that produced a residue which closely matches the hydrocarbon component of the Murchison carbonaceous chondrite. This experiment suggests that the parent material of the meteoritic component originated as polycyclic aromatic hydrocarbon species in carbon stars during their later stages of evolution. The experiments also indicate that the pathway from those formation sites to eventual incorporation into the meteorite parent body involved hydrogenation in a plasma in the solar nebula or in H II regions prior to the solar nebula. This model is consistent with what is known about the meteoritic hydrocarbon component including deuterium abundance, the observation of cosmic infrared emission bands best attributed to polycyclic aromatic hydrocarbon molecules, and the inherent stability of these molecules that allows their formation in stars and subsequent survival in the interstellar medium.     

Radiation pressure effects in the oscillations of compressible rotating homogeneous spheroids

Earlier models of compressible, rotating, and homogeneous ellipsoids with gas pressure are generalized to include the presence of radiation pressure. Under the assumptions of a linear velocity field of the fluid and a bounded ellipsoidal surface, the dynamical behaviour of these models can be described by ordinary differential equations. These equations are used to study the finite oscillations of massive radiative models with masses 10M and 30M in which the effects of radiation pressure are expected to be important.Models with two different degrees of equilibrium are chosen: an equilibrium (i.e., dynamically stable) model with an initial asymmetric inward velocity, and a nonequilibrium model with a nonequilibrium central temperature and which falls inwards from rest. For each of these two degrees of equilibrium, two initial configurations are considered: rotating spheroidal and nonrotating spherical models.From the numerical integration of the differential equations for these models, we obtain the time evolution of their principal semi-diametersa ₁ anda ₃, and of their central temperatures, which are graphically displayed by making plots of the trajectories in the (a ₁,a ₃) phase space, and of botha ₁ and the total central pressureP _c against time.It is found that in all the equilibrium radiative models (in which radiation pressure is taken into account), the periods of the oscillations of botha ₁ andP _c are longer than those of the corresponding nonradiative models, while the reverse is true for the nonequilibrium radiative models. The envelopes of thea ₁ oscillations of the equilibrium radiative models also have much longer periods; this result also holds for the nonequilibrium models whenever the envelope is well defined. Further, as compared to the nonradiative models, almost all the radiative models collapse to smaller volumes before rebouncing, with the more massive model undergoing a larger collapse and attaining a correspondingly larger peakP _c.When the mass is increased, the dynamical behavior of the radiative model generally becomes more nonperiodic. The ratio of the central radiation pressure to the central gas pressure, which is small for low mass models, increases with mass, and at the center of the more massive model, the radiation pressure can be comparable in magnitude to the gas pressure. In all the radiative models, the average periods as well as the average amplitudes of both thea ₁ andP _c oscillations also increase with mass.When either rotation or radiation pressure effects or both are included in the equilibrium nonradiative model, the period of the envelope of thea ₁ oscillations is increased. The presence of rotation in the equilibrium radiative model, however, decreases this period.Some astrophysical implications of this work are briefly discussed.     

Carbonatites from Eastern Paraguay and genetic relationships with potassic magmatism: C, O, Sr and Nd Isotopes

Summary Geochemical characteristics were systematically determined for Early Cretaceous samples of carbonatitic rocks from Eastern Paraguay (Rio Apa, Amambay and Central Provinces). The data show that all the occurrences have an enriched isotopic signature and that the carbonatites have negligible or absent crustal signature. A petrogenetic model (parent liquids, fractional crystallization, hydrothermal interactions and weathering) is proposed as a function of incompatible trace element, stable (O-C) and radiogenic (Sr-Nd) isotope variations with the aim to test the significance of carbonatitic complexes as a marker of the metasomatized subcontinental lithospheric mantle. The results indicate that the carbonatites and primary carbonates from eastern Paraguay, and those from the north eastern Paraná Basin (SE Brazil), were affected by metasomatic events distinct in time and composition.

---

Karbonatite aus Ost-Paraguay und ihre genetische Beziehung zu Kalium-Magmatismus: C O, Sr und Nd isotope

Zusammenfassung Die geochemischen Charakteristika von frühkretazischen Karbonatitproben aus Ostparaguay (Rio Alpa, Amambay und Zentrale Provinzen) wurden untersucht. Die Daten belegen, daß alle Vorkommen eine isotopische Anreicherungssignatur zeigen und daß ihnen eine entsprechende Krustensignatur fehlt. Ein Petrologisches Modell (Ausgangsschmelze, fraktionierte Kristallisation, hydrothermale Interaktion und Verwitterung) wird auf Grund der Verteilung der inkompatiblen Spurenelemente, der stabilen (C-O) und radiogenen (Sr-Nd) Isotope vorgeschlagen. Es versucht die Bedeutung der Karbonatitkomplexe als Markerhorizonte des metasomatischen subkontinentalen Mantels zu überprüfen. Die Ergebnisse zeigen, daß die Karbonatite und die primären Karbonate in Ostparaguay, und jene aus dem Paraná Becken SüdostBrasiliens durch zeitlich und zusammensetzungsmäßig unterschiedliche metasomatische Prozesse erfaßt wurden.

---

---

With 10 Figures     

Compressibility and crystal structure of sillimanite, Al2SiO5, at high pressure

The unit-cell dimensions and crystal structure of sillimanite at various pressures up to 5.29 GPa have been refined from single-crystal X-ray diffraction data. As pressure increases, a and b decrease linearly, whereas c decreases nonlinearly with a slightly positive curvature. The axial compression ratios at room pressure are β_a:β_b:β_c=1.22:1.63:1.00. Sillimanite exhibits the least compressibility along c, but the least thermal expansivity along a (Skinner et al. 1961; Winter and Ghose 1979). The bulk modulus of sillimanite is 171(1) GPa with K′=4 (3), larger than that of andalusite (151 GPa), but smaller than that of kyanite (193 GPa). The bulk moduli of the [Al1O₆], [Al2O₄], and [SiO₄] polyhedra are 162(8), 269(33), and 367(89) GPa, respectively. Comparison of high-pressure data for Al₂SiO₅ polymorphs reveals that the [SiO₄] tetrahedra are the most rigid units in all these polymorphic structures, whereas the [AlO₆] octahedra are most compressible. Furthermore, [AlO₆] octahedral compressibilities decrease from kyanite to sillimanite, to andalusite, the same order as their bulk moduli, suggesting that [AlO₆] octahedra control the compression of the Al₂SiO₅ polymorphs. The compression of the [Al1O₆] octahedron in sillimanite is anisotropic with the longest Al1-OD bond shortening by ～1.9% between room pressure and 5.29 GPa and the shortest Al1-OB bond by only 0.3%. The compression anisotropy of sillimanite is primarily a consequence of its topological anisotropy, coupled with the compression anisotropy of the Al-O bonds within the [Al1O₆] octahedron.     

Estimation for partially observed Markov processes

Many stochastic process models for environmental data sets assume a process of relatively simple structure which is in some sense partially observed. That is, there is an underlying process (X_n, n 0) or (X_t, t 0) for which the parameters are of interest and physically meaningful, and an observable process (Y_n, n 0) or (Y_t, t 0) which depends on the X process but not otherwise on those parameters. Examples are wide ranging: the Y process may be the X process with missing observations; the Y process may be the X process observed with a noise component; the X process might constitute a random environment for the Y process, as with hidden Markov models; the Y process might be a lower dimensional function or reduction of the X process. In principle, maximum likelihood estimation for the X process parameters can be carried out by some form of the EM algorithm applied to the Y process data. In the paper we review some current methods for exact and approximate maximum likelihood estimation. We illustrate some of the issues by considering how to estimate the parameters of a stochastic Nash cascade model for runoff. In the case of k reservoirs, the outputs of these reservoirs form a k dimensional vector Markov process, of which only the kth coordinate process is observed, usually at a discrete sample of time points.     

Crystallization processes in an artificial magma: variations in crystal shape,growth rate and composition with melt cooling history

A large (4.8 m³, 1.3x10⁷ g) artificial mafic melt with a bulk composition similar to that of a basalt (but with a high CaO content of 17 wt%) was generated during a demonstration of in situ vitrification and was allowed to cool naturally. During the melting process, convection was vigorous, resulting in a chemically and thermally homogeneous melt body. Once heating was complete, the cooling rate was rapid with the temperature dropping from 1500°C to 500°C in 6 days within the interior of the 3 m diameter, 1.5 m thick body. A 20 h period of constant temperature (1140°C) observed during colling was the result of latent heat released by widespread crystallization. The final crystalline assemblage consists of diopsidic to hedenbergitic pyroxene and anorthitic feldspar, with a subordinate amount of potassic feldspar, plus a small amount of evolved glass. The compositions and proportions of phases agree well with those predicted by the MELTS thermodynamic model. Thermal and textural evidence suggest that convection within the melt ceased coincident with formation of the first crystals. Textural investigation of core samples reveals large (up to 1 cm in length) acicular diopsidic pyroxenes in a matrix of smaller feldspar and zoned pyroxene crystals (500 m in length). Crystal shape and pyroxene composition vary as a function of position within the solidified body, as a function of cooling rate. Both crystal size and degree of crystallinity are highest in the central, most slowly-cooled parts of the rock. Crystal shape is characterized by tabular, equilibrium-growth forms in the slowly-cooled areas, grading to highly skeletal, dendritic forms at the rapidly-cooled edges of the body. The pyroxene crystals are dominantly homogeneous diopside, but crystals are characterized by thin Fe-rich hedenbergitic rims. These rims were deposited when Mg in the melt was depleted by diopside growth, and melt temperature had cooled sufficiently to allow Fe-rich pyroxene growth. Crystal growth rates can be calculated based on thermal behavior of the melt, reinforced by thermodynamic modelling, and are determined to be between 10^-7 and 10^-8 cm/s in the central part of the melt. These estimates agree well with growth rates in natural systems with similar cooling rates. Pyroxene crystals that formed at a higher cooling rates are characterized by higher Al and lower Mg contents relative to tabular equilibrium crystalline forms, presumably as a result of disequilibrium melt compositions at the crystal-melt interface.New Mexico Bureau of Mines and Mineral Resources, New Mexico Institute of Mining and Technology, Socorro, NM 87801, USA     

A study of bernalite,Fe(OH)3, using Mössbauer spectroscopy,optical spectroscopy and transmission electron microscopy

To study the crystal chemistry of bernalite, Fe(OH)₃, and the nature of the octahedral Fe³⁺ environment, Mössbauer spectra were recorded from 80 to 350 K, optical spectra were recorded at room temperature and a sample was studied using transmission electron microscopy. The Mössbauer spectrum of bernalite consists of a single six-line magnetic spectrum at 80 K. A broadened six-line magnetic spectrum with significantly less intensity is observed at higher temperatures, and is attributed to a small fraction of bernalite occurring as small particles. The variation of hyperfine magnetic field data for bulk bernalite with temperature is well described by the Weiss molecular field model with parameters of H ₀ = 55.7±0.3 T and T _N = 427±5K. The centre shift data were fitted to the Debye model with parameters ₀=0.482±0.005 mm/s (relative to -Fe) and _M=492±30 K. The quadrupole shift is near zero at 300 K, and does not vary significantly with temperature. Absorption spectra in the visible and near infrared range show three crystal field bands of Fe³⁺ at 11 300, 16000 and 23 200 cm^-1, giving a crystal field splitting of 14 570 cm^-1 and Racah parameters of B=629 cm^-1 and C=3381 cm^-1. Infrared reflection spectra show two distinct OH-stretching frequencies, which could correspond to two structurally different types of OH groups. A band was also observed at 2250 cm^-1, suggesting the presence of molecular CO₂ in the large cation site. Analytical transmission electron microscopy indicates that Si occurs within the bernalite structure as well as along domain boundaries. Electron diffraction and imaging show that bernalite is polysynthetically twinned along {100} planes with twin domains ranging from 3 to 20 nm in thickness. Results are discussed with respect to the nature of the octahedral Fe³⁺ site, and compared with values for other iron oxides and hydroxides.     

Episode 49 of the Pu'u 'O'o-Kupaianaha eruption of Kilauea volcano-breakdown of a steady-state eruptive era

The Pu'u 'O'o-Kupaianaha eruption (1983-present) is the longest lived rift eruption of either Kilauea or neighboring Mauna Loa in recorded history. The initial fissure opening in January 1983 was followed by three years of episodic fire fountaining at the Pu'u 'O'o vent on Kilauea's east rift zone 19km from the summit (episodes 4–47). These spectacular events gave way in July 1986 to five and a half years of nearcontinuous, low-level effusion from the Kupaianaha vent, 3km to the cast (episode 48). A 49th episode began in November 1991 with the opening of a new fissure between Pu'u 'O'o and Kupaianaha. this three week long outburst heralded an era of more erratic eruptive behavior characterized by the shut down of Kupaianaha in February 1992 and subsequent intermittent eruption from vents on the west flank of Pu'u 'O'o (episodes 50 and 51). The events occurring over this period are due to progressive shrinkage of the rift-zone reservoir beneath the eruption site, and had limited impact on eruption temperatures and lava composition.     

A search for ultra-high energy counterparts to gamma-ray bursts

A small air shower array operating over many years has been used to search for ultra-high energy (UHE) gamma radiation ( 50 TeV) associated with gamma-ray bursts (GRBs) detected by the BATSE instrument on the Compton Gamma-Ray Observatory (CGRO). Upper limits for a one minute interval after each burst are presented for seven GRBs located with zenith angles < 20°. A 4.3 excess over background was observed between 10 and 20 minutes following the onset of a GRB on 11 May 1991. The confidence level that this is due to a real effect and not a background fluctuation is 99.8%. If this effect is real then cosmological models are excluded for this burst because of absorption of UHE gamma rays by the intergalactic radiation fields.     

Search for low energyν e,μ,τ andv e,μ,τ in coincidence with BATSE gamma ray bursts

A search for low energy neutrinos of all flavours in correlation with 553 ray bursts detected by BATSE aboard the Compton Observatory has been performed by the LSD (Liquid Scintillator Detector) neutrino telescope. No excess ofe,, orv _e,, candidate has been detected by LSD during the time interval in which BATSE detected the 90% of the photon flux for any of the GRBs analyzed. Upper limits on the neutrino fluxes are given in the paper.