Metamorphism and PGE-Au content of chromitite from the Ipanema mafic/ultramafic Complex, Minas Gerais, Brazil

The Santa Cruz massif, which forms part of the Ipanema mafic/ultramafic Complex, Minas Gerais, Brazil, has an exposed upward sequence of metadunite, metaharzburgite (including three separate chromitite layers), metapyroxenite, metagabbro, and metaanorthosite. Primary igneous chromite grains in the main chromitite layer are poikiloblastic and tectonically fragmented, and have a narrow (10-20 mum) margin of chromian spinel. Cataclased chromite fragments are extensively replaced and mantled by chromian spinel; they have a composite margin comprised of an inner zone of more aluminous spinel and an euhedral outer zone of more Cr-rich spinel, representing granulite and amphibolite facies metamorphic events, respectively. The contents of platinum-group elements (PGE) and Au in chromite separates are relatively high (Os 45, Ir 23, Ru 136, Rh 19, Pt 98, Pd 63, and Au 83 ppb), and significantly enriched (similar to 4x) over whole rock values. Platinum-group minerals are not observed and micrometre-sized inclusions of sulfide minerals (chalcopyrite and pentlandite) in relict chromite are rare. However, comparison of mineral proportions in the separated chromite and whole rock shows that the precious metals are hosted predominantly in the relict igneous chromite grains, rather than the secondary chromian spinel and primary and secondary Mg-rich silicates. The major element composition and average chondrite-normalized PGE pattern of the separated chromite correspond to S-poor stratiform chromitite. We suggest that the precious metals accumulated with chromite during crystallization of a S-poor magma, and were not remobilized in the relict chromite during the subsequent high grade metamorphism.

The growth of large mafic intrusions: Comparing Niquelandia and Ivrea igneous complexes

The Niquelandia Complex, Brazil, is one of the world's largest mafic-ultramafic plutonic complexes. Like the Mafic Complex of the Ivrea-Verbano Zone, it is affected by a pervasive high-T foliation and shows hypersolidus deformation structures, contains significant inclusions of country-rock paragneiss, and is subdivided into a Lower and an Upper Complex. In this paper, we present new SHRIMP U-Pb zircon ages that provide compelling evidence that the Upper and the Lower Niquelandia Complexes formed during the same igneous event at ca. 790 Ma. Coexistence of syn-magmatic and high-T subsolidus deformation structures indicates that both complexes grew incrementally as large crystal mush bodies which were continuously stretched while fed by pulses of fresh magma. Syn-magmatic recrystallization during this deformation resulted in textures and structures which, although appearing metamorphic, are not ascribable to post-magmatic metamorphic event(s), but are instead characteristic of the growth process in huge and deep mafic intrusions such as both the Niquelandia and Ivrea Complexes. Melting of incorporated country-rock paragneiss continued producing hybrid rocks during the last, vanishing stages of magmatic crystallization. This resulted in the formation of minor, late-stage hybrid rocks, whose presence obscures the record of the main processes of interaction between mantle magmas and crustal components, which may be active at the peak of the igneous events and lead to the generation of eruptible hybrid magmas. (C) 2012 Elsevier B.V. All rights reserved.

Chemical composition of glass inclusions from ODP Holes 157-953C and 157-956B

We report S concentrations and relative proportions of (SO4)2- and S2- in OL- and CPX-hosted glass inclusions and in host glassy lapilli from Miocene basaltic hyaloclastites drilled north and south of Gran Canaria during ODP Leg 157. Compositions of glass inclusions and lapilli resemble those of subaerial Miocene shield basalts on Gran Canaria and comprise mafic to more evolved tholeiitic to alkali basalt and basanite (10.3-3.7 wt.% MgO, 44.5-56.9 wt.% SiO2). Glass inclusions fall into three groups based on their S concentrations: a high-sulfur group (1050 to 5810 ppm S), an intermediate-sulfur group (510 to 1740 ppm S), and a low-sulfur group (<500 ppm S). The most S-rich inclusions have the highest and nearly constant proportion of sulfur dissolved as sulfate determined by electron microprobe measurements of SKa peak shift. Their average S6+/S_total value is 0.75+/-0.09, unusually high for ocean island basalt magmas. The low-sulfur group inclusions have low S6+/S_total ratios (0.08+/-0.05), whereas intermediate sulfur group inclusions show a wide range of S6+/S_total (0.05-0.83). Glassy lapilli and their crystal-hosted glass inclusions with S concentrations of 50 to 1140 ppm S have very similar S6+/S_total ratios of 0.36+/-0.06 implying that sulfur degassing does not affect the proportion of (SO4)2- and S2- in the magma. The oxygen fugacities estimated from S6+/S_total ratios and from Fe3+/Fe2+ ratios in spinel inclusions range from NNO-1.1 to NNO+1.8. The origin of S-rich magmas is unclear. We discuss (1) partial melting of a mantle source at relatively oxidized fO2 conditions, and (2) magma contamination by seawater either directly or through magma interaction with seawater-altered Jurassic oceanic crust. The intermediate sulfur group inclusions represent undegassed or slightly degassed magmas similar to submarine OIB glasses, whereas the low-sulfur group inclusions are likely to have formed from magmas significantly degassed in near-surface reservoirs. Mixing of these degassed magmas with stored volatile-rich ones or volatile-rich magma replenishing the chamber filled by partially degassed magmas may produce hybrid melts with strongly varying S concentrations and S6+/S_total ratios.

Composition of melt inclusions and age of zircons of plagiogneisses from the Kola Superdeep Borehole

A comprehensive study of melt inclusions and SHRIMP dating of zircons from trondhjemite gneisses of the sequence VIII from the Kola Superdeep Borehole has revealed presence of old primary magmatic crystals with age up to 2887+/-15 Ma. This is not consistent with the previous view, according to which the oldest zircons from the Archean Complex in SG-3 are products of granulite metamorphism. Primary magmatic zircons of early generation (from 2887 to 2842 Ma) formed in deep-seated magma chambers during partial crystallization of CO2-saturated trondhjemite estimates on duration of generation of tonalite-trondhjemite-granite melts through partial melting of mafic rocks.

Investigation of inclusions trapped inside Libyan Desert glass by Raman microscopy

Several specimens of Libyan Desert Glass (LDG), an enigmatic natural glass from Egypt, were subjected to investigation by micro-Raman spectroscopy. The spectra of inclusions inside the LDG samples were successfully measured through the layers of glass and the mineral species were identified on this basis. The presence of cristobalite as typical for high-temperature melt products was confirmed, together with co-existing quartz. TiO2 was determined in two polymorphic species, rutile and anatase. Micro-Raman spectroscopy proved also the presence of minerals unusual for high-temperature glasses such as anhydrite and aragonite.

Timing and source constraints on the relationship between mafic and felsic intrusions in the Emeiashan large igneous province

Several I- and A-type granite, syenite plutons and spatially associated, giant Fe–Ti–V deposit-bearing mafic ultramafic layered intrusions occur in the Pan–Xi(Panzhihua–Xichang) area within the inner zone of the Emeishan large igneous province (ELIP). These complexes are interpreted to be related to the Emeishan mantle plume. We present LA-ICP-MS and SIMS zircon U–Pb ages and Hf–Nd isotopic compositions for the gabbros, syenites and granites from these complexes. The dating shows that the age of the felsic intrusive magmatism (256.2 ± 3.0–259.8 ± 1.6 Ma) is indistinguishable from that of the mafic intrusive magmatism (255.4 ± 3.1–259.5 ± 2.7 Ma) and represents the final phase of a continuous magmatic episode that lasted no more than 10 Myr. The upper gabbros in the mafic–ultramafic intrusions are generally more isotopically enriched (lower eNd and eHf) than the middle and lower gabbros, suggesting that the upper gabbros have experienced a higher level of crustal contamination than the lower gabbros. The significantly positive eHf(t) values of the A-type granites and syenites (+4.9 to +10.8) are higher than those of the upper gabbros of the associated mafic intrusion, which shows that they cannot be derived by fractional crystallization of these bodies. They are however identical to those of the mafic enclaves (+7.0 to +11.4) and middle and lower gabbros, implying that they are cogenetic. We suggest that they were generated by fractionation of large-volume, plume-related basaltic magmas that ponded deep in the crust. The deep-seated magma chamber erupted in two stages: the first near a density minimum in the basaltic fractionation trend and the second during the final stage of fractionation when the magma was a low density Fe-poor, Si-rich felsic magma. The basaltic magmas emplaced in the shallowlevel magma chambers differentiated to form mafic–ultramafic layered intrusions accompanied by a small amount of crustal assimilation through roof melting. Evolved A-type granites (synenites and syenodiorites) were produced dominantly by crystallization in the deep crustal magma chamber. In contrast, the I-type granites have negative eNd(t) [-6.3 to -7.5] and eHf(t) [-1.3 to -6.7] values, with the Nd model ages (T Nd DM2) of 1.63-1.67 Ga and Hf model ages (T Hf DM2) of 1.56-1.58 Ga, suggesting that they were mainly derived from partial melting of Mesoproterozoic crust. In combination with previous studies, this study also shows that plume activity not only gave rise to reworking of ancient crust, but also significant growth of juvenile crust in the center of the ELIP.

Lower mantle diamonds from Rio Soriso (Juina area, Mato Grosso, Brazil)

The morphology, colour, fluorescence, cathodoluminescence, nitrogen content and aggregation state, internal structure and mineral inclusions have been studied for 69 alluvial diamonds from the Rio Soriso (Juina area, Mato Grosso State, Brazil). Nitrogen in most diamonds (53%) is fully aggregated as B centres, but there is also a large proportion of N-free stones (38%). A strong positive correlation between nitrogen and IR-active hydrogen concentrations is observed. The diamonds contain (in order of decreasing abundance) ferropericlase, CaSi-perovskite, magnetite, MgSi-perovskite, pyrrhotite, 'olivine', SiO2, perovskite, tetragonal almandine-pyrope phase and some other minerals represented by single grains. The Rio Soriso diamond suite is subdivided into several subpopulations that originated in upper and lower mantle of ultramafic and mafic compositions, with the largest subgroup forming in the ultramafic lower mantle. Analysed ferropericlase grains are enriched in Fe (Mg#=0.43-0.89), which is ascribed to their origin in the lowermost mantle. The Juina kimberlites may be unique in sampling the material from depths below 1,700 km that ascended in a plume formed at the core-mantle boundary.

Two-scale computational modelling of water flow in unsaturated soils containing irregular-shaped inclusions

The focus of this paper is two-dimensional computational modelling of water flow in unsaturated soils consisting of weakly conductive disconnected inclusions embedded in a highly conductive connected matrix. When the inclusions are small, a two-scale Richards’ equation-based model has been proposed in the literature taking the form of an equation with effective parameters governing the macroscopic flow coupled with a microscopic equation, defined at each point in the macroscopic domain, governing the flow in the inclusions. This paper is devoted to a number of advances in the numerical implementation of this model. Namely, by treating the micro-scale as a two-dimensional problem, our solution approach based on a control volume finite element method can be applied to irregular inclusion geometries, and, if necessary, modified to account for additional phenomena (e.g. imposing the macroscopic gradient on the micro-scale via a linear approximation of the macroscopic variable along the microscopic boundary). This is achieved with the help of an exponential integrator for advancing the solution in time. This time integration method completely avoids generation of the Jacobian matrix of the system and hence eases the computation when solving the two-scale model in a completely coupled manner. Numerical simulations are presented for a two-dimensional infiltration problem.

Stresses around two equal circular elastic inclusions in a pressurised cylindrical shell

The stress problem of two equal circular elastic inclusions in a pressurised cylindrical shell has been solved by using single inclusion solutions together with Graf’s addition theorem. The effect of the inter-inclusion distance on the interface stresses in the shell as well as in the inclusion is studied. The results obtained for small values of curvature parameter fi @*=(a*/8Rt) [12(1-v*)]“*, a, R, t being inclusion radius and shell radius and thickness) when compared with the flat-plate results show good agreement. The results obtained in non-dimensional form are presented graphically.

Magmatic evolution of topaz-bearing granite stocks within the Wiborg rapakivi granite batholith

This study brings new insights into the magmatic evolution of natural F-enriched peraluminous granitic systems. The Artjärvi, Sääskjärvi and Kymi granite stocks within the 1.64 Ga Wiborg rapakivi granite batholith have been investigated by petrographic, geochemical, experimental and melt inclusion methods. These stocks represent late-stage leucocratic and weakly peraluminous intrusive phases typical of rapakivi granites worldwide. The Artjärvi and Sääskjärvi stocks are multiphase intrusions in which the most evolved phase is topaz granite. The Kymi stock contains topaz throughout and has a well-developed zoned structure, from the rim to the center: stockscheider pegmatite equigranular topaz granite porphyritic topaz granite. Geochemically the topaz granites are enriched in F, Li, Be, Ga, Rb, Sn and Nb and depleted in Mg, Fe, Ti, Ba, Sr, Zr and Eu. The anomalous geochemistry and mineralogy of the topaz granites are essentially magmatic in origin; postmagmatic reactions have only slightly modified the compositions. The Kymi equigranular topaz granite shows the most evolved character, and the topaz granites at Artjärvi and Sääskjärvi resemble the less evolved porphyritic topaz granite of the Kymi stock. Stockscheiders are found at the roof contacts of the Artjärvi and Kymi stocks. The stockscheider at Artjärvi is composed of biotite-rich schlieren and pegmatite layers parallel to the contact. The schlieren layering is considered to have formed by velocity-gradient sorting mechanism parallel to the flow, which led to the accumulation of mafic minerals along the upper contact of the topaz granite. Cooling and contraction of the topaz granite formed fractures parallel to the roof contact and residual pegmatite magmas were injected along the fractures and formed the pegmatite layers. The zoned structure of the Kymi stock is the result of intrusion of highly evolved residual melt from deeper parts of the magma chamber along the fractured contact between the porphyritic granite crystal mush and country rock. The equigranular topaz granite and marginal pegmatite (stockscheider) crystallized from this evolved melt. Phase relations of the Kymi equigranular topaz granite have been investigated utilizing crystallization experiments at 100 to 500 MPa as a function of water activity and F content. Fluorite and topaz can crystallize as liquidus phases in F-rich peraluminous systems, but the F content of the melt should exceed 2.5 - 3.0 wt % to facilitate crystallization of topaz. In peraluminous F-bearing melts containing more than 1 wt % F, topaz and muscovite are expected to be the first F-bearing phases to crystallize at high pressure, whereas fluorite and topaz should crystallize first at low pressure. Overall, the saturation of fluorite and topaz follows the reaction: CaAl2Si2O8 (plagioclase) + 2[AlF3]melt = CaF2 (fluorite) + 2Al2SiO4F2 (topaz). The obtained partition coefficient for F between biotite and glass D(F)Bt/glass is 1.89 to 0.80 (average 1.29) and can be used as an empirical fluormeter to determine the F content of coexisting melts. In order to study the magmatic evolution of the Kymi stock, crystallized melt inclusions in quartz and topaz grains in the porphyritic and the equigranular topaz granites and the marginal pegmatite were rehomogenized and analyzed. The homogenization conditions for the melt inclusions from the granites were 700 °C, 300 MPa, and 24 h, and for melt inclusions from the pegmatite, 700 °C, 100 MPa, and 24/96 h. The majority of the melt inclusions is chemically similar to the bulk rocks (excluding H2O content), but a few melt inclusions in the equigranular granite show clearly higher F and low K2O contents (on average 11.6 wt % F, 0.65 wt % K2O). The melt inclusion compositions indicate coexistence of two melt fractions, a prevailing peraluminous and a very volatile-rich, possibly peralkaline. Combined petrological, experimental and melt inclusion studies of the Kymi equigranular topaz granite indicate that plagioclase was the liquidus phase at nearly water-saturated (fluid-saturated) conditions and that the F content of the melt was at least 2 wt %. The early crystallization of biotite and the presence of muscovite in crystallization experiments at 200 MPa contrasts with the late-stage crystallization of biotite and the absence of muscovite in the equigranular granite, indicating that crystallization pressure may have been lower than 200 MPa for the granite.

Ductile deformation of single inclusions in simple shear with a finite-strain hyperelastoviscoplastic rheology

We examine the 2D plane-­strain deformation of initially round, matrix-­bonded, deformable single inclusions in isothermal simple shear using a recently introduced hyperelastoviscoplastic rheology. The broad parameter space spanned by the wide range of effective viscosities, yield stresses, relaxation times, and strain rates encountered in the ductile lithosphere is explored systematically for weak and strong inclusions, the effective viscosity of which varies with respect to the matrix. Most inclusion studies to date focused on elastic or purely viscous rheologies. Comparing our results with linear-­viscous inclusions in a linear-­viscous matrix, we observe significantly different shape evolution of weak and strong inclusions over most of the relevant parameter space. The evolution of inclusion inclination relative to the shear plane is more strongly affected by elastic and plastic contributions to rheology in the case of strong inclusions. In addition, we found that strong inclusions deform in the transient viscoelastic stress regime at high Weissenberg numbers (≥0.01) up to bulk shear strains larger than 3. Studies using the shapes of deformed objects for finite-­strain analysis or viscosity-­ratio estimation should establish carefully which rheology and loading conditions reflect material and deformation properties. We suggest that relatively strong, deformable clasts in shear zones retain stored energy up to fairly high shear strains. Hence, purely viscous models of clast deformation may overlook an important contribution to the energy budget, which may drive dissipation processes within and around natural inclusions.