Apollo 16 neutron stratigraphy

The Apollo 16 soils have the largest low energy neutron fluences (up to 10¹⁷ n/cm², E < 0.18eV) yet observed in lunar samples. Variations in the isotopic ratios ¹⁵⁸Gd/¹⁵⁷Gd and ¹⁵⁰Sm/¹⁴⁹Sm (up to 1.9% and 2.0% respectively) indicate that the low energy neutron fluence in the Apollo 16 drill stem increases with depth throughout the section sampled. Such a variation implies that accretion has been the dominant regolith “gardening” process at this location. The data may be fit by a model of continuous accretion of pre-irradiated material at a rate of ～70 g/(cm² · 10⁸yr) or by models involving as few as two slabs of material in which the first slab could have been deposited as long as 10⁹ yr ago.The ratio of the number of neutrons captured per atom by Sm to the number captured per atom by Gd is lower than in previously measured lunar samples, which implies a lower energy neutron spectrum at this site. The variation of this ratio with chemical composition is qualitatively similar to that predicted by Lingenfelter, Canfield and Hampel.Variations are observed in the ratio ¹⁵²Gd/¹⁶⁰Gd which are fluence correlated and probably result from neutron capture by¹⁵¹Eu.     

Ab initio study of MgSiO3 C2/c enstatite

We have used a newly developed ab initio constant-pressure molecular dynamics with variable cell shape technique to investigate the zero temperature behaviour of high pressure clinoenstatite (MgSiO₃-C2/c) from 0 up to 30 GPa. The optimum structure at 8 GPa, as well as structural trends under pressure, compare very well with experimental data. At this pressure, we find noticeable “fluctuations” in the chain configuration which suggests the structure is on the verge of a mechanical instability. Two distinct compressive behaviours then appear: one below and another above 8 GPa. This phenomenon may be related to the observed transition to a lower symmetry P2₁/c phase which involves a reconfiguration of the silicate chains, and suggests that the C2/c structure at low pressures found here, may be an artifact of the dynamical algorithm which preserves space group in the absence of symmetry breaking fluctuations. Comparison with calculations in other magnesium silicate phases, indicates that the size and shape of the silicate units (isolated and/or linked tetrahedra and octahedra) are generally well described by the local density approximation; however, the weaker linkages provided by the O-Mg-O bonds, are not as well described. This trend suggests that, as in the recently studied case of H₂O-ice, the structural properties of more inhomogeneous systems, like enstatite, may be improved by using gradientcorrected density functionals.     

Intercomparison of coastal area morphodynamic models

Five different coastal area morphodynamic models have been set up to run on the same offshore breakwater layout and an intercomparison carried out on the hydrodynamic and morphodynamic output produced by each scheme. In addition, the predicted morphodynamics was checked against available laboratory and field data.It is concluded that the models are capable of producing realistic estimates for the dominant morphodynamic features associated with offshore breakwaters. Coupling of the wave, current and sediment transport components of each scheme is shown to yield bathymetry which attains a state of equilibrium, unlike models which are based on the initial transport field only.     

Small-scale heterogeneity and soil-moisture variability in the unsaturated zone

The distribution of water within a soil profile can only be partly explained by the time distribution and rate of surface-water input. Observed differences in soil moisture within the unsaturated zone result from the interaction of surface-water inputs with spatially inhomogeneous soil characteristics. Water which initially percolates vertically is differentially impeded as a result of subtle textural changes in the soil, and is then preferentially retained in such zones of transition, causing large differences in soil-water content to occur. The scale of this vertical variability is of tenths of metres, whilst lateral variability of soil moisture reflects textural changes over a few metres. The observed influence of small-scale heterogeneity on soil-water content suggests that the conventional assumptions of isotropicity and homogeneity of the textural and hydraulic properties of porous media used in drainage basin and hillslope hydrological models need scrutiny, even for single stratigraphic units.     

Thermal expansion and crystal structure of FeSi between 4 and 1173 K determined by time-of-flight neutron powder diffraction

The thermal expansion and crystal structure of FeSi has been determined by neutron powder diffraction between 4 and 1173?K. No evidence was seen of any structural or magnetic transitions at low temperatures. The average volumetric thermal expansion coefficient above room temperature was found to be 4.85(5)?×?10^?5?K^?1. The cell volume was fitted over the complete temperature range using Grüneisen approximations to the zero pressure equation of state, with the internal energy calculated via a Debye model; a Grüneisen second-order approximation gave the following parameters: θ_D=445(11)?K, V ₀=89.596(8)?ų, K ₀′=4.4(4) and γ′=2.33(3), where θ_D is the Debye temperature, V ₀ is V at T=0?K, K ₀′ is the first derivative with respect to pressure of the incompressibility and γ′ is a Grüneisen parameter. The thermodynamic Grüneisen parameter, γ_th, has been calculated from experimental data in the range 4–400?K. The crystal structure was found to be almost invariant with temperature. The thermal vibrations of the Fe atoms are almost isotropic at all temperatures; those of the Si atoms become more anisotropic as the temperature increases.     

Computer simulations of the structural and physical properties of the olivine and spinel polymorphs of Mg2SiO4

The aim of the work presented is to develop a computer simulation technique which will predict the structure and physical properties of forsterite and ringwoodite, the major mantle-forming polymorphs of Mg₂SiO₄. The technique is based upon energy minimization, in which all structural parameters are varied until the configuration with the lowest energy is achieved. The lattice energy and physical properties (e.g. elasticity and dielectric constants) are calculated from interatomic potentials, which generally include electrostatic and short-range terms. We investigate several types of traditional potential models, and present a new type of model which includes partial ionic charges and a Morse potential to describe the effect of covalency on the Si-O bond. This new form of potential model is highly successful, and not only reproduces the zero-pressure structural, elastic and dielectric properties of forsterite and ringwoodite, but also accurately describes their pressure dependence.     

Palaeomonsoon and palaeoproductivity records of δ18O, δ{13}C and CaCO3 variations in the northern Indian Ocean sedimentsC and CaCO3 variations in the northern Indian Ocean sediments

δ¹⁸ O and δ¹³C of G.sacculifer have been measured in five cores from the northern Indian Ocean. In addition, high resolution analysis (1 to 2 cm) was performed on one core (SK-20-185) for both δ¹⁸O and gd¹³C in five species of planktonic foraminifera. CaCO₃ variation was measured in two cores. The results, presented here, show that

The lattice dynamics and thermodynamics of the Mg2SiO4 polymorphs

We use an approach based upon the Born model of solids, in which potential functions represent the interactions between atoms in a structure, to calculate the phonon dispersion of forsterite and the lattice dynamical behaviour of the beta-phase and spinel polymorphs of Mg₂SiO₄. The potential used (THB1) was derived largely empirically using data from simple binary oxides, and has previously been successfully used to model the infrared and Raman behaviour of forsterite. It includes ‘bond bending’ terms, that model the directionality of the Si-O bond, in addition to the pair-wise additive Coulombic and short range terms. The phonon dispersion relationships of the Mg₂SiO₄ polymorphs predicted by THB1 were used to calculate the heat capacities, entropies, thermal expansion coefficients and Gruneisen parameters of these phases. The predicted heat capacities and entropies are in outstandingly good agreement with those determined experimentally. The predicted thermodynamic data of these phases were used to construct a phase diagram for this system, which has Clausius-Clapeyron slopes in very close agreement with those found by experiment, but which has predicted transformation pressures that show less close agreement with those inferred from experiment. The overall success, however, that we have in predicting the lattice dynamical and thermodynamic properties of the Mg₂SiO₄ polymorphs shows that our potential THB1 represents a significant step towards finding the elusive quantitative link between the microscopic or atomistic behaviour of minerals and their macroscopic properties.     

Petrogenesis of the Greenhills Complex,Southland, New Zealand: magmatic differentiation and cumulate formation at the roots of a Permian island-arc volcano

A Permian (~265 Ma) intrusive complex which formed as a magmatic feeder reservoir to an immature island-arc volcano is fortuitously exposed in southern New Zealand. Known as the Greenhills Complex, this intrusion was emplaced at shallow crustal levels and consists of two layered bodies which were later intruded by a variety of dykes. Cumulates, which include dunite, olivine clinopyroxenite, olivine gabbro, and hornblende gabbro-norite, are related products of parent-magma fractionation. Both primary (magmatic) and secondary platinum-group minerals occur within dunite at one locality. Using the composition of cumulus minerals, mafic dykes and melt inclusions, we have determined that the parent magmas of the complex were hydrous, low-K island-arc tholeiites of ankaramitic affinities. Progressive magmatic differentiation of this parent magma generated fractionated melt of high-alumina basalt composition which is now preserved only as dykes which cut the Complex. Field evidence and cumulus mineral profiles reveal that the magma chambers experienced turbulent magmatic conditions during cumulate-rock formation. Recharge of the chambers by primitive magma is likely to have coincided with eruption of residual melt at the surface. Similar processes are inferred to account for volcanic-rock compositions in other parts of this arc terrane and in modern island-arc systems.