Sulfur and selenium systematics of the subcontinental lithospheric mantle: Inferences from the Massif Central xenolith suite (France)

Selenium has been analyzed in addition to S in 58 spinel peridotite xenoliths collected in Cenozoic alkali basalts from the Massif Central (France). The S concentration range now available for this suite, calculated from 123 samples, is the largest ever reported for alkali basalt-hosted xenoliths (<3-592 ppm). Likewise, the Se concentrations range between 0.2 and 67 ppb. No partial melting signature can be identified from the S and Se systematic. Half of the analyzed xenoliths have lost S during supergene weathering. By contrast, neither surficial alteration, nor loss of chalcophile elements during eruption can explain the regional-scale variations of S and Se concentrations. A first group of lherzolite xenoliths sampled in Southern Massif Central, from volcanic centers older and spatially unrelated to the Massif Central plume that triggered the Cenozoic volcanism, contains between 20 and 250 ppm S (with occasional S concentrations up to 592 ppm) and 12-67 ppb Se. It is clear that the highest S values, originally interpreted as representing S abundances in the primitive mantle, were in fact enriched by metasomatism. Highly variable S and Se contents (<5-360 ppm; 9-52 ppb) have also been observed in peridotite xenoliths collected in the Northern Massif Central, from volcanic centers mostly older than the plume. Like Group I xenoliths, these Group II xenoliths were strongly metasomatized by volatile-rich carbonated/silicated melts which precipitated Cu-rich sulfides. A third group of xenoliths from Plio-Quaternary basalts spatially related to the Massif Central Plume are uniformly poor in S (10-60 ppm) and Se (9-29 ppb). In this Group III, poikiloblastic textured xenoliths have lost most of their S and Se budget by peridotite-melt interactions at high melt/rock ratios. Taken as a whole, the Massif Central xenolith suite provides further evidence for strong heterogeneities in the S and Se budget of the sub-continental lithospheric mantle. However, the few LREE-depleted fertile lherzolites that escaped strong metasomatic alterations suggest a S- and Se-depleted primitive mantle reservoir compared to currently accepted primitive mantle estimates.     

Constraints from deuterium on the formation of icy bodies in the Jovian system and beyond

We consider the role of deuterium as a potential marker of location and ambient conditions during the formation of small bodies in our Solar system. We concentrate in particular on the formation of the regular icy satellites of Jupiter and the other giant planets, but include a discussion of the implications for the Trojan asteroids and the irregular satellites. We examine in detail the formation of regular planetary satellites within the paradigm of a circum-Jovian subnebula. Particular attention is paid to the two extreme potential subnebulae—“hot” and “cold”. In particular, we show that, for the case of the “hot” subnebula model, the D:H ratio in water ice measured from the regular satellites would be expected to be near-Solar. In contrast, satellites which formed in a “cold” subnebula would be expected to display a D:H ratio that is distinctly over-Solar. We then compare the results obtained with the enrichment regimes which could be expected for other families of icy small bodies in the outer Solar system—the Trojan asteroids and the irregular satellites. In doing so, we demonstrate how measurements by Laplace, the James Webb Space Telescope, HERSCHEL and ALMA will play an important role in determining the true formation locations and mechanisms of these objects.     

Numerical simulation of the propagation of P waves in fractured media

We study the propagation of P waves through media containing open fractures by performing numerical simulations. The important parameter in such problems is the ratio between crack length and incident wavelength. When the wavelength of the incident wavefield is close to or shorter than the crack length, the scattered waves are efficiently excited and the attenuation of the primary waves can be observed on synthetic seismograms. On the other hand, when the incident wavelength is greater than the crack length, we can simulate the anisotropic behaviour of fractured media resulting from the scattering of seismic waves by the cracks through the time delay of the arrival of the transmitted wave. The method of calculation used is a boundary element method in which the Green's functions are computed by the discrete wavenumber method. For simplicity, the 2-D elastodynamic diffraction problem is considered. The rock matrix is supposed to be elastic, isotropic and homogeneous, while the cracks are all empty and have the same length and strike direction. An iterative method of calculation of the diffracted wavefield is developed in the case where a large number of cracks are present in order to reduce the computation time. The attenuation factor Q ⁻¹ of the direct waves passing through a fractured zone is measured in several frequency bands. We observe that the attenuation factor Q ⁻¹ of the direct P wave peaks around kd = 2, where k is the incident wavenumber and d the crack length, and decreases proportionally to ( kd ) ⁻¹ in the high-wavenumber range. In the long-wavelength domain, the velocity of the direct P wave measured for two different crack realizations is very close to the value predicted by Hudson's theory on the overall elastic properties of fractured materials.     

Constraints in 3D gravity inversion

A three-dimensional (3D) inversion program is developed to interpret gravity data using a selection of constraints. This selection includes minimum distance, flatness, smoothness and compactness constraints, which can be combined using a Lagrangian formulation. A multigrid technique is also implemented to resolve separately large and short gravity wavelengths. The subsurface in the survey area is divided into rectangular prismatic blocks and the problem is solved by calculating the model parameters, i.e. the densities of each block. Weights are given to each block depending on depth, a priori information on density and the density range allowed for the region under investigation. The present computer code is tested on modelled data for a dipping dike and multiple bodies. Results combining different constraints and a weight depending on depth are shown for the dipping dike. The advantages and behaviour of each method are compared in the 3D reconstruction. Recovery of geometry (depth, size) and density distribution of the original model is dependent on the set of constraints used. From experimentation, the best combination of constraints for multiple bodies seems to be flatness and a minimum volume for the multiple bodies. The inversion method is tested on real gravity data from the Rouyn-Noranda (Quebec) mining camp. The 3D inversion model for the first 10 km is in agreement with the known major lithological contacts at the surface; it enables the determination of the geometry of plutons and intrusive rocks at depth.     

The catastrophe structure of thermohaline convection in a two-dimensional fluid model and a comparison with low-order box models

Abstract

We impose a surface forcing on the 2D, Boussinesq, thermohaline equations in a rectangular domain, in the form of equatorially symmetric cosine distributions of salinity flux and temperature. This system may be seen as an idealization of the ocean thermohaline circulation on the global scale over intervals of centuries or millenia. Multiple steady states are found numerically. They reflect the competition between the opposite signs of the temperature and salinity-driven equatorially symmetric circulations. There are also pole-to-pole, equatorially asymmetric circulations. In the control space of the temperature and salinity-flux forcing amplitudes, these equilibria form two cusp catastrophes, and transitions between stable equilibria occur through several distinct bifurcations. These catastrophes can be reproduced in simple box models connecting stirred reservoirs through capillary pipes. This steady-state analysis may provide a framework for a better understanding of climatic transitions between different stable regimes of the ocean-atmosphere system.     

An automated salt‐tracing gauge for flow‐velocity measurement

This article introduces the SVG (salt‐velocity gauge), a novel automated technique for measuring flow velocity by means of salt tracing. SVG allows a high measuring rate (up to one every 2 seconds), short control section length (down to 10 cm), high accuracy (+[sol ]?1·5 cm s^?1), and unbiased calculation of the mean velocity in experimental conditions with turbulent, supercritical flow. A few cubic centimetres of saturated salt solution (NaCl) are injected into the flow at regular time intervals using a programmable solenoid valve. The tracer successively passes two conductivity probes placed a short distance downstream. The transformation of the signal between the two probes is modelled as a one‐dimensional diffusion wave equation. Model calibration gives an estimation of the mean velocity and the diffusion for each salt plume. Two implementations of the SVG technique are described. The first was an outdoors simulated rainfall experiment in Senegal (conductivity probes at 40 cm apart, 8 Hz measurement rate, salt injections at 10 second intervals). Mean velocity was estimated to range between 0·1 and 0·3 m s^?1. The second was a laboratory‐based flume experiment (conductivity probes at 10 cm apart, 32 Hz, salt injections at 2 second intervals). Another SVG with probes at 34 cm apart was used for comparison. An acoustic Doppler velocimeter (ADV) was also used to give an independent assessment of velocity. Using the 10 cm salt gauge, estimated mean velocity ranged from 0·6 to 0·9 m s^?1 with a standard deviation of 1·5 cm s^?1. Comparisons between ADV, 10 cm SVG and 34 cm SVG were consistent and demonstrated that the salt‐tracing results were unbiased and independent of distance between probes. Most peaks were modelled with r² > 90 per cent. The SVG technology offers an alternative to the dye‐tracing technique, which has been severely criticized in the literature because of the wide interval of recommended values for the correction factor α to be applied to the timings. This article demonstrates that a fixed value of α is inappropriate, since the correction factor varies with velocity, diffusion and the length of the control section. Copyright © 2005 John Wiley & Sons, Ltd.     

Climate change impacts on the energy system: a review of trends and gaps

Major transformation of the global energy system is required for climate change mitigation. However, energy demand patterns and supply systems are themselves subject to climate change impacts. These impacts will variously help and hinder mitigation and adaptation efforts, so it is vital they are well understood and incorporated into models used to study energy system decarbonisation pathways. To assess the current state of understanding of this topic and identify research priorities, this paper critically reviews the literature on the impacts of climate change on the energy supply system, summarising the regional coverage of studies, trends in their results and sources of disagreement. We then examine the ways in which these impacts have been represented in integrated assessment models of the electricity or energy system.Studies tend to agree broadly on impacts for wind, solar and thermal power stations. Projections for impacts on hydropower and bioenergy resources are more varied. Key uncertainties and gaps remain due to the variation between climate projections, modelling limitations and the regional bias of research interests. Priorities for future research include the following: further regional impact studies for developing countries; studies examining impacts of the changing variability of renewable resources, extreme weather events and combined hazards; inclusion of multiple climate feedback mechanisms in IAMs, accounting for adaptation options and climate model uncertainty.     

Iron isotope fractionation in subterranean estuaries

Dissolved Fe concentrations in subterranean estuaries, like their river-seawater counterparts, are strongly controlled by non-conservative behavior during mixing of groundwater and seawater in coastal aquifers. Previous studies at a subterranean estuary of Waquoit Bay on Cape Cod, USA demonstrate extensive precipitation of groundwater-borne dissolved ferrous iron and subsequent accumulation of iron oxides onto subsurface sands. Waquoit Bay is thus an excellent natural laboratory to assess the mechanisms of Fe-isotope fractionation in redox-stratified environments and determine potential Fe-isotope signatures of groundwater sources to coastal seawater. Here, we report Fe isotope compositions of iron-coated sands and porewaters beneath the intertidal zone of Waquoit Bay. The distribution of pore water Fe shows two distinct sources of Fe: one residing in the upward rising plume of Fe-rich groundwater and the second in the salt-wedge zone of pore water. The groundwater source has high Fe(II) concentration consistent with anoxic conditions and yield δ⁵⁶Fe values between 0.3 and −1.3‰. In contrast, sediment porewaters occurring in the mixing zone of the subterranean estuary have very low δ⁵⁶Fe values down to −5‰. These low δ⁵⁶Fe values reflect Fe-redox cycling and result from the preferential retention of heavy Fe-isotopes onto newly formed Fe-oxyhydroxides. Analysis of Fe-oxides precipitated onto subsurface sands in two cores from the subterranean estuary revealed strong δ⁵⁶Fe and Fe concentration gradients over less than 2m, yielding an overall range of δ⁵⁶Fe values between −2 and 1.5‰. The relationship between Fe concentration and δ⁵⁶Fe of Fe-rich sands can be modeled by the progressive precipitation of Fe-oxides along fluid flow through the subterranean estuary. These results demonstrate that large-scale Fe isotope fractionation (up to 5‰) can occur in subterranean estuaries, which could lead to coastal seawater characterized by very low δ⁵⁶Fe values relative to river values.     

Organisation hydrochimique,biologique et sédimentologique d’un lac intracontinental à peuplements lagunaires : la Birket Karoun (Fayoum,Egypte)

Résumé

La Birket Karoun a été étudiée d’un point de vue biogéologique sur une période de 12 mois. Les données ainsi obtenues, ajoutées à celles de la littérature disponible, permettent de présenter les caractéristiques biogéologiques essentielles de ce bassin original.

Le lac n’est alimenté actuellement que par des eaux continentales issues du système nilotique ce qui entraîne des propriétés hydrochimiques assez différentes des corps d’eau paraliques influencés par la mer. Cependant, les caractéristiques hydrologiques, biologiques et sédimen-taires de la Birket Karoun la classent parmi les écosystèmes de type lagunaire. Il est même possible d’y reconnaître une organisation biogéologique zonale très semblable à celle des lagunes, la partie la plus « marine » (ou la moins confinée) se situant à l’opposé de la zone des apports d’eau douce. Une connexion entre le lac et la mer, actuelle ou passée, étant à peu près exclues, il faut admettre que les eaux continentales concentrées peuvent présenter, dans certains cas, des caractéristiques chimiques et biochimiques presque marines vis-à-vis des êtres vivants.

Dans cette optique, la Birket Karoun apparaît comme un modèle actuel susceptible de fournir une alternative à l’hypothèse d’une « Mer Saharienne » pour rendre compte de l’existence des lacs sahariens pleistocènes à faunes fossiles lagunaires.