Genetic significance of fluid inclusions in the CSA Cu–Pb–Zn deposit, Cobar, Australia

CSA mine exploits a ‘Cobar-type’ Cu–Pb–Zn±Au±Ag deposit within a cleaved and metamorphosed portion of the Cobar Supergroup, central New South Wales. The deposit comprises systems of ‘lenses’ that encompass veins, disseminations and semi-massive to massive Cu–Pb–Zn ores. The systems and contained lenses truncate bedding, are approximately coplanar with regional cleavage and similarly oriented shear zones and plunge parallel to the elongation lineation. Systems have extreme vertical continuity (>1000 m), short strike length (400 m) and narrow width (100 m), exhibit vertical and lateral ore-type variation and have alteration haloes. Models of ore formation include classical hydrothermalism, structurally controlled remobilisation and polymodal concepts; syntectonic emplacement now holds sway.Fluid inclusions were examined from quartz±sulphide veins adjacent to now-extracted ore, from coexisting quartz–sulphide within ore, and from vughs in barren quartz veins. Lack of early primary inclusions precluded direct determination of fluids associated with D₂–D₃ ore and vein emplacement. Similarly, decrepitation (by near-isobaric heating) of the two oldest secondary populations precluded direct determination of fluid phases immediately following D₂–D₃ ore and vein emplacement. Post-decrepitation outflow (late D₃ to early post-D₃) is recorded by monophase CH₄ inclusions. Entrained outflow of deeply circulated meteoric fluid modified the CH₄ system; modification is recorded by H₂O+CH₄ and H₂O+(trace CH₄) secondary populations and by an H₂O+(trace CH₄) primary population. The contractional tectonics (D₂–D₃) of ore emplacement was superseded by relaxational tectonics (D_4P) that facilitated meteoric water penetration and return flow.Under D₂ prograde metamorphism, entrapment temperatures (Tt) and pressures (Pt) for pre-decrepitation secondary inclusions are estimated as Tt300–330 °C and Pt1.5–2 kbar≈Plith (the lithostatic pressure). Decrepitation accompanied peak metamorphism (T350–380 °C) in mid- to late-D₃, while in late-D₃ to early post-D₃, essentially monophase CH₄ secondary inclusions were entrapped at Tt350 °C and Pt=1.5–2 kbar≈Plith. Subsequently, abundant CH₄ and entrained meteoric water were entrapped as H₂O+CH₄ secondaries under slowly decreasing temperature (Tt330–350 °C) and constant pressure (Pt1.5–2 kbar). Finally, with increasingly dominant meteoric outflow, H₂O+(trace CH₄) populations record decreasing temperatures (Tt>300 to <350 down to 275–300 °C) at pressures of Phydrostatic<Pt (1 kbar) <Plith (1.5 kbar).The populations of inclusions provide insight into fluid types, flow regimes and P–T conditions during parts of the deposit's evolution. They indirectly support the role of basin-derived CH₄ fluids in ore formation, but provide no insight into a basement-sourced ore-forming fluid. They fully support post-ore involvement of meteoric water. The poorly constrained entrapment history is believed to span 10 Ma from 395 to 385 Ma.     

Lu–Hf and U–Pb isotope systematics of zircons from the Storgangen intrusion, Rogaland Intrusive Complex, SW Norway: implications for the composition and evolution of Precambrian lower crust in the Baltic Shield

The Storgangen orebody is a concordantly layered, sill-like body of ilmenite-rich norite, intruding anorthosites of the Rogaland Intrusive Complex (RIC), SW Norway. 17 zircon grains were separated from ca. 5 kg of sand-size flotation waste collected from the on-site repository from ilmenite mining. These zircons were analysed for major and trace elements by electron microprobe, and for U–Pb and Lu–Hf isotopes by laser ablation microprobe plasma source mass spectrometry. Eight of the zircons define a well-constrained (MSWD=0.37) concordant population with an age of 949±7 Ma, which is significantly older than the 920–930 Ma ages previously reported for zircon inclusions in orthopyroxene megacrysts from the RIC. The remaining zircons, interpreted as inherited grains, show a range of ²⁰⁷Pb/²⁰⁶Pb ages up to 1407±14 Ma, with an upper intercept age at ca. 1520 Ma. The concordant zircons have similar trace element patterns, and a mean initial Hf isotope composition of ¹⁷⁶Hf/¹⁷⁷Hf_{949 Ma}=0.28223±5 (_Hf=+2±2). This is similar to the Hf-isotope composition of zircons in a range of post-tectonic Sveconorwegian granites from South Norway, and slightly more radiogenic than expected for mid-Proterozoic juvenile crust. The older, inherited zircons show Lu–Hf crustal residence ages in the range 1.85–2.04 Ga. One (undated) zircon plots well within the field of Hf isotope evolution of Paleoproterozoic rocks of the Baltic Shield. These findings indicate the presence of Paleoproterozoic components in the deep crust of the Rogaland area, but do not demonstrate that such rocks, or a Sveconorwegian mantle-derived component, contributed significantly to the petrogenesis of the RIC. If the parent magma was derived from a homogeneous, lower crustal mafic granulite source, the lower crustal protolith must be at least 1.5 Ga old, and it must have an elevated Rb/Sr ratio. This component would be indistinguishable in Sr, Nd and Hf isotopes from some intermediate mixtures between Sveconorwegian mantle and Paleoprotoerzoic felsic crust, but it cannot account for the initial ¹⁴³Nd/¹⁴⁴Nd of the most primitive, late Sveconorwegian granite in the region, without the addition of mantle-derived material.     

The geology of kimberlite pipes of the Ekati property, Northwest Territories, Canada

This paper reviews key characteristics of kimberlites on the Ekati property, NWT, Canada. To date 150 kimberlites have been discovered on the property, five of which are mined for diamonds. The kimberlites intrude Archean basement of the central Slave craton. Numerous Proterozoic diabase dykes intrude the area. The Precambrian rocks are overlain by Quaternary glacial sediments. No Phanerozoic rocks are present. However, mudstone xenoliths and disaggregated sediment within the kimberlites indicate that late-Cretaceous and Tertiary cover (likely <200 m) was present at the time of emplacement. The Ekati kimberlites range in age from 45 to 75 Ma. They are mostly small pipe-like bodies (surface area mostly <3 ha but up to 20 ha) that typically extend to projected depths of 400–600 m below current surface. Pipe morphologies are strongly controlled by joints and faults. The kimberlites consist primarily of variably bedded volcaniclastic kimberlite (VK). This is dominated by juvenile constituents (olivine and lesser kimberlitic ash) and variable amounts of exotic sediment (primarily mud), with minor amounts of xenolithic wall-rock material (generally <5%). Kimberlite types include: mud-rich resedimented VK (mRVK); olivine-rich VK (oVK); sedimentary kimberlite; primary VK (PVK); tuffisitic kimberlite (TK) and magmatic kimberlite (MK). The presence and arrangement of these rock types varies widely. The majority of bodies are dominated by oVK and mRVK, but PVK is prominent in the lower portions of certain kimberlites. TK is rare. MK occurs primarily as precursor dykes but, in a few cases, forms pipe-filling intrusions. The internal geology of the kimberlites ranges from simple single-phase pipes (RVK or MK), to complex bodies with multiple, distinct units of VK. The latter include pipes infilled with steep, irregular VK blocks/wedges and at least one case in which the pipe is occupied by well-defined sub-horizontal VK phases, including a unique, 100-m-thick graded sequence. The whole-rock compositions of VK samples suggest significant loss of kimberlitic fines during eruption followed by variable dilution by surface sediment and concurrent incorporation of kimberlitic ash. Diamond distribution within the kimberlites reflects the amount and nature of mantle material sampled by individual kimberlite phases, but is modified considerably by eruption and depositional processes. The characteristics of the Ekati kimberlites are consistent with a two-stage emplacement process: (1) explosive eruption/s causing vent clearing followed by formation of a significant tephra rim/cone of highly fragmented, olivine-enriched juvenile material with varying amounts of kimberlitic ash and surface sediments (predominantly mud); and (2) infilling of the vent by direct deposition from the eruption column and/or resedimentation of crater rim materials. The presence of less fragmented, juvenile-rich PVK in the lower portions of certain pipes and the intrusion of large volumes of MK to shallow levels in some bodies suggest emplacement of relatively volatile-depleted, less explosive kimberlite in the later stages of pipe formation and/or filling. Explosive devolatilisation of CO₂-rich kimberlite magma is interpreted to have been the dominant eruption mechanism, but phreatomagmatism is thought to have played a role and, in certain cases, may have been dominant.     

Inclusions in diamonds from the K14 and K10 kimberlites, Buffalo Hills, Alberta, Canada: diamond growth in a plume?

Analyses of mineral inclusions, carbon isotopes, nitrogen contents and nitrogen aggregation states in 29 diamonds from two Buffalo Hills kimberlites in northern Alberta, Canada were conducted. From 25 inclusion bearing diamonds, the following paragenetic abundances were found: peridotitic (48%), eclogitic (32%), eclogitic/websteritic (8%), websteritic (4%), ultradeep? (4%) and unknown (4%). Diamonds containing mineral inclusions of ferropericlase, and mixed eclogitic-asthenospheric-websteritic and eclogitic-websteritic mineral associations suggests the possibility of diamond growth over a range of depths and in a variety of mantle environments (lithosphere, asthenosphere and possibly lower mantle).

Eclogitic diamonds have a broad range of C-isotopic composition (δ¹³C=−21‰ to −5‰). Peridotitic, websteritic and ultradeep diamonds have typical mantle C-isotope values (δ¹³C=−4.9‰ av.), except for two ¹³C-depleted peridotitic (δ¹³C=−11.8‰, −14.6‰) and one ¹³C-depleted websteritic diamond (δ¹³C=−11.9‰). Infrared spectra from 29 diamonds identified two diamond groups: 75% are nitrogen-free (Type II) or have fully aggregated nitrogen defects (Type IaB) with platelet degradation and low to moderate nitrogen contents (av. 330 ppm-N); 25% have lower nitrogen aggregation states and higher nitrogen contents (30% IaB; <1600 ppm-N).

The combined evidence suggests two generations of diamond growth. Type II and Type IaB diamonds with ultradeep, peridotitic, eclogitic and websteritic inclusions crystallised from eclogitic and peridotitic rocks while moving in a dynamic environment from the asthenosphere and possibly the lower mantle to the base of the lithosphere. Mechanisms for diamond movement through the mantle could be by mantle convection, or an ascending plume. The interaction of partial melts with eclogitic and peridotitic lithologies may have produced the intermediate websteritic inclusion compositions, and can explain diamonds of mixed parageneses, and the overlap in C-isotope values between parageneses. Strong deformation and extremely high nitrogen aggregation states in some diamonds may indicate high mantle storage temperatures and strain in the diamond growth environment. A second diamond group, with Type IaA–IaB nitrogen aggregation and peridotitic inclusions, crystallised at the base of the cratonic lithosphere. All diamonds were subsequently sampled by kimberlites and transported to the Earth's surface.     

Genesis and evolution of the lithospheric mantle beneath the Buffalo Head Terrane, Alberta (Canada)

Mantle xenoliths and xenocrysts were retrieved from three of the 88–86 Ma Buffalo Hills kimberlites (K6, K11, K14) for a reconnaissance study of the subcontinental lithospheric mantle (SCLM) beneath the Buffalo Head Terrane (Alberta, Canada). The xenoliths include spinel lherzolites, one garnet spinel lherzolite, garnet harzburgites, one sheared garnet lherzolite and pyroxenites. Pyroxenitic and wehrlitic garnet xenocrysts are derived primarily from the shallow mantle and lherzolitic garnet xenocrysts from the deep mantle. Harzburgite with Ca-saturated garnets is concentrated in a layer between 135–165 km depth. Garnet xenocrysts define a model conductive paleogeotherm corresponding to a heat flow of 38–39 mW/m². The sheared garnet lherzolite lies on an inflection of this geotherm and may constrain the depth of the lithosphere–asthenosphere boundary (LAB) beneath this region to ca 180 km depth.

A loss of >20% partial melt is recorded by spinel lherzolites and up to 60% by the garnet harzburgites, which may be related to lithosphere formation. The mantle was subsequently modified during at least two metasomatic events. An older metasomatic event is evident in incompatible-element enrichments in homogeneous equilibrated garnet and clinopyroxene. Silicate melt metasomatism predominated in the deep lithosphere and led to enrichments in the HFSE with minor enrichments in LREE. Metasomatism by small-volume volatile-rich melts, such as carbonatite, appears to have been more important in the shallow lithosphere and led to enrichments in LREE with minor enrichments in HFSE. An intermediate metasomatic style, possibly a signature of volatile-rich silicate melts, is also recognised. These metasomatic styles may be related through modification of a single melt during progressive interaction with the mantle. This metasomatism is suggested to have occurred during Paleoproterozoic rifting of the Buffalo Head Terrane from the neighbouring Rae Province and may be responsible for the evolution of some samples toward unradiogenic Nd and Hf isotopic compositions.

Disturbed Re–Os isotope systematics, evident in implausible model ages, were obtained in situ for sulfides in several spinel lherzolites and suggest that many sulfides are secondary (metasomatic) or mixtures of primary and secondary sulfides. Sulfide in one peridotite has unradiogenic ¹⁸⁷Os/¹⁸⁸Os and gives a model age of 1.89±0.38 Ga. This age coincides with the inferred emplacement of mafic sheets in the crust and suggests that the melts parental to the intrusions interacted with the lithospheric mantle.

A younger metasomatic event is indicated by the occurrence of sulfide-rich melt patches, unequilibrated mineral compositions and overgrowths on spinel that are Ti-, Cr- and Fe-rich but Zn-poor. Subsequent cooling is recorded by fine exsolution lamellae in the pyroxenes and by arrested mineral reactions.

If the lithosphere beneath the Buffalo Head Terrane was formed in the Archaean, any unambiguous signatures of this ancient origin may have been obliterated during these multiple events.     

Inference and uncertainty of snow depth spatial distribution at the kilometre scale in the Colorado Rocky Mountains: the effects of sample size,random sampling,predictor quality,and validation procedures

Historically, observing snow depth over large areas has been difficult. When snow depth observations are sparse, regression models can be used to infer the snow depth over a given area. Data sparsity has also left many important questions about such inference unexamined. Improved inference, or estimation, of snow depth and its spatial distribution from a given set of observations can benefit a wide range of applications from water resource management, to ecological studies, to validation of satellite estimates of snow pack. The development of Light Detection and Ranging (LiDAR) technology has provided non‐sparse snow depth measurements, which we use in this study, to address fundamental questions about snow depth inference using both sparse and non‐sparse observations. For example, when are more data needed and when are data redundant? Results apply to both traditional and manual snow depth measurements and to LiDAR observations. Through sampling experiments on high‐resolution LiDAR snow depth observations at six separate 1.17‐km² sites in the Colorado Rocky Mountains, we provide novel perspectives on a variety of issues affecting the regression estimation of snow depth from sparse observations. We measure the effects of observation count, random selection of observations, quality of predictor variables, and cross‐validation procedures using three skill metrics: percent error in total snow volume, root mean squared error (RMSE), and R². Extremes of predictor quality are used to understand the range of its effect; how do predictors downloaded from internet perform against more accurate predictors measured by LiDAR? Whereas cross validation remains the only option for validating inference from sparse observations, in our experiments, the full set of LiDAR‐measured snow depths can be considered the ‘true’ spatial distribution and used to understand cross‐validation bias at the spatial scale of inference. We model at the 30‐m resolution of readily available predictors, which is a popular spatial resolution in the literature. Three regression models are also compared, and we briefly examine how sampling design affects model skill. Results quantify the primary dependence of each skill metric on observation count that ranges over three orders of magnitude, doubling at each step from 25 up to 3200. Whereas uncertainty (resulting from random selection of observations) in percent error of true total snow volume is typically well constrained by 100–200 observations, there is considerable uncertainty in the inferred spatial distribution (R²) even at medium observation counts (200–800). We show that percent error in total snow volume is not sensitive to predictor quality, although RMSE and R² (measures of spatial distribution) often depend critically on it. Inaccuracies of downloaded predictors (most often the vegetation predictors) can easily require a quadrupling of observation count to match RMSE and R² scores obtained by LiDAR‐measured predictors. Under cross validation, the RMSE and R² skill measures are consistently biased towards poorer results than their true validations. This is primarily a result of greater variance at the spatial scales of point observations used for cross validation than at the 30‐m resolution of the model. The magnitude of this bias depends on individual site characteristics, observation count (for our experimental design), and sampling design. Sampling designs that maximize independent information maximize cross‐validation bias but also maximize true R². The bagging tree model is found to generally outperform the other regression models in the study on several criteria. Finally, we discuss and recommend use of LiDAR in conjunction with regression modelling to advance understanding of snow depth spatial distribution at spatial scales of thousands of square kilometres. Copyright © 2012 John Wiley & Sons, Ltd.     

The impact of rural land management changes on soil hydraulic properties and runoff processes: results from experimental plots in upland UK

To develop an evidence base to help predict the impacts of land management change on flood generation, four experimental sites were established on improved grassland used for sheep grazing at the Pontbren catchment in upland Wales, UK. At each site, three plots were established where surface runoff was measured, supplemented by measurements of soil infiltration rates and soil and vegetation physical properties. Following baseline monitoring, treatments were applied to two of the plots: exclusion of sheep (ungrazed) and exclusion of sheep and planting with native broadleaf tree species (tree planted), with the third plot acting as a control (grazed pasture). Due to a particularly dry summer that occurred pre‐treatment, the soil hydrological responses were initially impacted by the effects of the climate on soil structure. Nevertheless, treatments did have a clear influence on soil hydrological response. On average, post‐treatment runoff volumes were reduced by 48% and 78% in ungrazed and tree‐planted plots relative to the control, although all results varied greatly over the sites. Five years following treatment application, near‐surface soil bulk density was reduced and median soil infiltration rates were 67 times greater in plots planted with trees compared to grazed pasture. The results illustrate the potential use of upland land management for ameliorating local‐scale flood generation but emphasise the need for long‐term monitoring to more clearly separate the effects of land management from those of climatic variability. Copyright © 2013 John Wiley & Sons, Ltd.