Title: Neoarchean large igneous provinces on the Kaapvaal Craton in southern Africa re-define the formation of the Ventersdorp Supergroup and its temporal equivalents

U-Pb geochronology on baddeleyite is a powerful technique that can be applied effectively to chronostratigraphy. In southern Africa, the Kaapvaal Craton hosts a wellpreserved Mesoarchean to Paleoproterozoic geological record, including the Neoarchean Ventersdorp Supergroup. It overlies the Witwatersrand Supergroup and its world-class gold deposits. The Ventersdorp Supergroup comprises the Klipriviersberg Group, Platberg Group, and Pniel Group. However, the exact timing of formation of the Ventersdorp Supergroup is controversial. Here we present 2789 ± 4 Ma and 2787 ± 2 Ma U-Pb isotope dilution-thermal ionization mass spectrometry (ID-TIMS) baddeleyite ages and geochemistry on mafic sills intruding the Witwatersrand Supergroup, and we interpret these sills as feeders to the overlying Klipriviersberg Group flood basalts. This constrains the age of the Witwatersrand Supergroup and gold mineralization to at least ca. 2.79 Ga. We also report 2729 ± 5 Ma and 2724 ± 7 Ma U-Pb ID-TIMS baddeleyite ages and geochemistry from a mafic sill intruding the Pongola Supergroup and on an east-northeast–trending mafic dike, respectively. These new ages distinguish two of the Ventersdorp Supergroup magmatic events: the Klipriviersberg and Platberg. The Ventersdorp Supergroup can now be shown to initiate and terminate with two large igneous provinces (LIPs), the Klipriviersberg and Allanridge, which are separated by Platberg volcanism and sedimentation. The age of the Klipriviersberg LIP is 2791–2779 Ma, and Platberg volcanism occurred at 2754– 2709 Ma. The Allanridge LIP occurred between 2709–2683 Ma. Klipriviersberg, Platberg, and Allanridge magmatism may be genetically related to mantle plume(s). Higher heat flow and crustal melting resulted as a mantle plume impinged below the Kaapvaal Craton lithosphere, and this was associated with rifting and the formation of LIPs.


INTRODUCTION
Major intraplate magmatic events, many of which generate large igneous provinces (LIPs), are globally becoming better understood through U-Pb dating of their magmatic feeders: mafic dike swarms and sill provinces. Although the definition of LIPs has become increasingly complex, a LIP is typically defined as a mainly mafic magmatic province that is >0.1 Mkm 2 in area, with volumes >0.1 Mkm 3 (Ernst, 2014). Usually, LIPs are emplaced within 1-5 m.y. or consist of multiple short pulses over a maximum of a few tens of millions of years (Ernst and Youbi, 2017) and geochemically have intraplate characteristics. Additionally, many mafic dike swarms and sill provinces connected with LIPs intrude supracrustal successions and can be crucial temporal markers in the stratigraphic record or help to constrain the timing of critical events in the Earth's history (e.g., Ernst and Youbi, 2017;Gumsley et al., 2017).
The Ventersdorp Supergroup (Fig. 1) on the Kaapvaal Craton in southern Africa is the largest and oldest predominantly volcanic supracrustal succession in the world and has been linked to the slightly lower volume Fortescue Group on the Pilbara Craton; both may have been contiguous in the hypothetical Vaalbara Supercraton (de Kock et al., 2009;Evans et al., 2017). The Ventersdorp Supergroup overlies the worldclass gold-and uranium-bearing conglomerate reefs of the Witwatersrand Supergroup. However, the exact timing of the formation of the Ventersdorp Supergroup has been a matter of controversy, with implications for the termination of Witwatersrand Supergroup and gold deposition (Armstrong et al., 1991;Wingate, 1998;de Kock et al., 2012;Cornell et al., 2017). Armstrong et al. (1991) reported 2714 ± 16 Ma and 2709 ± 8 Ma ages (Table 1) for the lower (Klipriviersberg Group) and middle (Makwassie Formation of the Platberg Group) part of the Ventersdorp Supergroup, respectively, making the whole volcanic province a single, but pulsed LIP, ∼10 m.y. long, connected with a mantle plume (e.g., Eriksson et al., 2002). However, several authors have argued for a longer and older period of formation between ca. 2.78 Ga and ca. 2.71 Ga (Wingate, 1998;de Kock et al., 2012;Cornell et al., 2017). For example, Cornell et al. (2017) re-dated the quartz porphyry of the Makwassie Formation at 2720 ± 2 Ma (Table 1), calling into question the age of ca. 2709 Ma obtained by Armstrong et al. (1991) for the same unit, as well as the ca. 2714 Ma age for underlying flood basalts of the Klipriviersberg Group. Similarly contradicting ages (i.e., older than ca. 2714 Ma) were determined from units that are lithostratigraphically similar to the Ventersdorp Supergroup from the northwestern and western regions of the Kaapvaal Craton (Table 1; e.g., Wingate, 1998;de Kock et al., 2012). The exact correlation of these units to the Ventersdorp Supergroup, however, has yet to be confirmed. It has also been argued that Ventersdorp Supergroup magmatism represent three pulses, or potentially three different LIP events (Ernst and Buchan, 2001;Ernst, 2014).
In this study, mafic sills and dikes intruding units older than the Ventersdorp Supergroup (i.e., the Witwatersrand Supergroup, the Pongola Supergroup, and the basement granitoids and granitic gneisses of the Kaapvaal Craton) provide new chronostratigraphic information for the Ventersdorp Supergroup. This information further constrains the timing and formation of the Ventersdorp Supergroup and temporally related units. For the first time, in conjunction with previous studies, we can now define LIPs as initiating, developing, and terminating the formation of the Ventersdorp Supergroup over a period of ∼100 m.y. (e.g., Ernst and Buchan, 2001). The Witwatersrand Supergroup and its world-class gold deposits can also now be shown to have formed by ca. 2.79 Ga.
Downloaded from https://pubs.geoscienceworld.org/gsa/gsabulletin/article-pdf/4913103/b35237.pdf by Uniwersytet Slaski user as the Ventersdorp contact reef or VCR) are sporadically developed (e.g., Pelletier, 1937;Chunnett, 1994). This formation gives way to the ca. 2714 Ma Klipriviersberg Group flood basalts (Winter, 1976;Armstrong et al., 1991). In the North West Province of South Africa, and southwestern Botswana (Fig. 1), a lithostratigraphically similar succession of flood basalts within the Derdepoort Formation was correlated to the Klipriviersberg Group (e.g., Tyler, 1979). The Derdepoort Formation, however, was inferred by Wingate (1998) Grobler and Walraven, 1993;Moore et al., 1993;Walraven et al., 1996), while the Modipe Complex gabbro and Gabarone Complex granite ( Fig. 1) have yielded ages between 2784 Ma and 2779 Ma, respectively (Grobler and Walraven, 1993;Moore et al., 1993;Mapeo et al., 2004;Denyszyn et al., 2013). The only indication that this dated series of coeval magmatic events may have been spatially more significant is an age of 2781 ± 5 Ma for a dacite in a drill core ∼480 km to the southwest, near Kimberley ( Fig. 1; Cornell et al., 2017). It has thus been speculated that the Klipriviersberg Group forms part of this ca. 2.78 Ga magmatic pulse (Wingate, 1998;de Kock et al., 2012), but prior to this study, such ages were not reported from the central region of the Kaapvaal Craton, which is the type area of the Klipriviersberg Group. An unconformity separates the Klipriviersberg Group from the Platberg Group (van der Westhuizen et al., 2006); the Platberg Group formed within numerous grabens and half grabens following the termination of Klipriviersberg Group volcanism (Winter, 1976;Stanistreet and McCarthy, 1991). The Platberg Group is best developed over the western Witwatersrand Supergroup in the central part of the Kaapvaal Craton and farther westward ( Fig. 1; e.g., Tankard et al., 1982). The typical Platberg Group graben-fill succession consists of basal, coarsegrained sedimentary rocks of the Kameeldoorns Formation, followed by predominantly more mafic volcanic rocks of the Goedgenoeg Formation and more felsic volcanic rocks of the Makwassie Formation (Winter, 1976;Armstrong et al., 1991;van der Westhuizen et al., 2006). In some instances, interbedded mafic volcanic and sedimentary rock (i.e., the Rietgat Formation) is developed overlying the typical successions, but its occurrence is limited to reactivated grabens (van der Westhuizen et al., 2006). Graben-fill is characterized by rapid facies changes and wedgeshaped geometries (de Kock et al., 2012), with possible caldera formation accounting for the excessively thick formations (Meintjes and van der Westhuizen, 2018a). The Makwassie Formation has been dated by Armstrong et al. (1991) at 2709 ± 8 Ma and more recently by Cornell et al. (2017) at 2720 ± 2 Ma at several localities (Table 1). Cornell et al. (2017) further presented an age for the underlying Goedgenoeg Formation of 2746 ± 9 Ma (Table 1). These ages by Cornell et al. (2017) for the Platberg Group are in agreement with age constraints from many lithological correlatives of the Ventersdorp Supergroup that have been described under various names due to non-continuous outcrop in isolated grabens ( Fig. 1; van der Westhuizen et al., 2006). These include the 2739 ± 10 Ma Sodium Group (Altermann and Lenhardt, 2012), the 2714 ± 3 Ma Zoetlief Group (Walraven et al., 1991), the 2729 ± 3 Ma Amalia Group (Poujol et al., 2005), and the 2733-2724 Ma Hartswater Group (de Kock et al., 2012) (see Table 1). It has been suggested that the age differences between Platberg Group equivalents is a result of diachronous graben development across the Kaapvaal Craton (de Kock et al., 2012). All these U-Pb ages are in disagreement with the younger ca. 2714 Ma age for the Klipriviersberg Group, which is stratigraphically beneath the Platberg Group.
The so-called Pniel Group overlies the Platberg Group with an unconformity (Winter, 1976) that followed the cessation of graben development (Stanistreet and McCarthy, 1991). The Pniel Group covers an extensive area in the central and western parts of the Kaapvaal Craton (van der Westhuizen et al., 2006). The base of the Pniel Group, the Bothaville Formation, is composed of an upward-fining succession of conglomerates and sandstones (Visser and Grobler, 1985). The basaltic andesites of the Allanridge Formation (Winter, 1976) are conformable with the underlying Bothaville Formation (van der Westhuizen et al., 2006). The Pniel Group has not been dated directly, but a maximum age constraint of 2720 ± 4 Ma is available for the lithostratigraphically equivalent Zeekoebaart Formation of the southwestern margin of the Kaapvaal Craton ( Fig. 1 and Table 1; Cornell et al., 2018). A minimum age constraint comes from a 2664 ± 1 Ma age for volcanic rocks at the base of the overlying Transvaal Supergroup in the Buffelsfontein Group (Table 1; Barton et al., 1995). The Allanridge Formation is ascribed to renewed rifting and thermal subsidence of the Kaapvaal Craton after the cessation of volcanism in the Platberg Group (Burke et al., 1985;Clendenin et al., 1988).
Mafic dike emplacement on the stabilized Kaapvaal Craton is summarized in de Kock et al.
(2019), and the first known dike swarm occurred between ca. 2980 Ma and ca. 2966 Ma, with the southeast-trending Badplaas Dike Swarm on the southeastern region of the Kaapvaal Craton Gumsley et al., 2015). These dikes are temporally and spatially associated with the Nsuze Group of the Pongola Supergroup and the Piet Retief Suite of the Usushwana Complex Gumsley et al., 2015). Further mafic magmatism was manifest later in a mafic sill province hosted within the Pongola Supergroup, which includes the 2866 ± 2 Ma Hlagothi Complex (Gumsley et al., 2013;. The Modipe Complex was then emplaced at 2784 ± 1 Ma ( Fig. 1 and Table 1; Denyszyn et al., 2013), which is temporally associated with the Gabarone Complex granite and volcanism within the nearby Derdepoort and Kanye Formations (Grobler and Walraven, 1993;Moore et al., 1993;Walraven et al., 1996;Wingate, 1998;Mapeo et al., 2004) on the northwestern region of the Kaapvaal Craton. Approximately 80 m.y. later, mafic magmatism was again documented between ca. 2701 Ma and ca. 2654 Ma on the eastern and southeastern regions of the Kaapvaal Craton ( Fig. 1 and Table 1; Olsson et al., 2010Olsson et al., , 2011Gumsley et al., 2016). This includes the 2701 ± 11 Ma to 2659 ± 13 Ma Rykoppies Dike Swarm, which appears to radiate out from beneath the 2056-2055 Ma Bushveld Complex (Olsson et al., 2011;Zeh et al., 2015) and the 2664-2654 Ma northeast-trending White Mfolozi Dike Swarm (Table 1; Gumsley et al., 2016). The Rykoppies and White Mfolozi Dike Swarms potentially mark the termination of Ventersdorp Supergroup-related magmatism at ca. 2.65 Ga, as these dike swarms are also coeval with 2664 ± 1 Ma proto-basinal volcanism to the Transvaal Supergroup in the Buffelsfontein Group (Barton et al., 1995), as well as other proto-basins. Lastly, there are the so-called Kopenang dikes, which were documented by Meier et al. (2009) on the central region of the Kaapvaal Craton within the Witwatersrand Supergroup, although they have random trends. Meier et al. (2009) linked these mafic dikes geochemically to the Klipriviersberg Group. However, their exact relationship remains to be confirmed.

SAMPLING
In this study, three samples of mafic sills and one sample from a mafic dike are examined (Table DR1 1 ). Two samples, named AM1 and VJ1, were collected from drillcore near one of the type areas of the Ventersdorp Supergroup in the central region of the Kaapvaal Craton in the Free State Province of South Africa (Fig. 1). Another two samples, PRGE and ENED08, were collected from outcrops in the southeastern (Mpumalanga Province) and southeasternmost (KwaZulu-Natal Province) regions of the Kaapvaal Craton in South Africa (Fig. 1).
Mafic intrusive rocks are commonly intersected in boreholes through Precambrian strata and especially within the Witwatersrand Supergroup on the central region of the Kaapvaal Craton. Many are sills that are over ∼100 m thick, including the studied sills, which are stratigraphically below the flood basalts of the Klipriviersberg Group. The sills are weakly metamorphosed and deformed along with the Witwatersrand Supergroup strata they intruded. The AM1 and VJ1 mafic sills are located east of Klerksdorp and Welkom but south of the Vredefort Dome. The sills intruded into the Government Subgroup of the upper West Rand Group from the Witwatersrand Supergroup (Fig. 2). Sample VJ1 was collected from drill core ∼10 km northeast of Edenville from a ∼200-m-thick mafic sill at ∼1400 m depth, which intruded into sandstones of the Palmietfontein Formation. Sample AM1 was collected from a drill core ∼10 km north of Steynsrus ( Fig. 2) and comes from a ∼110-mthick mafic sill at ∼1800 m depth that intruded into sandstones of the Afrikaner Formation.
Another mafic sill, PRGE, was sampled near the road between Piet Retief and Amsterdam in southeast Mpumalanga (Fig. 3A). The outcrop is on the border between the Evergreen and Redcliff Farms, with the mafic sill intruding into the Redcliff Formation. The Redcliff Formation, comprising shales and banded iron formation, and has been assigned to the Mozaan Group of the Pongola Supergroup. The mafic sill itself was previously ascribed to the 2990-2978 Ma Piet Retief Suite of the Usushwana Complex (Hammerbeck, 1982;Gumsley et al., 2015).
Sample ENED-08 is from an east-northeasttrending mafic dike that was sampled near the White Mfolozi River between Ulundi and Vryheid (Fig. 3B). The dike intrudes the 3254-3234 Ma granitoid gneiss basement in the White Mfolozi inlier of northern KwaZulu-Natal  and is part of a mafic dike swarm, herein termed the Ulundi Dike Swarm. This dike swarm is in turn cut by the 2664-2654 Ma plagioclase-megacrystic, northeast-trending White Mfolozi Dike Swarm (Gumsley et al., 2016). Further, the sampled dike cuts an undated southeast-trending mafic dike; both are in turn cut by a ca. 2423 Ma mafic sheet . Some geochemistry and paleomagnetism on these mafic units was conducted by Klausen et al. (2010) and Lubnina et al. (2010).

Petrography
Petrographic analysis of thin sections was undertaken at the Institute of Earth Sciences in the University of Silesia in Katowice using an Olympus BX-51 optical microscope. Mineral chemical analysis of the main rock-forming and accessory minerals was carried out at the Inter-Institutional Laboratory of Microanalyses of Minerals and Synthetic Substances, University of Warsaw, using a CAMECA SX-100 electron microprobe. The analytical conditions employed an accelerating voltage of 15 kV, a beam current of 20 nA, counting times of 4 s for the peaks and background, and a beam diameter of 1-5 µm.
Reference materials, analytical lines, and mean detection limits (in wt. %) can be found in Table  DR2 (see footnote 1). Diffracting crystals include LIF, PET, TAP, and LPET.

Whole-Rock Major and Trace Element Geochemistry
For whole-rock geochemistry, samples were selected from the most homogenous parts of the sampled mafic sills and dike after the weathered material was removed. They were then hand crushed and milled in a tungsten-carbide ring mill. The resulting material was then coned and quartered before being dispatched for analysis of major, minor, and trace elements at Bureau Veritas Laboratories in Vancouver, Canada. Analyses were made on glass beads prepared from the powdered samples with a sample-to-flux (lithium tetraborate) ratio of 1:10, and the resulting molten bead was rapidly digested in weak nitric acid solution. Inductively coupled plasma emission spectroscopy (ICPES) was used for major and minor elements, and inductively coupled plasma mass spectrometry (ICPMS) was used for trace elements, including rare-earth elements (REE). The volatile content of each sample was determined by loss on ignition (LOI). The data were plotted in diagrams using the GeoChemical Data toolkit (GCDkit; Janoušek et al., 2006).

U-Pb Baddeleyite Geochronology
The sampled units were processed for baddeleyite U-Pb dating by the isotope dilution-thermal ionization mass spectrometry (ID-TIMS) and laser ablation (LA)-ICPMS analytical methodologies. The LA-ICPMS methodology was deployed to test the reproducibility between the two independent methods at the LA-ICPMS Laboratory at Lund University. Samples were hand crushed and ground to a coarse-grained powder using a chrome-steel ring mill in the Department of Geology, Lund University. The water-based separation process developed by Söderlund and Johansson (2002) was used to extract baddeleyite grains for ID-TIMS; however, only VJ1 yielded suitably sized baddeleyite grains for LA-ICPMS (greater than 60 µm long and 20 µm wide).

ID-TIMS
The extracted baddeleyite grains were split into fractions that were transferred in ethanol to pre-cleaned Teflon capsules. The fractions were washed for several cycles in ultrapure HNO 3 and H 2 O. A small amount of a 236-233 U-205 Pb tracer solution and an ultrapure HF:HNO 3 (10:1) mixture was added to each Teflon capsule. The bad-deleyite fractions were completely dissolved in an oven after three days at ∼190 °C. The capsules were placed on a hotplate at ∼100 °C until the HF solution had evaporated. A mixture of ultrapure 0.25 M H 3 PO 4 and 6.2 M HCl was added to each capsule and dried again on a hotplate. The remaining sample droplet in the Teflon capsules was mixed with 2 µl of a prepared Si gel and placed on an outgassed Re filament. The samples on the outgassed Re filaments were heated until the H 3 PO 4 burnt off. The resultant filaments were then placed in a carousel inside a Finnigan TRITON thermal ionization mass spectrometer at the Department of Geosciences, Swedish Museum of National History in Stockholm. The mass spectrometer is equipped with a secondary electron multiplier and Faraday Cups.
The Re filaments and samples were heated to temperatures between ∼1210 °C and ∼1250 °C, where the intensities of the isotopic masses of 204 Pb, 205 Pb, 206 Pb, 207 Pb, and 208 Pb were measured in cycles of between 20 and 140 at gradually increasing temperatures. The Pb signal intensities were measured in either static mode with Faraday Cups or in dynamic mode with peak-switching using the secondary electron multipliers. The isotopic masses of 233 U, 236 U, and 238 U were measured as oxides at temperatures of ∼1300 °C to ∼1340 °C at increasing temperatures. The 235 U was calculated from 238 U, using 238 U/ 235 U = 137.818, following Hiess et al. (2012). No significant differences due to interferences from 18 O were detected, and therefore corrections were not made. The measurements were made in dynamic mode for between 40 and 80 cycles. Initial data reduction was made using an "in-house" Microsoft Excel program using algorithms from Ludwig (2003). However, the Microsoft Excel macro, Isoplot 4.15, from Ludwig (2012) using U decay constants from Jaffey et al. (1971), was used for final data calculations and concordia diagrams. Initial common Pb corrections were made using the isotopic compositions from the global common Pb evolution model of Stacey and Kramers (1975) at the age of the sample. Error propagations in dates are given at 2σ and do not include decay constant errors.

LA-ICPMS
Baddeleyite grains from sample VJ1 were placed on double-sided adhesive tape and mounted into epoxy resin in a 1-inch-wide Teflon ring and left to harden for three days. The hardened epoxy mount was then separated from the tape and polished for between 20 s and 140 s, with 9 µm, 3 µ, and 1 µm diamond paste using a Struers Rotopol-22 automated polisher. The mount was subsequently cleaned with ethanol. The reference material used was baddeleyite from the Phalaborwa Complex, which was dated to 2059.60 ± 0.35 Ma (Heaman, 2009).
The analyses were carried out at the LA-ICPMS Laboratory in the Department of Geology at Lund University. The LA-ICPMS uses a 193 nm Analyte G2 laser ablation unit with a two volume HelEx sample holder and an Aurora Elite quadrupole inductively coupled plasma mass spectrometer. Ablation was done with a 5 Hz repetition rate, a fluence of ∼4 j/cm 2 , using 20 µm spot diameters. The total acquisition time for each analysis was 55 s, of which the first 20 s were used to determine the Ar/He/N 2 gas blank, followed by 28 s of laser ablation and 3-4 s of wash-out delay. The He carrier gas was mixed downstream with Ar and N 2 before entering a "squid" signal smoothing device made up of several tubes that split the stream and re-joined it before entering the plasma. The mass spectrometer was calibrated using NIST612 glass to give stable 206 Pb, 207 Pb, and 238 U signals and low oxide production rates ( 238 U 16 O/ 238 U below 0.5%) and a 232 Th/ 238 U ratio of ∼1. The 202 Hg, 204 (Pb + Hg), 206 Pb, 207 Pb, 208 Pb, 232 Th, 235 U, and 238 U intensities were determined using secondary electron multipliers. The interference of 204 Hg on 204 Pb was monitored by measuring 202 Hg, assuming a 202 Hg/ 204 Hg ratio of 4.36 (natural abundance). Although 235 U was measured, the 207 Pb/ 235 U was calculated from 238 U, using 238 U/ 235 U = 137.818 following Hiess et al. (2012). The laser-induced elemental fractionation and instrumental mass bias on measured isotopic ratios were corrected through standard-sample bracketing using the Phalaborwa baddeleyite standard (Heaman 2009). A total of 36 analyses were made in the following order: six standards followed by seven unknowns (i.e., VJ1), four standards, six unknowns, four standards, five unknowns, and four standards. Data reduction was done using the Iolite software (Paton et al., 2011) with the U-Pb geochronology data reduction scheme routine of Paton et al. (2010). Correction routines for downhole fractionation and instrumental drift were applied in the Iolite software. Common-Pb correction was done in the Iolite software using VizaulAge by Petrus and Kamber (2012) from the terrestrial Pb-isotope composition of Stacey and Kramers (1975). All age data were calculated using Isoplot 4.15 (Ludwig, 2012), which was also used to plot and evaluate data. Error propagation in dates reported are given at 2σ and do not include decay constant errors.

Petrography
The four samples investigated have a mineral assemblage typical of mafic intrusions (Table  DR1). The primary rock-forming minerals are mostly zoned plagioclase feldspar (An 41 Ab 59 -An 18 Ab 82 ; Fig. 4A Table  DR3). Further evidence of metamorphic alteration is provided by secondary Al-bearing titanite rims on ilmenite and ilmenite-magnetite inter-growths (Fig. 4D). Locally, secondary biotite (Ti = 0.37-0.43 a.p.f.u.; #fm = 0.64-0.67; Table  DR3) was formed at the expense of opaque minerals (Fig. 4B). Late chlorite-epidote veins cut the rocks. The metamorphic alteration is typical of greenschist facies metamorphic conditions. Samples also show secondary, low-temperature alteration manifest by the sericitization of plagioclase feldspar and chloritization of mafic silicates and aluminosilicates. The least altered samples are PRGE and AM1, while the other two samples show pervasive alteration (Table  DR1). Sample PRGE displayed graphic intergrowths of quartz and alkali feldspar (Or 79-56 , Ab 36-17 An 1-5 ), forming a groundmass between the clusters of mafic phases (Figs. 4C and 4E).

Whole-Rock Major and Trace Geochemistry
Whole-rock major and trace element compositions for the studied mafic sills and dike are listed in Table DR4 (see footnote 1). Data from previous studies (Crow and Condie, 1988;Marsh et al., 1992;Nelson et al., 1992;Meier et al., 2009;Klausen et al., 2010;Gumsley et al., 2016;Meintjes and van der Westhuizen, 2018b) are used for comparison and are shown in Table  DR5 (see footnote 1). Data were filtered using an LOI of less than 10% to account for extensive alteration. A thorough review of the geochemistry and petrogenesis of the Ventersdorp Supergroup can be found in Humbert et al. (2019). The mafic sills in this study (AM1 and VJ1) straddle the basalt and the basaltic andesite/andesite fields in the Zr/Ti-Nb/Y classification diagram ( Fig. 5A; Pearce, 1996, after Winchester andFloyd, 1977), and are similar to the Klipriviersberg Group basalts and Kopenang dikes, with the lowest Nb/Y and Zr/Ti. The mafic sill PRGE and mafic dike ENED08 are within the basalt and basaltic andesite/andesite range of the Platberg Group volcanic rocks (and Rykoppies Dike Swarm), which are bimodal in Zr/Ti and show greater Zr/Ti in the dacites/rhyolites of the Makwassie Formation. The Allanridge Formation basaltic andesites have a slightly higher Nb/ Yb but are equivalent to the lower Zr/Ti grouping from the Platberg Group. The White Mfolozi Dike Swarm, however, is geochemically distinct from the other groupings, having lower Zr/Ti and Nb/Yb. In the AFM diagram (Irvine and Baragar, 1971), samples AM1 and VJ1 and the basalts of the Klipriviersberg Group and the Kopenang dikes plot in the field of tholeiitic basalts with slight calc-alkaline affinities (Fig. 5B). Samples PRGE and ENED08 again show more similarities to the Platberg Group mafic volcanic rocks, Rykoppies Dike Swarm, and Allanridge Formation basaltic andesites, which have a more calc-alkaline differentiation trend. Using the Th/Yb-Nb/Yb diagram (Pearce, 2008), AM1 and VJ1 are in close proximity to the Kopenang dikes and the Klipriviersberg Group basalts, whereas samples PRGE and ENED08 plot near the mafic volcanic rocks of the Platberg Group, Rykoppies Dike Swarm, and Allanridge Formation basaltic andesites (Fig. 5C)

U-Pb Baddeleyite Geochronology
The water-based separation process of Söderlund and Johansson (2002) yielded ∼40 grains of baddeleyite each for VJ1 and PRGE, and 30 grains for AM1, whereas only 20 grains of ENED08 were recovered. The best-quality grains analyzed in the VJ1, AM1, and ENED08 samples were light brown and generally clear with some slight to moderate frostiness with a length of between 50 µm and 60 µm for VJ1 and ∼30 µm for AM1. ENED08 baddeleyite grains were between 20 µm and 30 µm long. Sample PRGE had baddeleyite grains that were dark brown and generally clear with some slight frostiness with lengths of ∼40 µm. U and Pb isotopic measurements are listed in Tables DR6 and DR7 (see footnote 1), and data are plotted on Wetherill concordia diagrams in Figures 7 and 8.

ID-TIMS
Four fractions of baddeleyite from VJ1 were analyzed, which were composed of between two and three baddeleyite grains each (Fig. 7A). Regression of these fractions produced a precise upper intercept date of 2787 ± 2 Ma (mean square of weighted deviates [MSWD] = 1.5). Three fractions are between 2% to 1% discordant and a fourth fraction is 6% discordant, likely due to zircon alteration. The most discordant fraction constrains the lower intercept of 295 ± 98 Ma. A weighted mean 207 Pb/ 206 Pb date of the three least discordant fractions is 2785 ± 1 Ma (MSWD = 0.67). The fourth fraction is excluded because it is not within error of the 207 Pb/ 206 Pb dates of the other three fractions. Two fractions of AM1, composed of six baddeleyite grains each, produced a preliminary upper intercept date of 2789 ± 4 Ma (Fig. 7B), which is within error of the date of VJ1. Regression yields a preliminary lower free intercept date of 106 ± 190 Ma. The two fractions were 5% and 2% discordant, with 207 Pb/ 206 Pb dates between ca. 2788 Ma and ca. 2786 Ma. For sample ENED08, three baddeleyite fractions were analyzed, with three grains in each (Fig. 7C). Free regression yields upper and lower intercept dates of 2729 ± 5 Ma and 228 ± 57 Ma, respectively (MSWD <0.1). Two fractions plot between 6% and 5% discordant and have 207 Pb/ 206 Pb dates between ca. 2723 Ma and ca. 2722 Ma, and a third fraction is 16% discordant, likely as a result of zircon alteration, and constrains the lower intercept.
Sample PRGE yielded an upper intercept date of 2727 ± 3 Ma and a lower intercept date of 540 ± 280 Ma (MSWD = 0.71) from analysis of four fractions of baddeleyite (Fig. 7D). The discordance varies between 10% and 0%, with the three least discordant fractions being between 1% and 0% discordant. A weighted mean 207 Pb/ 206 Pb age calculated on these three least discordant dates is 2724 ± 7 Ma (MSWD = 3.3). The fourth fraction is excluded due to not being within error of the 207 Pb/ 206 Pb dates of the other three fractions.

LA-ICPMS
Analytical results from the Phalaborwa Complex reference material (Heaman, 2009) are presented in Table DR7 (see footnote 1). No significant difference was observed between 207 Pb/ 206 Pb ratios of VJ1 measured in the sequence, but the isotopic data show variable discordance (between ∼17% and ∼1%), forming a discordant array in a Wetherill concordia diagram (Fig. 8A). Using all data (Table DR8;  Difficulties in determining the true U-Pb dates of grains analyzed using LA-ICPMS might arise from difficulties with correction of common Pb. In the presented data, the average mass 204 background signal was relatively high, around 1900 counts per second, which increases the risk that small amounts of common Pb cannot be detected. When detecting isotopes close to the detection limit, it is not possible to quantify concentrations; however, by plotting the 206 Pb/ 204 Pb against common Pb corrected and uncorrected 207 Pb/ 206 Pb dates, we observe an increasing age difference for analyses with low 206 Pb/ 204 Pb, which mainly shows up in decreasing common Pb corrected dates. This supports the idea that some grains contain a small amount of common Pb but also that we seemingly overcorrect the data, as observed when 206 Pb/ 204 Pb ratios are below 9000. Using all uncorrected data, we obtain an upper intercept date of 2799 ± 15 Ma (MSWD = 3.4) and a weighted mean 207 Pb/ 206 Pb date of 2796 ± 6 Ma (MSWD = 3.2). Using only uncorrected data from analyses with 206 Pb /204 Pb above 9000, we   Fig. 8B). As observed, there is, within the limits of uncertainty, no significant difference between these age estimates. We, however, prefer to use the data that are arguably least affected by common Pb, and thus our best date estimate of VJ1 is 2798 ± 12 Ma, using the upper intercept, as with the ID-TIMS data.

Crystallization Ages by ID-TIMS and LA-ICPMS on Baddeleyite
Regression of the U-Pb ID-TIMS data for VJ1 yields an upper intercept date of 2787 ± 2 Ma (MSWD = 1.5) and is within analytical uncertainty of the LA-ICPMS result as well as the preliminary AM1 date of 2789 ± 4 Ma, both of which are also obtained by free regression. This is preferred age over a weighted mean 207 Pb/ 206 Pb age on the three least discordant analyses at 2785 ± 1 Ma (MSWD = 0.67), as no data overlap the concordia, making the weighted mean date a minimum age. Although the lower intercept is close to 0 Ma, which indicates recent Pb disturbance, it is still beyond uncertainty at 295 ± 98 Ma, and therefore an isotopic disturbance at an earlier time is possible. This is interpreted as likely belonging to the Karoo magmatic event at 183-179 Ma (Duncan et al., 1997). The relatively precise 2787 ± 2 Ma date is our preferred crystallization age, since ID-TIMS normally allows for robust common-Pb correction. Noteworthy is that the discordance of ID-TIMS fractions is ∼3%, i.e., much less than for the LA-ICPMS analyses. The crystallization age of ENED08 is also interpreted as 2729 ± 5 Ma (MSWD < 0.01) using a free regression, as no analyses overlap with the concordia, with Pb disturbance again likely attributed to Karoo magmatism. The upper intercept date of PRGE is 2727 ± 3 Ma (MSWD = 0.71). However, a weighted mean 207 Pb/ 206 Pb date on the three least discordant fractions, which overlap the concordia, is 2724 ± 7 Ma (MSWD = 3.3). The weighted mean 207 Pb/ 206 Pb date is thus our preferred crystallization age for PRGE.

Source of Magma
The incompatible element patterns of all the magmatic events presented in Figure 6 are different, but with the exception of the White Mfolozi Dike Swarm, they all show "subduction-related" or "arc-like" signatures (e.g., Pearce, 2008;Ernst, 2014;O'Neill and Jenner, 2016;Humbert et al., 2018; that are notably marked by negative Nb and Ta anomalies. These signatures are common in many Archean and Paleoproterozoic mafic rocks of the Kaapvaal Craton (e.g., O'Neill and Jenner 2016;Humbert et al., 2018; and in other intraplate events (including LIPs) in the rest of the world (e.g., Pearce, 2008;Ernst, 2014). These events also plot higher Th/ Yb and Th/Nb values than the MORB-OIB array (Fig. 5C). According to Pearce (2008), this preferential enrichment in Th over Nb may either indicate crustal contamination or assimilation of the magma and/or a lithospheric mantle source that was fluid-metasomatized by past subduction. Substantial crustal contamination actually has been inferred through whole-rock Sm-Nd isotopic data in the case of the Platberg Group (average εNd 2720 Ma = -3.4) and Allanridge Formation (average εNd 2700 Ma = -3.0), with an average εNd 2780 Ma = 0 for the source of the Klipriviersberg magmas probably also possibly contaminated to a lesser extent (Humbert et al., 2019). Note that no isotopic data exist for the other events so far. In addition, the quartz and alkali feldspar groundmass in the PRGE sample can be interpreted as the result of local mixing-mingling phenomena (e.g., Burda et al., 2011) from assimilation of silica-rich material. It can also simply represent the highly fractionated portions of rift-related intracontinental basaltic magma.
Samples from the White Mfolozi Dike Swarm are clearly geochemically distinct from any other unit on the Kaapvaal Craton (Figs. 5 and 6). They notably present lower Th/Yb, Nb/Yb, Th/ Ta, and La/Yb ratios, as well as "flat" and moderately enriched signatures compared to chondrite/primitive mantle in REE/trace-element diagrams. The White Mfolozi Dike Swarm samples plot close to primitive mantle in Figures 5C and   5D and might be from such a source. On the contrary, they also may have been generated from a depleted mantle source (i.e., NMORB-like) that later underwent minor crustal contamination.

The Klipriviersberg Large Igneous Province
The herein dated 2787 ± 2 Ma and 2789 ± 4 Ma mafic sills VJ1 and AM1 are the first documented occurrence of this magmatic event from the central region of the Kaapvaal Craton (Fig. 9). These sills are situated stratigraphically beneath the Klipriviersberg Group flood basalts of the Ventersdorp Supergroup in the Witwatersrand Supergroup and are interpreted as likely related to them. Further, these mafic sills closely resemble the flood basalts of the Klipriviersberg Group geochemically (Figs. 5 and 6; e.g., Winter, 1976;Marsh et al., 1992;Humbert et al., 2019). Meier et al. (2009) has already shown the similar geochemical  Table 1 and Table DR9 (see footnote 1). Data were filtered, and only ages with 2σ less than 20 Ma are plotted. Probability density distribution diagrams are shown for detrital zircon ages from the Postmasburg Group, the Pretoria Group, and Ventersdorp Supergroup. Data sources are given in Table  DR10 (see footnote 1). Only results with discordance less than 5% are plotted. Uncertainties in age are given at 2σ; n refers to the total number of data points analyzed.

A B
Downloaded from https://pubs.geoscienceworld.org/gsa/gsabulletin/article-pdf/4913103/b35237.pdf by Uniwersytet Slaski user characteristics of mafic feeders to the Klipriviersberg Group flood basalts in the form of the Kopenang dikes intruding the Witwatersrand Supergroup in the same area (Fig. 1). Therefore, it is reasonable to interpret that the mafic sills AM1 and VJ1 analyzed in this study belong to a single magmatic system represented by mafic intrusions that fed the overlying flood basalts of the Klipriviersberg Group.
As shown in Figure 9, we now recognize a widespread 2791-2779 Ma magmatic event across the central, western, and northwestern regions of the Kaapvaal Craton. The mafic sills from the central region of the Kaapvaal Craton in this study, AM1 and VJ1, are coeval with the bimodal volcanic Kanye Formation and its plutonic equivalent, the Gabarone Complex granites. The Kanye Formation has been considered a protobasinal or precursor phase to the Ventersdorp Supergroup in the Swartruggens Trough (Tyler, 1979;Ernst, 2014). The nearby Derdepoort Formation and Modipe Complex gabbro was dated to 2782 ± 5 Ma and 2784 ± 1 Ma, respectively (Table 1; Wingate, 1998;Denyszyn et al., 2013). These units were originally correlated with the Klipriviersberg Group flood basalts and other units of the Ventersdorp Supergroup (Tyler, 1979). Similarly aged units are documented in the Amalia-Kraaipan granite-greenstone terrane of the western and northwestern regions of the Kaapvaal Craton (Figs. 9A and 9B). These include the Mosita, Mmathete, and "red" granites, as well as the Turfloop and Rooibokvlei granites in the Makoppa Dome (e.g., Grobler and Walraven, 1993;Moore et al., 1993;Mapeo et al., 2004;de Kock et al., 2012;Denyszyn et al., 2013;Ernst, 2014). This magmatic event is also recorded in detrital zircon data from the Ventersdorp and Transvaal supergroups (Figs. 9A and 9B; Tables DR9 and DR10; see footnote 1).
Furthermore, the ages proposed for the Klipriviersberg Group flood basalts in this study do not agree with the 2714 ± 16 Ma age obtained by Armstrong et al. (1991) from the type locality of the Klipriviersberg Group flood basalts. This is the only volcanic event occurring in the region of the mafic sills dated in this study other than the ca. 2914 Ma Crown Formation lava of the Witwatersrand Supergroup (Armstrong et al., 1991;Fig. 9A). However, the age presented by Armstrong et al. (1991) for the Klipriviersberg Group flood basalts has been questioned by Wingate (1998), de Kock et al. (2012), Cornell et al. (2017), and Humbert et al. (2019. A recalculation of the Armstrong et al. (1991) isotopic data using the least discordant analyses at ca. 2.7 Ga and ca. 1.0 Ga results in an upper intercept date of 2757 ± 57 Ma and a lower intercept date of 1082 ± 140 Ma. This result, as originally shown by Wingate (1998), indicates that the magmatic/ metamorphic event that formed the younger ca. 1.0 Ga generation of zircons may have affected the ca. 2.7 Ga zircons and should be regarded as a minimum age of crystallization for the flood basalts. Recent geochronological studies on the Platberg Group (in the Makwassie Formation) and one of its stratigraphic equivalents (i.e., the Hartswater Group) also add temporal constraints on the age of the Klipriviersberg Group, indicating an older age (e.g., de Kock et al., 2012;Cornell et al., 2017). Furthermore, Barton et al. (1989) analyzed detrital zircons from the Venterspost Formation, which constitutes the base of the Klipriviersberg Group, and obtained a maximum age of 2780 ± 5 Ma for deposition of the formation from the youngest concordant zircon grains.
The magmatic event described above meets the criteria for a LIP, as specified by Ernst (2014). As illustrated in Figure 9, we can now redefine this event to encompass the Klipriviersberg Group flood basalts and coeval mafic sills, the Kanye and Derdepoort bimodal volcanic rocks, and the temporally related Gabarone Complex and Modipe Complex granites and gabbros as well as granites farther afield to the west (Fig. 9B). The LIP had an area of ∼250,000 km 2 , with the thickness of the flood basalt package up to 2 km. Further material from the magmatic province has likely eroded away, as shown by the extensive detrital zircon spectra peak at ca. 2.78 Ga. An age range for all units belonging to this newly defined Klipriviersberg LIP is 2791-2779 Ma, although through more modern methodology, this age range can likely be narrowed further. Additionally, the lack of any intervening sedimentary layers or unconformities in the flood basalts provides further evidence of the short-lived duration of the magmatic event.

The Platberg Volcanic Province
Following a magmatic hiatus from the Klipriviersberg LIP, large-scale rifting, sedimentation, and volcanism dominated the western and central Kaapvaal Craton between 2754 Ma and 2709 Ma (de Kock et al., 2012;Cornell et al., 2017). In the eastern Kaapvaal Craton, numerous granites were emplaced simultaneously (e.g., Hofmann et al., 2015), which Taylor et al. (2010) argued were likely related to the extension and magmatism that deposited the Ventersdorp Supergroup (Figs. 9A and 9B;Tables DR9 and DR10). In this study, the mafic sill PRGE from the southeastern part of the Kaapvaal Craton intruding into the Pongola Supergroup was dated to 2724 ± 7 Ma, and thus it is coeval with this widespread magmatic event. Geochemically, this sill resembles volcanic rocks within the Platberg Group, specifically the Makwassie Formation (Figs. 5 and 6). Near the sill, the Am-sterdam Formation is the only unit with which it may be linked. The Amsterdam Formation unconformably overlies the Pongola Supergroup and is a volcanic succession composed of the Gobosha Member dacites near Amsterdam and the Vaalkop Member rhyolites near Piet Retief (Hammerbeck, 1982). The sill geologically and geochemically resembles the Gobosha dacites of the Amsterdam Formation (Hammerbeck, 1982), and therefore it is hypothesized as a likely part of a feeder system to the volcanic rocks of the Amsterdam Formation. This interpretation implies that the Amsterdam Formation is the easternmost equivalent of the Platberg Group (Fig. 9B). Additionally, farther to the east in the Mkondo Suite, leucosomes dated to 2730 ± 5 Ma have been reported to date granulite facies metamorphism (Taylor et al., 2010).
To the south of the dated sill, our age of 2729 ± 5 Ma for ENED08 on a coeval eastnortheast-trending dike establishes a newly identified dike swarm on the southeasternmost part of the Kaapvaal Craton. This dike swarm, which we name the "Ulundi Dike Swarm" (Figs. 1 and 3), is the first recognized dike swarm associated with the Platberg Group of the Ventersdorp Supergroup and extends this magmatic event further to the southeast. The ENED08 dike is cut by plagioclase megacrystic northeast-trending dikes of the 2664-2654 Ma White Mfolozi Dike Swarm (Gumsley et al., 2016) as well as by a mafic sheet dated at 2420 ± 3 Ma .
Using the new age constraints presented in this study of 2729 ± 5 Ma for ENED08 and 2724 ± 7 Ma for PRGE, and age constraints of 2733-2724 Ma in de Kock et al. (2012) and 2720 ± 2 Ma in Cornell et al. (2017), the Platberg volcanic province (or possibly a LIP pulse of a combined Platberg-Allanridge LIP) is constrained to between ca. 2754 Ma and ca. 2709 Ma. For the Makwassie Formation, the Cornell et al. (2017) age of 2720 ± 2 Ma replaces the original 2709 ± 8 Ma age presented by Armstrong et al. (1991). A detrital zircon-spectra peak again provides evidence of magmatism and erosion at ca. 2.71 Ga.

The Allanridge Large Igneous Province
These new data for samples ENED08 and PRGE help to refine the minimum and maximum age constraints on the Allanridge Formation basaltic andesites, which unconformably overlie the Platberg Group along with the Bothaville Formation (both part of the Pniel Group). The Allanridge Formation basaltic andesites are sufficiently voluminous to be defined as a LIP according to Ernst (2014), being over 200,000 km 2 in area and up to 1 km thick. Age constraints on the Allanridge LIP bracket its deposition between ca. 2709 Ma, which is the age of the underlying Makwassie Formation of the Platberg Group (Cornell et al., 2017), and ca. 2664 Ma, which corresponds to the age of the Buffelsfontein Group in the overlying proto-basins of the Transvaal Supergroup (Barton et al., 1995).
However, the Rykoppies Dike Swarm, emplaced between ca. 2701 Ma and ca. 2664 Ma, may also be linked to the Allanridge LIP and appears geochemically similar (Figs. 5 and 6). This mafic magmatism documents peaks at 2701-2692 Ma, 2686-2683 Ma, and 2664-2659 Ma (Gumsley et al., 2016). Reinterpretation of the complex discordant U-Pb isotopic data and Pb loss histories of these dike swarms appears to indicate two magmatic events at 2701-2683 Ma and 2664-2654 Ma using only the least discordant U-Pb baddeleyite isotopic data from the Rykoppies and White Mfolozi dike swarms. If this is correct, the age of the Allanridge LIP may reflect this older pulse at 2701-2683 Ma in the Rykoppies Dike Swarm, with further magmatism at 2664-2654 Ma represented by the geochemically distinct White Mfolozi Dike Swarm, which is coeval to parts of the Rykoppies Dike Swarm (despite being geochemically different) and volcanism and renewed rifting in proto-basinal fills to the Transvaal Supergroup, such as the bimodal volcanic rocks of the Buffelsfontein Formation at 2664 ± 1 Ma (Barton et al., 1995). The White Mfolozi Dike Swarm may then be regarded as a separate magmatic event, potentially another LIP, especially when including coeval parts of the Rykoppies Dike Swarm. The White Mfolozi Dike Swarm has a unique geochemical signature but has an approximately orthogonal linear trend to the Rykoppies Dike Swarm, which is radiating (Fig. 1). While the White Mfolozi Dike Swarm could represent part of a separate linear or radiating swarm, we speculate that the White Mfolozi Dike Swarm could alternatively be evidence for a giant circumferential swarm Ernst 2018, 2019) that circumscribes the plume center at the focus of the Rykoppies Dike Swarm.

LIPs and the Ventersdorp Supergroup
In summary, Ventersdorp magmatism can be divided into distinct volcanic/magmatic events, some of which are of LIP scale, including the 2791-2779 Ma Klipriviersberg LIP and the 2701-2683 Ma Allanridge LIP, which are separated from the 2664-2654 Ma White Mfolozi/ Rykoppies magmatic province, which may also be a LIP. Despite the unconformity between the Platberg Group and the Pniel Group, the Allanridge LIP was closely preceded by the 2754-2709 Ma Platberg volcanic province, which is composed of a mixed succession of sedimentary and bimodal volcanic rocks. The Platberg Group was, therefore, likely genetically related to a mantle plume responsible for the emplacement of the Allanridge LIP beneath the lithosphere, making the plume distinctly pulsed (with fractionation or arc-related magmatism, with a crustal contribution, as seen from geochemistry; Humbert et al., 2019). The mantle plume center for the Allanridge LIP is well defined by the radiating Rykoppies Dike Swarm. In addition, the associated Ventersdorp rift zone trends approximately from this plume center, suggesting that this rifting was in part coeval with the Allanridge LIP. The location of a plume center for the 2791-2779 Ma Klipriviersberg LIP and 2664-2654 Ma White Mfolozi volcanic province remains unconstrained.

Implications for the Age of Gold Deposition within the Witwatersrand Supergroup
The 2787 ± 2 Ma and 2789 ± 4 Ma mafic sills reported in this study intrude and crosscut the gold-bearing Witwatersrand Supergroup successions, indicating that the placer gold deposits must predate the mafic dikes and sills as suggested by Meier et al. (2009). This makes the minimum depositional age for placer gold and the Witwatersrand Supergroup greater than ca. 2.79 Ga and is in agreement with the ca. 2.78 Ga and ca. 2.76 Ga diagenetic ages reported from xenotime (England et al., 2001;Kositcin et al., 2003) as well as a ca. 2.78 Ga detrital zircon age for the Venterspost Formation of the Klipriviersberg Group, which conformably overlies the Witwatersrand Supergroup (Barton et al., 1989). This further confirms the approximate minimum and maximum age constraints for the cessation of the Witwatersrand Supergroup deposition and the initiation of the Ventersdorp Supergroup. The Venterspost Formation was deposited conformably onto the Witwatersrand Supergroup, with the sedimentary rocks still unconsolidated with the eruption of the Klipriviersberg Group flood basalts (Chunnett, 1994). Diagenesis occurred between ca. 2.78 Ga and ca. 2.76 Ga (England et al., 2001;Kositcin et al., 2003). This has implications for gold metallogeny, indicating that the Witwatersrand Supergroup sedimentary rocks were deposited within a maximum timespan of ∼180 m.y., as opposed to the previously suggested maximum timespan of ∼256 m.y., dating gold metallogeny to ca. 2.78 Ga (Armstrong et al., 1991;Meier et al., 2009).

Tectonic Model
Using the chronostratigraphic constraints above (Figs. 9A and 9B), at 2791-2779 Ma, large amounts of magmatic material were preserved in the central and northwestern regions of the Kaapvaal Craton, such as the Klipriviersberg LIP. This LIP incorporates the Klipriviersberg flood basalts that flowed across the central Kaapvaal Craton and the Kanye and Derdepoort basalts and rhyolites on the northwestern Kaapvaal Craton. This was accompanied by the emplacement of anorogenic granites, which were likely produced from crustal melting associated with an underlying mantle plume responsible for the Klipriviersberg LIP (Grobler and Walraven, 1993). These granites were emplaced along the Colesberg and Thabazimbi-Murchison lineaments, which likely reflected suture zones of thin, weaker crust susceptible to higher heat flows (e.g., Ernst, 2014). A mantle plume is shown in LIP-scale magmatism, incorporating mantle-derived komatiites and high-MgO basalts at the base of the Klipriviersberg Group as well as uplift and erosion of unconsolidated Witwatersrand Supergroup and Venterspost Formation sediments before deposition of the Klipriviersberg Group flood basalts.
As extensive, short-lived magmatism of the Klipriviersberg LIP terminated, the impact of the same or potentially a different mantle plume some ∼30 m.y. later led to extension and graben development associated with volcanism and sedimentation across the central and western Kaapvaal Craton; this is now preserved as the Platberg Group between ca. 2754 and ca. 2709 Ma. In the southeastern Kaapvaal Craton, this magmatism is also now documented in the 2729 ± 5 Ma east-northeast-trending Ulundi Dike Swarm and possibly in the Amsterdam Formation, which accompanied metamorphism, core complex emplacement, and granite formation in Swaziland and adjacent areas of South Africa. Mafic magmatism in this area is also evident in sills in the Pongola Supergroup dated here to 2724 ± 7 Ma.
As the grabens were filled with the Platberg Group, sedimentation and volcanism continued after a small hiatus in the Pniel Group sometime after ca. 2709 Ma, with basaltic andesites preserved as the Allanridge LIP likely erupted between ca. 2709 Ma and ca. 2683 Ma, which was possibly related to the same mantle plume responsible for the Platberg volcanic province. Once again, the LIP itself, uplift and erosion of the Platberg Group rocks, and emplacement of komatiites and high-MgO basalts at the base of the Allanridge Formation basaltic andesites provide evidence for a mantle plume. These volcanic rocks blanketed the underlying western regions of the Kaapvaal Craton, while mafic dike swarms were emplaced at 2701-2683 Ma on the eastern region of the Kaapvaal Craton in the Rykoppies Dike Swarm. The Rykoppies Dike Swarm radiates out from a focus underneath the eastern Bushveld Complex and was likely produced by a mantle plume at 2701-2683 Ma (Olsson et al., 2011), if the three pulses of magmatism proposed by Gumsley et al. (2016) are reinterpreted and reduced to two (Olsson et al., , 2011 at 2701-2683 Ma and at 2664-2654 Ma, with the older pulse being responsible for the Allanridge LIP. The 2664-2654 Ma pulse was possibly associated with another LIP and renewed rifting in the so-called proto-basins to the Transvaal Supergroup and the White Mfolozi Dike Swarm, with some magma injected along the same trends as those of the Rykoppies Dike Swarm, possibly under the same stress regime or lines of weakness in the crust. This terminated volcanism on the Kaapvaal Craton during the proto-basinal phases to the Transvaal Supergroup, which included failed rifting.

CONCLUSION
Baddeleyite from two mafic sills interpreted as feeders of Klipriviersberg magmatism, hosted within the Witwatersrand Supergroup on the central Kaapvaal Craton, yields U-Pb ID-TIMS crystallization ages of 2787 ± 2 Ma and 2789 ± 4 Ma. Complementary U-Pb LA-ICPMS analysis of baddeleyite grains from one of the sills yielded an similar crystallization age of 2793 ± 7 Ma. All dates above overlap within error, and a 2791-2779 Ma age is interpreted for emplacement of these sills and the larger Klipriviersberg LIP. Baddeleyite from a mafic sill in the Pongola Supergroup on the southeastern Kaapvaal Craton yielded a crystallization age of 2724 ± 7 Ma. This extends magmatism associated with the Platberg Group farther east to the southeastern Kaapvaal Craton. Additionally, an east-northeast-trending mafic dike from the southeasternmost Kaapvaal Craton, part of the newly defined Ulundi Dike Swarm, produced a U-Pb baddeleyite crystallization age of 2729 ± 7 Ma, which makes it part of the first dike swarm directly related to the Ventersdorp Supergroup (specifically the Platberg Group). Our results extend the 2791-2779 Ma magmatic event into the central Kaapvaal Craton and 2754-2709 Ma magmatism into the southeastern Kaapvaal Craton.
These new data, combined with previous data associated with the 2701-2654 Ma Rykoppies and White Mfolozi magmatism, allow us to distinguish two separate LIPs within the Ventersdorp Supergroup and shortly thereafter: at 2791-2779 Ma in the Klipriviersberg LIP, at 2701-2683 Ma in the Allanridge LIP, and potentially at 2664-2654 Ma in the White Mfolozi Dike Swarm after the formation of the Ventersdorp Supergroup but before formation of the Transvaal Supergroup. The first two LIPs initiated and terminated Ventersdorp Supergroup development across the entire Kaapvaal Craton, which can be associated with magmatism of the Klipriviersberg Group and the Pniel Group, respectively. These LIPs are separated by rifting and graben development with associated sedimentation and bimodal volcanism during the formation of the 2754-2709 Ma Platberg Group during the arrival of the Allanridge mantle plume at the base of the lithosphere. These new ages also allow us to reinterpret the timing of gold deposition within the Witwatersrand Supergroup, with gold deposition in the conglomerates occurring before ca. 2.79 Ga.