Karoo Large Igneous Province Research Articles

The timing and duration of large-scale carbon release in the Early Jurassic

The Toarcian oceanic anoxic event (T-OAE, ca. 183 Ma) in the Early Jurassic was one of the most significant warming events of the Phanerozoic, associated with large-scale carbon emissions, mass extinction, and perturbations to hydrology and ocean chemistry. However, the age and duration of this hyperthermal have long been uncertain, hindering our understanding of the timing and pace of carbon release and the associated environmental and biotic changes. We present high-precision radioisotopic dates bracketing a biostratigraphically constrained record of the T-OAE in Japan. Our geochronology reveals an unexpectedly short T-OAE duration of ∼300 k.y. and a temporal coincidence with Ferrar large igneous province (LIP) magmatism. Our dates refute previous work linking the T-OAE to the earlier Karoo LIP, which was coincident with an earlier carbon cycle perturbation at the Pliensbachian−Toarcian boundary. Our results suggest both events were driven by extremely rapid (perhaps sub-millennial scale) thermogenic carbon degassing during LIP sill intrusion.

Plume‐Driven Subduction Termination in 3‐D Mantle Convection Models

AbstractThe effect of mantle plumes is secondary to that of subducting slabs for modern plate tectonics when considering plate driving forces. However, the impact of plumes on tectonics and planetary surface evolution may nonetheless have been significant. We use numerical mantle convection models in a 3‐D spherical chunk geometry with damage rheology to study some of the dynamics of plume‐slab interactions. Substantiating our earlier 2‐D results, we observe a range of interaction scenarios, and that the plume‐driven subduction terminations we had identified earlier persist in more realistic convective flow. We analyze the dynamics of plume affected subduction, including in terms of their geometry, frequency, and the overall effect of plumes on surface dynamics as a function of the fraction of internal to bottom heating. Some versions of such plume‐slab interplay may be relevant for geologic events, for example, for the inferred ∼183 Ma Karoo large igneous province formation and associated slab disruption. More recent examples may include the impingement of the Afar plume underneath Africa leading to disruption of the Hellenic slab, and the current complex structure imaged for the subduction of the Nazca plate under South America. Our results imply that plumes may play a significant role not just in kick‐starting plate tectonics, but also in major modifications of slab‐driven plate motions, including for the present‐day mantle.

Release of mercury during contact metamorphism of shale: Implications for understanding the impacts of large igneous province volcanism

Elevated mercury (Hg) in sedimentary strata are a widely used tracer for assessing the relationship between large igneous province (LIP) activity and global environmental change. A key unknown in applying this proxy is the extent to which Hg was sourced from contact metamorphism of sedimentary rocks during sill intrusions versus gaseous emissions of the magmas themselves. Here, we investigate Hg behaviour during contact metamorphism of shales. We show loss of 80–99% of the sedimentary Hg in contact aureoles in four case studies covering the interactions around dykes, sills and plutons associated the High Arctic LIP (Sverdrup Basin, Canada), the Karoo LIP (South Africa) and the Skagerrak-centred LIP (Oslo Rift, Norway). A combination of geochemical data and thermal modelling around a dyke from the High Arctic LIP shows 33% Hg volatilization in the aureole at 265–300°C. The other cases show similar behaviours with significant lowering of organic-bound Hg, more significantly in the innermost 60% of the contact aureoles. We hypothesize that gaseous Hg is transported out of aureoles during metamorphism, together with CH4 and CO2. Furthermore, we estimate the thermogenic Hg mobilization from Karoo LIP aureoles as 72–192 t per km3 of aureole, which is between 1–3 times the estimated volumetric Hg release from Karoo magmas. When scaling our results to the size of the shale portions of the Karoo Basin affected by the LIP and a timescale of 100 kyr of sill emplacement, the average Hg flux is calculated to have been 78–207 t/y with maximum values up to ∼300 t/y. The pulsed nature of intrusive volcanism suggests that this thermogenic Hg flux could have dominated LIP Hg emissions during periods of their life span. Our results demonstrate that the global Hg cycle can be significantly perturbed following LIP-scale sill emplacement into organic-rich sedimentary rocks and our quantification of the emissions based on source-rock analysis provides important information for independent interpretation of the sedimentary Hg record.

The nature of the Karoo mantle source region from the perspective of olivine in the Luenha picrites: an oxygen isotopic and chemical study

The Karoo large igneous province has been divided into rift zone and basin-related groups, with picrites from the Luenha river, Mozambique, representing an end-member of the latter. New O isotope, major and trace element data for olivine have been combined with MELTS crystallisation modelling to deconvolve compositional diversity associated with magma differentiation from source-derived heterogeneity. Three olivines populations have been discerned as follows: (1) the ‘main trend’, which records crystallisation from a variety of magma compositions; (2) the ‘low Cu trend’, which is inferred to represent xenocrysts or antecrysts; and (3) the ‘high CaO’ olivines, which record polybaric crystallisation of a primitive, little fractionated magma. The trace element variability in olivine phenocrysts relates partially to sampling of different parts of the same overall magma transport and storage systems, and partly to heterogeneity of parental magmas and their mantle sources. When the measured δ18Oolivine values have been converted into δ18Omelt values, the mean δ18Omelt values for the ‘main trend’ and ‘low Cu’ groups are indistinguishable from each other (5.7 ± 0.1‰, 2σ); however, the mean δ18Omelt value of 6.1 ± 0.1‰ for the ‘high CaO’ group is distinctly enriched. These data record source heterogeneity and suggest contributions from two mantle sources, one with elevated δ18O, and another with more ‘typical’ mantle δ18O. Combining these data with previously reported trace element and Nd and Sr isotope data support derivation from a mantle source similar to non-chondritic bulk silicate earth, but with minor contributions (1–2% for the enriched magmas) from a recycled sedimentary component. This points to the importance of a primitive mantle source for the basin-related successions in Karoo.

Geochemistry of dolerite intrusions of the southeastern Main Karoo Basin, South Africa: Magma evolution, evidence for thermogenic gas sequestration, and potential implications for the early Toarcian Oceanic Anoxic Event

The igneous intrusions of the Karoo Large Igneous Province (LIP) provide critical information on the magma plumbing system, which served as source rock maturation agent for organic-rich sedimentary rocks of the Main Karoo Basin as well as pathway and barrier to subsurface fluid flow. Here, we present results from the research borehole KWV-01 drilled in the Eastern Cape Province of South Africa. The intersected succession, spanning the Dwyka to the lower Beaufort Group, provides representative samples of the Karoo LIP magma plumbing system.Using whole-rock and Sr-Nd isotope analyses, our data suggest that primary melts were derived from a sub-lithospheric mantle source and acquired a subduction signature through interaction with a metasomatised lithospheric mantle. Three magma types identified within this study show a geochemical overlap with previously reported data from the Golden Valley Sill Complex, ca. 220 km to the NW of the analysed borehole. The results of MELTS and EC-RAFC modelling on the three magma types reveal that their compositional differences are related to variable degrees of fractional crystallisation and contamination by lower crustal rocks and/or the sedimentary host rocks of the Main Karoo Basin.Our results indicate that at least one of the sills acted as fractured reservoir, thereby providing evidence that igneous intrusions are important for CO2 sequestration. We quantified the implications of this process for the Karoo LIP-related Toarcian Oceanic Anoxic Event using the thermogenic gas budget produced by state-of-the-art computer simulations and tested various end-member scenarios for partial thermogenic gas sequestration. In agreement with recent studies, we show that all the flat parts of the Karoo sills acted as sequestration reservoirs. Our results have general implications for LIP volcanism in organic-rich basins, suggesting that a good first-order approximation of thermogenic gas emission is half the amount mobilised in the contact aureoles of LIP plumbing systems.

Heavy Rare Earth Elements and the Sources of Continental Flood Basalts

AbstractHeavy rare earth elements (HREEs) in mafic and ultramafic volcanic rocks are useful recorders of mantle source processes because their ratios are not easily modified by differentiation. Here we utilize REEBOX PRO, a simulator of adiabatic decompression melting of the mantle, to study the behavior of HREEs in the formation of continental flood basalt (CFB) parental magmas in the mantle. We simulate partial melting of depleted peridotite, pyrolitic peridotite, pyroxenite, and peridotite-pyroxenite mixtures at mantle potential temperatures of 1350–1650°C and lithospheric thicknesses of 50–150 km, and compare the results to natural data. Many large igneous provinces are typified by low-Ti and high-Ti CFBs with contrasting HREE patterns. Our results show that low-Ti CFBs originate mainly from peridotitic sources. Flat mid-ocean ridge basalt-like HREE patterns typical of low-Ti CFBs can be generated beneath thick lithosphere (~100 km), given that mantle potential temperatures are high (>1500°C) and garnet is completely consumed from the source. We thus challenge the common interpretation that flat HREE patterns always indicate shallow sources for CFB parental magmas. High-Ti CFBs require pyroxenite-bearing sources (≥10%). Contrary to a common view, their steep oceanic island basalt-like HREE patterns can be generated beneath quite a thin lithosphere (~50 km), which is due to increased garnet stability in pyroxenite sources. When applied to CFBs of the Karoo large igneous province, the results are compatible with a model where a mantle plume penetrates a progressively thinning Gondwana lithosphere.

40Ar/39Ar geochronology of the Drakensberg continental flood basalts: Understanding large argon isotopic variations in mafic groundmass and plagioclase size fractions

The extrusive component of the Karoo Large Igneous Province, the Drakensberg Continental Flood basalts of South Africa and Lesotho, is mainly composed of altered tholeiitic basalts without uraniferous mineral phases that are suitable for high-precision dating (e.g. baddeleyite or zircon), and thus dating is restricted to 40Ar/39Ar groundmass and plagioclase. Previous 40Ar/39Ar dates of plagioclase span from 178.9 ± 4.2 Ma to 189.2 ± 3.3 Ma, while most dates range between 180 and 183 Ma (Jourdan et al., 2007b; Moulin et al., 2011; Moulin et al., 2017). Dating mafic rocks can be hindered by the presence of plagioclase xenocrysts (inherited 40Ar), and the ubiquitous presence of alteration phases. Separation protocols and acid leaching are used to minimize the effect of these issues, although previous 40Ar/39Ar studies of the Drakensburg lavas did not interpret their dates within a petrochronological context, and thus they are open to inaccuracy. We have used cathodoluminescence imaging, QEMSCAN analysis and optical petrology of thin sections to show that two distinct generations of plagioclase occur within the lavas, and at least two different phases of K-rich (adularia and K-rich clays) and Ca- alteration have occurred. 40Ar/39Ar dates have been obtained from some acid-leached and unleached aliquots of groundmass and both plagioclase generations. Our results show that i) large, cluster-forming plagioclase yield older 40Ar/39Ar dates and are interpreted to be xenocrystic, ii) apparently unaltered plagioclase in thin sections hosts secondary Ar reservoirs related to alteration, and iii) alteration and leaching can enhance the effects of 37,39Ar recoil in plagioclase. We conclude that the most accurate date for the Drakensburg lavas that applies to the entire thickness of the lava stacks is 184.4 ± 1.2 Ma, which was obtained from the matrix-forming plagioclase that contains the highest proportion of non-inherited radiogenic 40Ar, is least affected by secondary alteration and did not apparently suffer the effects of Ar isotope recoil.

Paleomagnetism and Rock Magnetism of Permian sedimentary rocks and Early Jurassic Karoo (Large Igneous Province) sills intersected in the KWV-1 borecore, Eastern Cape Province, South Africa

The nearly 100 Ma record of almost continuous shallow marine to non-marine sedimentation from the Late Carboniferous to the Early Jurassic recorded in the Karoo Basin, South Africa, includes coal deposits and potentially significant shale gas occurrences. As such, numerous studies have addressed the thermal history of the Karoo Basin, using a range of approaches. We use paleomagnetic and rock magnetic data to evaluate the paleotemperature history of the host Permian strata (Ecca Group and lowermost Beaufort Group) intersected in the ∼2320 m deep KARIN KWV -1 drilling project core, obtained in 2015. Where cored, the Karoo sedimentary sequence has been intruded by several Early Jurassic mafic sills, belonging to the Karoo Large Igneous Province, of thickness ranging from less than 10 cm to ∼150 m. The continuous core provides a unique natural laboratory for studying the effects of Karoo magmatism on the basin's sedimentary sequence. The KWV-1 drilling project was spudded ∼10 km east of Willowvale town, Eastern Cape Province.The spatial variation in thermal effects in host sedimentary strata, relative to proximity to Karoo sill contacts, is well-defined on the basis of anisotropy of magnetic susceptibility (AMS) data, paleomagnetic demagnetization response, and thermal alteration index (A40) values based on temperature vs. magnetic susceptibility tests. A zone of pervasive remagnetization of host sedimentary strata characterizes the contact aureoles, the spatial extents of which are to a first order dictated by the thickness of each Karoo sill. A40 data are interpreted to show that ambient paleotemperature estimates increase from ∼150 °C in lower Beaufort/upper Ecca strata to ∼285° in Lower Ecca strata that are sufficiently removed from Karoo sills. However, lower Ecca Group rocks in close proximity to Karoo sills were likely heated to ∼500 °C. Natural remanent magnetization (NRM) intensities of rocks in proximity to sills range from ∼17 to ∼ 438 mA/m, and susceptibilities range from 1.3 to 3.2 × 10−3 SI volume. In contrast, for strata sampled at distances greater than the width of the nearest sill, NRM intensities range from 0.056 to ∼4 mA/m, and susceptibility values range between 193 and 488 X 10−6 SI volume. For host sedimentary rocks, remanence is completely unblocked at temperatures ranging from ∼350 °C to ∼630 °C and acquisition and backfield DC demagnetization of isothermal remanent magnetization (IRM) indicate that specimens are fully saturated in fields of ∼200–600 mT, with coercivity of remanence of ∼30–100 mT. All of the eleven Karoo sills examined in this study yield remanence inclinations that are moderate to steep negative, and are inferred to have recorded a normal polarity field, possibly during a single chron. Our results are interpreted to suggest that substantial thermal effects of Karoo sills on the host sedimentary strata are largely confined to the contact aureoles, with a width that is dependent on sill thickness. We hypothesize that parts of the Ecca sedimentary sequence have remained at ambient, basinal burial temperatures below about 150 °C, and thus that an early-acquired Permian remanence may possibly be preserved in these sedimentary rocks at sufficient distance from thermal contact aureoles. Alternatively, much of the sedimentary section may have been subject to a combination of thermoviscous/chemical remagnetization.

Lithostratigraphy of the Bumbeni Complex and its associated subdivisions, South Africa

Abstract The Bumbeni Complex represents an Early Cretaceous volcanic complex that developed at the southern termination of the Lebombo mountain range in northern KwaZulu-Natal, South Africa. The volcano-sedimentary complex intrudes and disconformably overlies the Lebombo Group of the Karoo Large Igneous Province. The volcano-stratigraphy is defined by an assemblage of volcanic and volcaniclastic rocks as well as late stage quartz-syenite intrusions. Based on current geological and geophysical studies the complex is considered to represent the remnants of a collapsed cauldron structure. The stratigraphic succession is intricate and laterally discontinuous, comprising mafic and felsic lava flows, welded and non-welded tuffs, ignimbrites, lahar deposits, rhyolite domes and dykes and coarse-grained syenitic rocks. Rhyolite dome formation and collapse, and associated plinian eruptions and column collapse, were largely responsible for the variety of volcaniclastic rocks, including classic air-fall tuff deposits. The complex forms part of a northeast trending volcanic basement lineament known as the Bumbeni Ridge defined by aeromagnetic data.

Periodicity of Karoo rift zone magmatism inferred from zircon ages of silicic rocks: Implications for the origin and environmental impact of the large igneous province

New U-Pb ages for zircons constrain the duration of silicic magmatism and timing of coeval mafic magmatism across the main rift zone of the Karoo large igneous province in Mozambique. Our 190 ± 2 Ma, 188.4 ± 0.9 Ma, 181.7 ± 1.0 Ma, 180 ± 3 Ma, 178 ± 2 Ma and 172 ± 2 Ma ages support periodicity of Karoo magmatism previously inferred from 40Ar/39Ar age data. The ∼ 190–188 Ma ages confirm early onset of magmatism and the ∼ 182–178 Ma ages correlate the bimodal volcanic successions of the Lower Zambezi and the Movene Formation with widespread silicic magmatism across the rift zone. The ∼ 172 Ma age corresponds to waning magmatic activity. The age range and Hf isotopic compositions of zircons indicate up to ∼9 Ma lifespan for the Jurassic silicic magma chambers and suggest that the ∼2700–400 Ma xenocrysts represent crustal sources of the host rocks. The available chronological data indicate that the ∼183 Ma main phase magmatism was largely confined within the main Karoo and Kalahari basins and that the preceding and subsequent phases were mainly associated with the Karoo rift zone. Judging from geochemical literature, different kinds of magmas were erupted during the successive magmatic phases. We calculate from published geochemical data that the mafic main phase magmas were relatively poor in CO2 and SO2 and the lava stacking patterns point to low eruption rates, which suggests that degassing of sedimentary wall-rocks of intrusions probably triggered the coeval Pliensbachian-Toarcian extinction. In contrast, the mafic late phase magmas were rich in CO2 and SO2 and at least some of the lavas indicate high eruption rates. We propose that efficient degassing from widespread mafic magmatism and explosive eruption of over 30,000 km3 of silicic magmas in the Karoo rift zone linked the ∼182–178 Ma late phase magmatism with contemporaneous global biosphere crises.

Constraining the timescales of mafic magmatism of the Central Karoo Large Igneous Province using high precision U-Pb zircon geochronology

AbstractRecent U-Pb high-precision geochronological studies have shown rapid emplacement of the intrusive doleritic component of the Karoo Large Igneous Province (KLIP) in Southern Africa. However, these studies focused on a relatively small geographic and altitudinal region of the KLIP. Additionally, the timing of initiation of extrusive volcanism, preserved in the Drakensberg-Lesotho highlands and its relationship to the intrusive suite, has only been imprecisely constrained by Ar-Ar dates. Here, we present new high-resolution U-Pb zircon ages on dolerite sills and dykes from across the central eastern Karoo Basin (South Africa) at elevations between mean sea level and 1 560 m, as well as U-Pb detrital zircon data that can be used to estimate the maximum age of volcaniclastic deposition near the base of the extrusive component of the KLIP. Dolerite samples were taken across two areas: (1) thick dykes exposed along the coast of the Indian Ocean to ~1 600 m flanking the Drakensberg Escarpment in the Eastern Cape; and (2) sills between 20 and 220 m below surface, in a borehole core within the interior of the Karoo Basin, 400 km hinterland from the coastline. Our estimated dolerite emplacement ages span a range of ca. 80 thousand years (Kyr), between 183.122 ± 0.029/-0.061 and 183.042 ± 0.042/-0.072 million years ago (Ma), and fall within the 331 +60/-54 Kyr age range previously established for magmatism related to the KLIP, despite the marked increase in sampling coverage in terms of area and altitude in this study. Therefore, KLIP geochronology is consistent with other LIPS such as the Siberian and Deccan Traps that supports the hypothesis of rapid emplacement timescales (<1 Myr). Additionally, these data are consistent with, but better delineate that the KLIP in southern Africa appears to be ca. 500 Kyr older than the main phase of magmatism in the Ferrar LIP of Antarctica. Detrital zircons from the basal volcanic sequence of the Drakensberg Group exhibit age peaks at ca. 1 and 0.5 Ga, typical of the surrounding Namaqua-Natal and Pan-African basement rocks, as well as younger peaks at ca. 260 and 200 Ma that likely relate to source provenances from south-western Gondwana and reworking of the Karoo Supergroup sedimentary rocks. High-precision U-Pb dates of the youngest zircon grains result in a maximum depositional age for the basal pyroclastics of 185.25 ± 0.25 Ma, allowing for a ca. 2 Myr offset with the intrusive Karoo dolerite suite.

Local melt contamination and global climate impact: Dating the emplacement of Karoo LIP sills into organic-rich shale

The emplacement of Large Igneous Provinces (LIPs) is commonly correlated with global climate change, and environmental and biological crises. To establish this complex causative link, chemical proxies from marine sedimentary sections must be temporally tied to LIP activity through high-precision geochronology. The temporal link is often ambiguous, especially as periods of environmental change are commonly shorter than the full duration of magmatic and volcanic activity in a LIP. Most importantly, temporal relations need to be established through precise and accurate U-Pb geochronology; however, mafic rocks are often undersaturated in zircon or baddeleyite. New U-Pb age determinations of pegmatitic pods from Karoo LIP dolerite sills into the Ecca Group of the Karoo Basin yield single crystal dates with uncertainties below 50 ka, which allow for a new, detailed assessment of this intrusive event. These samples yield complicated age spectra and Hf geochemistry, which indicate that localized assimilation of sedimentary rocks enabled zircon saturation within these mafic melts. Zircon crystallization age interpretations from these data indicate that sill emplacement within the Ecca Group was coeval across the basin within uncertainty of U-Pb dating, and therefore represents a period of particularly high magma flux. Since the Ecca Group is locally exceptionally rich in organic matter, the widespread and rapid sill-related contact metamorphism generated a significant volume of thermogenic volatiles. Given the temporal overlap of the emplacement of sills between 183.147 ± 0.059 Ma and 183.187 ± 0.133 Ma and the timing of the global Toarcian Oceanic Anoxic Event (TOAE), our data support a causal link between this discrete period of high flux emplacement of Karoo sills into the Ecca Group and the global climate change that took place during the Early Jurassic.

Machine learning-based prediction of trace element concentrations using data from the Karoo large igneous province and its application in prospectivity mapping

In this study, we present a machine learning-based method to predict trace element concentrations from major and minor element concentration data using a legacy lithogeochemical database of magmatic rocks from the Karoo large igneous province (Gondwana Supercontinent). Wedemonstrate that a variety of trace elements, including most of the lanthanides, chalcophile, lithophile, and siderophile elements, can be predicted with excellent accuracy. This finding reveals that there are reliable, high-dimensional elemental associations that can be used to predict trace elements in a range of plutonic and volcanic rocks. Since the major and minor elements are used as predictors, prediction performance can be used as a direct proxy for geochemical anomalies. As such, our proposed method is suitable for prospective exploration by identifying anomalous trace element concentrations. Compared to multivariate compositional data analysis methods, the new method does not rely on assumptions of stoichiometric combinations of elements in the data to discover geochemical anomalies. Because we do not use multivariate compositional data analysis techniques (e.g. principal component analysis and combined use of major, minor and trace elements data), we also show that log-ratio transforms do not increase the performance of the proposed approach and are unnecessary for algorithms that are not spatially aware in the feature space. Therefore, we demonstrate that high-dimensional elemental associations can be modelled in an automated manner through a data-driven approach and without assumptions of stoichiometry within the data. The approach proposed in this study can be used as a replacement method to the multivariate compositional data analysis technique that is used for prospectivity mapping, or be used as a pre-processor to reduce the detection of false geochemical anomalies, particularly where the data is of variable quality. • Machine learning can predict trace element concentrations from major and minor element geochemical data. • Various machine learning algorithms can accurately predict trace element concentrations. • Log-ratio transformation of geochemical data is not always necessary when using machine learning. • Prediction residuals (difference between predicted and actual quantities) can be used to map geochemical anomalies.

Variable magma flow in sills: Can a magma source be constrained?

It has been increasingly recognized that emplacement of mafic sills plays an important role in magma transport throughout the crust. However, further understanding of sill propagation and merging of flow within sill complexes in large igneous provinces is needed. This study focuses on outcrop along the KwaZulu-Natal North Coast (South Africa), where dolerite sills of the Karoo Large Igneous Province intruded shales, siltstones, and sandstones of the Permian-aged Vryheid Formation. These sills intruded at high crustal levels and are associated with significant vesiculation, host rock fluidization and brecciation. An integrated approach is used to examine sill morphology, intrusive relationships, and inherent structures in the sills in relation to the magma flow processes. Mesoscale magma flow indicators preserved in the sills include bridge structures, intrusive steps, deformed vesicles, and magma lobes. The predominantly prolate magnetic fabrics determined by anisotropy of magnetic susceptibility techniques indicate magma flow in the sills was oriented NW, NE–SW and SSE, and largely locally coaxial with field evidence. The mesoscale magma flow indicators developed from linkage between numerous sill segments and variation in the magma propagation direction within the plane of intrusion and show directional variation within (and between) sills. The local flow dynamics within each sill can be revealed in this way; however, the possible variation in flow in the intrusion must be considered when interpreting the gross magma flow direction of the entire sill network. The sills occur in close vertical proximity and are inferred to have intruded through lateral and upward magma flow within merging sill segments creating an interconnected sill network, fed from an eastward local magma source. Interpretation of magma flow in sills should be done with care, to avoid oversimplification regarding regional flow, and highly variable flow in sills is best reconciled at the field scale.

Assessing the importance of thermogenic degassing from the Karoo Large Igneous Province (LIP) in driving Toarcian carbon cycle perturbations

The emplacement of the Karoo Large Igneous Province (LIP) occurred synchronously with the Toarcian crisis (ca. 183 Ma), which is characterized by major carbon cycle perturbations. A marked increase in the atmospheric concentration of CO2 (pCO2) attests to significant input of carbon, while negative carbon isotope excursions (CIEs) in marine and terrestrial records suggest the involvement of a 12C-enriched source. Here we explore the effects of pulsed carbon release from the Karoo LIP on atmospheric pCO2 and δ13C of marine sediments, using the GEOCLIM carbon cycle model. We show that a total of 20,500 Gt C replicates the Toarcian pCO2 and δ13C proxy data, and that thermogenic carbon (δ13C of −36 ‰) represents a plausible source for the observed negative CIEs. Importantly, an extremely isotopically depleted carbon source, such as methane clathrates, is not required in order to replicate the negative CIEs. Although exact values of individual degassing pulses represent estimates, we consider our emission scenario realistic as it incorporates the available geological knowledge of the Karoo LIP and a representative framework for Earth system processes during the Toarcian.

West Antarctic mantle deduced from mafic magmatism

Abstract Distinct mantle compositions recorded in primitive West Antarctic magmatic rocks vary by tectonic setting and time. Deep asthenospheric mantle-plume sources or shallow metasomatized mantle sources may operate either coincidently or independently to supply melts for magmatism. For example, contemporaneous subduction and plume dynamics produced the Ferrar–Karoo Large Igneous Province; subduction-related melting followed by slab-rollback or melting of slab-hosted pyroxenite explains Antarctic Peninsula volcanism through time; Marie Byrd Land magmatism results from plume materials variably mixed with subduction-modified mantle; while magmatism in Victoria Land and western Ross Sea is best explained by plate dynamics and melting of asthenospheric and metasomatized lithospheric sources, and not by an upwelling plume. Element and isotopic ratios show a fundamental change between Marie Byrd Land and Victoria Land mantle domains. Specifically, Pb isotopes indicate that Victoria Land magmatism sources have a stronger focal zone (FOZO) mantle component, while Marie Byrd Land magmatism possesses more of the high μ = high 238 U/ 204 Pb (HIMU) mantle component that leads to high 206 Pb/ 204 Pb over time. The chemical and isotopic heterogeneity of relatively unfractionated igneous rocks in West Antarctica reflects fundamental differences in mantle domains and melting conditions. This mantle variability coincides with changes in crustal structure and composition, and has a geophysical signature that is manifest in seismic data and tomographic models.

Sr‐Enriched Glassy Picrites From the Karoo Large Igneous Province are Evolved, not Primitive Magmatic Rocks

AbstractMagmatic products of the Karoo Large Igneous Province can be divided into a volumetrically dominant, compositionally uniform low‐Ti tholeiitic suite, and a subordinate, geographically restricted, compositionally diverse, incompatible‐rich high‐Ti suite. High‐Ti picrites contain up to 2,400 ppm Sr, 1,900 ppm Ba, and 550 ppm Zr, which seems unusual for olivine‐enriched rocks. We studied six Karoo picrites to determine the phase(s) in which Sr resides. Samples consist of 10–30% olivine phenocrysts in a groundmass of brown glass, augite, feldspar, ilmenite, and apatite. Glass compositions vary, but are generally evolved, ranging from basaltic trachyandesite to dacite. X‐ray intensity maps demonstrate that most of the Sr resides in the glasses, and to a lesser extent, in feldspars, if present. Some samples contain two texturally and compositionally distinct glasses, best modeled in terms of magma hybridization: variably accumulated olivine phenocrysts surrounded by evolved Sr‐rich (to 9,470 ppm) Type 2 melts formed by extensive olivine fractionation were infiltrated by chemically distinct Type 1 melts. Upon eruption, Type 2 melt quenched to minor glass around olivine phenocrysts, and the dominant Type 1 glass acquired its evolved composition by quench crystallization of groundmass mineral phases. Both glasses are rich in H2O (up to 3.8 wt. %) but are nearly devoid of CO2. Calculated parental melts have much higher K2O and incompatible trace elements (e.g., Sr or Ba >1,200 ppm) relative to low‐Ti tholeiites. The unusual parental melt compositions imply derivation by small degrees of partial melting from SCLM mantle sources enriched in Sr and other incompatibles.

Serial interaction of primitive magmas with felsic and mafic crust recorded by gabbroic dikes from the Antarctic extension of the Karoo large igneous province

Two subvertical gabbroic dikes with widths of ~ 350 m (East-Muren) and ≥ 500 m (West-Muren) crosscut continental flood basalts in the Antarctic extension of the ~ 180 Ma Karoo large igneous province (LIP) in Vestfjella, western Dronning Maud Land. The dikes exhibit unusual geochemical profiles; most significantly, initial (at 180 Ma) εNd values increase from the dike interiors towards the hornfelsed wallrock basalts (from − 15.3 to − 7.8 in East-Muren and more gradually from − 9.0 to − 5.5 in West-Muren). In this study, we utilize models of partial melting and energy-constrained assimilation‒fractional crystallization in deciphering the magmatic evolution of the dikes and their contact aureoles. The modeling indicates that both gabbroic dikes acquired the distinctly negative εNd values recorded by their central parts by varying degrees of assimilation of Archean crust at depth. This first phase of deep contamination was followed by a second event at or close to the emplacement level and is related to the interaction of the magmas with the wallrock basalts. These basalts belong to a distinct Karoo LIP magma type having initial εNd from − 2.1 to + 2.5, which provides a stark contrast to the εNd composition of the dike parental magmas (− 15.3 for East-Muren, − 9.0 for West-Muren) previously contaminated by Archean crust. For East-Muren, the distal hornfelses represent partially melted wallrock basalts and the proximal contact zones represent hybrids of such residues with differentiated melts from the intrusion; the magmas that were contaminated by the partial melts of the wallrock basalts were likely transported away from the currently exposed parts of the conduit before the magma–wallrock contact was sealed and further assimilation prevented. In contrast, for West-Muren, the assimilation of the wallrock basalt partial melts is recorded by the gradually increasing εNd of the presently exposed gabbroic rocks towards the roof contact with the basalts. Our study shows that primitive LIP magmas release enough sensible and latent heat to partially melt and potentially assimilate wallrocks in multiple stages. This type of multi-stage assimilation is difficult to detect in general, especially if the associated wallrocks show broad compositional similarity with the intruding magmas. Notably, trace element and isotopic heterogeneity in LIP magmas can be homogenized by such processes (basaltic cannibalism). If similar processes work at larger scales, they may affect the geochemical evolution of the crust and influence the generation of, for example, massif-type anorthosites and “ghost plagioclase” geochemical signature.

Spatial variations of sills and implications for magma dispersal across the Karoo basin

Abstract Dolerite sill complexes of the Karoo Large Igneous Province (ca. 183 Ma) show systematic variations in emplacement style and size throughout the Karoo basin. These variations are explained in terms of three main, interrelated factors, namely the overburden thickness or emplacement depth, variations in host rock rigidities as a result of sedimentary facies changes in the Karoo basin, and proximity to magma feeders. In the northern parts of the thinner (<500 m) and more coarse-clastic Karoo stratigraphy, sills intrude preferentially below more rigid sandstone horizons that acted as stress barriers causing the arrest of magma ascent and lateral spreading below sandstone beds. The low overburden promotes roof uplift above sills and associated brittle faulting can initiate the formation of inclined sheets that limits the lateral propagation path of inner sills. Roof uplift is further promoted by the proximity to magma feeders in the basement and resulting variations in magma pressure that control the spreading rate and inflation of sills. Localised dyke networks spaced at regular intervals and rooted in underlying sills reflect the stretching of roof rocks above inflating sills. The combination of these effects results in relatively small (<10 km) diameters of sills in the northern parts of the basin. Sills emplaced at intermediate depths (ca. 700 m) in the central Karoo basin are marked by larger diameters (>30 km) and thicknesses of up to 100 m. This reflects the higher overburden pressures and the delay of roof failure and subsequent formation of inclined sheets. Dyke networks in the roof of these sills become more irregular and non-systematic at these greater depths. At even greater depths of up to 2 km in the southern parts of the Karoo basin, mega-sills reach diameters of 50 to 80 km, but thicknesses of only up to 35 m. Thick shale-rich sequences in the southern Karoo basin facilitate sill emplacement through internal host-rock deformation and ductile flow. The thicker overburden and different host rock rigidity delay or suppress roof failure and formation of inclined sheet, thus allowing for the lateral propagation of sills. The deeper-seated sills are typically not associated with local dyke networks.

New high precision U-Pb ages and Hf isotope data from the Karoo large igneous province; implications for pulsed magmatism and early Toarcian environmental perturbations

To better constrain the age and duration of the magmatism associated with the Karoo large igneous province (LIP), we present new U-Pb ID-TIMS dates and εHf values from baddeleyite and zircon grains from Karoo basin mafic sills and from felsic samples from the Lebombo and Mwenezi monoclines, together with an 40Ar/39Ar age database of Karoo rocks that has been filtered for true plateau ages with >70% of 39Ar released and in which all 40Ar/39Ar ages were recalculated using the current best estimates for the decay constants. Zircon and baddeleyite ages from three Karoo basin sills range from 183.36 ± 0.17/0.27 to 183.06 ± 0.07/0.21 Ma, where the two uncertainties reflect the analytical error and the additional error associated with decay constant uncertainty. Zircon from the Mutandawhe pluton are dated to 176.84 ± 0.06/0.20 Ma, which represents the first high-precision U-Pb age of the late stage Karoo-LIP magmatism in the northern Lebombo-Mwenezi region. Initial hafnium isotopes are close to chondritic for the Karoo basin and central Lebombo samples (εHf from −2 to +3), but more negative for zircon grains from the Mutandawhe pluton (−11.3 ± 1.1, 2SD). In combination with previous studies and in agreement with the updated 40Ar/39Ar ages, we show that the sill complex that intruded the Karoo basin was short-lived at ~320 ± 180 ka and that it pre-dated the magmatism of the Ferrar-LIP by around 460 ka, whereas the entire Karoo-LIP was emplaced over a period of ca. 6.5 Ma. Based on high-precision U-Pb geochronology, Karoo-LIP magmatism occurred after 183.36 ± 0.17 Ma and therefore postdated the extinction pulses of the late Pliensbachian and likely the Pliensbachian-Toarcian boundary. However, we support previous conclusions that the start of the Karoo-LIP activity agrees with the onset of the Toarcian oceanic anoxic event and the early Toarcian warming, indicating that these environmental changes were likely a response to the magmatic activity of the Karoo-LIP.