A search for radio pulsars around low-mass white dwarfs

Low-mass white dwarfs can be produced either in low-mass X-ray binaries by stable mass transfer to a neutron star, or in a common envelope phase with a heavier white dwarf companion. We have searched eight low-mass white dwarf candidates recently identified in the Sloan Digital Sky Survey for radio pulsations from pulsar companions, using the Green Bank Telescope at 340 MHz. We have found no pulsations down to flux densities of 0.6–0.8 mJy kpc⁻² and conclude that a given low-mass helium-core white dwarf has a probability of  <0.18 ± 0.05  of being in a binary with a radio pulsar.     

Two-dimensional modelling of subduction zone anisotropy with application to southwestern Japan

We present a series of 2-D numerical models of viscous flow in the mantle wedge induced by a subducting lithospheric plate. We use a kinematically defined slab geometry approximating the subduction of the Philippine Sea plate beneath Eurasia. Through finite element modelling we explore the effects of different rheological and thermal constraints (e.g. a low-viscosity region in the wedge corner, power law versus Newtonian rheology, the inclusion of thermal buoyancy forces and a temperature-dependent viscosity law) on the velocity and finite strain field in the mantle wedge. From the numerical flow models we construct models of anisotropy in the wedge by calculating the evolution of the finite strain ellipse and combining its geometry with appropriate elastic constants for effective transversely isotropic mantle material. We then predict shear wave splitting for stations located above the model domain using expressions derived from anisotropic perturbation theory, and compare the predictions to ∼500 previously published shear wave splitting measurements from seventeen stations of the broad-band F-net array located in southwestern Japan. Although the use of different model parameters can have a substantial effect on the character of the finite strain field, the effect on the average predicted splitting parameters is small. However, the variations with backazimuth and ray parameter of individual splitting intensity measurements at a given station for different models are often different, and rigorous analysis of details in the splitting patterns allows us to discriminate among different rheological models for flow in the mantle wedge. The splitting observed in southwestern Japan agrees well with the predictions of trench-perpendicular flow in the mantle wedge along with B-type olivine fabric dominating in a region from the wedge corner to about 125 km from the trench.     

Quantifying garnet-melt trace element partitioning using lattice-strain theory: assessment of statistically significant controls and a new predictive model

As a complement to our efforts to update and revise the thermodynamic basis for predicting garnet-melt trace element partitioning using lattice-strain theory (van Westrenen and Draper in Contrib Mineral Petrol, this issue), we have performed detailed statistical evaluations of possible correlations between intensive and extensive variables and experimentally determined garnet-melt partitioning values for trivalent cations (rare earth elements, Y, and Sc) entering the dodecahedral garnet X-site. We applied these evaluations to a database containing over 300 partition coefficient determinations, compiled both from literature values and from our own work designed in part to expand that database. Available data include partitioning measurements in ultramafic to basaltic to intermediate bulk compositions, and recent studies in Fe-rich systems relevant to extraterrestrial petrogenesis, at pressures sufficiently high such that a significant component of majorite, the high-pressure form of garnet, is present. Through the application of lattice-strain theory, we obtained best-fit values for the ideal ionic radius of the dodecahedral garnet X-site, r ₀(3+), its apparent Young’s modulus E(3+), and the strain-free partition coefficient D ₀(3+) for a fictive REE element J of ionic radius r ₀(3+). Resulting values of E, D ₀, and r ₀ were used in multiple linear regressions involving sixteen variables that reflect the possible influence of garnet composition and stoichiometry, melt composition and structure, major-element partitioning, pressure, and temperature. We find no statistically significant correlations between fitted r ₀ and E values and any combination of variables. However, a highly robust correlation between fitted D ₀ and garnet-melt Fe–Mg exchange and D _Mg is identified. The identification of more explicit melt-compositional influence is a first for this type of predictive modeling. We combine this statistically-derived expression for predicting D ₀ with the new expressions for predicting E and r ₀ outlined in the first of our pair of companion papers into an updated set of formulae that use easy-to-measure quantities (e.g. garnet composition, pressure, temperature) to predict variations in E, r ₀, and D ₀. These values are used in turn to calculate D values for those garnets. The updated model substantially improves upon a previous model (van Westrenen et al. in Contrib Mineral Petrol 142:219–234, 2001), and accounts well for trivalent cation partitioning in nominally anhydrous systems up to at least 15 GPa, including for eclogitic bulk compositions and for Fe-rich systems appropriate to magmagenesis on the Moon and Mars. The new model is slightly less successful in predicting partitioning with strongly majoritic garnets, although the mismatch is much less than with the original 2001 model. Although it also improves upon the 2001 model in predicting partitioning in hydrous systems, the mismatch between model and observation is still unacceptably large. The same statistical tools were applied in an attempt to predict tetravalent partitioning as well, because lattice-strain based techniques are not applicable to such partitioning. However, no statistically significant predictive relationships emerged from that effort. Our analyses show that future efforts should focus on filling the gap in partitioning data between ∼10 and 25 GPa to evaluate more closely the gradual transition of garnet to majorite, and on systematically expanding the hydrous partitioning database to allow extension of our model to water-bearing systems.     

Effects of orientation on the diffusive properties of fluid-filled grain boundaries during pressure solution

An unresolved issue in the study of pressure solution in rock materials is the dependence of grain boundary structure and diffusive properties on the mutual orientation of neighbouring grain lattices. We report electrical measurements yielding the diffusivity of differently oriented halite–glass and halite–halite contacts loaded in the presence of brine. The halite–glass contact experiments show pressure solution of the halite and an effect of halite lattice orientation on grain boundary transport. Post-mortem observations show an orientation-dependent grain boundary texture controlled by the periodic bond chains in the halite structure. It is inferred that this texture determines the internal grain boundary structure and properties during pressure solution. In the halite–halite experiments neck-growth occurred, its rate depending on twist-misorientation. The results imply that deformation by pressure solution may lead to lattice-preferred orientation development, and that polymineralic rocks may deform faster at lower stresses than monomineralic rocks.     

First evidence for the Cenomanian–Turonian oceanic anoxic event (OAE2, ‘Bonarelli’ event) from the Ionian Zone,western continental Greece

Integrated biostratigraphic (planktonic foraminifera, calcareous nannofossils), chemostratigraphic (bulk C and O isotopes) and compound-specific organic geochemical studies of a mid-Cretaceous pelagic carbonate—black shale succession of the Ionian Zone (western Greece), provide the first evidence for the Cenomanian–Turonian oceanic anoxic event (OAE2, ‘Bonarelli’ event) in mainland Greece. The event is manifested by the occurrence of a relatively thin (35 cm), yet exceptionally organic carbon-rich (44.5 wt% TOC), carbonate-free black shale, near the Cenomanian–Turonian boundary within the Vigla limestone formation (Berriasian–Turonian). Compared to the ‘Bonarelli’ black-shale interval from the type locality of OAE2 in Marche–Umbria, Italy, this black shale exhibits greatly reduced stratigraphic thickness, coupled with a considerable relative enrichment in TOC. Isotopically, enriched δ¹³C values for both bulk organic matter (−22.2‰) and specific organic compounds are up to 5‰ higher than those of underlying organic-rich strata of the Aptian-lower Albian Vigla Shale member, and thus compare very well with similar values of Cenomanian–Turonian black shale occurrences elsewhere. The relative predominance of bacterial hopanoids in the saturated, apolar lipid fraction of the OAE2 black shale of the Ionian Zone supports recent findings suggesting the abundance of N₂-fixing cyanobacteria in Cretaceous oceans during the Cenomanian–Turonian and early Aptian oceanic anoxic events.     

Problems in predicting the mobility of slow-moving landslides

This paper discusses problems in predicting the mobility of slow-moving landslides. Three case studies are presented here where research has been carried out by the Utrecht University: The La Valette landslide complex in the French Alps, the La Mure landslide in the French pre-Alps near Grenoble and the Hau landslide in Switzerland.To predict field velocities of these slow-moving landslides the viscosity parameters of the material of these landslides were determined by strain-controlled tests in a ring shear apparatus based on Bishop's design at Utrecht University. The viscosity parameters from the laboratory proved to be 10 to 1000 times lower than viscosities obtained from back analyses on the observed velocities in the field. This discrepancy may be explained by the development of negative pore pressures when the plastic material slides over a rigid, wavy slip surface and/or by convergent flow effects. The associated gain in strength results in a higher apparent viscosity.A more detailed analysis is made for the movements of the La Valette landslide. Observed velocities at the La Valette landslide are difficult to describe by one parameter set as the response to a change in groundwater level is not the same during a rise or fall in the piezometric level. This deviation may be explained by rapid changes in total stresses and consequently changes in pore pressure under (partly) undrained conditions. The emerging hysteresis with local and temporal variations in pore pressure makes it difficult to predict in detail the moving pattern of landslides.     

On the spreading mechanism of shallow groundwater in the Hinterland of the Dutch Dune hill area

The low lying Western part of the Netherlands is protected from the sea by a 5 km wide stretch of dunes rising to some 50 m of height. The fresh water pocket in the dunes overlies saline groundwater and a brackish transition zone. There was during a century an extraction of fresh groundwater for drinking water, supported by artificial infiltration. This has been stopped some 30 years ago. The consequent wetting of the valuable farm area (flower cultures) behind the dunes is stronger and more extensive than could be expected from mere replenishment of the fresh water zones in the dunes. It is shown in this paper that the lateral shear flows in the brackish and saline groundwater area have displaced (and are displacing) the interfaces vertically downward. The effect of more fresh and less saline groundwater in an arbitrary groundwater column is an (extra) rise of the groundwater head of the upper fresh water part. The described slow process will continue for decades until a new equilibrium has been established. In the mean time the inner dune farm area will have to cope with a surprisingly strong and extensive water level rise.