Kimberlite petrogenesis: Insights from clinopyroxene-melt partitioning experiments at 6 GPa in the CaO-MgO-Al2O3-SiO2-CO2 system

In this experimental study, we examine the mineral-melt partitioning of major and trace elements between clinopyroxene and CO₂-rich kimberlitic melts at a pressure of 6 GPa and temperatures of 1410°C and 1430°C. The melts produced contain ∼ 28 wt% dissolved CO₂, and are saturated with olivine and clinopyroxene. To assess the effects of temperature, crystal and melt compositions on trace element partitioning, experiments were performed in the model CaO-MgO-Al₂O₃-SiO₂-CO₂ system. Our results reveal that all the elements studied, except Al, Mg, Si, and Ga, are incompatible in clinopyroxene. Partition coefficients show a considerable range in magnitude, from ∼ 10⁻³ for D_U and D_Ba to ∼ 2.5 for D_Si. The two experimental runs show similar overall partitioning patterns with the D values being lower at 1430°C. Rare earth elements display a wide range of partition coefficients, D_La (0.012-0.026) being approximately one order of magnitude lower than D_Lu (0.18-0.23). Partition coefficients for the 2+ and 3+ cations entering the M2-site exhibit a near-parabolic dependence on radius of the incorporated cations as predicted from the lattice strain model. This underlines the contribution made by the crystal structure toward controlling the distribution of trace elements. Using data obtained in this study combined with that in the published literature, we also discuss the effects that other important parameters, namely, melt composition, pressure, and temperature, could have on partitioning.Our partition coefficients have been used to model the generation of the Group I (GI) kimberlites from South Africa. The numerical modeling shows that kimberlitic melts can be produced by ∼0.5% melting of a MORB-type depleted source that has been enriched by small-degree melts originating from a similar depleted source. This result suggests that the source of GI kimberlites may be located at the lithosphere-asthenosphere transition. Percolation of small degree melts from the asthenosphere would essentially create a metasomatic horizon near the bottom of the non-convecting sublithospheric mantle. Accumulation of such small degree melts together with the presence of volatiles and conductive heating would trigger melting of the ambient mantle and subsequently lead to eruption of kimberlitic melts. Additionally, our model shows that the GI source can be generated by metasomatism of a 2 Ga old MORB source ca. 1 Ga ago. Assuming that MORB-type mantle is the most depleted source of magmas on earth, then this is the oldest age at which the GI source could have existed. However, this age most likely reflects the average age of a series of metasomatic events than that of a single event.     

The lunar rock and mineral characterization consortium: Deconstruction and integrated mineralogical,petrologic, and spectroscopic analyses of mare basalts

The lunar rock and mineral characterization consortium (LRMCC) has conducted coordinated mineralogy/petrography/spectroscopy analyses of a suite of pristine lunar basalts. Four basalt slabs (two low‐Ti, two high‐Ti) and paired thin sections were analyzed. Thin sections were analyzed for mineralogy/petrography, while the slabs were used to prepare particulate separates of major mineral phases and bulk samples. Mineral separates and particulate bulk samples were crushed to controlled grain sizes and their reflectance spectra measured in the NASA RELAB at Brown University. The resulting data set provides an essential foundation for spectral mixing models, offers valuable endmember constraints for space weathering analyses, and represents critical new ground truth results for lunar science and exploration efforts.     

Anatomy of individual agglutinates from a lunar highland soil

Abstract— We report results of our investigation of the relationship between values of I_s/FeO (relative concentration of nanophase Fe⁰ divided by total FeO content), glass abundance, total Fe content, and degree of digestion of <20 μm clasts for 22 individual agglutinates (250–1000 μm) from the mature Apollo 16 soil 61181 (I_s/FeO = 82 units in the <250 μm fraction). Agglutinates are important products of space weathering on the Moon, and they influence spectral observations at visible and near-IR wavelengths. Values of I_s/FeO for individual agglutinates (250–1000 μm) within this single soil span a range from 3 to 262 units which is larger than the range observed for all Apollo 16 bulk soils (～0 to 110 units). No correlation was observed between I_s/FeO and glass abundance and FeO concentrations for either agglutinitic glass or whole agglutinate particles under investigation. Our results suggest that the variation in I_s/FeO for agglutinates from a single soil may be in part a consequence of natural mixing processes on the Moon that produce highly-variable environments (with respect to surface exposure) for agglutinate formation and in part to variable kinetics of reactions in an agglutinate melt, which are influenced by a variety of factors including melt composition, temperature, impactor velocity, and quench rate. We cannot exclude but do not see evidence for other processes including addition of exotic agglutinates, micrometeoritic bombardment into compositionally-diverse microtargets, recycling of agglutinates, preferential melting of very fine soil particles, and production of nanophase Fe⁰ in amorphous rims of very fine irradiated lunar grains contributing to the observed variation of I_s/FeO.     

Self-regulated gravitational accretion in protostellar discs

We present a numerical model for the evolution of a protostellar disc that has formed self-consistently from the collapse of a molecular cloud core. The global evolution of the disc is followed for several million years after its formation. The capture of a wide range of spatial and temporal scales is made possible by use of the thin-disc approximation. We focus on the role of gravitational torques in transporting mass inward and angular momentum outward during different evolutionary phases of a protostellar disc with disc-to-star mass ratio of order 0.1. In the early phase, when the infall of matter from the surrounding envelope is substantial, mass is transported inward by the gravitational torques from spiral arms that are a manifestation of the envelope-induced gravitational instability in the disc. In the late phase, when the gas reservoir of the envelope is depleted, the distinct spiral structure is replaced by ongoing irregular non-axisymmetric density perturbations. The amplitude of these density perturbations decreases with time, though this process is moderated by swing amplification aided by the existence of the disc's sharp outer edge. Our global modelling of the protostellar disc reveals that there is typically a residual non-zero gravitational torque from these density perturbations, i.e. their effects do not exactly cancel out in each region. In particular, the net gravitational torque in the inner disc tends to be negative during first several million years of the evolution, while the outer disc has a net positive gravitational torque. Our global model of a self-consistently formed disc shows that it is also self-regulated in the late phase, so that it is near the Toomre stability limit, with a near-uniform Toomre parameter Q ≈ 1.5–2.0. Since the disc also has near-Keplerian rotation, and comparatively weak temperature variation, it maintains a near-power-law surface density profile proportional to r ^−3/2.     

Magnetically-regulated fragmentation induced by nonlinear flows and ambipolar diffusion

We present a parameter study of simulations of fragmentation regulated by gravity, magnetic fields, ambipolar diffusion, and nonlinear flows. The thin-sheet approximation is employed with periodic lateral boundary conditions, and the nonlinear flow field (“turbulence”) is allowed to freely decay. In agreement with previous results in the literature, our results show that the onset of runaway collapse (formation of the first star) in subcritical clouds is significantly accelerated by nonlinear flows in which a large-scale wave mode dominates the power spectrum. In addition, we find that a power spectrum with equal energy on all scales also accelerates collapse, but by a lesser amount. For a highly super-Alfvénic initial velocity field with most power on the largest scales, the collapse occurs promptly during the initial compression wave. However, for trans-Alfvénic perturbations, a subcritical magnetic field causes a rebound from the initial compression, and the system undergoes several oscillations before runaway collapse occurs. Models that undergo prompt collapse have highly supersonic infall motions at the core boundaries. Cores in magnetically subcritical models with trans-Alfvénic initial perturbations also pick up significant systematic speeds by inheriting motions associated with magnetically-driven oscillations. Core mass distributions are much broader than in models with small-amplitude initial perturbations, although the disturbed structure of cores that form due to nonlinear flows does not guarantee subsequent monolithic collapse. Our simulations also demonstrate that significant power (if present initially) can be maintained with negligible dissipation in large-scale compressive modes of a magnetic thin sheet, in the limit of perfect flux freezing.     

Five dimensional Bianchi type-V cosmological models described by a binary mixture of perfect fluid and dark energy with Λ-term

Some recent experimental observations have been shown that inclusion of electron collisions damping in inertial Alfvén wave (IAW) dynamics may be important for laboratory as well as space plasmas. This paper presents the numerical simulation of model equation governing the nonlinear dynamics of IAW in low-beta plasmas. When the nonlinearity arises due to the ponderomotive force and Joule heating driven density perturbations, the model equation turns out to be a modified nonlinear Schr?dinger equation (MNLS). The electron collisions are introduced only in the electron momentum equation. The damped localized structures of IAW with sidebands are obtained. Also, the effect of collisional damping on power spectra of magnetic fluctuations with different scaling laws has been studied. These turbulent structures may be responsible for particle acceleration in laboratory and space plasmas.     

Analysis of buried pipelines subjected to reverse fault motion

Presently available simplified analytical methods and semi-empirical methods for the analysis of buried pipelines subjected to fault motion are suitable only for the strike-slip and the normal-slip type fault motions, and cannot be used for the reverse fault crossing case. A simple finite element model, which uses beam elements for the pipeline and discrete nonlinear springs for the soil, has been proposed to analyse buried pipeline subjected to reverse fault motion. The material nonlinearities associated with pipe-material and soil, and geometric nonlinearity associated with large deformations were incorporated in the analysis. Complex reverse fault motion was simulated using suitable constraints between pipe-nodes and ground ends of the soil spring. Results of the parametric study suggest that the pipeline's capacity to accommodate reverse fault offset can be increased significantly by choosing a near-parallel orientation in plan with respect to the fault line. Further improvement in the response of the pipeline is possible by adopting loose backfill, smooth and hard surface coating, and shallow burial depth in the fault crossing region. For normal or near normal orientations, pipeline is expected to fail due to beam buckling at very small fault offsets.