Precambrian Research 224 (2013) 11–22
Contents lists available at SciVerse ScienceDirect
Precambrian Research
journa l homepage: www.elsevier.com/locate/precamres
Trading partners: Tectonic ancestry of southern Africa and western Australia, in
Archean supercratons Vaalbara and Zimgarn
a,b,∗ c d,e f g
Aleksey V. Smirnov , David A.D. Evans , Richard E. Ernst , Ulf Söderlund , Zheng-Xiang Li
a
Department of Geological and Mining Engineering and Sciences, Michigan Technological University, Houghton, MI 49931, USA
b
Department of Physics, Michigan Technological University, Houghton, MI 49931, USA
c
Department of Geology and Geophysics, Yale University, New Haven, CT 06520, USA
d
Ernst Geosciences, Ottawa K1T 3Y2, Canada
e
Carleton University, Ottawa K1S 5B6, Canada
f
Department of Earth and Ecosystem Sciences, Division of Geology, Lund University, SE 223 62 Lund, Sweden
g
Center of Excellence for Core to Crust Fluid Systems, Department of Applied Geology, Curtin University, Perth, WA 6845, Australia
a r t i c l e i n f o
a b s t r a c t
Article history: Original connections among the world’s extant Archean cratons are becoming tractable by the use of
Received 26 April 2012
integrated paleomagnetic and geochronologic studies on Paleoproterozoic mafic dyke swarms. Here we
Received in revised form
∼
report new high-quality paleomagnetic data from the 2.41 Ga Widgiemooltha dyke swarm of the Yil-
19 September 2012
garn craton in western Australia, confirming earlier results from that unit, in which the primary origin
Accepted 21 September 2012
of characteristic remanent magnetization is now confirmed by baked-contact tests. The correspond-
Available online xxx ◦ ◦ ◦
ing paleomagnetic pole (10.2 S, 159.2 E, A95 = 7.5 ), in combination with newly available ages on dykes
from Zimbabwe, allow for a direct connection between the Zimbabwe and Yilgarn cratons at 2.41 Ga,
Keywords:
Paleomagnetism with implied connections as early as their cratonization intervals at 2.7–2.6 Ga. The proposed “Zimgarn”
Paleoproterozoic supercraton was likely distinct from Vaalbara (Kaapvaal plus Pilbara) at 2.4 Ga, but both of those entities
Yilgarn independently fragmented at ca. 2.1–2.0 Ga, reassembling into the Kalahari and West Australian cratons
Pilbara by 1.95–1.8 Ga.
Zimbabve craton © 2012 Elsevier B.V. All rights reserved.
Kaapvaal craton
1. Introduction The Widgiemooltha mafic dykes of the Yilgarn craton (West-
ern Australia) constitute one of the most prominent dyke swarms
Fundamental changes across the entire Earth system during of the world’s Precambrian shields (e.g., Sofoulis, 1966; Parker
the Archean–Paleoproterozoic transition (ca. 3000–2000 Ma) are et al., 1987). The dykes generally trend east-west and outcrop
reflected in global dynamics that became more similar to modern- over the entire craton (Fig. 1). Most of the dykes are of picrite,
style plate tectonics through that interval (Condie and Kröner, olivine–dolerite or quartz–dolerite composition and have not been
2008; Reddy and Evans, 2009), although the kinematic record from affected by significant metamorphic events since their formation
that era is currently limited by a dearth of high-quality paleomag- (e.g., Hallberg, 1987).
netic data (Evans and Pisarevsky, 2008). However, recent work The age of the Widgiemooltha dyke swarm is constrained by
using a combined paleomagnetic and geochronologic approach on isotopic dating of its three largest dykes. Nemchin and Pidgeon
mafic dyke swarms within several Archean cratons, for example (1998) reported a baddeleyite U–Pb age of 2418 ± 3 Ma for the
Slave (Buchan et al., 2009, 2012), Superior (Evans and Halls, 2010; westernmost extension of the ∼600-km-long Binneringie Dyke
defining member of the Superia supercraton of Bleeker, 2003), and (Fig. 1). Electron microprobe U–Pb analyses of baddeleyites from
Vaalbara (de Kock et al., 2009), has started to change this situation, an eastern location of the Binneringie Dyke yielded a slightly
enabling us to test and quantify the original connections that may younger age of 2410.3 ± 2.1 Ma (French et al., 2002). A baddeleyite
−
have existed between them. U–Pb age of 2410.6 + 2.1/ 1.6 Ma was also obtained for the Cele-
bration Dyke (Doehler and Heaman, 1998). Fletcher et al. (1987)
reported Rb–Sr and Sm–Nd isochron age of 2411 ± 38 Ma for the
ultramafic core of the ∼180-km-long Jimberlana intrusion, located
∗
Corresponding author at: Department of Geological and Mining Engineering and
south of the Binneringie Dyke (Fig. 1). These data broadly con-
Sciences, Michigan Technological University, 630 Dow ECE Building, 1400 Townsend
firm the earliest geochronological study of the Widgiemooltha
Drive, Houghton, MI 49931, USA. Tel.: +1 906 487 2365.
dyke swarm (Turek, 1966), which yielded a 2420 ± 30 Ma Rb–Sr
E-mail address: [email protected] (A.V. Smirnov).
0301-9268/$ – see front matter © 2012 Elsevier B.V. All rights reserved. http://dx.doi.org/10.1016/j.precamres.2012.09.020
12 A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22
-20o processing techniques, and it lacked a rigorous field test to confirm
the primary origin of dyke magnetization. Here we report new high-
Pilbara
quality paleomagnetic data from the Widgiemooltha dyke swarm,
with implications for tectonic style of craton amalgamation and
separation during the Archean–Paleoproterozoic transition.
-24o
2. Methods 23
A 24 Yilgarn
We sampled 27 sites representing 14 separate Widgiemooltha
10
11,12 19,20 dykes spanning most of the Yilgarn craton (Fig. 1). Five to thirty-
9
two field-drilled oriented core samples were collected from each
-28o C,I 22 18
26 21 site. Orientation was done with both solar and magnetic compasses.
17
16 D,G Low-field thermomagnetic analyses, remanence unblocking tem-
5 7
15 perature and magnetic hysteresis measurements indicate the
H E 14
F 13 presence of pseudosingle-domain magnetite or low-Ti titanomag-
CD
RD netite in most dykes (Fig. 2).
BD Magnetic remanence measurements were conducted at Yale
o Perth 4 JD
-32 University and Michigan Technological University using automated
1 three-axis DC SQUID 2G rock magnetometers housed in magnet-
27 2 3
ically shielded environments. After measurement of the natural
remanent magnetization (NRM), the samples were cycled through
∼
the Verwey transition at 120 K (Verwey, 1939) by immersing
them into liquid nitrogen for about two hours in order to reduce
o a viscous component carried by larger magnetite grains (Schmidt,
-36
116o 120o 124o
1993). Next, a low alternating field (AF) pre-treatment (to 100 mT in
20 mT steps) was applied to remove any remaining low-coercivity
Fig. 1. Locations of sampling sites of this study (circles) and of Evans’ (1968) (dia-
viscous or isothermal remanence. Finally, 15–20 thermal demagne-
monds) sampling sites. Open (or grey) symbols show the sites used (or not used)
tization steps were performed in an inert (nitrogen) atmosphere.
to calculate the mean paleomagnetic directions. Shaded areas show the Yilgarn and
Pilbara cratons. Thin gray lines indicate the Widgiemooltha dyke swarm; solid lines Progressive demagnetization was carried out until the magnetic
show the Binneringie (BD), Jimberlana (JD), Celebration (CD), and Randalls (RD) intensity of the samples fell below noise level or until the measured
dykes. Inset map shows the area studied by Evans (1968). ◦
directions became erratic and unstable (typically at 580–590 C).
The characteristic remanent magnetization (ChRM) for samples
− −
11 1 displaying nearly linear demagnetization trajectories was isolated
age (using a decay constant of 1.39 × 10 yr ); and the consis-
using principal-component analysis (PCA) (Kirschvink, 1980). The
tency of isotopic ages attests to the excellent preservation of the
dykes. best-fit line was used if defined by at least three consecutive demag-
netization steps that trended toward the origin and had a maximum
Magnetic surveys show both positive and negative anoma-
◦
angle of deviation (MAD) less than 20 . For the majority of dykes,
lies associated with the Widgiemooltha dykes (e.g., Tucker and
◦
the ChRM was isolated within a narrow (10–25 C) temperature
Boyd, 1987), suggesting the presence of dual-polarity remanent
◦
range above 520 C (Figs. 3–5). The mean directions were calculated
magnetization. The only previous paleomagnetic study of the
using Fisher statistics (Fisher, 1953). A site mean was accepted for
Widgiemooltha dyke swarm (Evans, 1968) found remanence direc-
further calculations if it was obtained from three or more samples,
tions of both polarities from four independently cooled dykes.
◦
the confidence circle (˛ ) was smaller than 15 , and the precision
However, that study was confined to a limited area in the southeast- 95
parameter k was greater than 50.
ern Yilgarn (Fig. 1), it did not employ modern demagnetization and a b 900 15
Mrs/Ms = 0.143
Hc = 159 Oe
600 Hc = 408 Oe
0
(SI units) 300 κ
0 Magnetic Moment (memu) -15
-200 0 200 400 600 -0.5 0.0 0.5 Temperature (oC) Magnetic field (kOe)
◦ ◦
Fig. 2. (a) A typical low-field thermomagnetic curve ( (T)) measured in argon from −192 C to 700 C using an AGICO KLY-4S magnetic susceptibility meter at Yale University.
◦
The (T) curve indicates the presence of a single magnetic phase with a Curie temperature between 570–580 C (low-Ti titanomagnetite or magnetite). The presence of nearly
− ◦
stoichiometric magnetite is further supported by a characteristic peak observed at about 153 C, associated with the Verwey (1939) transition. (b) A typical magnetic
hysteresis loop indicating pseudo-single domain behavior of the sample. Abbreviations are Hc, coercivity; Hcr, coercivity of remanence; Mrs, saturation remanence; Ms,
saturation magnetization.
A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22 13 a b
510
540 590
585 c d
525 545
575 585
e
Fig. 3. Paleomagnetic results from the Binneringie dyke. (a–d) Typical orthogonal vector plots of thermal demagnetization (vertical/horizontal projections shown by
open/filled symbols). Equal-area plots show the accepted paleomagnetic directions (open squares) and their mean (solid circle) with the 95% confidence circle (˛95). Numbers
◦
indicate temperature steps (in C). (e) The site-mean directions from this study (circles) and from Evans (1968) (triangles) used to calculate the mean direction (star) and its
˛95 circle (solid line) for the Binneringie dyke (Table 1). Open square shows the direction from Site K (Evans, 1968) collected from baked granite about five meters from the
Binneringie dyke. Although based on only two independently oriented samples (Table 1), the site K direction hints at the primary origin of the dyke magnetization.
3. Results geomagnetic field at the Archean–Proterozoic boundary (Smirnov
et al., 2011).
Eighteen sites from thirteen separate dykes yielded well- Although Evans (1968) found consistent directions between
◦
defined (˛95 < 15 ) site mean directions, including twelve sites dykes and baked granite host rocks at two sites, those are not
(from eight different dykes) that yielded directions consistent strictly positive baked-contact tests because they lack stable, dif-
with those reported by Evans (1968) (Fig. 6a, Table 1), and the ferent results from unbaked host rocks. We augment the aggregate
◦ ◦
combined dataset yields a paleomagnetic pole at 10.2 S, 159.2 E Widgiemooltha dataset with positive baked-contact tests con-
◦
(A95 = 7.5 ). The angular dispersion of virtual geomagnetic poles firming the primary age of magnetization in these dykes. First is the
◦
(VGPs) S = (11.7 ± 1.9) from the Widgiemooltha dykes is lower inverse baked contact test from Narrogin quarry (Site 1) near the
◦
than that of 15–17 observed at the comparable latitude band for western end of the Binneringie Dyke, where the dyke is cross-cut
the last 5 million years (e.g., Johnson et al., 2008). However, the by a thin, younger, east-west dyke bearing a remanence direction
low S value is consistent with the systematically low VGP disper- identical to those of the Marnda Moorn large igneous province
sion values measured from other reliable paleomagnetic datasets of (LIP), also known in the Narrogin region as Wheatbelt or Fraser
the Neoarchean and Paleoproterozoic age, suggesting more stable dykes (Pisarevsky et al., 2003; Wingate and Pidgeon, 2005) (Table 2;
14 A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22 ) for the
◦
( is
95 9.9 7.5 N used
10.9 A 14.8
sites
of
longitude.
K 38 35 13.1 28 27 48 and
number
the
is
latitude
B
E)
◦ ( site
p distribution.
340.7 340.1 ˚ 342.1 339.9 353.3 336.8 336.3 327.9 332.2 338.3 341.6 324.8 357.3 338.0 339.8 346.6 347.1 342.4 333.7 335.3 344.1 335.9 325.6 336.0 339.2 the
VGP
are directions.
s for
˚
,
s N)
. ◦ ( p 1.3 1.0 7.0 2.8 5.2 9.3 8.3 8.6 6.2 6.1 8.2 08.6 20.0 10.7 10.2 16.3 22.1 –8.2 17.9 26.7 14.7 13.9 17.9 19.9 21.3 − − − parameter
(1968)
paleomagnetic
for
Evans
) ◦
( e e e e e e e e e from 95 1.7 1.3 4.0 7.0 5.2 9.2 2.8 4.3 4.6 2.7 4.0 4.0 2.4 5.8 2.5 5.8 4.8
– – – – – – – – 10.9 ˛ 12.5 − concentration
parameter
sites
the
the
and
9 9571 9.5 26 25 75 7527 8.9 11.9 71 76 51 24 25 86 92 500 160 731 137 166 119 166 562 146 145 189 797 172 866 117 circle denote
1607 1212 1.6
k concentration
the
letters
and
confidence
study, )
◦ 95% 60.0 40.2 64.3 62.3 59.8 64.1 57.6 75.5 68.0 65.4 62.7 67.1 42.8 69.5 61.0 66.5 68.0 67.8 65.2 56.4 68.8 81.4 66.1 54.4 61.2 67.6 68.1 67 65.5 66.6 48.5 73.5 (
circle I
− − − − − − − − − − − − − − − − − − − − − − − − − − − − − this
the
for
are
K
confidence )
◦ and (
collected
95 260.5 D 248.1 236.2 225.1 233.5 245.8 217.3 278.5 243.0 252.5 251.6 262.7 264.8 354.6 299.0 352.4 228.3 356.9 238.1 266.1 239.0 245.3 247.5 95%
A
sites the
directions. direction.
are (VGP).
mean.
k
B 7 5 247.8 5 8 9 original the mean pole
and scattered for the 95
˛ total directions.
used (I);
the
489 44 2 258.9 272.3 7 75 246.0 3 3 244.5 3 6 130.4 8 5 7 78 6 5 7 3 7 335.5 51.0 denote yielded 10 N 11 19 16 242.5 of
Not
geomagnetic
). scattered
which
1968 inclination 8,
numbers
virtual
E) calculation yielded
◦ ID: and and the (
s for Evans, 4,
of
(
120.52 120.59 120.59 ˚ 116.95 117.11 118.54 121.42 122.32 121.7 121.7 122.1 121.74 121.31 121.8 122.0 121.7 121.8 121.59 121.97 121.93 121.33 119.92 119.45 121.97 122.38 122.0 121.8
(D) Site
2,
which
used
Dyke 16,
sites
dykes. not
longitude
and
S)
calculation. declination was ◦
( 13, and
s 30.90 30.90 30.9 31.53 31.7 31.6 32.16 32.18 31.32 29.46 28.63 28.06 26.38 26.36 27.83 28.04 and 12,
Binneringie
mean Binneringie
23
the
latitude the Sites
total directions. Widgiemooltha
Site b
g
only for for the
from the
paleomagnetic the
our
only) mean are
for Dyke to
P the
study
away from
˚ the
, Hill polarities m Dyke p
study
(1968) of
calculated calculated This 2 5
, Dyke – the
inverted
results (this calculate
of b was was
equivalent Evans
to
is was (1968) granite
mean
in
)
King Dyke Dyke direction.
Dyke Celebration Randalls
calculation d d d d d 25
used d,f a,d
d
27 33.0 11 19 20 21 23 18 7 D 32.1 E 32.1 13269, CI K 5 32.94 32.49 F 31.80 G 31.7 1415 17 32.2 22 32.1 24 A H 30.18 32.1 1968 direction direction Evans for
baked
mean
direction
provided c
plus paleomagnetic H
used
from
mean mean
Evans, samples
of ( not Binneringie
of dikes not A Site – – Binneringie – Jimberlana
is
Widgiemooltha site site
study ID 1
95 Sample No No Site ˛ Site The f c a e Site Binneringie Mean Mean Total Jimberlana This Other Mean g b d between-site
Table Summary number a
A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22 15
a Site 7 c N
0 500
-45
E D
575 -90 G W E F 7 S
-45 b
-45 0 180
0
Fig. 4. Paleomagnetic results from the Jimberlana dyke. (a) A typical orthogonal vector plot of thermal demagnetization (vertical/horizontal projections shown by open/filled
◦
symbols) from Site 7. Numbers indicate temperature steps (in C). (b) Accepted paleomagnetic directions (open squares) and their mean (solid circle) with the 95% confidence
circle (˛95) for Site 7 (equal-area projection). (c) The mean direction for the Jimberlana dyke (open square) and its ˛95 area (solid circle) calculated from Site 7 (this study)
and Sites D, E, F, and G (Evans, 1968) (Table 1).
Fig. 6b). Although the younger Narrogin dyke is undated, its orien- extension of the Muggamurra swarm, a part of the Mardna Moorn
tation and paleomagnetic directions suggest a Marnda Moorn age LIP (Myers, 1990; Wingate and Pidgeon, 2005) into the northern
that would imply the Binneringie remanence pre-dates 1.2 Ga. Yilgarn. Alternatively, these sites may represent Widgiemooltha
An additional positive baked contact test that was conducted dykes that recorded transitional field directions. One of Jimberlana
at Booylgoo Hills (Site 21) where a Widgiemooltha dyke intrudes dyke sites (Site 5) resulted in a mean direction that is significantly
an Archean lava sequence provides more direct evidence of a different from the consistent directions obtained from the other
primary magnetization in the Widgiemooltha swarm (Table 2; Jimberlana sites (Table 1; Fig. 4c). We speculate that this direc-
Fig. 6c). These data indicate that the Booylgoo dyke carries a pri- tion may represent a local remagnetization event. Finally, Sites 19
mary magnetic remanence, whose direction is identical to other and 22 yielded well-defined northerly directions with intermediate
Widgiemooltha dykes and thus confirms that the Widgiemooltha upward inclination possibly reflecting a recent field overprint.
grand mean pole is representative of the Yilgarn craton at 2.41 Ga.
The combined result from our study and that of Evans (1968) mer- 4. Discussion and conclusions
its a perfect 7 point classification on Van der Voo’s (1990) Q-scale
of paleomagnetic reliability. A moderate to high paleolatitude for Yilgarn is directly com-
Six additional sites yielded site-mean directions significantly parable to that of the Zimbabwe craton, for which recent U–Pb
different from the Widgiemooltha grand mean (Table 1). The mean dating allows reconstruction at 2.41 Ga according to the published
directions from three dykes (Sites 17, 18, and 20) are close to the Sebanga Poort Dyke virtual geomagnetic pole (Mushayandebvu
Fraser dyke mean (Pisarevsky et al., 2003) and to the directions et al., 1995; see Appendix A). The moderate uncertainties in both
yielded by the intruding thin dyke at Narrogin (Table 2; Fig. 6b) thus poles permit a variety of allowable juxtapositions between the two
offering a possibility that these dykes may represent an eastward cratons, which we propose were joined together in a supercraton
Table 2
Summary of paleomagnetic directions used for the baked contact tests shown in Figure 6b,c. N is the number of samples used to calculate the paleomagnetic declination (D)
and inclination (I); ˛95 and k are the 95% confidence circle and the concentration parameter for paleomagnetic directions.
◦ ◦ ◦
Site ID N D ( ) I ( ) k ˛95 ( )
Total Widgiemooltha mean (Table 1) 9 247.5 −66.6 117 4.8
Narrogin inverse baked contact test
Site 1 (>5 m from the intruding dyke) 11 248.1 −62.3 731 1.7
−
Site 1 (<5 cm from the intruding dyke) 1 322.8 72.5 – –
1 325.7 −73.2 – –
−
Samples within the intruding dyke 1 326.3 61.6 – –
1 331.1 −63.4 – –
1 355.7 −58.9 – –
Fraser Dyke mean (Pisarevsky et al., 2003) 23 332.6 −74.1 110 2.9
Booylgoo Hills baked contact test
Widgiemooltha dyke 3 265.4 −60.0 44 18.7
Archean lavas at 3 m form the contact 7 266.1 −64.8 54 8.3
Archean lavas at >9 m from the contact 15 333.8 −40.0 25 7.2
16 A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22
Site 11 Site 14 a b -90 -90
N N 555 -45 -45
535 575 W E
600 0 W E 0 180 180 S S
Site 15 Site 21 c d -90 -90
N 500 N 535 W E -45 -45 585
0
565 0 180 W E 180 S S
Site 23 Site 24 e f -90 -90
N N 475 525 -45 -45
0 585 0
W 565 E 180 W E 180
S S
Fig. 5. Paleomagnetic results from the other dykes used to calculate the grand Widgiemooltha direction (Table 1). Typical orthogonal vector plots of thermal demagnetization
◦
(vertical/horizontal projections shown by open/filled symbols). Numbers indicate temperature steps (in C). Equal-area plots show the accepted paleomagnetic directions
˛
(open squares) and their mean (solid circle) with the 95% confidence circle ( 95) for each site.
“Zimgarn” (Fig. 7). We emphasize that until additional paleomag- the sweeping pattern of younger craton stabilization from west to
netic data are obtained from the two cratons, a precise model for east across Zimbabwe (Jelsma and Dirks, 2002). Also, sinistral off-
Zimgarn will remain elusive; nonetheless, we tentatively choose set on major terrane-bounding shear zones in the eastern Yilgarn
a preferred juxtaposition guided by several additional geometric (Czarnota et al., 2010) can be restored to a proto-Zimgarn configu-
constraints. First, we align the Sebanga Poort Dyke approximately ration (Fig. 7b) that nearly aligns the proposed Kurnalpi magmatic
parallel to the Widgiemooltha dykes. Such parallelism is allowed arc at ca. 2700 Ma with another proposed arc crossing Zimbabwe
paleomagnetically because both dyke suites have remanence decli- at the same time interval (Jelsma and Dirks, 2002).
nations roughly in line with the dyke strikes. However, even after Postulation of a Zimgarn supercraton presents a novel dilemma
satisfying this constraint, within the uncertainties of the poles, any for Neoarchean reconstructions: are dyke swarms or basement
relative position of the two cratons is permitted (i.e., they could provinces more diagnostic of original connections between extant
be translated around each other). Second, we restrict such trans- cratons? Bleeker and Ernst (2006) presented the simple geomet-
lation so that the 2.57 Ga Great Dyke of Zimbabwe projects off the ric argument that as vertical features, dykes provide piercing
margin of Yilgarn, because it has not yet been identified in Western points whose positions are not sensitive to current exposure level;
Australia. This constraint restricts the cratons to be either north in addition, analogy with more completely preserved radiating
or south of each other in paleogeographic coordinates at 2.41 Ga dyke swarms (e.g., Central Atlantic Magmatic Province) shows
(Fig. 7a). Third, we choose the northern option for Zimbabwe rel- that the rift-related dykes can show nearest-neighbor connections
ative to Yilgarn, because the more juvenile Neoarchean character even among suture-bounded provinces of different age. This rea-
of the eastern Yilgarn (Czarnota et al., 2010) is a better match for soning provides a rationale for juxtaposing cratons with distinct
A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22 17
0o 0o a b Narrogin inverse baked contact test
intruding dyke
< 5 cm from intruding dyke
o
o o Fraser dyke mean
90 90û
-90
-90
> 0.5 m from 0o intruding dyke
Booylgoo Hills baked contact test
Archean lava 180û > 9 meters 180û c from contact
Dyke
o
o
3 meters 90
-90 from contact
180o
Fig. 6. (a) Stereonet of accepted site-mean paleomagnetic directions (diamonds – this study; triangles – Evans, 1968). Thick and thin circles (also in (b and c)) show the mean
˛
between-site direction and its 95 cone calculated from all accepted directions (Table 1). (b) The inverse baked contact test at the Narrogin quarry (Site 1; see text). Diamond
shows the mean ChRM direction for Site 1 (the ˛95 is not shown). Squares show the ChRM directions for the samples taken from the Binneringie Dyke within 5 cm from the
contact with the intruding dyke. Triangles show the ChRM directions yielded by the intruding dyke. The mean Fraser dyke direction (star) is calculated for Site 1 based on
the paleomagnetic pole from Pisarevsky et al. (2003). (c) The Booylgoo baked contact test (see text). Square and triangle show the mean direction from the Archean lavas
sampled at 3 m and >9 m from the contact, respectively.
cratonization histories, such as Yilgarn and Zimbabwe, which were cratonic stabilization (Dirks and Jelsma, 1998; Frei et al., 1999).
previously grouped into separate supercratons (Bleeker, 2003). Second, Zimgarn and Vaalbara magmatic “barcodes” share no
The distinctive 2.41-Ga age match between Widgiemooltha and known matches through the interval 2.6–2.0 Ga interval (e.g., Ernst
Sebanga Poort dykes (Söderlund et al., 2010) may be diagnostic and Bleeker, 2010). Third, paleomagnetic latitudes are somewhat
of their original connection despite the mismatch of cratonization distinct between the two supercratons at 2.69 Ga (although uncer-
ages for Yilgarn and Zimbabwe. tainties are large; Table 3).
The Zimgarn-like assemblage might have been joined to Superia Zimgarn probably fragmented into separately drifting
from Neoarchean cratonization to early Paleoproterozoic separa- Zimbabwe and Yilgarn cratons sometime between 2.2 and
tion, as proposed by Söderlund et al. (2010), but precisely coeval 2.0 Ga (Figs. 7c and 8a). Two rift- to passive-margin successions
paleomagnetic poles at 2.41 Ga are not yet available from Supe- are evident: the Magondi Supergroup in Zimbabwe, which is older
rior or adjoining cratons to test that hypothesis. In the Scotland than 2.0 Ga (McCourt et al., 2001) and may indicate rifting as
sector of the North Atlantic craton, the older Scourie dykes are early as 2.26 Ga (Manyeruke et al., 2004); and the Yerrida Basin in
−
nearly the same age (2418 + 7/ 4 Ma; Heaman and Tarney, 1989) northern Yilgarn, that extends from as early as ca. 2.2–2.1 Ga in a
as the Widgiemooltha and Sebanga dykes, but Scourie paleo- sag phase (Windplain Group) to a rift stage (Mooloogool Group)
magnetic remanences have been reset by Laxfordian (ca. 1.7 Ga) that is only constrained to be older than 1.84 Ga (Rasmussen and
metamorphism (Beckmann, 1976; Piper, 1979). Superia lay at low Fletcher, 2002; Pirajno et al., 2004). Despite the uncertainty of
paleolatitudes at 2.45 Ga, as well documented by data from the Mat- associating preserved epicratonic strata with a northern Yilgarn rift
achewan dykes (summarized by Evans and Halls, 2010), but both margin, that margin must have been established prior to accretion
Yilgarn and Zimbabwe lack paleomagnetic results from that older of the Narracoota/Dalgaringa arc (perhaps also the Pilbara craton)
age. in the early stages of the 2.0–1.95 Ga Glenburgh orogeny (Pirajno
The Vaalbara supercraton has recently been quantified for the et al., 2004; Johnson et al., 2011). Both the Yerrida and Magondi
2.8–2.7 Ga interval (de Kock et al., 2009). It is likely that Vaal- successions contain carbonate–evaporite platforms characterized
bara was paleogeographically separated from Zimgarn, for several by the Lomagundi carbon-isotopic excursion; according to our
reasons. First, basement ages and cover successions of Vaalbara model, they are correlative successions along a common cratonic
substantially differ from those of Zimgarn, including extensive margin (Fig. 7c), representing rifting of Zimgarn from other,
carbonate platform development at 2.6 Ga (reviewed by Knoll originally conjoined cratons yet to be identified (possibly Superia,
and Beukes, 2009), when Zimgarn was still undergoing final according to Söderlund et al., 2010). By 1.88 Ga, Zimbabwe and
18 A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22
a Sebanga VGP
Widgie- moo ltha Zimbabwe
6 Zim babwe
Yilgarn Yilgar n
3
Paleo-E quator 2410 Ma 2410 Ma
b v Zimbabwe c v v v Zimbabwe v
v v Yilgarn v 2690 Ma v
v
2100 Ma
Yilgarn paleolatitude and orientation unknown
◦ ◦ ◦
Fig. 7. (a) Zimgarn supercraton at 2410 Ma, achieved by rotation of Yilgarn to Zimbabwe according to Euler parameters (–47 , 077 , +157 ), and paleogeographic restoration
according to the Widgiemooltha dyke paleomagnetic pole determined herein. Greenstone belts are shaded light gray. GD = Great Dyke (paralleled by the Umvimeela and
East Dykes), SPD = Sebanga Poort Dyke, BD = Binneringie Dyke, JD = Jimberlana Dyke. Younger marginal basins and orogens are labeled for reference, despite anachronism.
Present north is indicated for each craton. In the right-side panel, poles are shown with their 95% uncertainty regions. VGP = virtual geomagnetic pole. (b) Tectonic elements
of Zimgarn at 2690 Ma. Paleolatitudes and orientation of Zimbabwe conform to the Reliance pole (Table 3), with polarity chosen to minimize motion between 2690, 2575, and
2410 Ma. Yilgarn is displaced to restore post-2690-Ma sinistral offset on major terrane-bounding shear zones (simplified in this illustration by the generalized shear couple),
thus bringing proposed 2690-Ma volcanic arcs (Jelsma and Dirks, 2002; Czarnota et al., 2010) nearer to alignment. (c) Proposed continuous passive margin extending from
the Yerrida margin (Yilgarn) to the Magondi margin (Zimbabwe) at ca. 2100 Ma. Paleolatitude and orientation of Zimgarn are not constrained at that age.
Yilgarn were widely separated from each other, as indicated by is so widespread across Kaapvaal at 2.06 Ga (Reischmann, 1995;
the moderate versus low paleolatitudes, respectively implied by Cawthorn et al., 2006; Mapeo and Wingate, 2009; Prendergast,
the Mashonaland and Frere poles (Table 3). 2012), has not yet been recognized in Pilbara, suggesting that
Kaapvaal and Pilbara share similar tectonostratigraphic records Vaalbara could have split prior to 2.06 Ga. If so, then the 2.2 Ga
through much of the early Paleoproterozoic Era, such simi- mafic events in Kaapvaal (Cornell et al., 1996) and Pilbara (Müller
larity providing the original basis for the Vaalbara hypothesis et al., 2005) could document their initial separation. Unfortu-
(Cheney, 1996). However, the Bushveld-related magmatism that nately, among the four cratons discussed herein, only Kaapvaal has Vaalbara a SPD b Kalahari breakup Zimbabwe 50° craton Zimbabwe assembly SPD Kaapvaal Pilbara 40° Kaapvaal Yilgarn
BD
JD 30°
West Australian craton
20° BDJD assembly Zimgarn
Yilgarn
breakup 10°
2050 Ma Pilbara 1800 Ma
Fig. 8. Kinematic sequence of the trading cratonic partners. (a) Vaalbara and Zimgarn in the early stages of breakup, 2050 Ma. Only Kaapvaal is constrained paleomagnetically
at this age (Table 3), but other cratons are positioned for consistency in the kinematic sequence. (b) Kalahari and West Australian cratons after final assembly, 1800 Ma;
paleolatitudes and orientations of both cratons are consistent with the ca. 1800 Ma poles listed in Table 3.
A.V. Smirnov et al. / Precambrian Research 224 (2013) 11–22 19
Table 3
Paleomagnetic poles bearing on Zimgarn, Vaalbara, and their transition toward united Kalahari and West Australian cratons.
◦ ◦ ◦ ◦ ◦ ◦
Craton/Paleomagnetic pole Age (Ma) p ( N), p ( E) r ( N), r ( E) A95 ( ) 1234567 Q Test ( ) Reference
Zimgarn: Zimbabwe Ref.
Reliance middle-temp. 2690 44, 128 44, 128 17 1011101 5 f 21 Yoshihara and Hamano (2004)
Great Dyke grand mean 2575 24, 057 24, 057 9 1111100 5 C,c 37 See Appendix A
Sebanga Poort Dyke (VGP) 2410 17, 006 17, 006 14 1001100 3 C 45 See Appendix A
Widgiemooltha dikes 2410 −10, 159 06, 013 8 1111111 7 C,c 50 This study
Vaalbara: Kaapvaal Ref.
Allanridge basalt combined 2685 −70, 346 −70, 346 6 1111101 6 G 36 de Kock et al. (2009)
−
−
Rykoppies E-W dikes only 2685 63, 339 63, 339 7 1111101 6C 37 Lubnina et al. (2010)
Woongarra/Weeli Wolli 2450 −50, 220 −78, 235 5 1111100 5 g 11 Evans (2007) (abstract only)
Ongeluk lava 2220 −01, 101 −01, 101 5 1111101 6 G 17 Evans et al. (1997)
Waterberg UBS-I 2055 37, 051 37, 051 11 1111100 5 F,G 28 de Kock et al. (2006)
Bushveld complex 2050? 19, 031 19, 031 6 1111110 6 f 49 Letts et al. (2009)
Transitional poles: Kaapvaal Ref.
Waterberg UBS-II ∼1950? −11, 330 −11, 330 10 0111110 5 c,g 33 de Kock et al. (2006)
∼
Black Hills dikes 1900? 09, 352 09, 352 5 0111110 5 C 42 Lubnina et al. (2010)
Mashonaland + Mazowe 1880 07, 338 18, 016 7 1111110 6 C 48 Hanson et al. (2011)
Waterberg sills & Soutp. lavas 1875 16, 017 16, 017 9 1111110 6 C 50 Hanson et al. (2004)
−
Frere iron formation (Yilgarn) 1890 45, 220 2 1110111 6 (f) 10 Williams et al. (2004)
Kalahari craton:
Soutpansberg UBS-II ∼1800? 25, 008 25, 008 8 0111110 5 M 39 de Kock (2007)
West Australian craton:
Hamersley overprint HP3 ∼1800? −35, 212 −35, 212 3 0110101 4 – 12 Li et al. (2000)