Geochemistry and Petrogenesis of Mafic and Ultramafic Inclusions in Sublayer and Offset Dikes, Sudbury Igneous Complex, Canada

Abstract

Abstract The c. 1·85 Ga Sudbury Igneous Complex (SIC) is the igneous remnant of one of the oldest, largest and best-preserved impact structures on Earth and contains some of the world’s largest magmatic Ni–Cu–PGE sulfide deposits. Most of the mineralization occurs in Sublayer, Footwall Breccia and inclusion-bearing quartz diorite (IQD), all of which contain significant (Sublayer and IQD) to minor (Footwall Breccia) amounts of olivine-bearing mafic–ultramafic inclusions. These inclusions have only rare equivalents in the country rocks and are closely associated with the Ni–Cu–PGE sulfide mineralization. They can be divided into three groups on the basis of petrography and geochemical characteristics. Group I (n = 47) includes igneous-textured olivine melanorite and olivine melagabbronorite inclusions in the Whistle and Levack embayments on the North Range with Zr/Y, Zr/Nb, Nb/U and Zr/Hf similar to igneous-textured Sublayer matrix. Group I inclusions are interpreted to be anteliths that crystallized from a mixture of SIC impact melt and a more mafic melt, probably derived by melting of ultramafic footwall rocks. Group II includes Group IIA (n = 17) shock metamorphosed wehrlite and olivine clinopyroxenite inclusions in the Levack embayment and Group IIB (n = 2) shock metamorphosed olivine melanorite inclusions in the Foy Offset on the North Range. Group II inclusions have similar trace element patterns [e.g. negative Th–U, Nb–Ta–(Ti), Sr and Zr–Hf anomalies] and overlapping Nb/U to a layered mafic–ultramafic intrusion in the footwall of the Levack and Fraser deposits, which together with their limited distribution suggests that Group II inclusions are locally-derived xenoliths. Group III (n = 21) includes phlogopite lherzolite and feldspar lherzolite inclusions with igneous, recrystallized and shock-metamorphic textures in the Trill, Levack and Bowell embayments, and the Foy Offset dike on the North Range. They have no equivalents in the exposed country rocks. The calculated parental magma is similar to continental arc basalt formed by approximately 5% partial melting of garnet peridotite. Ol–Cpx–Pl thermobarometry of several Group III inclusions indicate equilibration at 900–1120 ºC and 210 ± 166 MPa to 300 ± 178 MPa, suggesting crystallization in the upper-middle crust (7·7 ± 6·6 to 10·9 ± 6·5 km), prior to being incorporated into the lower parts of the impact melt sheet during impact excavation. The exotic xenoliths provide information about the depth of impact and composition of upper-middle crust in the Sudbury region at 1850 Ma, the local xenoliths provide information about the thermomechanical erosion process that followed generation of the impact melt, the anteliths provide information about the early crystallization history of the SIC, and all of the inclusions provide constraints on the genesis of Sublayer, IQD, footwall breccia, and associated Ni–Cu–PGE mineralization.