Sulphur isotopes in carbonatites and associated silicate rocks from the Superior Province, Canada

Thirty-five S isotope analyses obtained from six carbonatite complexes from the Superior Province, Canadian Shield, ranging in age from 1,897 Ma to 1,093 Ma, have δ³⁴S_CDT values of between ?4.5‰ and +3.4‰. Pyrrhotite, chalcopyrite and pyrite mineral separates were used. Each complex possesses its own distinct range and mean S isotope composition. The range for Schryburt Lake is: ?4.5‰ to ?3.4‰ ( mean?=??3.9‰), for Big Beaver House: ?3.6‰ to ?1.5‰ (mean?=??2.2‰), for Cargill: ?1.5‰–+0.5‰ (mean?=??0.7‰), for Spanish River: ?0.1‰–+0.1‰ (mean?=?0.0‰), and for Firesand River: +1.3‰–+3.4‰ (mean?=?+1.7‰). A single sample from Carb Lake yielded a δ³⁴S_CDT value of +2.8‰. Differences in isotope compositions can be related to isotope effects brought about during melt generation and emplacment, such as variations in fo₂ and temperature. The different S and C isotope data for most complexes, however, suggest that the parental melts could have been generated from a heterogeneous mantle source, although process-driven changes cannot be completely ruled out.     

Source of parental melts to carbonatites–critical isotopic constraints

The restriction of most carbonatites to continental areas raises questions as to whether the parental liquids to carbonatites are generated within the continental lithosphere, or whether they are derived from deeper parts of the mantle with the lithosphere playing a subsidiary, but important, role in trapping volatile-rich melts/fluids. The constraints imposed by both radiogenic and stable isotopic data from carbonatites world-wide are consistent with a sub-lithospheric source for the parental melts, associated with either asthenospheric ‘upwellings’ or more deep-seated, plume-related activity. Crucial evidences that support the generation of carbonated melts from sub-lithospheric mantle are: the petrogenetic and temporal association of carbonatites with large igneous provinces (LIPs; e.g., Deccan, Parana), carbonatites with primitive noble gas isotopic signatures, radiogenic isotope ratios similar to OIBs (i.e. involving HIMU, EM1 and FOZO mantle components), and the uniform, time-integrated Rb/Sr and Sm/Nd development lines for Sr and Nd isotopic data for carbonatites from the Superior Province, Canada, and the Kola Peninsula, Russia. Such findings are difficult to reconcile with a lithospheric origin. Even if a metasomatized lithospheric mantle is considered to be the sole source of carbonated melts, an unknown mantle process is required, large enough to generate widespread, synchronous metasomatism spanning regions in excess of 1?×?10⁶?km².     

Anisotropic diffusion in a finite cylinder, with geochemical applications

Atomic diffusion in minerals may not be well represented by solutions to the diffusion equation for a sphere with a single-valued diffusivity, either because they have platy or elongated habits or because the energetics of diffusion is sensitive to crystallographic direction. In many cases, a cylinder having characteristic radial and axial diffusivities is arguably a better model, but rigorous solutions to the anisotropic diffusion equation for a finite cylinder have not been available. Here we develop general analytical solutions that capture both the internal distribution of diffusant as a function of time, C(r, z, t), and the fraction, F, of diffusant lost during a specified thermal history. These solutions are shown to conform with existing analytical expressions for limiting cases of diffusion in a slab or infinite cylinder. We present, in addition, a simple numerical (finite difference) approach that not only reproduces the results of our analytical expressions but also enables us to move beyond some of the limitations of the equations to simulate complex natural scenarios involving non-zero and time-dependent boundary conditions, arbitrary initial distribution of diffusant within the cylinder and simultaneous diffusion and radiogenic ingrowth. The complementary nature of the two approaches is emphasized and several illustrative applications to ‘real-world’ problems are described, including noble-gas thermochronometry and halogen-hydroxyl interdiffusion in apatite.     

Thermal constraints on the early history of the H-chondrite parent body reconsidered

Reconstructions of the early thermal history of the H-chondrite parent body have focused on two competing hypotheses. The first posits an undisturbed thermal evolution in which the degree of metamorphism increases with depth, yielding an “onion-shell” structure. The second posits an early fragmentation-reassembly event that interrupted this orderly cooling process. Here, we test these hypotheses by collecting a large number of previously published closure age and cooling rate data and comparing them to a suite of numerical models of thermal evolution in an idealized parent body. We find that the onion-shell hypothesis, when applied to a parent body of radius 75-130 km with a thermally insulating regolith, is able to explain 20 of the 21 closure age data and 62 of the 71 cooling rates. Furthermore, six of the eight meteorites for which multiple data (at different temperatures) are available, can be accounted for by onion-shell thermal histories. We therefore conclude that no catastrophic disruption of the H-chondrite parent body occurred during its early thermal history. The relatively small number of data not explained by the onion-shell hypothesis may indicate the formation of impact craters on the parent body which, while large enough to excavate all petrologic types, were small enough to leave the parent body largely intact. Impact events fulfilling these requirements would likely have produced transient crater diameters at least 30% of the parent body diameter.     

The stress field at Soultz-sous-Forêts from focal mechanisms of induced seismic events: Cases of the wells GPK2 and GPK3

The stress field at the EGS geothermal site of Soultz-sous-Forêts has been the subject of many studies, because it largely controls the response of the reservoir to fluid injection. The analysis of borehole logging data, especially breakouts and drilling-induced tension fractures, in the four geothermal wells define an average SHmax orientation ranging between 170° and 180° down to 5 km depth. It also reveals strong heterogeneities in several depth intervals. In this paper, the inversion of double-couple source mechanisms of seismic events induced during GPK2 and GPK3 stimulation tests is performed to retrieve the orientation and shape factor of the stress tensor, using the Slickenside Analysis Package of Michael, 1984, Michael, 1987a, Michael, 1987b. The results indicate a well-determined orientation of Shmin in GPK2 and in GPK3; in GPK3 Shmin is clockwise rotated by about 10°. The stress tensor defines an uniaxial extension. The results from both methods, analysis of borehole logging data and inversion of focal mechanisms, are only slightly different; the discrepancy may be due to the larger reservoir volume covered by the focal mechanisms, which can include strong stress heterogeneities.     

Contribution of the exploration of deep crystalline fractured reservoir of Soultz to the knowledge of enhanced geothermal systems (EGS)

Over the past 20 years, the Soultz experimental geothermal site in Alsace, France, has been explored in detail by the drilling of five boreholes, three of which extend to 5 km depth. Data on geology, fluid geochemistry, temperature, microseismicity, hydraulics and geomechanics have been collected and interpreted by the various teams from the participating European countries and their international collaborators. Two reservoirs have been developed within granite at depths of 3.5 and 5 km. The reservoir at 3.5 km was formed from two wells, 450 m apart, both of which were subjected to hydraulic stimulation injections. The system was circulated continuously for 4 months at 25 kg/s in 1997 using a downhole pump, and yielded results that were extremely encouraging. The impedance reduced to 0.1 MPa/l/s, the first time this long-standing target had been attained. Construction of a deeper system began shortly afterwards with the drilling of 3 deviated wells to 5 km true vertical depth, where the temperature was 200 °C. The wells were drilled in a line, 600 m apart at reservoir depth, and all were hydraulically stimulated and subjected to acidization injections. The 3-well system was circulated under buoyancy drive for 5 months in 2005 with injection in the central well, GPK-3, and production from the two outer wells, GPK-2 and GPK-4. This showed good linkage between one doublet pair, but not the other. Further acidization operations on the low-productivity well led to its productivity increasing to almost the same level as the other wells. Construction of a power plant at the site was completed in 2008 and a trial circulation with a production pump in one well and the other shut-in was conducted with power production. Downhole pumps are now installed in both production wells in preparation for long-term circulation of the system. In this article we present an overview of the principal accomplishments at Soultz over the past two decades, and highlight the main results, issues identified, and lessons learnt.     

The impacts of hysteresis on variably saturated hydrologic response and slope failure

This investigation employs 3D, variably saturated subsurface flow simulation to examine hysteretic effects upon the hydrologic response used to drive unsaturated slope stability assessments at the Coos Bay 1 (CB1) experimental catchment in the Oregon Coast Range, USA. Slope stability is evaluated using the relatively simple infinite slope model for unsaturated soils driven by simulated pore-water pressures for an intense storm that triggered a slope failure at CB1 on 18 November 1996. Simulations employing both hysteretic and non-hysteretic soil–water retention curves indicate that using either the drying soil–water retention curve or an intermediate soil–water retention curve that attempts to average the wetting and drying retention curves underestimates the near-surface hydrologic response and subsequently the potential for slope failure. If hysteresis cannot be considered in the hydrologic simulation, the wetting soil–water retention curve, which is seldom measured, should be used for more physically based slope stability assessment. Without considering hysteresis or using the wetting soil–water retention curve, the potential for landsliding in unsaturated materials may be underestimated and a slope failure could occur when simulations predict stability.     

A shallow subsurface controlled release facility in Bozeman, Montana, USA, for testing near surface CO2 detection techniques and transport models

A controlled field pilot has been developed in Bozeman, Montana, USA, to study near surface CO₂ transport and detection technologies. A slotted horizontal well divided into six zones was installed in the shallow subsurface. The scale and CO₂ release rates were chosen to be relevant to developing monitoring strategies for geological carbon storage. The field site was characterized before injection, and CO₂ transport and concentrations in saturated soil and the vadose zone were modeled. Controlled releases of CO₂ from the horizontal well were performed in the summers of 2007 and 2008, and collaborators from six national labs, three universities, and the U.S. Geological Survey investigated movement of CO₂ through the soil, water, plants, and air with a wide range of near surface detection techniques. An overview of these results will be presented.     

Scaling Properties of Temperature Spectra and Heat-Flux Cospectra in the Surface Friction Layer Beneath an Unstable Outer Layer

Temperature variance and temperature power spectra in the unstable surface layer have always presented a problem to the standard Monin-Obukhov similarity model. Recently that problem has intensified with the demonstration by Smedman et al. (2007, Q J Roy Meteorol Soc 133: 37–51) that temperature spectra and heat-flux cospectra can have two distinct peaks in slightly unstable conditions, and by McNaughton et al. (2007, Nonlinear Process Geophys 14: 257–271) who showed that the wavenumber of the peak of temperature spectra in a convective boundary layer (CBL), closely above the surface friction layer (SFL), can be sensitive to the CBL depth, z _i. Neither the two-peak form at slight instability nor the dependence of peak position on z _i at large instability is compatible with the Monin-Obukhov model. Here we examine the properties of temperature spectra and heat-flux cospectra from between these extremes, i.e. from within the unstable SFL, in two experiments. The analysis is based on McNaughton’s model of the turbulence structure in the SFL. According to this model, heat is transported through most of the SFL by sheet plumes, created by the action of impinging outer eddies. The smallest and most effective of these outer eddies have sizes that scale on SFL depth, z _s. The z _s-scale eddies and plumes are organised within the overall convection pattern in the CBL, and in turn they organise the motion of smaller eddies within the SFL, whose sizes scale on height, z. The main experimental results are: (1) the peak amplitudes of the temperature spectra in the SFL are collapsed with a scaling factor (z_sz)^1/3e_o^2/3{(z_{\rm s}z)^{1/3}\varepsilon_{\rm o}^{2/3}} divided by the square of the surface temperature flux, where e_o{\varepsilon_{\rm o}} is the dissipation rate of turbulent energy in the outer CBL (above the SFL); (2) the peak wavenumbers of the temperature spectra are collapsed with the mixed length scale (z _i z _s)^1/2; (3) the peak wavenumbers of the heat-flux cospectra are collapsed with the doubly-mixed length scale (z _i z _s)^1/4 z ^1/2; (4) for z/z _s < 0.03, the peak in the cospectrum is replaced by another peak at a wavenumber about a magnitude larger. This peak’s position scales on z; (5) all these findings are consistent with the observations of Smedman et al.