A subduction wedge origin for Paleoarchean peridotitic diamonds and harzburgites from the Panda kimberlite, Slave craton: evidence from Re–Os isotope systematics

An extensive study of peridotitic sulfide inclusion bearing diamonds and their prospective harzburgitic host rocks from the 53 Ma Panda kimberlite pipe, Ekati Mine, NWT Canada, has been undertaken with the Re–Os system to establish their age and petrogenesis. Diamonds with peridotitic sulfide inclusions have poorly aggregated nitrogen (<30% N as B centers) at N contents of 200–800 ppm which differs from that of chromite and silicate bearing diamonds and indicates residence in the cooler portion of the Slave craton lithospheric mantle. For most of the sulfide inclusions, relatively low Re contents (average 0.457 ppm) and high Os contents (average 339 ppm) lead to extremely low ¹⁸⁷Re/¹⁸⁸Os, typically << 0.05. An age of 3.52 ± 0.17 Ga (MSWD = 0.46) and a precise initial ¹⁸⁷Os/¹⁸⁸Os of 0.1093 ± 0.0001 are given by a single regression of 11 inclusions from five diamonds that individually provide coincident internal isochrons. This initial Os isotopic composition is 6% enriched in ¹⁸⁷Os over 3.5 Ga chondritic or primitive mantle. Sulfide inclusions with less radiogenic initial Os isotopic compositions reflect isotopic heterogeneity in diamond forming fluids. The harzburgites have even lower initial ¹⁸⁷Os/¹⁸⁸Os than the sulfide inclusions, some approaching the isotopic composition of 3.5 Ga chondritic mantle. In several cases isotopically distinct sulfides occur in different growth zones of the same diamond. This supports a model where C–O–H–S fluids carrying a radiogenic Os signature were introduced into depleted harzburgite and produced diamonds containing sulfides conforming to the 3.5 Ga isochron. Reaction of this fluid with harzburgite led to diamonds with less radiogenic inclusions while elevating the Os isotope ratios of some harzburgites. Subduction is a viable way of introducing such fluids. This implies a role for subduction in creating early continental nuclei at 3.5 Ga and generating peridotitic diamonds.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Origin of the Nchanga copper–cobalt deposits of the Zambian Copperbelt

The Zambian Copperbelt forms the southeastern part of the 900-km-long Neoproterozoic Lufilian Arc and contains one of the world’s largest accumulations of sediment-hosted stratiform copper mineralization. The Nchanga deposit is one of the most significant ore systems in the Zambian Copperbelt and contains two major economic concentrations of copper and cobalt, hosted within the Lower Roan Group of the Katangan Supergroup. A Lower Orebody (copper only) and Upper Orebody (copper and cobalt) occur towards the top of arkosic units and within the base of overlying shales. The sulfide mineralogy includes pyrite, bornite, chalcopyrite, and chalcocite, although in the Lower Orebody, sulfide phases are partially or completely replaced by malachite and copper oxides. Carrollite is the major cobalt-bearing phase and is restricted to fault-propagation fold zones within a feldspathic arenite. Hydrothermal alteration minerals include dolomite, phlogophite, sericite, rutile, quartz, tourmaline, and chlorite. Quartz veins from the mine sequence show halite-saturated fluid inclusions, ranging from ~31 to 38 wt% equivalent NaCl, with homogenisation temperatures (Th_TOT) ranging between 140 and 180°C. Diagenetic pyrites in the lower orebody show distinct, relatively low δ ³⁴S, ranging from −1 to −17‰ whereas arenite- and shale-hosted copper and cobalt sulfides reveal distinctly different δ ³⁴S from −1 to +12‰ for the Lower Orebody and +5 to +18‰ for the Upper Orebody. There is also a clear distinction between the δ ³⁴S mean of +12.1±3.3‰ (n=65) for the Upper Orebody compared with +5.2±3.6‰ (n=23) for the Lower Orebody. The δ ¹³C of dolomites from units above the Upper Orebody give δ ¹³C values of +1.4 to +2.5‰ consistent with marine carbon. However, dolomite from the shear-zones and the alteration assemblages within the Upper Orebody show more negative δ ¹³C values: −2.9 to −4.0‰ and −5.6 to −8.3‰, respectively. Similarly, shear zone and Upper Orebody dolomites give a δ ¹⁸O of +11.7 to +16.9‰ compared to Lower Roan Dolomites, which show δ ¹⁸O of +22.4 to +23.0‰_. Two distinct structural regimes are recognized in the Nchanga area: a weakly deformed zone consisting of basement and overlying footwall siliciclastics, and a moderate to tightly folded zone of meta-sediments of the Katangan succession. The fold geometry of the Lower Roan package is controlled by internal thrust fault-propagation folds, which detach at the top of the lowermost arkose or within the base of the overlying stratigraphy and show vergence towards the NE. Faulting and folding are considered to be synchronous, as folding predominantly occurred at the tips of propagating thrust faults, with local thrust breakthrough. The data from Nchanga suggests a strong link between ore formation and the development of structures during basin inversion as part of the Lufilian Orogeny. Sulfides tend to be concentrated within arenites or coarser-grained layers within shale units, suggesting that host-rock porosity and possibly permeability played a role in ore formation. However, sulfides are also commonly orientated along, but not deformed by, a tectonic fabric or hosted within small fractures that suggest a significant role for deformation in the development of the mineralization. The ore mineralogy, hydrothermal alteration, and stable isotope data lend support to models consistent with the thermochemical reduction of a sulfate- (and metal) enriched hydrothermal fluid, at the site of mineralization. There is no evidence at Nchanga for a contribution of bacteriogenic sulfide, produced during sedimentation or early diagenesis, to the ores.Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.Editorial handling: H. Frimmel     

The evolution of groundwater in the Tyrrell catchment, south-central Murray Basin, Victoria, Australia

The Tyrell catchment lies on the western margin of the Riverine Province in the south-central Murray Basin, one of Australia’s most important groundwater resources. Groundwater from the shallow, unconfined Pliocene Sands aquifer and the underlying Renmark Group aquifer is saline (total dissolved solids up to 150,000 mg/L) and is Na-Cl-Mg type. There is no systematic change in salinity along hydraulic gradients implying that the aquifers are hydraulically connected and mixing during vertical flow is important. Stable isotopes (¹⁸O+²H) and Cl/Br ratios indicate that groundwater is entirely of meteoric origin and salts in this system have largely been derived by evapotranspiration of rainfall with only minor halite dissolution, rock weathering (mainly feldspar dissolution), and ion exchange between Na and Mg on clays. Similarity in chemistry of all groundwater in the catchment implies relative consistency in processes over time, independent of any climatic variation. Groundwater in both the Pliocene Sands and Renmark Group aquifers yield ages of up to 25 ka. The Tyrrell Catchment is arid to semi-arid and has low topography. This has resulted in relatively low recharge rates and hydraulic gradients that have resulted in long groundwater residence times.     

Computer microtomography of zircon: A new approach to the selection of targets for U-Pb geochronologic studies

The paper reports results obtained by the complex studying of zircon crystals from rocks of various genesis. Zircon is one of the minerals most often used as geochronometers. It also provides genetic information on superimposed processes that is “recorded” in the external and internal morphology of its crystals. This mineral is thoroughly examined to select its unaltered crystals for U-Pb dating by the single-grain method. Zircon grains are selected with the use of optical and electron microscopy and cathodoluminescence. This publication presents the first results obtained by examining zircon by computer microtomography (μCT) and the results of the studying of the external and internal structure by conventional techniques ([optical microscopy and SEM (SE and CL)]. The paper demonstrates the advantages of the application of the μCT techniques to the selection of targets for U-Pb zircon dating: there is no need for the special preparation of the samples and no material is destructed. However, this technique may be not informative enough if the zircon contains inherited core whose density does not differ from the density of the surrounding mineral.     

Geochronological constraints on the Legs Lake shear zone with implications for regional exhumation of lower continental crust, western Churchill Province, Canadian Shield

The Legs Lake shear zone marks the southeastern boundary of an extensive region (>20,000 km²) of high-pressure (0.8–1.5+ GPa) granulite-facies rocks in the western Churchill Province, Canada. The shear zone is one of the largest exhumation-related structures in the Canadian Shield and coincides with the central segment of the ∼2,800 km long Snowbird tectonic zone. The movement history of this shear zone is critical for the development of models for the exhumation history of the high-pressure region. We used electron microprobe U–Th–Pb dating of monazite with supplemental ID-TIMS U–Pb geochronology to place constraints on the timing of shear zone activity. Combining these and other data, we suggest that regional exhumation occurred during at least three distinct phases over an ∼150 million year period. The first phase involved high temperature decompression from ∼1.0 to 0.8–0.7 GPa shortly following 1.9 Ga peak metamorphism, possibly under an extensional regime. The second phase involved rock uplift and decompression of the hanging wall to 0.5–0.4 GPa during east-vergent thrusting across the Legs Lake shear zone at ca. 1.85 Ga. This phase was likely driven by early collision-related convergence in the Trans-Hudson orogen. The final phase of regional exhumation, involving the removal of 15–20 km of overburden from both footwall and hanging wall, likely occurred after ∼1.78 Ga and may have been related to regional extensional faulting.Electronic Supplementary Material Supplementary material is available to authorised users in the online version of this article at .     

Timescales and significance of high-pressure, high-temperature metamorphism and mafic dike anatexis, Snowbird tectonic zone, Canada

New geochronological, isotopic and geochemical data for a spectacular swarm of deep crustal migmatitic mafic dikes offer important insight into processes operative during 1.9 Ga high pressure, high temperature metamorphism along the Snowbird tectonic zone in northern Saskatchewan. High-precision U–Pb zircon dates reveal anatexis of Chipman mafic dikes at 1,896.2 ± 0.3 Ma during syntectonic and synmetamorphic intrusion at conditions of 1.0–1.2 GPa, >750°C. U–Pb zircon dates of 1,894–1,891 Ma for cross-cutting pegmatites place a lower bound on major metamorphism and deformation at the currently exposed crustal levels. The persistence of elevated temperatures for ~14 m.y. following peak conditions is implied by younger U–Pb titanite dates, and by Sm–Nd whole rock isotopic data that suggest the derivation of the pegmatites by melting of a mafic source. Limited melting of the host felsic gneiss at 1.9 Ga despite high temperature is consistent with evidence for their previous dehydration by granulite facies metamorphism in the Archean. Spatial heterogeneity in patterns of mafic dike and tonalitic gneiss anatexis can be attributed to lateral peak temperature and compositional variability. We correlate 1,896 Ma Chipman mafic dike emplacement and metamorphism with substantial 1.9 Ga mafic magmatism over a minimum along-strike extent of 1,200 km of the Snowbird tectonic zone. This suggests a significant, continent-wide period of asthenospheric upwelling that induced incipient continental rifting. Extension was subsequently terminated by hinterland contraction associated with Trans-Hudson accretion and orogenesis. Little activity in the lower crust for ca. 650 m.y. prior to Proterozoic metamorphism and mafic magmatism implies an extended interval of cratonic stability that was disrupted at 1.9 Ga. This episode of destabilization contrasts with the record of long-term stability in most preserved cratons, and is important for understanding the lithospheric characteristics and tectonic circumstances that control the destruction or survival of continents.     

Site Investigation for the Effects of Vegetation on Ground Stability

The procedure for geotechnical site investigation is well established but little attention is currently given to investigating the potential of vegetation to assist with ground stability. This paper describes how routine investigation procedures may be adapted to consider the effects of the vegetation. It is recommended that the major part of the vegetation investigation is carried out, at relatively low cost, during the preliminary (desk) study phase of the investigation when there is maximum flexibility to take account of findings in the proposed design and construction. The techniques available for investigation of the effects of vegetation are reviewed and references provided for further consideration. As for general geotechnical investigation work, it is important that a balance of effort is maintained in the vegetation investigation between (a) site characterisation (defining and identifying the existing and proposed vegetation to suit the site and ground conditions), (b) testing (in-situ and laboratory testing of the vegetation and root systems to provide design parameters) and (c) modelling (to analyse the vegetation effects).     

Multi-reference evaluation of uncertainty in earth orientation parameter measurements

Uncertainties in polar motion and length-of-day measurements are evaluated empirically using several data series from the space-geodetic techniques of the global positioning system (GPS), satellite laser ranging (SLR), and very long baseline interferometry (VLBI) during 1997–2002. In the evaluation procedure employed here, known as the three-corner hat (TCH) technique, the signal common to each series is eliminated by forming pair-wise differences between the series, thus requiring no assumed values for the “truth” signal. From the variances of the differenced series, the uncertainty of each series can be recovered when reasonable assumptions are made about the correlations between the series. In order to form the pair-wise differences, the series data must be given at the same epoch. All measurement data sets studied here were sampled at noon (UTC); except for the VLBI series, whose data are interpolated to noon and whose UT1 values are also numerically differentiated to obtain LOD. The numerical error introduced to the VLBI values by the interpolation and differentiation is shown to be comparable in magnitude to the values determined by the TCH method for the uncertainties of the VLBI series. The TCH estimates for the VLBI series are corrupted by such numerical errors mostly as a result of the relatively large data intervals. Of the remaining data sets studied here, it is found that the IGS Final combined series has the smallest polar motion and length-of-day uncertainties.     

Outlier detection algorithms and their performance in GOCE gravity field processing

The satellite missions CHAMP, GRACE, and GOCE mark the beginning of a new era in gravity field determination and modeling. They provide unique models of the global stationary gravity field and its variation in time. Due to inevitable measurement errors, sophisticated pre-processing steps have to be applied before further use of the satellite measurements. In the framework of the GOCE mission, this includes outlier detection, absolute calibration and validation of the SGG (satellite gravity gradiometry) measurements, and removal of temporal effects. In general, outliers are defined as observations that appear to be inconsistent with the remainder of the data set. One goal is to evaluate the effect of additive, innovative and bulk outliers on the estimates of the spherical harmonic coefficients. It can be shown that even a small number of undetected outliers (<0.2 of all data points) can have an adverse effect on the coefficient estimates. Consequently, concepts for the identification and removal of outliers have to be developed. Novel outlier detection algorithms are derived and statistical methods are presented that may be used for this purpose. The methods aim at high outlier identification rates as well as small failure rates. A combined algorithm, based on wavelets and a statistical method, shows best performance with an identification rate of about 99%. To further reduce the influence of undetected outliers, an outlier detection algorithm is implemented inside the gravity field solver (the Quick-Look Gravity Field Analysis tool was used). This results in spherical harmonic coefficient estimates that are of similar quality to those obtained without outliers in the input data.