High-precision Nd/Nd measurements in terrestrial rocks: Constraints on the early differentiation of the Earth’s mantle

We present new ultra-high precision ¹⁴²Nd/¹⁴⁴Nd measurements of early Archaean rocks using the new generation thermal ionization mass spectrometer Triton. Repeated measurements of the Ames Nd standard demonstrate that the ¹⁴²Nd/¹⁴⁴Nd ratio can be determined with external precision of 2 ppm (2σ), allowing confident resolution of anomalies as small as 5 ppm. A major analytical improvement lies in the elimination of the double normalization procedure required to correct our former measurements from a secondary mass fractionation effect. Our new results indicate that metasediments, metabasalts, and orthogneisses from the 3.6 to 3.8 Ga West Greenland craton display positive ¹⁴²Nd anomalies ranging from 8 to 15 ppm. Using a simple two-stage model with an initial ε¹⁴³Nd value of 1.9 ± 0.6 ε-units, coupled ¹⁴⁷Sm-¹⁴³Nd and ¹⁴⁶Sm-¹⁴²Nd chronometry constrains mantle differentiation to 50-200 Ma after formation of the solar system. This chronological constraint is consistent with differentiation of the Earth’s mantle during the late stage of crystallization of a magma ocean. We have developed a two-box model describing ¹⁴²Nd and ¹⁴³Nd isotopic evolution of depleted mantle during the subsequent evolution of the crust-mantle system. Our results indicate that early terrestrial protocrust had a lifetime of ca. 0.7-1 Ga in order to produce the observed Nd isotope signature of Archaean rocks. In the context of this two box mantle-crust system, we model the evolution of isotopic and chemical heterogeneity of depleted mantle as a function of the mantle stirring time. Using the dispersion of ¹⁴²Nd/¹⁴⁴Nd and ¹⁴³Nd/¹⁴⁴Nd ratios observed in early Archaean rocks, we constrain the stirring time of early Earth’s mantle to 100-250 Ma, a factor of 5 shorter than the stirring time inferred from modern oceanic basalts.     

Study of the spatial and historical distribution of sediment inorganic contamination in the Toulon bay (France)

This study presents trace elements levels in surface and deep sediments of the Toulon bay (SE France) subjected to anthropogenic inputs (navy base, harbors, etc.). The studied elements (As, Cd, Cr, Cu, Hg, Ni, Pb, Zn) are defined as priority contaminants in aquatic systems. Fifty-five points scattered on the entire bay were sampled, allowing the determination of contaminants distribution with a high resolution. Several approaches were used to assess the degree of contamination and the potential toxicity of the Toulon bay sediments: comparison to the French legislation, surface-weighted average metal concentrations, enrichment factors (EF), geoaccumulation indices (Igeo), trace element stock calculation and comparison to sediment quality guidelines. A principal component analysis was performed to reveal common behavior of the studied contaminants. Results demonstrated the very high contamination of the small bay, especially in Hg (EF up to 1500), Cu, Pb and Zn, with export to the large bay further governed by hydrodynamics.     

Optimal removal of topological artefacts in microtomographic images of porous materials

The analysis of flow at the pore scale in porous media has been facilitated with the use of microtomography. A powerful tool for quantifying the fluid structure using these tomographic 3D reconstructions is skeletonisation, but the significant disadvantage of this method is its sensitivity to noise, resulting in artefacts in the skeleton. A pre-processing of the 3D image is therefore required, but no method has yet proven to completely solve this problem. By developing a new procedure that, by construction, directly identifies the voxels and only those that are responsible for topological artefacts in the skeleton, we are able to remove all artefacts, and furthermore can prove that we do so by modifying a minimal amount of voxels in the segmented 3D image (i.e. the tomographic image in which each voxel has been assigned to either the porous or the solid phase). This is possible by identifying the three fundamental types of artefacts that can arise in a 3D skeleton, and dealing with each appropriately. Application to a microtomographic image of a sintered glass powder is presented. Impact of the different processing methods on the flow within its porosity is measured through the computed permeability deviations.     

Numerical models of the thermomechanical evolution of planetesimals: Application to the acapulcoite‐lodranite parent body

The acapulcoite‐lodranite meteorites are members of the primitive achondrite class. The observation of partial melting and resulting partial removal of Fe‐FeS indicates that this meteorite group could be an important link between achondrite and iron meteorites, on the one hand, and chondrite meteorites, on the other. Thus, a better understanding of the thermomechanical evolution of the parent body of this meteorite group can help improve our understanding of the evolution of early planetesimals. Here, we use 2‐D and 3‐D finite‐difference numerical models to determine the formation time, initial radius of the parent body of the acapulcoite‐lodranite meteorites, and their formation depth inside the body by applying available geochronological, thermal, and textural constraints to our numerical data. Our results indicate that the best fit to the data can be obtained for a parent body with 25–65 km radius, which formed around 1.3 Ma after calcium‐aluminum‐rich inclusions. The 2‐D and 3‐D results considering various initial temperatures and the effect of porosity indicate possible formation depths of the acapulcoite‐lodranite meteorites of 9–19 and 14–25 km, respectively. Our data also suggest that other meteorite classes could form at different depths inside the same parent body, supporting recently proposed models (Elkins‐Tanton et al. 2011 ; Weiss and Elkins‐Tanton 2013 ).     

Looking for the “missing endmember” in South Atlantic Ocean mantle around Ascension Island

Summary Published bathymetric, geophysical and geochemical data suggest the presence of a hotspot in the South Atlantic, between the Bode Verde and Ascension fracture zones, in the vicinity of Ascension Island which has been named Ascension hotspot. We provide here new isotope and trace element data on Grattan (a near-ridge seamount) basalts which suggest that this edifice was produced by the plume presently feeding the adjacent Mid-Atlantic Ridge (MAR) segment and that it could represent the geochemical “missing end-member” described by Hanan et al. (1986) in the South Atlantic Ocean. Its isotopic signature differs slightly from that of Ascension Island and is more radiogenic in ²⁰⁸Pb. Trace element data are intermediate between those of MORB and characteristic OIB. We suggest that near Grattan seamount there exists a weak, perhaps intermittent, and slightly chemically heterogeneous plume.

---

Résumé De nombreuses données bathymétriques, géophysiques et géochimiques suggérent qu’il existe un point chaud dans l’Atlantique Sud entre les zones de fracture d’Ascencion et de Bode Verde, dans les environs de l’?le d’Ascencion et qui a été nommé point chaud d’Ascencion. Nous publions ici de nouvelles données géochimiques et isotopiques sur des basaltes du mont sous-marin Grattan qui est trés proche de la dorsale. Ces données montrent que cet édifice a été produit par le panache mantellique qui influence également la source des basaltes de la dorsale, et qu’il pourrait s’agir du “p?le manquant” décrit par Hanan et al. (1986) dans l’Atlantique Sud. Sa signature isotopique diffère légèrement de celle de l’?le d’Ascencion et présente notamment des valeurs plus radiogéniques en ²⁰⁸Pb. Les données en éléments en trace sont intermédiaires entre celle d’un MORB et d’un OIB caractéristique. Nous suggérons donc qu’il existe preès du mont sous-marin Grattan un panache mantellique faiblement actif, peut-être intermittent et chimiquement hétérogène.

---

---

Received January 15, 2000; revised version accepted April 17, 2000     

50 Myr recovery from the largest negative δ13C excursion in the Ediacaran ocean

Sedimentary rocks deposited during the Ediacaran period (～630–542 Ma) contain carbonates whose carbon isotopic ratios show a marked negative excursion consisting of a precipitous drop from +5‰ to ?12‰, followed by a sub‐linear recovery to positive δ¹³C values. Isotopic ages (U/Pb) and thermal subsidence modelling are combined to constrain the excursion in time and indicate an onset at ～600 Ma, and duration of recovery of approximately 50 Myr. The excursion is widely recognized in Oman and has potential correlatives in Ediacaran strata elsewhere, and may thus represent a characteristic feature of the Ediacaran period. The amplitude of this carbon isotope excursion far exceeds those of other Neoproterozoic anomalies. The isotopic trend of negative excursion and long‐term recovery spanned at least one short‐lived glacial episode (at 580 Ma), but appears unrelated to glaciation, which indicates that negative anomalies in the Neoproterozoic marine carbon isotope record are not directly or uniquely linked to ice ages.     

Effects of Tidal-Forcing Variations on Tidal Properties Along a Narrow Convergent Estuary

A 1D analytical framework is implemented in a narrow convergent estuary that is 78 km in length (the Guadiana, Southern Iberia) to evaluate the tidal dynamics along the channel, including the effects of neap-spring amplitude variations at the mouth. The close match between the observations (damping from the mouth to ～ 30 km, shoaling upstream) and outputs from semi-closed channel solutions indicates that the M₂ tide is reflected at the estuary head. The model is used to determine the contribution of reflection to the dynamics of the propagating wave. This contribution is mainly confined to the upper one third of the estuary. The relatively constant mean wave height along the channel (<?10% variations) partly results from reflection effects that also modify significantly the wave celerity and the phase difference between tidal velocity and elevation (contradicting the definition of an “ideal” estuary). Furthermore, from the mouth to ～ 50 km, the variable friction experienced by the incident wave at neap and spring tides produces wave shoaling and damping, respectively. As a result, the wave celerity is largest at neap tide along this lower reach, although the mean water level is highest in spring. Overall, the presented analytical framework is useful for describing the main tidal properties along estuaries considering various forcings (amplitude, period) at the estuary mouth and the proposed method could be applicable to other estuaries with small tidal amplitude to depth ratio and negligible river discharge.     

Accounting for Aquifer Heterogeneity from Geological Data to Management Tools

A nested workflow of multiple‐point geostatistics (MPG) and sequential Gaussian simulation (SGS) was tested on a study area of 6 km² located about 20 km northwest of Quebec City, Canada. In order to assess its geological and hydrogeological parameter heterogeneity and to provide tools to evaluate uncertainties in aquifer management, direct and indirect field measurements are used as inputs in the geostatistical simulations to reproduce large and small‐scale heterogeneities. To do so, the lithological information is first associated to equivalent hydrogeological facies (hydrofacies) according to hydraulic properties measured at several wells. Then, heterogeneous hydrofacies (HF) realizations are generated using a prior geological model as training image (TI) with the MPG algorithm. The hydraulic conductivity (K) heterogeneity modeling within each HF is finally computed using SGS algorithm. Different K models are integrated in a finite‐element hydrogeological model to calculate multiple transport simulations. Different scenarios exhibit variations in mass transport path and dispersion associated with the large‐ and small‐scale heterogeneity respectively. Three‐dimensional maps showing the probability of overpassing different thresholds are presented as examples of management tools.