Osmium-strontium-neodymium-lead isotopic covariations in mid-ocean ridge basalt glasses and the heterogeneity of the upper mantle

Osmium, strontium, neodymium, and lead isotopic data have been obtained for 30 hand picked samples of basaltic glass from the Pacific, Atlantic and Indian mid-oceanic ridges. Large variations in Os isotopic ratios exist in the glasses, from abyssal peridotite-like values to radiogenic compositions similar to oceanic island basalts (¹⁸⁷Os/¹⁸⁶Os and ¹⁸⁷Os/¹⁸⁸Os ratios range from 1.06 to 1.36 and from 0.128 to 0.163, respectively). Os isotopic and elemental data suggest the existence of mixing correlations. This relationship might be ascribed to secondary contamination processes; however, such a hypothesis cannot account for the negative correlation observed between Os and Nd isotopes and the existence of complementary covariations between Os and SrPb isotopes. In this case, OsSrNdPb isotopic variations are unrelated to late post-eruption or shallow level contamination. These relationships provide strong evidence that the Os isotopic composition of the samples are derived from the mantle and thus implies a global chemical heterogeneity of the oceanic upper mantle. The results are consistent with the presence of recycled oceanic crust in the mantle sources of mid-ocean ridge basalts, and indicate that the unique composition of the upper mantle below the Indian ocean results from its contamination by a mantle component characterized by radiogenic Os and particularly unradiogenic Nd and Pb isotopic compositions.     

Simulation of normal distributed smooth fields by Karhunen-Loéve expansion in combination with kriging

Simulation of multigaussian stochastic fields can be made after a Karhunen-Loéve expansion of a given covariance function. This method is also called simulation by Empirical Orthogonal Functions. The simulations are made by drawing stochastic coefficients from a random generator. These numbers are multiplied with eigenfunctions and eigenvalues derived from the predefined covariance model. The number of eigenfunctions necessary to reproduce the stochastic process within a predefined variance error, turns out to be a cardinal question. Some ordinary analytical covariance functions are used to evaluate how quickly the series of eigenfunctions can be truncated. This analysis demonstrates extremely quick convergence to 99.5% of total variance for the 2nd order exponential (‘gaussian’) covariance function, while the opposite is true for the 1st order exponential covariance function. Due to these convergence characteristics, the Karhunen-Loéve method is most suitable for simulating smooth fields with ‘gaussian’ shaped covariance functions. Practical applications of Karhunen-Loéve simulations can be improved by spatial interpolation of the eigenfunctions. In this paper, we suggest interpolation by kriging and limits for reproduction of the predefined covariance functions are evaluated.     

Crystal morphologies in pillow basalts: implications for mid-ocean ridge processes

We present the results of a detailed petrological study of a sparsely phyric basalt (MAPCO CH98-DR11) dredged along the Mid-Atlantic Ridge (30°41′N). The sample contains microphenocrysts of olivine that display four different rapid-growth morphologies. Comparison of these morphologies with those obtained in dynamic crystallization experiments allows us to constrain the thermal history of the sample. The dendritic morphology (swallowtail, chain and lattice olivine) is directly related to the final quenching during magma–seawater interaction. In contrast, the three other morphologies, namely the complex polyhedral crystal, the closed hopper and the complex swallowtail morphology result from several cycles of cooling–heating (corresponding to a maximum degree of undercooling of 20–25°C) during crystal growth. These thermal variations occurred before eruption and are interpreted to be the result of turbulent convection in a small magmatic body beneath the ridge. The results suggest that the Mid-Atlantic Ridge is underlain by a mush zone that releases batches of liquid during tectonic segregation. Aphyric basalts are emitted during eruptions controlled by the tectonic activity, whereas phyric basalts correspond to small fractions of magma from the mush zone mobilized by reinjections of primitive magmas.     

Influence of the zone of weakness on dip angle and shear heating of subducted slabs

We study the importance of the zones of weakness and the pattern of downgoing flow in steady-state models of subducting lithosphere, which interacts mechanically and thermally with the ambient mantle. The non-linear system of governing equations consists of (i) the momentum equation in stream function formulation and (ii) the steady-state heat transfer equation including conduction and advection of heat and dissipation. A finite element method has been applied to this system. We consider the viscosity to be a non-linear function of both the temperature and the stream function. In steady-state two-dimensional (2D) flow, the stream function isolines follow material trajectories. They are used to follow the top of the subducting slab, which because of its possible increase in water content, is assumed to have a lower viscosity. The zone of weakness has been thus obtained in the self-consistent fashion since the stream function as well as the temperature are the output from our modeling and no a priori assumptions about the shape of the bending lithosphere are taken into account. It was shown that several orders decrease of viscosity in the zone of weakness is required to obtain the dip angle of about 45°. If the decrease of viscosity is not sufficient enough, the subducted slab either sinks almost vertically or does not exhibit a plate-like behavior. We have also demonstrated that shear heating can unrealistically increase at the zone of weakness for fast subductions if decrease of viscosity is underestimated.     

A Modeling Study of the Effects of Anomalous Snow Cover over the Tibetan Plateau upon the South Asian Summer Monsoon

The effect of anomalous snow cover over the Tibetan Plateau upon the South Asian summer monsoon is investigated by numerical simulations using the NCAR regional climate model (RegCM2) into which gravity wave drag has been introduced. The simulations adopt relatively realistic snow mass forcings based on Scanning Multi-channel Microwave Radiometer (SNINIR) pentad snow depth data. The physical mechanism and spatial structure of the sensitivity of the South Asian early summer monsoon to snow cover anomaly over the Tibetan Plateau are revealed. The main results are summarized as follows. The heavier than normal snow cover over the Plateau can obviously reduce the shortwave radiation absorbed by surface through the albedo effect, which is compensated by weaker upward sensible heat flux associated with colder surface temperature, whereas the effects of snow melting and evaporation are relatively smaller.The anomalies of surface heat fluxes can last until June and become unobvious in July. The decrease of the Plateau surface temperature caused by heavier snow cover reaches its maximum value from late April to early May. The atmospheric cooling in the mid-upper troposphere over the Plateau and its surrounding areas is most obvious in May and can keep a fairly strong intensity in June. In contrast, there is warming to the south of the Plateau in the mid-lower troposphere from April to June with a maximum value in May.The heavier snow cover over the Plateau can reduce the intensity of the South Asian summer monsoon and rainfall to some extent, but this influence is only obvious in early summer and almost disappears in later stages.     

Fitting the Linear Model of Coregionalization by Generalized Least Squares

In geostatistical studies, the fitting of the linear model of coregionalization (LMC) to direct and cross experimental semivariograms is usually performed with a weighted least-squares (WLS) procedure based on the number of pairs of observations at each lag. So far, no study has investigated the efficiency of other least-squares procedures, such as ordinary least squares (OLS), generalized least squares (GLS), and WLS with other weighing functions, in the context of the LMC. In this article, we compare the statistical properties of the sill estimators obtained with eight least-squares procedures for fitting the LMC: OLS, four WLS, and three GLS. The WLS procedures are based on approximations of the variance of semivariogram estimates at each distance lag. The GLS procedures use a variance–covariance matrix of semivariogram estimates that is (i) estimated using the fourth-order moments with sill estimates (GLS₁), (ii) calculated using the fourth-order moments with the theoretical sills (GLS₂), and (iii) based on an approximation using the correlation between semivariogram estimates in the case of spatial independence of the observations (GLS₃). The current algorithm for fitting the LMC by WLS while ensuring the positive semidefiniteness of sill matrix estimates is modified to include any least-squares procedure. A Monte Carlo study is performed for 16 scenarios corresponding to different combinations of the number of variables, number of spatial structures, values of ranges, and scale dependence of the correlations among variables. Simulation results show that the mean square error is accounted for mostly by the variance of the sill estimators instead of their squared bias. Overall, the estimated GLS₁ and theoretical GLS₂ are the most efficient, followed by the WLS procedure that is based on the number of pairs of observations and the average distance at each lag. On that basis, GLS₁ can be recommended for future studies using the LMC.     

Petrology and geochemistry of a diamondiferous lherzolite from the Premier diamond mine, South Africa

This paper reports on the petrology and geochemistry of a diamondiferous peridotite xenolith from the Premier diamond mine in South Africa.

The xenolith is altered with pervasive serpentinisation of olivine and orthopyroxene. Garnets are in an advanced state of kelyphitisation but partly fresh. Electron microprobe analyses of the garnets are consistent with a lherzolitic paragenesis (8.5 wt.% Cr₂O₃ and 6.6 wt.% CaO). The garnets show limited variation in trace element composition, with generally low concentrations of most trace elements, e.g. Y (<11 ppm), Zr (<18 ppm) and Sr (<0.5 ppm). Garnet rare earth element concentrations, when normalised against the C1 chondrite of McDonough and Sun (Chem. Geol. 120 (1995) 223), are characterised by a rare earth element pattern similar to garnet from fertile lherzolite.

All diamonds recovered are colourless. Most crystals are sharp-edged octahedra, some with minor development of the dodecahedral form. A number of crystals are twinned octahedral macles, while aggregates of two or more octahedra are also common. Mineral inclusions are rare. Where present they are predominantly small black rosettes believed to consist of sulfide. In one instance a polymineralic (presumably lherzolitic) assemblage of reddish garnet, green clinopyroxene and a colourless mineral is recognised.

Infrared analysis of the xenolith diamonds show nitrogen contents generally lower than 500 ppm and variable nitrogen aggregation state, from 20% to 80% of the ‘B’ form. When plotted on a nitrogen aggregation diagram a well defined trend of increasing nitrogen aggregation state with increasing nitrogen content is observed. Carbon isotopic compositions range from −3.6 ‰ to −1.3 ‰. These are broadly correlated with diamond nitrogen content as determined by infrared spectroscopy, with the most negative C-isotopic compositions correlating with the lowest nitrogen contents.

Xenolith mantle equilibration temperatures, calculated from nitrogen aggregation systematics as well as the Ni in garnet thermometer are on the order of 1100 to 1200 °C.

It is concluded that the xenolith is a fertile lherzolite, and that the lherzolitic character may have resulted from the total metasomatic overprinting of pre-existing harzburgite. Metasomatism occurred prior to, or accompanied, diamond growth.     

Modélisation morphodynamique cross-shore d'un estran vaseux

Numerical experiments were performed to simulate the profile evolution of an intertidal mudflat with a 1D cross-shore morphodynamical model. First, the hydrodynamical forcing is a cross-shore tidal current due to semi-diurnal variations of the free surface elevation at the open boundary. Further, considering the conservation of the action density of surface gravity waves, a wave height (and resulting bottom shear stress) calculation is added to the morphodynamical model. Results of the numerical experiments show that the shape of the profile reaches equilibrium. The mudflat progrades continually when the forcing is tide only, whereas it can be steady under the simultaneous action of tide and waves. To cite this article: B. Waeles et al., C. R. Geoscience 336 (2004).