Emplacement dynamics and geochemistry of the ca. 2.08 Ga Devarabanda dyke swarm, Dharwar craton: exploring sub-swarm variations
Srinjoy Datta1, Sayandeep Banerjee1, Amiya K. Samal1*, Rajesh K. Srivastava1, Richard E. Ernst2, Aishwarya Mohan1
1Department of Geology, Banaras Hindu University, Varanasi- 221005, India
2Department of Earth Sciences, Carleton University, Ottawa, Canada K1S 5B6
Extracted and updated from:
Datta, S., Banerjee, S., Samal, A.K., Srivastava, R.K., Ernst, R.E., Mohan, A., 2024. Insights into the magma dynamics of the multi-pulsed ca. 2.08 Ga Devarabanda dyke swarm, eastern Dharwar craton: Constraints from integrated geochemical, magnetic fabric and emplacement studies. Precambr. Res., 411, p.107507. https://doi.org/10.1016/j.precamres.2024.107507
Introduction:
Large Igneous Provinces (LIPs) essentially comprise volcanic flood basalt provinces, underlying plumbing architecture and magmatic underplating (Ernst, 2014). The flood basalt exposures of Precambrian LIPs are generally lost due to denudation, leaving behind the plumbing architecture as the only representative exposures (Ernst, 2014; Ernst et al., 2019; Srivastava et al., 2022). Therefore, only plumbing system of a LIP (dyke swarms and sill complexes) is available for its identification and additional studies. Further, it is accepted that dyke exposures close to the plume centre are generally emplaced under the stronger influence of the primary magmatic stress field, whereas those at greater distances are controlled by the regional stress regimes (e.g., Baragar et al., 1996; Ernst et al., 2001; Hou et al., 2010). However, emplacement mechanisms of dyke swarms can be more complicated than previously perceived (cf. Macdonald et al., 2010). Linear dyke swarms are generally emplaced through rifting and are perpendicular to the spreading direction. However, the stress fields associated with a radiating dyke swarm are usually more complex. The nature of influence, and the interplay between the local magmatic stress fields, regional stress regimes, and pre-existing structures in the emplacement of a radiating dyke swarm is still poorly constrained.
The Dharwar craton hosts several generations of Paleoproterozoic mafic dyke swarms that span in age from ca. 2.37 to ca. 1.79 Ga (Söderlund et al., 2019). Among these, the ca. 2.08 Ga Devarabanda dyke swarm stands out due to its remarkable radiating geometry, which is evident in three distinct sub-swarms: SS-1 (N- to NNE-trending), SS-2 (NW-trending), and SS-3 (NE-trending) (Fig. 1). This swarm is exposed around the Proterozoic Cudappah basin, and is suggested to have a proximal plume centre located beneath the basin (Kumar et al., 2015; Parashuramulu et al., 2022). The unique geometry of the swarm, and its proximity to the plume centre, deemed it as an appropriate candidate for investigating the influence of the regional and local (magmatic) stress regimes in the dyke emplacement. Further, the chemistry and emplacement mechanisms may often vary within a large radiating dyke swarm (Baragar et al., 1996). Therefore, investigating the three sub-swarms of the ca. 2.08 Ga Devarabanda dyke swarm offers an opportunity to explore potential variability in geochemistry and emplacement systems. We further integrated constraints from geochemistry, geophysics and anisotropy of magnetic susceptibility (AMS) to develop a comprehensive understanding about this dyke swarm. Geochemistry was primarily used to characterize the source and to verify the genetic relationship between the three sub-swarms. AMS characteristics were employed to study the magma flow dynamics and infer the emplacement direction. Further geophysical evidence was used to identify the sub-surface features which might have played a role in the dyke emplacement processes. Our study indicates the presence of multiple magma pulses and independent emplacement events for the three sub-swarms. This dyke swarm is also linked with opening of the Proterozoic Cudappah basin, further enhancing its tectonic significance. In the present work, the new geochemical and AMS data were also compared with available geochemical data (Samal et al., 2021; Parashuramulu et al., 2022) and magnetic fabric data (Kumar et al., 2015), substantially increasing the confidence in our inferences.
Figure 1: Regional geological map (after, Geological Survey of India, 1981; Halls, 1982; Nagaraja Rao et al., 1987; Murty et al., 1987; Rao et al., 1990; Pandey et al., 1997; Radhakrishna and Joseph, 1998; Moyen et al., 2003; French and Heaman, 2010; Samal et al., 2019) displaying exposures of distinct ca. 2.08 Ga Devarabanda dyke swarms (after Kumar et al., 2015, Söderlund et al., 2019, Samal et al., 2021, Parashuramulu et al., 2022), eastern Dharwar craton. The black dashed rectangles demarcate the spatial extent of the tree sub-swarms (SS-1, SS-2, and SS-3). The red star denotes the approximate location of the plume centre for the ca. 2.08 Ga dyke swarm.
Geochemical variation and mantle source characteristics
The initial observations were made from petrographic investigation, where the SS-3 samples hosted substantial primary olivine phenocrysts, and a coarser average grain size than the other two sub-swarms (Fig. 2). The primitive mantle normalized multielement spider diagram and chondrite normalized rare earth element (REE) plot suggests that the SS-3 samples were more depleted than the SS-1 and SS-2 dykes (Fig. 3). Furthermore, the SS-3 dykes had Mg# (50 - 68.4), Ni (150 – 298 ppm) and Cr (159 – 316 ppm) concentrations, which indicate their more primitive nature. This was also substantiated by available geochronological data, which suggest SS-3 to have comparatively older age (2087 ± 5 Ma) (Söderlund et al., 2019), than the other two subswarms (2080 ± 6 to 2083 ± 1 Ma) (Kumar et al., 2015; Söderlund et al., 2019). The differences in the major and trace element compositions between the three sub-swarms is shown in Table 1.
Figure 2: Photomicrographs of the mafic dyke samples from the three sub-swarms of ca. 2.08 Ga dyke swarm. (a and b) Photomicrographs with characteristic subophitic to ophitic relationship between of clinopyroxene (cpx) and plagioclase (pl) under plane polarized light (PPL), and between crossed Nicols (BCN), respectively. (c) Photomicrographs displaying the SS-3 samples with substantial modal concentration of primary olivine phenocrysts under BCN.
Figure 3: (a) Chondrite - normalized (Sun and Mcdonough, 1989) REE plot of the ca.2.08 Ga dyke swarm samples. The SS-1 samples were taken from Samal et al. (2021). One sample from SS-3 was plotted after Parashuramulu et al. (2022). (b) Primitive mantle - normalized (Sun and Mcdonough, 1989) multi element spider diagram of the ca. 2.08 Ga Devarabanda dyke swarm.
Table 1: Comparison of geochemistry and AMS parameters between the three sub-swarms of the ca. 2.08 Ga Devarabanda dyke swarm. |
||||
Subswarm |
SS-1 |
SS-2 |
SS-3 |
Unit |
---|---|---|---|---|
TiO2 |
1 - 1.7 |
0.8 - 2.7 |
0.5 - 1 |
wt% |
Al2O3 |
14.2 - 14.5 |
12.7 - 15.7 |
14.5 - 17 |
wt% |
Fe2O3 |
13 - 15.5 |
9.8 - 16.8 |
1.5 - 4 |
wt% |
MgO |
5.4 - 6.5 |
2.9 - 5.3 |
7.5 - 9.9 |
wt% |
SiO2 |
47.3 - 51.1 |
48.1 - 55.6 |
47.2 - 49.5 |
wt% |
CaO |
7.7 - 10 |
6.4 - 8.7 |
9.2 - 10.6 |
wt% |
Total Alkali |
2.4 - 4.6 |
2 - 5.4 |
1.5 - 3.9 |
wt% |
Mg# |
40.5 - 50.9 |
25.7 - 45.5 |
50.0 - 68.4 |
- |
Ni |
80 - 130 |
31 - 80 |
150 - 298 |
ppm |
Cr |
80 - 210 |
28 - 98 |
159 - 316 |
ppm |
Zr |
63 - 142 |
45.6 - 205 |
29.5 - 132 |
ppm |
Hf |
1.7 - 3.5 |
1.2 - 5.5 |
0.8 - 3.3 |
ppm |
Nb |
3.0 - 8.0 |
2.4 - 20 |
1.7 - 9.2 |
ppm |
Ta |
0.2 - 0.5 |
0.2 - 0.8 |
1.2 - 0.7 |
ppm |
Ce |
17.3 - 42.5 |
9.2 - 51.6 |
7 - 25.3 |
ppm |
Nd |
10.5 - 24 |
6.0 - 26.6 |
4.1 - 15.6 |
ppm |
Rb |
21 - 101 |
10.4 - 146.6 |
4.5 - 69.5 |
ppm |
Ba |
368 - 223 |
223 - 368 |
91.4 - 301.4 |
ppm |
Th |
0.7 - 2.2 |
0.7 - 2.2 |
0.3 - 2.2 |
ppm |
Sr |
145 - 309 |
197.1 - 371.9 |
280.6 - 672.3 |
ppm |
(La/Yb)CN |
2.2 - 4.7 |
2.4 - 7.2 |
2.2 - 3.1 |
- |
(Dy/Yb)CN |
1.0 - 1.3 |
1.1 - 1.5 |
1.2 - 1.4 |
- |
(Gd/Yb)CN |
1.0 - 1.6 |
1.2 - 1.9 |
1.4 - 1.6 |
- |
(La/Sm)CN |
1.6 - 2.3 |
1.38 - 3.07 |
1.59 - 1.69 |
- |
Km |
904-132830 |
585 - 19210 |
4236-36060 |
μSI |
K1Inc |
vertical -inclined |
vertical -inclined |
inclined - horizontal |
|
T |
-0.5 - 0.3 |
-0.2 - 0.6 |
-0.7 - 0.8 |
|
Pj |
1.007 - 1.08 |
1.007 – 1.08 |
1.02 – 1.05 |
|
Total Alkali: Na2O (wt%) + K2O (wt%) Mg#: molar Mg/(Mg+Fe)*100 Km, K1 Inc, and T: AMS parameters see S1 for more detials SS-1 samples after Samal et al., (2021) SS-3 One sample (DG 56) after Parashuramulu et al., (2022) |
Although there were some subtle geochemical differences in the major and trace element geochemistry, all the three sub-swarms displayed similar trends in the multielement and REE plots (Fig. 3). All the samples had inclined REE slope, with distinct LREE enrichment and flat HREE pattern (Fig. 3a). This suggested a limited contribution of garnet in the mantle source. Further, the multielement spider diagram displayed strong Nb-Ta-Ti negative anomaly, flat Zr-Hf pattern and positive Eu and Sr anomaly, which were dignified as distinct crustal signatures (Fig. 3b). All the three sub-swarms followed a similar differentiation trend of the major oxides and had overlapping HFSE values (Fig. 4 a-d) suggesting a common parental melt for all three sub-swarms. Figure 4e and f suggest that all samples were evolved under influence of olivine, clinopyroxene and Cr-spinel fractionation. However, the varying fractionation trends suggest differences in the degree of magma evolution in the magma chambers. The low Gd/Yb values and TiO2/Yb values for all the three sub-swarms suggested their derivation from a shallow spinel dominated lherzolite source. This was better displayed graphically in the (La/Sm)CN vs. (Gd/Yb)CN plot, and the Nb/Yb vs. TiO2/Yb plot (Fig. 5 a, b). Further, the robust LIP printing diagram (Pearce et al., 2021) attributed observed crustal signatures in the studied samples to subduction related metasomatism of the mantle source (Fig. 5c). Corroborating the geochemical evidences, we suggest that the three subswarms were derived from a common shallow mantle that was probably metasomatised by an earlier subduction event, prior to the dyke emplacement. Thereafter, they underwent magmatic differentiation in crustal magma chambers. The SS-3 was then emplaced as an earlier primitive pulse, and the other two subswarms were emplaced as later differentiated pulses.
Figure 4: (a, and b) Variance of major oxide concentration against MgO (wt.%). (c, and d) Variation of HFSE contents against MgO (wt.%). (e and f) Trace element variation against MgO (wt.%). The differentiation trends are demarcated by the black arrows, and the orange dashed arrows represent the overlapping concentration of the HFSEs in the three sub-swarms. The symbols for the three sub-swarms are same as Figure 3.
Figure 5: Trace element ratio plots estimating the mantle source, (a) (La/Sm)N vs. (Gd/Yb)N plot (after Adhikari et al., 2021), (b) Nb/Yb vs. TiO2/Yb plot (after, Pearce, 2008), (c) LIP printing diagram (after, Pearce et al., 2021). The sample symbols are same as Figure 3.
Magnetic fabrics and structural framework
The three sub-swarms also show distinct differences in their magnetic fabrics. The SS-3 samples were characterised by sub-horizontal to inclined magnetic lineation suggesting horizontal magma injection, whereas the SS-1 and SS-2 samples dominantly displayed sub-vertical to inclined magnetic lineation indicating vertical magma injection (Fig. 6; Table 1). Further, the AMS fabrics in the host granites were compared with the dyke geometry to test the influence of the pre-existing strain fabrics in the dyke emplacement events (Fig. 7). The SS-1 and SS-2 dyke geometries have a low angle relationship with the AMS fabrics of the associated granite samples, indicating a favorable influence of the pre-existing strain fabrics in the dyke emplacement systems. On the contrary, the SS-3 dykes have a high angle relationship with the granite fabrics suggesting its emplacement under a stronger influence of magmatic stress regime related to the mantle plume, probably through primary fracture opening. Additionally, we analysed the gravity anomaly map of the Dharwar craton (Fig. 8; reference), to identify any sub-surface structures that might have affected the dyke emplacement events. We observed good spatial correlation of a linear gravity low zone (GL-1) with the SS-2 dykes. It was inferred as a pre-existing shear zone which might have facilitated the SS-2 emplacement. Earlier workers have related this zone as the Kurduwadi Lineament (Kale and Peshwa, 1988), which is also suggested as a major tectonic boundary in recent studies (Vasanthi and Santosh, 2021). Moreover, a gravity high zone (GH-1) was identified around the SS-1 exposures which might be related to probable crustal thinning and active rifting during the dyke emplacement. On the contrary, no such sub-surface feature was associated with the SS-3 dykes, indicating minor or no influence of pre-existing structural framework in its emplacement.
Figure 6: Lower hemisphere equal area projection of principal anisotropy axes of the three sub-swarms of ca. 2.08 Ga dyke swarm. The K1 axes are denoted by red squares, K2 axes are marked with green triangles, and K3 axes are marked by blue circles. The empty symbols mark the mean susceptibility axes. The magnetic foliation plane (K1-K2) is marked by the red great circle, and the dyke plane is marked by the grey thick line. ‘ζ’ denotes the angle between the dyke plane and the magnetic foliation plane, and ‘n’ is the number of cylindrical core specimens analysed. The SS-1 samples were taken from Kumar et al., (2015), and the SS-2 and SS-3 samples belong to the present study.
Figure 7: Comparison of AMS fabrics of the host granites and the associated mafic dyke samples. The SS-1 and SS-2 dyke planes show low angle variation with the granite fabrics, and the SS-3 dyke plane displays high angle relation with the granite fabrics.
Figure 8: Enhanced gravity anomaly map (GMSI 2006). GH-1 represents a gravity high located in the SW part of Cudappah Basin, that is inferred as a possible paleo plume centre. A linear SE trending gravity low (GL-1) zone located to the west of the Cudappah basin, is inferred as a pre-existing shear, and a long wave gravity high (GH-2) located north of the Cudappah basin, is attributed to crustal thinning and active rifting, probably associated with the emplacement of SS-1.
Emplacement model
Corroborating our evidence, an emplacement model (Fig. 9) has been proposed for this dyke swarm. We suggest that a plume head located underneath the present-day Cudappah basin, might have resulted in a domal uplift, that was accommodated by concentric fractures and normal faults. The magma generated from the metasomatised lithospheric mantle was guided by these normal faults and stored in crustal magma chambers of different configuration, surrounding the central domal uplift. The SS-3 was emplaced in the pre-rift phase, as a primitive pulse by lateral magma injection into radial fractures related to the domal uplift. The SS-2 and SS-1 were emplaced as subsequent magma pulses in the later rift phase primarily by vertical magma injection. The SS-2 emplaced along a pre-existing Kurduwadi lineament, probably reclaiming it as a rift zone, whereas the SS-1 was emplaced by active rifting and crustal thinning. The spatial variation and multiple pulses in the dyke swarm indicate diverse emplacement systems having a combined influence of pre-existing strain fabrics and the magmatic stress regime related to the plume centre. The emplacement model supports the apical graben concept (Baragar et al., 1996), suggesting subsidence of the domal uplift as a probable precursor to the opening of the Proterozoic Cudappah basin.
Figure 9: 3-D emplacement model of the ca. 2.08 Ga dyke swarm. (a) The pre-rift phase, where SS-3 dykes are emplaced as an earlier pulse from its feeder magma chamber via lateral magma flow and primary opening of radiating fractures (green dashed lines). (b) The rift phase, where greater magma effusion results in active rifting and vertical emplacement of the SS-1, whereas the SS-2 dykes are guided by the NW trending pre-existing shear. The direction of magma flow is denoted by red arrows.
Conclusions
The studied dyke swarm display distinct intra-swarm geochemical variation; the SS-3 having more primitive composition as compared to the SS-1 and SS-2 sub-swarms. All the three sub-swarms were derived from a common parental melt generated from a metasomatised spinel lherzolite mantle source.
AMS studies indicate the SS-3 to be emplaced by lateral magma flow, whereas the SS-1 and SS-2 were emplaced by vertical magma injection. Comparison with the granite fabrics indicate a favorable relationship of the pre-existing strain fabrics in the host granites in the emplacement of the SS-1 and SS-2 samples, but the SS-3 was emplaced under the influence of the primary magmatic stress regime.
An emplacement model was proposed, which suggest the magma to have mobilized to crustal magma chambers surrounding a central domal uplift due to the underlying plume, where the magma evolved independently. The SS-3 was emplaced as a primitive magma pulse by lateral magma injection, during an earlier pre-rift stage. Whereas the SS-1 and SS-2 were emplaced by vertical magma injection in a later rift stage. The SS-1 was emplaced by active rifting and crustal thinning, and the SS-2 was emplaced along a pre-existing Kurduwadi lineament. Eventually the domal uplift subsided to form an apical graben that acted as a precursor to the Cudappah basin.
References
Baragar, W.R.A., Ernst, R.E., Hulbert, L., Peterson, T., 1996. Longitudinal petrochemical variation in the Mackenzie dyke swarm, northwestern Canadian Shield. Journal of Petrology. 37(2), 317-359.
Ernst, R.E., 2014. Large igneous provinces. Cambridge University Press.
Ernst, R.E., Liikane, D.A., Jowitt, S.M., Buchan, K.L., Blanchard, J.A., 2019. A new plumbing system framework for mantle plume-related continental Large Igneous Provinces and their mafic-ultramafic intrusions. Journal of Volcanology and Geothermal Research, 384, 75-84.
Ernst, R.E., Grosfils, E.B., Mege, D., 2001. Giant dike swarms: Earth, venus, and mars. Annual Review, Earth and Planetary Science, 29 (1), 489–534.
French, J.E., Heaman, L.M., 2010. Precise U-Pb dating of Paleoproterozoic mafic dyke swarms of the Dharwar craton, India: Implications for the existence of the Neoarchean supercraton Sclavia. Precambrian Research 183, 416–441.
Geological Survey of India, 1981. Geological and Mineral Map of Karnataka and Goa, 1:500,000 Scale Hyderabad.
GMSI., 2006. Gravity map series of India, scale 1:2000,000. Hyderabad: NGRI Publication.
Halls, H.C., 1982. The importance and potential of mafic dyke swarms in studies of geodynamic processes. Geoscience Canada. 82(03), 142-154.
Hou, G., Kusky, T.M., Wang, C., Wang, Y., 2010. Mechanics of the giant radiating Mackenzie dyke swarm: a paleostress field modeling. Journal of Geophysical Research, Solid Earth. 115, B02402.
Kale, V.S., Peshwa, V.V., 1988. Manifestations of the Kurduwadi lineament based onremote sensor data interpretation. Current Science. 57 (23), 1293–1295.
Kumar, A., Parashuramulu, V., Nagaraju, E., 2015. A 2082 Ma radiating dyke swarm in the Eastern Dharwar Craton, southern India and its implications to Cuddapah basin formation. Precambrian Research. 266, 490-505.
Macdonald, R., Baginski, B., Upton, B.G.J., Pinkerton, H., MacInnes, D.A., MacGillivray, J.C., 2010. The Mull Palaeogene dyke swarm: Insights into the evolution of the Mull igneous centre and dyke-emplacement mechanisms. Mineralogical Magazine 74 (4), 601–622.
Moyen, J.-F., Nédélec, A., Martin, H., Jayananda, M., 2003. Syntectonic granite emplacement at different structural levels: the Closepet granite, South India Journal of Structural Geology. 25, 611–631.
Murthy, N.G.K., 1987. Mafic dyke swarms of the Indian shield. In: Halls, H.C., Fahrig, W.F. (eds.), Mafic Dyke Swarms. Geological Association of Canada Special Paper 34, 393–400.
Nagaraja Rao, B.K., Rajurkar, S.T., Ramalingaswamy, G., Ravindra Babu, B., 1987. Stratigraphy, Structure and Evolution of the Cuddapah Basin. In: Radhakrishna, B.P. (ed.), Purana Basins of Peninsular India (Middle to Late Proterozoic). Geological Society of India Memoir 6, pp. 33–86.
Pandey, B.K., Gupta, J.N., Sarma, K.J., Sastry, C.A., 1997. Sm–Nd, Pb–Pb and Rb–Sr geochronology and petrogenesis of the mafic dyke swarm of Mahbubnagar, South India: implications for Paleoproterozoic crustal evolution of the Eastern Dharwar Craton. Precambrian Research 84, 181–196.
Parashuramulu, V., Yadav, P., Sarma, D.S., Babu, N.R., 2022. Geochemistry of~ 2.08 Ga radiating mafic dyke swarm from the Dharwar Craton, India, and their implications on initiation of the Cuddapah Basin. Journal of Earth System Science 132(1), 4
Pearce, J.A., 2008. Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust. Lithos, 100(1-4), 14-48.
Pearce, J.A., Ernst, R.E., Peate, D.W., Rogers, C., 2021. LIP printing: Use of immobile element proxies to characterize Large Igneous Provinces in the geologic record. Lithos, 392. 106068.
Radhakrishna, T., Joseph, M., 1998. Geochemistry and petrogenesis of the Proterozoic dykes in Tamil Nadu, southern India: another example of the Archaean lithospheric mantle source. International Journal of Earth Sciences 87, 268–282.
Rao, J.M., Rao, G.V.S.P., Patil, S.K., 1990. Geochemical and paleomagnetic studies on the Middle Proterozoic Karimnagar mafic dyke swarm, India. In: Parker, Rickwood, Tucker (eds.), Mafic Dykes and Emplacement Mechanisms, 373–382.
Samal, A.K., Srivastava, R.K., Rahaman, W., 2021. Sr-Nd isotope geochemistry and petrogenesis of ca. 2.26–2.25 Ga and ca. 2.08 Ga mafic dyke swarms from the Dharwar craton, India: Insights into their mantle sources and geodynamic implications. Lithos, 406, 106503.
Söderlund, U., Bleeker, W., Demirer, K., Srivastava, R.K., Hamilton, M., Nilsson, M., Pesonen, L.J., Samal, A.K., Jayananda, M., Ernst, R.E., 2019. Emplacement ages of Paleoproterozoic mafic dyke swarms in eastern Dharwar craton, India: Implications for paleoreconstructions and support for a~ 30◦ change in dyke trends from south to north. Precambrian Research, 329, 26–43.
Srivastava, R.K., Ernst, R.E., Buchan, K.L., de Kock, M., 2022a. An overview of the plumbing systems of large igneous provinces and their significance. Geological Society of London. Special Paper. 518, 1-16.
Sun, S.S., McDonough, W.F., 1989. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. Geological Society of London, Special Publication 42(1), 313-345.
Vasanthi, A., Santosh, M., 2021. Lithospheric architecture and geodynamics of the Archean Dharwar craton and surrounding terranes: New insights from satellite gravity investigation. Gondwana Research. 95, 14–28.