Dilution of non-reactive tracers in variably saturated sandy structures

Based on a dye tracer experiment in a sand tank we addressed the problem of local dispersion of conservative tracers in the unsaturated zone. The sand bedding was designed to have a defined spatial heterogeneity with a strong anisotropy. We estimated the parameters that characterize the local dispersion and dilution from concentration maps of a high spatial and temporal resolution obtained by image analysis. The plume spreading and mixing behavior was quantified on the basis of the coefficient of variation of the concentration and of the dilution index. The heterogeneous structure modified the flow pattern depending on water saturation. The shape of the tracer plumes revealed the structural signature of the sand bedding at low saturation only. In this case pronounced preferential flow was observed. At higher flow rates the structure remained hidden by a spatially almost homogeneous behavior of the plumes. In this context, we mainly discuss the mechanism of re-distributing a finite mass of inert solutes over a large volume, due to macro- and micro-heterogeneities of the structure. (C) 2001 Elsevier Science Ltd. AU rights reserved.

Upscaling of anisotropy in unsaturated Miller-similar porous media

Geological and pedological processes rarely form isotropic media as is usually assumed in transport studies. Anisotropy at the Darcy or field scale may be detected directly by measuring flow parameters or may become indirectly evident from movement and shape of solute plumes. Anisotropic behavior of a soil at one scale may, in many cases, be related to the presence of lower-scale directional structures. Miller similitude with different pore-scale geometries of the basic element is used to model macroscopic flow and transport behavior. Analytical expressions for the anisotropic conductivity tensor are derived based on the dynamic law that governs the flow problem at the pore scale. The effects of anisotropy on transport parameters are estimated by numerical modeling.

Tracer applications of noble gas radionuclides in the geosciences

Noble gas radionuclides, including 81Kr (t1/2 = 229,000 years), 85Kr (t1/2 = 10.8 years), and 39Ar (t1/2 = 269 years), possess nearly ideal chemical and physical properties for studies of earth and environmental processes. Recent advances in Atom Trap Trace Analysis (ATTA), a laser-based atom counting method, have enabled routine measurements of the radiokrypton isotopes, as well as the demonstration of the ability to measure 39Ar in environmental samples. Here we provide an overview of the ATTA technique, and a survey of recent progress made in several laboratories worldwide.We review the application of noble gas radionuclides in the geosciences and discuss how ATTA can help advance these fields, specifically: determination of groundwater residence times using 81Kr, 85Kr, and 39Ar; dating old glacial ice using 81Kr; and an 39Ar survey of the main water masses of the oceans, to study circulation pathways and estimate mean residence times. Other scientific questions involving a deeper circulation of fluids in the Earth's crust and mantle are also within the scope of future applications. We conclude that the geoscience community would greatly benefit from an ATTA facility dedicated to this field, with instrumentation for routine measurements, as well as for research on further development of ATTA methods.

Revisiting the Nd-142 deficits in the 1.48 Ga Khariar alkaline rocks, India

The 146Sm–142Nd system plays a central role in tracing the silicate differentiation of the Earth prior to 4.1 Ga. After this time, given its initial abundance, the 146Sm can be considered to be extinct. Upadhyay et al. (2009) reported unexpected negative 142Nd anomalies in 1.48 Ga rocks of the Khariar nepheline syenite complex (India) and inferred that an early enriched, low-Sm/Nd reservoir must have contributed to the mantle source rocks of the Khariar complex. As 146Sm had been effectively extinct for about 2.6 billion years before the crystallisation of the Khariar samples, this Nd signature should have remained isolated from the convective mantle for at least that long. It was thus suggested that the source rock of Khariar samples had been sequestered in the lithospheric root of the Indian craton. Using a different chemical separation method, and a different Thermal Ionization Mass Spectrometry (TIMS) analysis protocol, the present study attempted to replicate these negative 142Nd anomalies, but none were found. To determine which data set is correct, we investigated three possible sources of bias between them: imperfect cancellation of Faraday collector efficiencies during multidynamic TIMS analysis, rapid sample fractionation between the sequential measurement of 146Nd/144Nd and 142Nd/144Nd, and non-exponential law behaviour resulting from so-called “domain mixing.” Incomplete cancellation of collector efficiencies was found unlikely to cause resolvable biases at the estimated level of variation among collector efficiencies. Even in the case of highly variable efficiency and resolvable biases, there is no reason to suspect that they would reproducibly affect only four rocks out of 10 analysed by Upadhyay et al. (2009). Although domain mixing may explain apparent “reverse” fractionation trends observed in some TIMS analyses, it cannot be the cause of the apparent negative anomalies in the study of Upadhyay et al. (2009). It was determined that rapid mass fractionation during the course of a multidynamic TIMS analysis can bias all measured Nd ratios. After applying an approximate correction for this effect, only one rock from Upadhyay et al. (2009) retained an apparent negative 142Nd anomaly. This, in conjunction with our new, anomaly-free data set measured at fractionation rates too low to cause bias, leads to the conclusion that the anomalies reported by Upadhyay et al. (2009) are a subtle and reproducible analytical artefact. The absence of negative 142Nd anomalies in these rocks relaxes the need for a mechanism (other than crust formation) that can isolate a Nd reservoir from the convective mantle for billions of years.

Continental degassing of 4He by surficial discharge of deep groundwater

Radiogenic He is produced by the decay of uranium and thorium in the Earth’s mantle and crust. From here, it is degassed to the atmosphere and eventually escapes to space. Assuming that all of the 4He produced is degassed, about 70% of the total He degassed from Earth comes from the continental crust. However, the outgoing flux of crustal He has not been directly measured at the Earth’s surface and the migration pathways are poorly understood. Here we present measurements of helium isotopes and the long-lived cosmogenic radio-isotope Kr in the deep, continental-scale Guarani aquifer in Brazil and show that crustal He reaches the atmosphere primarily by the surficial discharge of deep groundwater. We estimate that He in Guarani groundwater discharge accounts for about 20% of the assumed global flux from continental crust, and that other large aquifers may account for about 33%. Old groundwater ages suggest that He in the Guarani aquifer accumulates over half- to one-millionyear timescales. We conclude that He degassing from the continents is regulated by groundwater discharge, rather than episodic tectonic events, and suggest that the assumed steady state between crustal production and degassing of He, and its resulting atmospheric residence time, should be re-examined.

Anisotropy of magnetic susceptibility in natural olivine single crystals

Mantle flow dynamics can cause preferential alignment of olivine crystals that results in anisotropy of physical properties. To interpret anisotropy in mantle rocks, it is necessary to understand the anisotropy of olivine single crystals. We determined anisotropy of magnetic susceptibility (AMS) for natural olivine crystals. High-field AMS allows for the isolation of the anisotropy due to olivine alone. The orientations of the principal susceptibility axes are related to the olivine’s crystallographic structure as soon as it contains >3 wt % FeO. The maximum susceptibility is parallel to the c axis both at room temperature (RT) and at 77 K. The orientation of the minimum axis at RT depends on iron content; it is generally parallel to the a axis in crystals with 3–5 wt % FeO, and along b in samples with 6–10 wt % FeO. The AMS ellipsoid is prolate and the standard deviatoric susceptibility, k0, is on the order of 8*10210 m3/kg for the samples with <1wt % FeO, and ranges from 3.1*1029 m3/kg to 5.7*1029 m3/kg for samples with 3–10 wt % FeO. At 77 K, the minimum susceptibility is along b, independent of iron content. The shape of the AMS ellipsoid is prolate for samples with <5 wt % FeO, but can be prolate or oblate for higher iron content. The degree of anisotropy increases at 77 K with p0 7757.160.5. The results from this study will allow AMS fabrics to be used as a proxy for olivine texture in ultramafic rocks with high olivine content.

Lithium isotope constraints on crust–mantle interactions and surface processes on Mars

Lithium abundances and isotope compositions are reported for a suite of martian meteorites that span the range of petrological and geochemical types recognized to date for Mars. Samples include twenty-one bulk-rock enriched, intermediate and depleted shergottites, six nakhlites, two chassignites, the orthopyroxenite Allan Hills (ALH) 84001 and the polymict breccia Northwest Africa (NWA) 7034. Shergottites unaffected by terrestrial weathering exhibit a range in δ7Li from 2.1 to 6.2‰, similar to that reported for pristine terrestrial peridotites and unaltered mid-ocean ridge and ocean island basalts. Two chassignites have δ7Li values (4.0‰) intermediate to the shergottite range, and combined, these meteorites provide the most robust current constraints on δ7Li of the martian mantle. The polymict breccia NWA 7034 has the lowest δ7Li (−0.2‰) of all terrestrially unaltered martian meteorites measured to date and may represent an isotopically light surface end-member. The new data for NWA 7034 imply that martian crustal surface materials had both a lighter Li isotope composition and elevated Li abundance compared with their associated mantle. These findings are supported by Li data for olivine-phyric shergotitte NWA 1068, a black glass phase isolated from the Tissint meteorite fall, and some nakhlites, which all show evidence for assimilation of a low-δ7Li crustal component. The range in δ7Li for nakhlites (1.8 to 5.2‰), and co-variations with chlorine abundance, suggests crustal contamination by Cl-rich brines. The differences in Li isotope composition and abundance between the martian mantle and estimated crust are not as large as the fractionations observed for terrestrial continental crust and mantle, suggesting a difference in the styles of alteration and weathering between water-dominated processes on Earth versus possibly Cl–S-rich brines on Mars. Using high-MgO shergottites (>15 wt.% MgO) it is possible to estimate the δ7Li of Bulk Silicate Mars (BSM) to be 4.2 ± 0.9‰ (2σ). This value is at the higher end of estimates for the Bulk Silicate Earth (BSE; 3.5 ± 1.0‰, 2σ), but overlaps within uncertainty.

The tectonometamorphic evolution of the Sesia–Dent Blanche nappes (internal Western Alps): review and synthesis

This study reviews and synthesizes the present knowledge on the Sesia–Dent Blanche nappes, the highest tectonic elements in the Western Alps (Switzerland and Italy), which comprise pieces of pre-Alpine basement and Mesozoic cover. All of the available data are integrated in a crustal-scale kinematic model with the aim to reconstruct the Alpine tectono-metamorphic evolution of the Sesia–Dent Blanche nappes. Although major uncertainties remain in the pre-Alpine geometry, the basement and cover sequences of the Sesia–Dent Blanche nappes are seen as part of a thinned continental crust derived from the Adriatic margin. The earliest stages of the Alpine evolution are interpreted as recording late Cretaceous subduction of the Adria-derived Sesia–Dent Blanche nappes below the South-Alpine domain. During this subduction, several sheets of crustal material were stacked and separated by shear zones that rework remnants of their Mesozoic cover. The recently described Roisan-Cignana Shear Zone of the Dent Blanche Tectonic System represents such a shear zone, indicating that the Sesia–Dent Blanche nappes represent a stack of several individual nappes. During the subsequent subduction of the Piemonte–Liguria Ocean large-scale folding of the nappe stack (including the Roisan-Cignana Shear Zone) took place under greenschist facies conditions, which indicates partial exhumation of the Dent Blanche Tectonic System. The entrance of the Briançonnais micro-continent within the subduction zone led to a drastic change in the deformation pattern of the Alpine belt, with rapid exhumation of the eclogite-facies ophiolite bearing units and thrust propagation towards the foreland. Slab breakoff probably was responsible for allowing partial melting in the mantle and Oligocene intrusions into the most internal parts of the Sesia–Dent Blanche nappes. Finally, indentation of the Adriatic plate into the orogenic wedge resulted in the formation of the Vanzone back-fold, which marks the end of the pervasive ductile deformation within the Sesia–Dent Blanche nappes during the earliest Miocene.

A deep crust–mantle boundary in the asteroid 4 Vesta

The asteroid 4 Vesta was recently found to have two large impact craters near its south pole, exposing subsurface material. Modelling suggested that surface material in the northern hemisphere of Vesta came from a depth of about 20 kilometres, whereas the exposed southern material comes from a depth of 60 to 100 kilometres. Large amounts of olivine from the mantle were not seen, suggesting that the outer 100 kilometres or so is mainly igneous crust. Here we analyse the data on Vesta and conclude that the crust–mantle boundary (or Moho) is deeper than 80 kilometres.

The behaviour of incompatible elements during hydrous melting of metasomatized peridotite at 4–6 GPa and 1000 °C–1200 °C

Trace element behavior during hydrous melting of a metasomatized garnet–peridotite was examined at pressures of 4–6 GPa and temperatures of 1000 °C–1200 °C, conditions appropriate for fluid penetrating the mantle wedge atop the subducting slab. Experiments were performed in a rocking multi-anvil apparatus using a diamond-trap setup. The compositions of the fluid and melt phases were measured using the cryogenic LA-ICP-MS technique. The water-saturated solidus of the K-lherzolite composition is located between 900 °C and 1000 °C at 4 GPa and between 1000 °C and 1100 °C at 5 and 6 GPa. The partition coefficients between fluid or melt and clinopyroxene reveal an asymmetric MREE trough with a minimum at Dy. The clinopyroxene in equilibrium with aqueous fluids is characterized by DUfluid–cpx > DThfluid–cpx while DUmelt–cpx tends to be similar to DThmelt–cpx. The partition coefficients between fluid or melt and garnet reveal very strong light to heavy REE fractionation, DLa/DLu from 95 (hydrous melt) to 1600 (aqueous fluid). The LILE are highly incompatible with partition coefficients > 50. The behavior of HFSE are decoupled, with DZr,Hf close to 1 while DNb,Ta > 10. Garnet is characterized by DUmelt/fluid–garnet < DThmelt/fluid–garnet. A comparison of our experimental partitioning results for trivalent cations as well as the results from the literature and the calculations carried out using the lattice strain model adapted to the presence of water in the bulk system indicates that H2O in the fluid or melt phase has a prominent effect on trace element partitioning. Garnet in mantle rocks in equilibrium with an aqueous fluid is characterized by significantly higher Do(3 +) for REE in the X site of the garnet compared with the partitioning values of the optimal cation in garnet in equilibrium with hydrous melts. Our data show for the first time that the change in the nature of the mobile phase (fluid vs. melt) does affect the affinities of trace elements into the garnet crystal at conditions below the second critical endpoint of the system. The same also applies for clinopyroxene, although this is less clear. Consequently, our new data allow for refinements in predictive modeling of element transfer from the slab to the mantle wedge and of possible compositions of metasomatized mantle that sources OIB magmatism.

Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.

Complex geomorphologic assemblage of terrains in association with the banded terrain in Hellas basin, Mars

Hellas basin acts as a major sink for the southern highlands of Mars and is likely to have recorded several episodes of sedimentation and erosion. The north-western part of the basin displays a potentially unique Amazonian landscape domain in the deepest part of Hellas, called “banded terrain”, which is a deposit characterized by an alternation of narrow band shapes and inter-bands displaying a sinuous and relatively smooth surface texture suggesting a viscous flow origin. Here we use high-resolution (HiRISE and CTX) images to assess the geomorphological interaction of the banded terrain with the surrounding geomorphologic domains in the NW interior of Hellas to gain a better understanding of the geological evolution of the region as a whole. Our analysis reveals that the banded terrain is associated with six geomorphologic domains: a central plateau named Alpheus Colles, plain deposits (P1 and P2), reticulate (RT1 and RT2) and honeycomb terrains. Based on the analysis of the geomorphology of these domains and their cross-cutting relationships, we show that no widespread deposition post-dates the formation of the banded terrain, which implies that this domain is the youngest and latest deposit of the interior of Hellas. Therefore, the level of geologic activity in the NW Hellas during the Amazonian appears to have been relatively low and restricted to modification of the landscape through mechanical weathering, aeolian and periglacial processes. Thermophysical data and cross-cutting relationships support hypotheses of modification of the honeycomb terrain via vertical rise of diapirs such as ice diapirism, and the formation of the plain deposits through deposition and remobilization of an ice-rich mantle deposit. Finally, the observed gradual transition between honeycomb and banded terrain suggests that the banded terrain may have covered a larger area of the NW interior of Hellas in the past than previously thought. This has implications on the understanding of the evolution of the deepest part of Hellas.

Eclogitic diamonds from variable crustal protoliths in the northeastern Siberian craton: trace elements and coupled δ13C – δ18O signatures in diamonds and garnet inclusions

Diamonds of eclogitic assemblages are dominant in the placer diamond deposits of the northeastern Siberian platform. In this study we present new trace elements and stable isotopes (δ13C and δ18O) data for alluvial diamonds and their garnet inclusions from this locality. Cr-rich garnets of peridotitic affinity in the studied diamonds have a narrow range of δ18O values from 5.7‰ to 6.2‰, which is largely overlapping with the accepted mantle range. This narrow range suggests that the garnet inclusions showing different REE patterns and little variations in oxygen isotopes may have formed by different processes involving fluid/melts that, however, were in oxygen isotopic equilibrium with the mantle. The trace element composition of the eclogitic garnet inclusions supports a crustal origin for at least the high-Ca garnets, which show flat HREE patterns and in some cases a positive Eu-anomaly. High-Ca eclogitic garnets generally show heavier oxygen isotope compositions (δ18O 6.5–9.6‰) than what is observed in low-Ca garnets (δ18O 5.7–7.4‰). The variability in oxygen isotopes and trace elements is suggested to be inherited from contrasting crustal protoliths. The relationship between the high δ18O values of inclusions and the low δ13C values of the host diamonds implies that the high-Ca garnet inclusions were derived from intensely hydrated (e.g., δ18O > 7‰) and typically oxidised basaltic rock close to the seawater interface, and that the carbon for diamonds was closely associated with this protolith.

Sublimation of ice–tholins mixtures: A morphological and spectro-photometric study

Sublimation, the direct transition from solid to gas phase, is a process responsible for shaping and changing the reflectance properties of many Solar System surfaces. In this study, we have characterized the evolution of the structure/texture and of the visible and near-infrared (VIS–NIR) spectral reflectance of surfaces made of water ice mixed with analogues of complex extraterrestrial organic matter, named tholins, under low temperature (<-70° C) and pressure (10-⁵mbar) conditions. The experiments were carried out in the SCITEAS simulation setup recently built as part of the Laboratory for Outflow Studies of Sublimating Materials (LOSSy) at the University of Bern (Pommerol, A. et al. [2015a]. Planet. Space Sci. 109–110, 106–122). As the water ice sublimated, we observed in situ the formation of a sublimation lag deposit made of a water-free porous (>90% porosity) network of organic filaments on top of the ice. The temporal evolution of the tholins and water ice spectral features (reflectance at the absorption bands wavelengths, red slope, from 0.40 to 1.90lm) are analyzed throughout the sublimation of the samples. We studied how different mixtures of tholins with water (0.1 wt.% tholins as coating or inclusions within the water particles), and different ice particle sizes (4.5 ± 2.5 or 67 ± 31lm) influence the morphological and spectral evolutions of the samples. The sublimation of the ice below the mantle produces a gas flow responsible for the ejection of mm to cm-sized fragments of the deposit in outbursts-like events. The results show remarkable differences between these samples in term of mantle structure, speed of mantle building, rates and surface area of mantle ejections. These data provide useful references for interpreting remote-sensing observations of icy Solar System surfaces, in particular the activity of comet nuclei where sublimation of organic-rich ices and deposition of organic-dust particles likely play a major role. Consequently, the data presented here could be of high interest for the interpretation of Rosetta, and also New Horizons, observations.

Sublimation of water ice mixed with silicates and tholins: Evolution of surface texture and reflectance spectra, with implications for comets

The surfaces of many objects in the Solar System comprise substantial quantities of water ice sometimes mixed with minerals and/or organic molecules. The sublimation of the ice changes the structural and optical properties of these objects. We present laboratory data on the evolution of the structure and the visible and near-infrared spectral reflectance of icy surface analogues of cometary ices, made of water ice, complex organic matter (tholins) and silicates, as they undergo sublimation under low temperature (<-70°C) and pressure (10-⁵mbar) conditions inside the SCITEAS simulation chamber. As the water ice sublimated, we observed in situ the formation of a porous sublimation lag deposit, or sublimation mantle, at the top of the ice. This mantle is a network of filaments made of the non-volatile particles. Organics or phyllosilicates grains, able to interact via stronger inter-particulate forces than olivine grains, can form a foam-like structure having internal cohesiveness, holding olivine grains together. As this mantle builds-up, the band depths of the sub-surface water ice are attenuated until complete extinction under only few millimeters of mantle. Optically thick sublimation mantles are mainly featureless in the near infrared. The absorption bands of the minerals present in the mantle are weak, or even totally absent if minerals are mixed with organics which largely dominate the VIS–NIR reflectance spectrum. During sublimation, ejections of large fragments of mantle, triggered by the gas flow, expose ice particles to the surface. The contrast of brightness between mantled and ice-exposed areas depends on the wavelength range and the dust/ice ratio considered. We describe how the chemical nature of the non-volatiles, the size of their particles, the way they are mixed with the ice and the dust/ice mass ratio influence the texture, activity and spectro-photometric properties of the sublimation mantles. These data provide useful references for interpreting remote-sensing observations of comets and also icy satellites or trans-neptunian objects.