Copper-gold endoskarns and high-Mg monzodiorite–tonalite intrusions at Mt. Shea,Kalgoorlie, Australia: implications for the origin of gold–pyrite–tennantite mineralization in the Golden Mile

Five Cu–Au epidote skarns are associated with the Mt. Shea intrusive complex, located in the 2.7–2.6 Ga Eastern Goldfields Province of the Archean Yilgarn craton, in greenstones bounded by the Boulder Lefroy and Golden Mile strike-slip faults, which control the Golden Mile (1,435 t Au) at Kalgoorlie and smaller “orogenic” gold deposits at Kambalda. The Cu–Au deposits studied are oxidized endoskarns replacing faulted and fractured quartz monzodiorite–granodiorite. The orebodies are up to 140 m long and 40 m thick. Typical grades are 0.5% Cu and 0.3 g/t Au although parts are richer in gold (1.5–4.5 g/t). At the Hannan South mine, the skarns consist of epidote, calcite, chlorite, magnetite (5–15%), and minor quartz, muscovite, and microcline. Gangue and magnetite are in equilibrium contact with pyrite and chalcopyrite. The As–Co–Ni-bearing pyrite contains inclusions of hematite, gold, and electrum and is intergrown with cobaltite and Cu–Pb–Bi sulfides. At the Shea prospect, massive, net-textured, and breccia skarns are composed of multistage epidote, actinolite, albite, magnetite (5%), and minor biotite, calcite, and quartz. Gangue and magnetite are in equilibrium with Co–Ni pyrite and chalcopyrite. Mineral-pair thermometry, mass-balance calculations, and stable-isotope data (pyrite δ³⁴S_CDT = 2.5‰, calcite δ¹³C_PDB = −5.3‰, and δ¹⁸O_SMOW = 12.9‰) indicate that the Cu–Au skarns formed at 500 ± 50°C by intense Ca–Fe–CO₂–S metasomatism from fluids marked by an igneous isotope signature. The Mt. Shea stock–dike–sill complex postdates the regional D1 folding and metamorphism and the main phase of D2 strike-slip faulting. The suite is calc-akaline and comprises hornblende–plagioclase monzodiorite, quartz monzodiorite, granodiorite, and quartz–plagioclase tonalite porphyry. The intrusions display a wide range in silica content (53–73 wt% SiO₂), in ratio (0.37–0.89), and in ratio (0.02–0.31). Chromium (62–345 ppm), Ni (23–158), Sr (311–1361 ppm), and Ba (250–2,581 ppm) contents are high, Sr/Y ratios are high (24–278, mostly >50), and the rare earth element patterns are fractionated . These features and a negative niobium anomaly relative to the normal mid-ocean ridge basalt indicate that the suite formed by hornblende fractionation from a subduction-related monzodiorite magma sourced from metasomatized peridotite in the upper mantle. The magnesian composition of many intrusions was enhanced due to hornblende crystallization under oxidizing hydrous conditions and during the subsequent destruction of igneous magnetite by subsolidus actinolite–albite alteration. At the Shea prospect, main-stage Cu–Au epidote skarn is cut by biotite–albite–dolomite schist and by red biotite–albite replacement bands. Post-skarn alteration includes 20-m-thick zones of sericite–chlorite–ankerite schist confined to two D3 reverse faults. The schists are mineralized with magnetite + pyrite + chalcopyrite (up to 0.62% Cu, 1.6 g/t Au) and are linked to skarn formation by shared Ca–Fe–CO₂ metasomatism. Red sericitic alteration, marked by magnetite + hematite + pyrite, occurs in fractured porphyry. The biotite/sericite alteration and oxidized ore assemblages at the Shea prospect are mineralogically identical to magnetite–hematite-bearing gold lodes at Kambalda and in the Golden Mile. Published fluid inclusion data suggest that a “high-pressure”, oxidized magmatic fluid (2–9 wt% NaCl equivalent, , 200–400 MPa) was responsible for gold mineralization in structural sites of the Boulder Lefroy and Golden Mile faults. The sericite–alkerite lodes in the Golden Mile share the assemblages pyrite + tennantite + chalcopyrite and bornite + pyrite, and accessory high-sulfidation enargite with late-stage sericitic alteration zones developed above porphyry copper deposits.     

Brecciation and chemical heterogeneities of CI chondrites

Fragments in the size range from 40 μm to several hundred μm in the CI chondrites Orgueil, Ivuna, Alais, and Tonk show a wide range of chemical compositions with variations in major elements such as iron (10.4-42.4 wt% FeO), silicon (12.7-42.2 wt% SiO₂), and sulfur (1.01-15.8 wt% SO₃), but also important minor elements such as phosphorous (up to 5.2 wt% P₂O₅) or calcium (up to 6.6 wt% CaO). These variations are the result of the varying mineralogical compositions of these fragments. The distribution of phyllosilicates, magnetites, and possibly ferrihydrite, in particular, control the abundances of these elements. High REE contents—up to 150 times the solar abundances—were observed in phosphates, while matrix and sulfates are REE-depleted. The studied 113 fragments were subdivided into eight lithologies with similar mineralogical and thus chemical properties. The most common is the CGA lithology, consisting of a groundmass of Mg-rich, coarse-grained phyllosilicates and varying abundances of inclusions such as magnetite. The second most abundant lithology is the FGA lithology, consisting of a groundmass of fine-grained Fe-rich phyllosilicates. A rare, but important lithology consists of fragments with high contents of phosphates and other minerals. The proposed model for the evolution of these lithologies is based on a closed system alteration, where mineralogical differences in the lithologies reflect heterogeneities in the starting material. Comparison of our results with literature data indicates a general similarity of the four CI chondrites analyzed. Further comparison of bulk analyses suggests that the mass ‘threshold’ for chemical heterogeneities in CI chondrite samples is smaller than ∼1-2 g.     

TTG-type plutonic rocks formed in a modern arc batholith by hydrous fractionation in the lower arc crust

We present the geochemistry and intrusion pressures of granitoids from the Kohistan batholith, which represents, together with the intruded volcanic and sedimentary units, the middle and upper arc crust of the Kohistan paleo-island arc. Based on Al-in-hornblende barometry, the batholith records intrusion pressures from ~0.2 GPa in the north (where the volcano-sedimentary cover is intruded) to max. ~0.9 GPa in the southeast. The Al-in-hornblende barometry demonstrates that the Kohistan batholith represents a complete cross section across an arc batholith, reaching from the top at ~8–9 km depth (north) to its bottom at 25–35 km (south-central to southeast). Despite the complete outcropping and accessibility of the entire batholith, there is no observable compositional stratification across the batholith. The geochemical characteristics of the granitoids define three groups. Group 1 is characterized by strongly enriched incompatible elements and unfractionated middle rare earth elements (MREE)/heavy rare earth element patterns (HREE); Group 2 has enriched incompatible element concentrations similar to Group 1 but strongly fractionated MREE/HREE. Group 3 is characterized by only a limited incompatible element enrichment and unfractionated MREE/HREE. The origin of the different groups can be modeled through a relatively hydrous (Group 1 and 2) and of a less hydrous (Group 3) fractional crystallization line from a primitive basaltic parent at different pressures. Appropriate mafic/ultramafic cumulates that explain the chemical characteristics of each group are preserved at the base of the arc. The Kohistan batholith strengthens the conclusion that hydrous fractionation is the most important mechanism to form volumetrically significant amounts of granitoids in arcs. The Kohistan Group 2 granitoids have essentially identical trace element characteristics as Archean tonalite–trondhjemite–granodiorite (TTG) suites. Based on these observations, it is most likely that similar to the Group 2 rocks in the Kohistan arc, TTG gneisses were to a large part formed by hydrous high-pressure differentiation of primitive arc magmas in subduction zones.     

Damage phenomenon and petrophysical properties of sandstones at the Phnom Bakheng Temple (Angkor,Cambodia): first investigations and possible conservation measures

Sandstones, clay in the form of bricks and laterite are the building materials used by the Khmer to construct the imposing and magnificent temples in Southeast Asia. Many of these monuments suffer from fracturing, sanding, contour scaling, crust formation and salt weathering. The affinity to weathering is closely connected to the type of material. Two sandstone types classified as feldspathic arenite and quartz arenite of Angkor as well as two arkosic sandstones from Thailand are described and investigated in this study. Important petrophysical properties determined for the different sandstones consist of hydric expansion, thermal expansion, pore radii distribution and ultrasonic velocity. Different investigations such as capillary water uptake, surface hardness, hygroscopic water sorption, and salt resistance tests were undertaken in the laboratory to characterize the various rock types. Observations and quantified damage mapping were done onsite at the Phnom Bakheng Temple. Contour scaling in the form of weathering crusts is one of the main deterioration features observable at the Angkor monuments. Comparisons are made between the building stone, the crust material from the Phnom Bakheng Temple and fresh stone material used for restoration. Significant differences in hydric and especially in thermal expansion of the crust and sandstone have been determined. The results seem to indicate that extensional processes occur, which can be considered a force for detachment (i.e., contour scaling, flaking). In an experimental trial, the hydric and thermal expansion of the weathering crust and the building stone was significantly reduced by using a weak acid for the crust and a swelling inhibitor for the original building stone.     

The Namuskluft and Dreigratberg sections in southern Namibia (Kalahari Craton,Gariep Belt): a geological history of Neoproterozoic rifting and recycling of cratonic crust during the dispersal of Rodinia until the amalgamation of Gondwana

This paper presents combined U/Pb, Th/U and Hf isotope analyses on detrital and magmatic zircon grains together with whole-rock geochemical analyses of two basement and eight sedimentary rock samples from the Namuskluft and the Dreigratberg in southern Namibia (Gariep Belt). The sedimentary sections evolved during the Cryogenian on the SW part of the Kalahari Craton and where therefore deposited in an active rift setting during the break-up of Rodinia. Due to insufficient palaeomagnetic data, the position of the Kalahari Craton within Rodinia is still under discussion. There are possibilities to locate Kalahari along the western side of Australia/Mawsonland (Pisarevski et al. in Proterozoic East Gondwana: supercontinent assembly and break-up, Geological Society, London, 2003; Evans in Ancient Orogens and modern analogues. Geological Society, London, 2009; and others) or together with the Congo-Sao Francisco and Rio de la Plata Cratons (Li et al. in Prec Res 45: 203–2014, 2008; Frimmel et al. in Int J Earth Sci (Geol Rundsch) 100: 323–354, 2011; and others). It is sill unclear which craton rifted away from the Kalahari Craton during the Cryogenian. Although Middle to Upper Cryogenian magmatic activity is known for the SE Kalahari Craton (our working area) (Richtersveld Suite, Rosh Pinah Fm), all the presented samples show no U/Pb zircon ages younger than ca. 1.0 Ga and non-older than 2.06 Ga. The obtained U/Pb ages fit very well to the exposed basement of the Kalahari Craton (1.0–1.4 Ga Namaqua Province, 1.7–2.0 Ga Vioolsdrif Granite Suite and Orange River Group) and allow no correlation with a foreign craton such as the Rio de la Plata or Australia/Mawsonland. Lu–Hf isotopic signatures of detrital zircon point to the recycling of mainly Palaeoproterozoic and to a smaller amount of Archean crust in the source areas. εHf(t) signatures range between ?24 and +14.8, which relate to T_DM model ages between 1.05 and 3.1 Ga. Only few detrital zircon grains derived from magmas generated from Mesoproterozoic crustal material show more juvenile εHf(t) signatures of +14, +8 to +4 with T_DM model ages of 1.05–1.6 Ga. During Neoproterozoic deposition, only old cratonic crust with an inherited continental arc signature was available in the source area clearly demonstrated by Hf isotope composition of detrital zircon and geochemical bulk analysis of sedimentary rocks. The granodiorites of the Palaeoproterozoic basement underlying Namuskluft section are ca. 1.9 Ga old and show εHf(t) signatures of ?3 to ?5.5 with T_DM model ages of 2.4–2.7 Ga. These basement rocks demonstrate the extreme uplift and deep erosion of the underlying Kalahari Craton at its western margin before general subsidence during Cryogenian and Ediacaran time. The sedimentary sequence of the two examined sections (Namuskluft and Dreigratberg) proposes the presence of a basin and an increasing subsidence at the SW part of the Kalahari Craton during the Cryogenian. Therefore, we propose the initial formation of an intra-cratonic sag basin during the Lower Cryogenian that evolved later to a rift basin at the cratonic margin due to increasing crustal tension and rifting together with the opening of the Adamastor Ocean. As the zircons of the sedimentary rocks filling this basin show neither rift-related U/Pb ages nor an exotic craton as a possible source area, the only plausible sedimentary transport direction providing the found U/Pb ages would be from the E or the SE, directly from the heart of the Kalahari Craton. Due to subsidence and ongoing sedimentation from E/SE directions, the rift-related magmatic rocks were simply covered by the input of old intra-cratonic material that explains the absence of Neoproterozoic zircon grains in our samples. The geochemical analyses show the erosion of a continental arc and related sedimentary rocks with an overall felsic provenance. The source area was a deeply eroded and incised magmatic arc that evolved on continental crust, without any evidence for a passive margin. All of this can be explained by the erosion of rocks related to the Namaqua Belt, which represents one of the two major peaks of zircon U–Pb ages in all analysed samples. Therefore, the Namaqua Belt was well exposed during the Cryogenian, available to erosion and apart from the also well-exposed Palaeoproterozoic basement of the Kalahari Craton one potential source area for the sedimentary rocks in the investigated areas.     

Larval dispersal and vertical migration behaviour – a simulation study for short dispersal times

Current speed often varies with depth, so vertical movements of larvae are expected to have profound effects on dispersal velocity and therefore dispersal potential. Systematic behaviours are expected to have strong effects on dispersal. However, reliable information on the presence of vertical migrations in larvae is scarce, but the few well investigated empirical examples justify a detailed simulation study and an analysis of potential effects. We present a spatially explicit 3D hydrodynamic model that incorporates biological information in the form of active particles advected in a Lagrangian fashion. The set‐up is designed to analyze the sensitivity of dispersal distances to variation in vertical behaviour of larvae. We simulated short (4 days) pelagic larval durations (PLDs) to determine whether behaviour might be important over short dispersal periods. We found that sinusoidal behaviours (slow vertical migration) in or out of phase with tides did not significantly change the dispersal patterns compared to those of larvae that remained at the surface. By contrast, a quadratic pattern of behaviour resulting in rapid vertical migration, in or out of phase with tides, had dramatic effects on both distance and direction of dispersal. The resulting dispersal kernels were found to be multimodal due to the interaction between tidal and meteorological components in flow. Incorporating biological information on larval migrations in Lagrangian simulation of dispersal will be important in estimates of connectivity and forecasting marine reserve networks.     

On the effect of topography on surface wave propagation in the ambient noise frequency range

Due to the increasing popularity of analyzing empirical Green’s functions obtained from ambient seismic noise, more and more regional tomographical studies based on short-period surface waves are published. Results could potentially be biased in mountainous regions where topography is not small compared to the wavelength and penetration depth of the considered waves. We investigate the effect of topography on the propagation of short-period Rayleigh waves empirically by means of synthetic data using a spectral element code and a 3-D model with real topography. We show that topography along a profile through the studied area can result in an underestimation of phase velocities of up to about 0.7% at the shortest investigated period (3 s). Contrary to the expectation that this bias results from the increased surface distance along topography, we find that this error can be estimated by local topographic contrasts in the vicinity of the receiver alone. We discuss and generalize our results by considering topographic profiles through other mountain ranges and find that southern Norway is a good proxy to assess the topography effect. Nevertheless, topographic bias on phase velocity measurements is in general not large enough to significantly affect recovered velocity variations in the ambient noise frequency range.     

Bedload transport measurements at the Erlenbach stream with geophones and automated basket samplers

In the Erlenbach stream, a pre‐alpine steep channel in Switzerland, sediment transport has been monitored for more than 25 years. Near the confluence with the main valley river, stream flow is monitored and sediment is collected in a retention basin with a capacity of about 2000 m³. The basin is surveyed at regular intervals and after large flood events. In addition, sediment transport has been continuously monitored with piezoelectric bedload impact and geophone sensors since 1986. In 2008–2009, the measuring system in the Erlenbach stream was enhanced by installing an automatic system to obtain bedload samples. Movable metal baskets are mounted on a rail at the downstream wall of the large check dam above the retention basin, and they can be moved automatically into the flow to take bedload transport samples. The wire mesh of the baskets has a spacing of 10 mm to sample all sediment particles coarser than this size (which is about the limiting grain size detected by the geophones). The upgraded measuring system permits to obtain bedload samples over short sampling periods and to measure the grain size distribution of the transported material and its variation over time and with discharge. The analysis of calibration relationships for the geophone measuring system confirms findings from very similar measurements which were performed until 1999 with piezoelectric bedload impact sensors; there is a linear relationship between impulse counts and bedload mass passing over the sensors. Findings from flume experiments are used to discuss the most important factors which affect the calibration of the geophone signal. The bedload transport rates as measured by the moving baskets are among the highest measured in natural streams, with values of the order of several kilograms per meter per second. Copyright © 2012 John Wiley & Sons, Ltd.     

Optimal spherical spline filters for the analysis and comparison of regional-scale tomographic models

Advances in seismic tomography lead to increasingly detailed models of the Earth that are often represented on irregular and resolution-adaptive grids. To take full advantage of such models, their assessment must progress beyond a purely visual analysis, and tools must become available for their quantitative comparison.We present a method for the spectral analysis and comparison of multi-scale tomographic models. The method is applicable to irregular grids on the sphere, and is more efficient that filters based on spherical-harmonic expansions or convolution integrals. The combination of a spherical spline representation of tomographic information with Abel-Poisson scaling enables the construction of targetted spatial filters by solving a nonlinear inverse problem for appropriate weighting coefficients. This can be readily achieved with a simulated annealing approach for the limited number of weights. Once suitable filters have been generated they can be employed to address issues such as the patterns of small-scale heterogeneity, transitional structures and comparison of independent models from a region.We illustrate our method in a series of applications where we use different bandpass filters to detect differences in the distribution of small-scale heterogeneity beneath central and eastern Europe, and to compare several recent tomographic models of the Australian region.     

Magnetic signature of hydrocarbon-contaminated soils and sediments at the former oil field H?nigsen, Germany

Magnetic properties of hydrocarbon (HC) containing soils and sediments from two sites (Site A and B) of the former oil-field H?nigsen were analyzed in order to determine whether magnetic methods can be employed to delineate HC contamination of soils and sediments. Magnetic parameters such as magnetic susceptibility and induced isothermal remanent magnetizations, as well as soil and sediment properties such as pH, iron content and water content, HC content and most probable number counts of iron-metabolizing microorganisms were determined. The magnetic concentration-dependent parameters for HC contaminated samples were 25 times higher in soils from Site A than in sediment samples from Site B. However, at Site B the magnetic susceptibility was still four times higher in comparison to lithologically similar non-contaminated sediment samples from a third Site C. Newly formed magnetite containing mainly single domain particles was responsible for the magnetic enhancement, whereas superparamagnetic grains represented only a minor component. Site A had an acidic pH compared to neutral pH at Site B, and a higher crystalline and bioavailable total iron content. Nevertheless, Site B samples contained significant numbers of both iron(II)-oxidizing and iron(III)-reducing microorganisms indicating that microbial iron cycling might have taken place at this site and potentially played a role for iron mineral transformation, including magnetite (trans)formation. The content of total non-polar hydrocarbons (TNPH) at Site A was one order of magnitude higher than at Site B. Only at Site A magnetic susceptibility correlated well with TNPH. Our results demonstrate that HC contaminated samples had an enhanced magnetite content compared to non-contaminated soils and sediments. Therefore, magnetic methods may provide a quick and cost-effective way to assess HC contamination in soils and sediments. However, more field sites and laboratory investigations are needed to reveal the complex nature of the processes involved.