Next Article in Journal
The Immobility of Uranium (U) in Metamorphic Fluids Explained by the Predominance of Aqueous U(IV)
Next Article in Special Issue
Reconstruction of the Magma Transport Patterns in the Permian-Triassic Siberian Traps from the Northwestern Siberian Platform on the Basis of Anisotropy of Magnetic Susceptibility Data
Previous Article in Journal
Potential Future Alternative Resources for Rare Earth Elements: Opportunities and Challenges
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Thermobarometry of the Rajmahal Continental Flood Basalts and Their Primary Magmas: Implications for the Magmatic Plumbing System

by
Nilanjan Chatterjee
1,* and
Naresh C. Ghose
2,†
1
Department of Earth, Atmospheric and Planetary Sciences, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
2
Department of Geology, Patna University, Patna 800005, India
*
Author to whom correspondence should be addressed.
Current address: G/608, Raheja Residency, Koramangala, Bangalore 560034, India.
Minerals 2023, 13(3), 426; https://doi.org/10.3390/min13030426
Submission received: 1 February 2023 / Revised: 2 March 2023 / Accepted: 16 March 2023 / Published: 17 March 2023
(This article belongs to the Special Issue Large Igneous Provinces: Research Frontiers)

Abstract

:
The Late Aptian Rajmahal Traps originated through Kerguelen-Plume-related volcanism at the eastern margin of the Indian Shield. Clinopyroxene and whole-rock thermobarometry reveals that the Rajmahal magmas crystallized at P-T conditions of ≤~5 kbar/~1100–1200 °C. These pressures correspond to upper crustal depths (≤~19 km). Modeling shows that the Rajmahal primary magmas were last in equilibrium with mantle at P-T conditions of ~9 kbar/~1280 °C. The corresponding depths (~33 km) are consistent with gravity data that indicate a high-density layer at lower crustal depths below an upwarped Moho. Thus, the high-density layer probably represents anomalous mantle. It is likely that the mantle-derived magmas accumulated below the upwarped Moho and were subsequently transported via trans-crustal faults/fractures to the upper crust where they evolved by fractional crystallization in small staging chambers before eruption. In the lower part of the Rajmahal plumbing system, buoyant melts from the Kerguelen Plume may have moved laterally and upward along the base of the lithosphere to accumulate and erode the eastern Indian lithospheric root. The Rajmahal plumbing system was probably shaped by tectonic forces related to the breakup of Gondwana.

1. Introduction

Flood basalt eruptions have important environmental effects, and magma volumes and eruption rates of flood basalts are directly related to the structure of the plumbing system through which mantle-derived magmas reach the Earth’s surface [1,2,3]. The structure of the plumbing system in continental flood basalt provinces may depend on the regional tectonics and proximity to a mantle plume [4,5]. Exposures of dikes, sills and other magma bodies, and geophysical surveys shed light on the structure of the plumbing system [6,7,8,9]. The plumbing system of plume-related continental large igneous provinces (LIP) can be divided into four parts corresponding to four depth levels [3]. These are the asthenospheric mantle level where melt is generated (level 1), the lithospheric mantle level through which magmas ascend and accumulate below the Moho (level 2), the crustal level where magma is transported via dikes and sills into and out of small staging chambers where magma differentiation occurs (level 3), and the surface level where volcanoes develop (level 4).
An LIP related to the activity of the Kerguelen Plume has been recognized in eastern India, southwestern Australia (Bunbury), the eastern Antarctic margin, the Kerguelen Plateau, the Ninetyeast Ridge, the Broken Ridge, and the Naturaliste Plateau [10,11,12,13,14,15] (Figure 1). The oldest basalts (~132 Ma) erupted in southwestern Australia. However, there are some reports of even older Kerguelen-Plume-related volcanic activity at ~145–130 Ma in the eastern Himalayas [16]. The eastern Indian basalts and most of the basalts of the Kerguelen Plateau erupted during the peak of Kerguelen Plume activity at ~120–95 Ma [15]. In eastern India, the Late Aptian (~118–114 Ma) Rajmahal–Bengal Basin–Sylhet Traps (RBST) and the associated dolerite dikes and alkalic-carbonatitic-ultramafic intrusives comprise a large flood basalt province with an area of ~1 million km2 [17,18,19,20,21,22,23,24,25,26]. The Rajmahal and Sylhet Traps are located in the western and eastern parts of the province (Figure 2), and drill-core data indicate that the two are continuous under the Gangetic alluvium of the Bengal Basin [27]. Dolerite dikes with NW to NNW-trend and similar age (118–109 Ma) intrude the Precambrian basement to the southwest and west of the Rajmahal outcrop [28] (Figure 2). Geochemical, isotopic, and textural studies of the RBST basalts have provided clues to the processes of fractional crystallization and crustal contamination in their origin [10,19,21,23,25,26,29,30,31,32]. However, the crystallization and mantle melting conditions of the basalts are yet to be established, and the structure of the pathways through which the basalts erupted is poorly understood. Although geophysical surveys have helped to clarify subsurface structures [33,34], and two- and three-dimensional modeling of gravity data have identified possible magmatic underplating below the lower crust [35,36], the structure of the magma plumbing system in the upper mantle and crust remains unclear. Thermobarometry based on the chemical composition of minerals and whole rock provides the temperatures and depths of magma equilibration, which may be used to complement geophysical studies to understand the structure of the magma plumbing system. This study aims to understand the structure of the Rajmahal plumbing system through thermobarometry based on clinopyroxene compositions analyzed here and whole-rock compositions available from the literature. The P-T conditions of crystallization obtained from thermobarometry are then used to backtrack along the fractionation paths of the basalts and calculate compositions of primary magmas and their equilibration depths.

2. Geological and Geochemical Background

The Rajmahal Traps (area: ~4300 km2) are exposed along the north–south-oriented Rajmahal Hills in the northern part of the eastern margin of the Proterozoic Chotanagpur Gneissic Complex (CGC) (Figure 2). The CGC forms the basement of the Gondwana Supergroup, the uppermost part of which contains the Rajmahal basalts [37,38]. The Rajmahal flows are sub-horizontal with the dip increasing to ~5° eastward at the eastern flank of the hills, and the total thickness of the flows increases from ~230 m in surface exposure to >332 m subsurface in the Bengal Basin [27,35]. The thickness of individual flows varies from <1 m to 85 m [39,40,41]. Thin beds of shale, black shale, mudstone, siltstone, cross-bedded sandstone, and oolite associated with bentonite lenses and volcaniclastic rocks occur between lava flows in the lower one-third of the volcanic sequence, indicating sub-aqueous eruptions during the early phase of volcanism [32,41,42,43]. The upper part of the sequence is devoid of sedimentary and volcaniclastic rocks, and consists of subaerially erupted lava flows [23,32]. The volcanics are dominated by tholeiitic basalt and basaltic andesite with minor trachyandesite, andesite, dacite and rhyolite, and rare orthopyroxene-bearing basalt and andesite [21,23,32]. The tholeiites have been classified into two geochemical groups [10] (Figure 3). The Group I samples have higher Ti/Zr ratios and lower Zr/Y ratios than the Group II samples. According to the total alkali-silica classification (TAS diagram, [44]), the Group I samples are basalt and basaltic andesite, whereas the Group II samples are all basaltic andesite. Isotopic and trace element data suggest that the Group I samples are uncontaminated, whereas the Group II samples with incompatible element patterns showing large-ion-lithophile element (LILE) enrichment and negative Nb-Ta anomalies are variably contaminated with the continental crust (Figure 3c) [10,18,21,23,26,31,32]. The chemical distinctions between the two groups of basalt also correlate with their textures, the early erupting Group I basalts in the lower part of the sequence being phyric with plagioclase and/or clinopyroxene phenocrysts, and the late erupting Group II basalts in the upper part of the sequence being mostly aphyric [23].
The eastern margin of the CGC is characterized by north-trending, east-dipping asymmetric folds that formed during high-grade metamorphism and deformation along the north–south-oriented, mid-Neoproterozoic Eastern Indian Tectonic Zone, a major linear orogen that extends from eastern India through Kerguelen Plateau to eastern Antarctica [46,47,48]. The eastern and western margins of the Rajmahal outcrop are bounded by two major N–S-oriented faults, the Rajmahal and Saithia–Brahmani faults (Figure 2), and the basement near the western margin of the Bengal Basin, adjacent to the Rajmahal Traps, also contains N–S-oriented, E-dipping step-faults that presumably formed during the pre-Rajmahal break-up of Gondwana [27,35,49]. Furthermore, there is a broad and elongated positive gravity anomaly with its axis along the western boundary, and several +10 mGal to +25 mGal gravity highs along either side of the southwestern boundary (Kathikund to Deocha) of the Rajmahal outcrop [35,36] (Figure 2). Multichannel reflection data and elastic thickness structure of the lithosphere also indicate the presence of a north–south-oriented pseudofault under the Rajmahal Traps whose formation is attributed to the activity of the Kerguelen plume [50,51]. Thus, the Rajmahal eruptions may have been guided by faults generated through reactivation of a major, pre-existing N–S-oriented mid-Neoproterozoic lineament during the break-up of Gondwana.

3. Sample Locations and Bulk Compositions

For the purpose of mineral thermobarometry, 40 basalt samples containing phenocrysts were analyzed from the northeastern (Sahibganj, Moti Jharna, Maharajpur, Taljhari, Tinpahar, Kherwa, Berhait), northwestern (Dariachak), and central (Litipara, Amrapara, Kundpahar) sectors of the Rajmahal outcrop (Figure 2, Table 1). The bulk compositions of four of the samples are available from [10]. According to the TAS diagram [44], samples TT3 and LH1 from near Litipara are basalts, and samples LA2 (near Litipara) and RB88-35 (Tinpahar) are basaltic andesites (Figure 3a). According to the Ti/Zr vs. Zr/Y plot [10], samples LA2, TT3, and LH1 belong to the uncontaminated Group I, and sample RB88-35 belongs to the contaminated Group II Rajmahal basalts (Figure 3b). The bulk analyses of these samples show totals higher than ~99 wt% [10], and LOI values determined by [21] for the Rajmahal basalts are <1 wt.%. This indicates that the Rajmahal magmas were essentially anhydrous.

4. Analytical Methods

Textural studies and mineral analyses were performed on a JEOL JXA-8200 Superprobe electron probe microanalyzer (EPMA) at Massachusetts Institute of Technology, Cambridge, MA, USA operating with a 15 kV accelerating voltage, a 10 nA beam current, and 1 μm beam diameter. Typical counting times were 40 s per element that yielded accumulated counts with 1σ standard deviations of 0.3%–1.0% for major elements and 1%–5% for minor elements from counting statistics. The EPMA was calibrated using a set of synthetic and natural standards (DJ35 diopside-jadeite, ALP7 aluminous orthopyroxene, Synthetic anorthite, Amelia albite, Marjalahti olivine, rutile, hematite). The raw data were corrected for matrix effects with the CITZAF package [52].

5. Petrography and Mineral Chemistry

The samples are all of porphyritic basalt with phenocrysts (~2–4 mm) and microphenocrysts (~100–500 μm) surrounded by a fine-grained groundmass (Figure 4). The mineral assemblage in each sample is provided in Table 1. The phenocrysts are dominated by tabular clinopyroxene (augite) and lath-shaped plagioclase. Pigeonite occurs in some samples, notably from Berhait and Kherwa. Olivine is present in only one sample from Tinpahar. Rare K-feldspar occurs in one sample from Kundpahar and two samples from Amrapara. The average phenocryst content is ~5 vol.% with subequal amounts of augite and plagioclase. The groundmass consists of clinopyroxene and plagioclase. Ilmenite and hematite are ubiquitous in the groundmass and two samples contain maghemite. Some of the samples contain small spherules filled with zeolites.
The composition of the phenocrysts and microphenocrysts are summarized in Table 2 and presented in detail in Table S1. Clinopyroxene does not show compositional zoning with the exception of two samples from Dariachak in the northwest. It is dominantly augite, and its composition range in samples from the northeast, northwest and central sectors are En45–56Fs8–21Wo29–39, En48–54Fs9–15Wo36–40, and En47–54Fs7–19Wo31–41, respectively (Table 2). Augite from the central sector are slightly higher in the wollastonite component compared to the northeastern sector (Figure 5a). The latter is also slightly higher in jadeite, though the jadeite and aegirine contents are low (≤2 mole%, Figure 5b). The composition range of pigeonite is En42–74Fs18–46Wo4–12 considering all samples.
Plagioclase is mildly zoned with higher albite contents near the rim, especially in samples from the northwest and central sectors (Figure 5c). Its composition ranges from labradorite to bytownite, and rarely andesine (Table 2, An57–82Ab18–42Or0–4, An45–71Ab28–54Or0–1, and An51–81Ab19–48Or0–1, in the northeast, northwest, and central sectors).
The olivine in a Tinpahar (northeast) sample has a composition of Fo71–76Fa24–30. The composition ranges of ilmenite and hematite are Hem0–7Ilm85–97Pph1Gk1–9, and Hem43–92Ilm6–53Pph0–12Gk0–6, respectively (Table 2, Figure 5d), considering all samples.

6. Thermobarometry

The pressure–temperature (P-T) of crystallization of the Rajmahal basalts were calculated using clinopyroxene compositions, clinopyroxene-bulk equilibria, and whole-rock (bulk) compositions.

6.1. Clinopyroxene Thermobarometry

This method uses only clinopyroxene composition to determine the P-T of crystallization. The T-dependent barometric expression (Equation (32a)) and P-dependent thermometric expression (Equation (32d)) of [54] were solved simultaneously to obtain P-T. These equations are based on multiple regression of clinopyroxene compositions (854 analyses for (32a), 910 analyses for (32d)) obtained from anhydrous partial melting experiments on basalts in the P-T range of 1 bar–75 kbar/800–2200 °C. The equations use the enstatite–ferrosilite and diopside–hedenbergite components and cation proportions of clinopyroxene calculated on the basis of six oxygen atoms. The quoted uncertainties are ±3.1 kbar and ±58 °C for clinopyroxene crystallizing from anhydrous melts [54]. In addition, the P-T conditions were also calculated with the random forest machine learning-based algorithms of [55,56] that provided independent estimates of the uncertainties. The two studies [55,56] use the same methodology, but [56] uses an expanded dataset that includes alkalic liquids. These thermobarometers are also based on clinopyroxene compositions obtained from partial melting experiments on basalts that cover a P-T range of 0.002–30 kbar/750–1250 °C.
Using the augite compositions in Table S1 and Equations (32a) and (32d) of [54], the P-T of clinopyroxene crystallization are 1 bar–5.4 kbar (±3.1 kbar) and 1134–1195 °C (±58 °C) (Table 3, negative pressure values >−2 kbar are considered as 1 bar). Samples from Tinpahar and Moti Jharna in the northeastern sector show the highest average pressures of 5.4 kbar and 3.2 kbar, respectively. All other samples show pressures of <2 kbar. The method of [55] also yields low pressures (2.0–4.6 kbar, highest pressures at Tinpahar and Moti Jharna) but with lower uncertainties (±0.6–2.6 kbar) compared to [54] (Table 3). With the method of [56], the pressures are 0.0–1.8 kbar with uncertainties of less than ±4.5 kbar. The temperature estimates are similar with the methods of [55] (1093–1151 °C, ±19–62 °C) and [56] (1119–1157 °C, ±21–48 °C) compared to [54] (Table 3).

6.2. Clinopyroxene-Bulk Thermobarometry

The P-T conditions of crystallization were determined with the clinopyroxene-anhydrous liquid thermobarometer (Equations (P1) and (T1)) of [57] with quoted uncertainties of ±1.4 kbar and ±27 °C. The Cpx-bulk equilibrium was assessed in terms of the Fe2+-Mg distribution, the equilibrium value of the Cpx-bulk KD(Fe2+-Mg) being 0.28 ± 0.08 [54]. A knowledge of the bulk FeO content is necessary to apply this thermobarometer. The oxidation state of the bulk can be determined from the olivine-liquid [58] and magnetite–ilmenite [59] equilibria. Equation (8) of [58] provides a method to calculate the bulk Fe3+/ΣFe ratio from the olivine-melt FeT-Mg distribution (where FeT is total Fe) in samples containing equilibrium olivine. The bulk Fe3+/ΣFe ratio can also be calculated from the magnetite–ilmenite equilibrium that yields oxygen fugacity [59], which can be converted to Fe3+/ΣFe using Equation (6b) of [60]. Unfortunately, bulk composition of the olivine-bearing sample RB88-39 is not available, and magnetite is absent in all of the samples analyzed. So, the bulk Fe3+/ΣFe ratio could not be determined, and assumed values were used to assess Cpx-bulk equilibrium.
Assuming a bulk Fe3+/ΣFe ratio of 0.1, augite in samples LA2, TT3 and LH1 shows equilibrium Cpx-bulk Fe2+-Mg distribution (KD(Fe2+-Mg) values between 0.29 and 0.32, Table 4). Assuming bulk Fe3+/ΣFe ratios of zero and 0.2 change the Cpx-bulk KD(Fe2+-Mg) to 0.26–0.29 and 0.32–0.35, respectively, which are also within the variability of the equilibrium value [54]. Hence, P-T of the abovementioned samples were calculated using the clinopyroxene-liquid equilibrium with the formulations of [57]. Using the augite compositions in Table S1 and bulk compositions in [10], Equations (P1) and (T1) of [57] yielded P-T of 1 bar–4.6 kbar (±1.4 kbar) and 1137–1185 °C (±27 °C) for samples LA2, TT3, and LH1. These P-T results are similar to the P-T obtained from augite composition for the same samples (0.1–3.0 kbar, 1155–1177 °C) with the formulations of [54] (Table 4).

6.3. Whole Rock Thermobarometry

This method uses the predicted pressure-dependent fractionation path of a basalt to determine its pressure of crystallization. During fractionation, the basaltic melt is in equilibrium with one or several different mineral phases. At a multiple saturation point (MSP), the melt is in equilibrium with multiple phases. Parameterized expressions of experimental data [61] predict the composition of a basaltic liquid at its plagioclase lherzolite (Ol + Pl + Cpx + Opx) multiple saturation point (PL-MSP) as a function of pressure and bulk composition. The PL-MSP provides an anchor point for the fractional crystallization path that is defined by the olivine control line, and the Ol-Pl and Ol-Pl-Cpx cotectic boundaries whose intersection is determined by the method of [62]. The bulk composition of the sample and its PL-MSPs at different pressures are plotted in the Ol-Pl-Cpx (from Qz) and Ol-Cpx-Qz (from Pl) pseudoternary projections of the basalt tetrahedron according to the methods of [63,64]. If the sample shows displacement from its fractionation path, it contains excess crystal accumulation. Subtracting the excess crystals from the bulk results in a corrected bulk composition, and the sample plots on its fractionation path at a specific pressure. The corrected bulk composition represents the composition of the melt derived directly from its primitive parental magma by fractional crystallization at the specific pressure. The estimated pressure is accurate within ±2.5 kbar [61]. Using the corrected composition of the melt, its temperature is then calculated with Equation (16) of [54], which has a quoted uncertainty of ±27 °C.
The bulk compositions of the uncontaminated Rajmahal Group I basalts (excluding basaltic andesites) [10,21] and their PL-MSPs at different pressures were plotted in the Ol-Pl-Cpx (from Qz) and Ol-Cpx-Qz (from Pl) pseudoternary projections of the basalt tetrahedron (Figure 6). The basalts show displacement from their projected PL-MSPs toward the plagioclase apex in the Ol-Pl-Cpx diagram (Figure 6a), indicating cumulus enrichment of plagioclase. Subtracting 4%–10% equilibrium plagioclase from the bulk results in corrected bulk compositions of the samples that plot on their PL-MSPs between 1 bar and 4 kbar pressures in Figure 6a. Note that the Ol + Pl + Cpx cotectic boundaries shown in the Ol-Cpx-Qz diagram (Figure 6b) approximately project to the same point as the PL-MSP in Figure 6a, and it is not possible to determine whether a sample plots on its Ol + Pl + Cpx cotectic or on its PL-MSP only from Figure 6a. However, Figure 6b clearly shows that the samples plot on their Ol + Pl + Cpx cotectics at 1 bar–4 kbar (±2.5 kbar) pressures. Note that because Figure 6b is projected from plagioclase, subtracting plagioclase from the bulk has no effect on the sample’s position. The corrected bulk compositions are the compositions of the Ol-Pl-Cpx saturated melts that are derived directly from their primitive parental magmas through fractional crystallization. The temperatures of these melts calculated with Equation (16) of [54] are 1147–1188 °C (±27 °C). Notably, the calculated P-T of samples TT3, LH1, and KP6 (1 bar–2 kbar, 1157–1169 °C) are similar to the P-T determined from augite composition (1 bar–1.7 kbar, 1125–1166 °C) with the formulations of [54] (Table 4).

7. Primary Magma Modeling

Modeling was performed only on the Rajmahal Group I basalts [10,21], whose isotopic compositions indicate that they are uncontaminated by continental crust [10,18,21,26,31]. Because these basalts crystallized at pressures of ≤5 kbar, they evolved by fractional crystallization of olivine, followed by Ol + Pl and Ol + Pl + Cpx from their primitive parental magmas [61,65]. Hence, their primary magmas were modeled through low-pressure reverse fractional crystallization involving addition of Ol + Pl + Cpx (stage 1) and Ol + Pl (stage 2, the olivine-only stage was not necessary) to the bulk [66,67,68,69] (Figure 6). The phase assemblages were added in small steps (step size < 0.5%) with phase proportions and compositions shown in Table 5 and Table S2. Equilibrium Fe2+-Mg distribution between olivine-liquid (KD(Fe2+-Mg) = 0.3, [70]) and Cpx-liquid (KD(Fe2+-Mg) = 0.25), and equilibrium Ca-Na distribution between plagioclase-liquid [65] were maintained at each step of the calculation. In the Ol-Cpx-Qz projection (Figure 6b), the melt moved toward the Ol-Cpx sidebar in stage 1, and toward the olivine apex in stage 2. In the Ol-Pl-Cpx projection (Figure 6a), the melt remained approximately stationary in stage 1, and moved toward the Ol-Pl sidebar in stage 2. The phase proportions and the switching point between stage 1 and stage 2 with constraints from [62] were adjusted so that the melt moved toward its spinel lherzolite MSP (SL-MSP) at high pressures predicted by parameterized expressions of experimental data [66] (uncertainties in P-T at the SL-MSP: ±1.5 kbar, ±11 °C). At the end of the calculation, the melt was in equilibrium with mantle olivine (Fo90–91.5, for different samples), and it plotted exactly on its SL-MSP at a specific high pressure (Figure 6). The result was unique, as any deviation from the phase proportions and switching point between stages 1–2 would result in a melt not on its lherzolite MSP at any pressure, though it may show equilibrium with mantle olivine. In stage 1, MgO, Ni, CaO, and Al2O3 increased, and SiO2, FeOT (total FeO), and Na2O decreased (Figure 7). In stage 2, the oxides and Ni variations followed the same trends except for CaO, which decreased slightly (Figure 7).
The calculations show that the basalts are the products of 47%–61% fractionation of their primary magmas, and the primary magmas were last equilibrated with mantle at a pressure of ~9 kbar and temperatures of 1275–1286 °C (excluding one sample, Table 5 and Table S2). The uncertainties in these results arise from the uncertainties in mineral-melt KD(Fe2+-Mg) and the mantle olivine composition that may vary between Fo88 and Fo92. Assuming a ±10% uncertainty in the Ol-melt and Cpx-melt KD(Fe2+-Mg), and Fo88–92 mantle olivine compositions, the uncertainties in the MgO and FeO contents of the model primary magma are about ±6%, and the uncertainties in P-T are ±15% and ±3% (e.g., 9.0 ± 1.4 kbar, 1280 ± 40 °C). The compositional uncertainties are depicted as error bars in Figure 7.

8. Discussion

The lowermost level (level 1, [3]) of the Rajmahal plumbing system is poorly constrained. It has been suggested on the basis of isotopic data that the Rajmahal basalts originated by melting of a MORB-type mantle source with only heat supplied by the Kerguelen Plume [21,31]. According to [21], viscous drag in a steady-state Kerguelen Plume conduit at the rifted eastern Indian margin may have helped in asthenospheric upwelling and melting of the MORB-type mantle. From inversion of rare-earth element data, [21] calculated a mantle potential temperature of ~1350 °C for the origin of the Rajmahal and Sylhet magmas, indicating a moderately hot mantle at the northern edge of the Kerguelen Plume. The P-T of last equilibration of the Rajmahal primary magmas with the mantle (~9 kbar and ~1280 °C, potential temperatures up to ~130 °C higher [69]) calculated in this study are consistent with the results of [21] and melting in the spinel lherzolite field concluded by [26]. However, [25,26] presented trace element and isotopic evidence to show that the source of the least contaminated Rajmahal and Sylhet basalts was the primitive Kerguelen Plume, as previously suggested by [10], and a MORB-type source was not involved. Using Nd-Sr isotopic data (εNd(I) between −8.6 and +3.2, 87Sr/86Sr(I) between 0.70347 and 0.70965), [25,26] also modeled the origin of the contaminated Rajmahal–Sylhet basalts by mixing between lherzolite-derived melts and the Eastern Ghats granulites, and speculated that the basalts were contaminated through interaction and erosion of the Indian lithospheric root by the Kerguelen Plume. Considering that the axis of the Kerguelen Plume was south of the Rajmahal–Sylhet eruption site [21], the buoyant plume material may have moved laterally and upward through sublithospheric corridors [71,72,73] to accumulate and interact with the eastern Indian lithosphere.
Geophysical studies and primary magma modeling provide clues to the structure of levels 2 and 3 [3] of the Rajmahal plumbing system. Bouguer gravity data indicate the presence of a broad and elongated positive gravity anomaly with its axis along the western boundary of the Rajmahal Traps [35,36] (Figure 2). Three-dimensional gravity modeling and integration with seismic data delineated a high-density (3 g/cm3) layer, 16–18 km thick under the Rajmahal Traps and ~12 km thick in the region to the south, above the 36–38 km deep regional Moho [36]. This indicates that the high-density layer under Rajmahal Traps is at lower crustal depths below an upwarped Moho that may be at a depth of ~20 km. The Rajmahal primary magmas were last in equilibrium with the mantle at a pressure of ~9 kbar that corresponds with a depth of ~33 km considering an average crustal density of 2.8 g/cm3. Thus, the Rajmahal primary magmas were last in equilibrium with the mantle within the high-density layer, which probably represents anomalous mantle within the lower crust.
The level 3 [3] of the Rajmahal plumbing system is characterized by the anomalous mantle at lower crustal depths, dikes and trans-crustal fractures/faults, and upper crustal magma staging chambers. Dikes are rare in the main Rajmahal outcrop [23,40], though NW to NNW-trending basaltic andesite dikes of similar age (118–109 Ma) to the Rajmahal basalts are common in the Raniganj–Giridih–Koderma region to the southwest and west of the Rajmahal outcrop [28]. The north–south oriented Eastern Indian Tectonic Zone (EITZ) is a trans-crustal orogen, as indicated by the high metamorphic pressures (~10 kbar) of the granulites at the western contact of the Rajmahal Traps [46,48]. The subvertical faults associated with the EITZ, the Rajmahal and Saithia–Brahmani boundary faults along the eastern and western margins of the Rajmahal outcrop (Figure 2), and the faults in the basement of the Bengal Basin adjacent to the outcrop [27,35,49] may have acted as pathways for upward ascent of the Rajmahal primary magmas [74,75]. Thermobarometry in this study shows that the Rajmahal magmas crystallized in the upper crust. Robust estimates of the P-T conditions for magma crystallization by three different methods are between 1 bar and ~5 kbar, and 1093–1195 °C (Table 3 and Table 4), indicating a maximum depth of ~19 km considering an upper crustal density of 2.7 g/cm3. Crystallization probably occurred in small, upper crustal magma chambers distributed throughout the Rajmahal province [8,9]. The Rajmahal magmas differentiated through fractional crystallization and upper crustal assimilation [32] in the small, near-surface magma staging chambers before erupting on the surface.
Based on plate tectonic reconstructions [76], and analogy with basin evolution at the southwest Australian rifted margin [21,77] concluded that rifting associated with the breakup of Gondwana at the eastern Indian margin preceded or occurred synchronously with the eruption of the Rajmahal–Sylhet basalts. The recently dated 118–109 Ma old dikes of the Raniganj–Giridih–Koderma area [28] indicate that dike intrusion was synchronous with the emplacement of the Rajmahal Traps, and possibly occurred by exploiting fractures created by Gondwana breakup. Thus, plate tectonics may have played an important role in shaping the structure of the Rajmahal plumbing system [4,5].

9. Conclusions

Thermobarometry shows that the Rajmahal basalts crystallized in the upper crust (P-T of ≤5 kbar and 1093–1195 °C), and the Rajmahal primary magmas last equilibrated with the mantle near the Moho (P-T of ~9 kbar and ~1280 °C). A high-density layer below an upwarped Moho previously discovered through modeling of gravity data is probably anomalous mantle at shallow, lower crustal depths. The P-T results of this study complements the results of gravity modeling and provides a clearer picture of the upper levels of the Rajmahal plumbing system. It is likely that the mantle-derived primary magmas accumulated below the Moho and were subsequently transported via trans-crustal faults/fractures to the upper crust where the magma evolved by fractional crystallization in small staging chambers before erupting on the surface.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/min13030426/s1, Table S1. Chemical composition of minerals in the Rajmahal basalts; Table S2. Bulk composition of the Rajmahal Group I basalts, their primary magmas, and P-T of equilibration. Figure S1. (a) Core and (b) rim compositions of augite in the Rajmahal basalts (En—enstatite, Di—diopside, Hd—hedenbergite, Wo—wollastonite, Quad—quadrilateral components, Jd—jadeite, and Ae—aegirine) plotted according to [53]. Figure S2. (a) Core and (b) rim compositions of plagioclase in the Rajmahal basalts (Ab—albite, An—anorthite, Or—orthoclase). Figure S3. Chemical composition of Fe-Ti oxides in the Rajmahal basalts.

Author Contributions

Conceptualization, N.C.; methodology, N.C.; validation, N.C. and N.C.G.; formal analysis, N.C.; investigation, N.C. and N.C.G.; resources, N.C. and N.C.G.; data curation, N.C.; writing—original draft preparation, N.C.; writing—review and editing, N.C. and N.C.G.; visualization, N.C.; supervision, N.C. and N.C.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data are included in the main text and online Supplementary Materials.

Acknowledgments

We are grateful for the constructive comments of two anonymous reviewers during the peer review process that substantially improved the presentation of the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bryan, S.E.; Ukstins Peate, I.; Peate, D.W.; Self, S.; Jerram, D.A.; Mawby, M.R.; Marsh, J.S.; Miller, J.A. The largest volcanic eruptions on Earth. Earth Sci. Rev. 2010, 102, 207–229. [Google Scholar] [CrossRef] [Green Version]
  2. Burchardt, S. (Ed.) Introduction to volcanic and igneous plumbing systems—Developing a discipline and common concepts. In Volcanic and Igneous Plumbing Systems: Understanding Magma Transport, Storage, and Evolution in the Earth’s Crust; Elsevier: Amsterdam, The Netherlands, 2018; pp. 1–12. [Google Scholar]
  3. Ernst, R.E.; Liikane, D.A.; Jowitt, S.M.; Buchan, K.L.; Blanchard, J.A. A new plumbing system framework for mantle plume-related continental Large Igneous Provinces and their mafic-ultramafic intrusions. J. Volcanol. Geotherm. Res. 2019, 384, 75–84. [Google Scholar] [CrossRef]
  4. Tibaldi, A. Structure of volcano plumbing systems: A review of multi-parametric effects. J. Volcanol. Geotherm. Res. 2015, 298, 85–135. [Google Scholar] [CrossRef]
  5. Van Wyk de Vries, B.; van Wyk de Vries, M. Tectonics and volcanic and igneous plumbing systems. In Volcanic and Igneous Plumbing Systems: Understanding Magma Transport, Storage, and Evolution in the Earth’s Crust; Burchardt, S., Ed.; Elsevier: Amsterdam, The Netherlands, 2018; pp. 167–189. [Google Scholar]
  6. Burchardt, S.; Walter, T.R.; Tuffen, H. Growth of a volcanic edifice through plumbing system processes—Volcanic rift zones, magmatic sheet-intrusion swarms and long lived conduits. In Volcanic and Igneous Plumbing Systems: Understanding Magma Transport, Storage, and Evolution in the Earth’s Crust; Burchardt, S., Ed.; Elsevier: Amsterdam, The Netherlands, 2018; pp. 89–112. [Google Scholar]
  7. Galland, O.; Bertelsen, H.S.; Eide, C.H.; Guldstrand, F.; Haug, Ø.T.; Leanza, H.A.; Mair, K.; Palma, O.; Planke, S.; Rabbel, O.; et al. Storage and transport of magma in the layered crust—Formation of sills and related flat-lying intrusions. In Volcanic and Igneous Plumbing Systems: Understanding Magma Transport, Storage, and Evolution in the Earth’s Crust; Burchardt, S., Ed.; Elsevier: Amsterdam, The Netherlands, 2018; pp. 113–138. [Google Scholar]
  8. Mittal, T.; Richards, M.A. The magmatic architecture of continental flood basalts: 2. A new conceptual model. J. Geophys. Res. Solid Earth 2021, 126, e2021JB021807. [Google Scholar] [CrossRef]
  9. Mittal, T.; Richards, M.A.; Fendley, I.M. The magmatic architecture of continental flood basalts 1: Observations from the Deccan Traps. J. Geophys. Res. Solid Earth 2021, 126, e2021JB021808. [Google Scholar] [CrossRef]
  10. Storey, M.; Kent, R.W.; Saunders, A.D.; Hergt, J.; Salters, V.J.M.; Whitechurch, H.; Sevigny, J.H.; Thirlwall, N.F.; Leat, P.; Ghose, N.C.; et al. Lower Cretaceous rocks on continental margins and their relationship to the Kerguelen plateau. Proc. Ocean Drill. Prog. Sci. Res. 1992, 120, 33–53. [Google Scholar]
  11. Mahoney, J.J.; Jones, W.B.; Frey, F.A.; Salters, V.J.M.; Pyle, D.G.; Davies, H.L. Geochemical characteristics of lavas from Broken Ridge, the Naturaliste Plateau and southernmost Kerguelen Plateau: Cretaceous plateau volcanism in the southeast Indian Ocean. Chem. Geol. 1995, 120, 315–345. [Google Scholar] [CrossRef]
  12. Frey, F.A.; Coffin, M.F.; Wallace, P.J.; Weis, D.; Zhao, X.; Wise, S.W. Origin and evolution of a submarine large igneous province: The Kerguelen Plateau and Broken Ridge, southern Indian Ocean. Earth Planet. Sci. Lett. 2000, 176, 73–89. [Google Scholar] [CrossRef] [Green Version]
  13. Frey, F.A.; McNaughton, N.J.; Nelson, D.R.; Delaeter, J.R.; Duncan, R.A. Petrogenesis of the Bunbury basalt, western Australia: Interaction between the Kerguelen plume and Gondwana lithosphere? Earth Planet. Sci. Lett. 1996, 144, 163–183. [Google Scholar] [CrossRef]
  14. Weis, D.; Frey, F.A. Isotope geochemistry of Ninetyeast Ridge basement basalts: Sr, Nd and Pb evidence for the involvement of the Kerguelen hotspot. Proc. Ocean Drill. Prog. Sci. Res. 1991, 121, 591–610. [Google Scholar]
  15. Coffin, M.F.; Pringle, M.S.; Duncan, R.A.; Gladczenko, T.P.; Storey, R.D.; Müller, R.D.; Gahagan, L.A. Kerguelen hotspot magma output since 130 Ma. J. Petrol. 2002, 43, 1121–1140. [Google Scholar] [CrossRef] [Green Version]
  16. Srivastava, R.K. Early Cretaceous Greater Kerguelen large igneous province and its plumbing systems: A contemplation on concurrent magmatic records of the eastern Indian Shield and adjoining regions. Geol. J. 2022, 57, 681–693. [Google Scholar] [CrossRef]
  17. Talukdar, S.C.; Murthy, M.V.N. The Sylhet traps, their tectonic history and their bearing on problems of Indian flood basalt provinces. Bull. Volcanol. 1970, 35, 602–618. [Google Scholar] [CrossRef]
  18. Baksi, A.K. Petrogenesis and timing of volcanism in the Rajmahal flood basalt province, Northern India. Chem. Geol. 1995, 121, 73–90. [Google Scholar] [CrossRef]
  19. Baksi, A.K.; Barman, T.R.; Paul, D.K.; Ferar, E. Widespread early Cretaceous flood basal volcanism in eastern India: Geochemical data from the Rajmahal-Bengal-Sylhet traps. Chem. Geol. 1987, 63, 133–141. [Google Scholar] [CrossRef]
  20. Kent, R.W.; Saunders, A.D.; Storey, M.; Ghose, N.C. Petrology of the Early Cretaceous flood basalts and dykes along the rifted volcanic margin of eastern India. J. Southeast Asian Earth Sci. 1996, 13, 95–111. [Google Scholar] [CrossRef]
  21. Kent, R.W.; Saunders, A.D.; Kempton, P.D.; Ghose, N.C. Rajmahal basalts, Eastern India: Mantle sources and melt distribution at a volcanic rifted margin. In Large Igneous Provinces: Continental, Oceanic, and Planetary Flood Volcanism; Mahoney, J.J., Coffin, M.F., Eds.; American Geophysical Union: Washington, DC, USA, 1997; Volume 100, pp. 145–182. [Google Scholar]
  22. Kent, R.W.; Pringle, M.S.; Muller, R.D.; Saunders, A.D.; Ghose, N.C. 40Ar/39Ar geochronology of the Rajmahal basalts, India and their relationship to the Kerguelen Plateau. J. Petrol. 2002, 43, 1141–1153. [Google Scholar] [CrossRef] [Green Version]
  23. Ghose, N.C.; Kent, R.W. The Rajmahal basalts: A review of their geology, composition and petrogenesis. Geol. Soc. India Mem. 2003, 53, 167–196. [Google Scholar]
  24. Ray, J.S.; Pattanayak, S.K.; Pande, K. Rapid emplacement of the Kerguelen plume-related Sylhet Traps, eastern India: Evidence from 40Ar-39Ar geochronology. Geophys. Res. Lett. 2005, 32, 1–4. [Google Scholar] [CrossRef]
  25. Ghatak, A.; Basu, A.R. The Sylhet Traps: Vestiges of the Kerguelen plume in northeastern India. Earth Planet. Sci. Lett. 2011, 308, 52–64. [Google Scholar]
  26. Ghatak, A.; Basu, A.R. Isotopic and trace element geochemistry of alkalic–mafic–ultramafic–carbonatitic complexes and flood basalts in NE India: Origin in a heterogeneous Kerguelen plume. Geochim. Cosmochim. Acta 2013, 115, 46–72. [Google Scholar] [CrossRef]
  27. Sengupta, S. Geological and geophysical studies in western part of Bengal basin, India. Amer. Assoc. Petrol. Geol. Bull. 1966, 50, 1001–1017. [Google Scholar]
  28. Srivastava, R.K.; Wang, F.; Shi, W.; Ernst, R.E. Early Cretaceous mafic dykes from the Chhota Nagpur Gneissic Terrane, eastern India: Evidence of multiple magma pulses for the main stage of the Greater Kerguelen mantle plume. J. Asian Earth Sci. 2023, 241, 105464. [Google Scholar] [CrossRef]
  29. Agrawal, J.K.; Rama, F.A. Chronology of Mesozoic volcanic of India. Proc. Indian Acad. Sci. 1976, 84A, 157–179. [Google Scholar] [CrossRef]
  30. Acharya, S.K. Chemical Behavior of the Volcanic Rocks of South Rajmahal, Bihar. Ph.D. Thesis, Patna University, Patna, India, 1988. [Google Scholar]
  31. Mahoney, J.J.; Macdougall, J.D.; Lugmair, G.W.; Gopalan, K. Kerguelen hotspot source for Rajmahal and Ninetyeast Ridge? Nature 1983, 303, 385–389. [Google Scholar] [CrossRef]
  32. Ghose, N.C.; Chatterjee, N.; Windley, B.F. Subaqueous early eruptive phase of the late Aptian Rajmahal volcanism, India: Evidence from volcaniclastic rocks, bentonite, black shales, and oolite. Geosci. Front. 2017, 8, 809–822. [Google Scholar] [CrossRef] [Green Version]
  33. NGRI. Gravity Maps of India Scale 1:5,000,000. NGRI/GPH-1 to 5; National Geophysical Research Institute (NGRI): Hyderabad, India, 1978. [Google Scholar]
  34. Kaila, K.L.; Reddy, P.R.; Mall, D.M.; Venkateswarlu, N.; Krishna, V.G.; Prasad, A.S.S.S.R.S. Crustal structure of the West Bengal basin, India, from deep seismic sounding investigations. Geophys. J. Int. 1992, 111, 45–66. [Google Scholar] [CrossRef]
  35. Mukhopadhyay, M.; Verma, R.K.; Ashraf, M.H. Gravity field and structures of the Rajmahal Hills: Examples of the Paleo-Mesozoic continental margin in eastern India. Tectonophysics 1986, 131, 353–367. [Google Scholar] [CrossRef]
  36. Singh, A.P.; Kumar, N.; Singh, B. Magmatic underplating beneath the Rajmahal Trap: Gravity signature and derived 3D configuration. Proc. Indian Acad. Sci. 2004, 113, 759–769. [Google Scholar] [CrossRef]
  37. Ghose, N.C. Geology, tectonics and evolution of the Chhotanagpur granite gneiss complex, eastern India. In Recent Research in Geology 10; Hindustan Publish. Co.: Delhi, India, 1983; pp. 211–247. [Google Scholar]
  38. Chatterjee, N.; Ghose, N.C. Extensive Early Neoproterozoic high-grade metamorphism in north Chotanagpur Gneissic Complex of the Central Indian Tectonic Zone. Gond. Res. 2011, 20, 362–379. [Google Scholar] [CrossRef]
  39. Raja Rao, R.C.S.; Purushottam, A. Pitchstone flows in the Rajmahal Hills, Santhal Paraganas, Bihar. Geol. Surv. India Rec. 1963, 91, 341–348. [Google Scholar]
  40. Sarbadhikari, T.K. Petrology of the northeastern portion of the Rajmahal Traps. Quar. J. Geol. Min. Met. Soc. India 1968, 60, 151–171. [Google Scholar]
  41. Ghose, N.C.; Singh, S.P.; Singh, R.N.; Mukherjee, D. Flow stratigraphy of a selected sections of Rajmahal basalts, eastern India. J. Southeast Asian Earth Sci. 1996, 13, 83–93. [Google Scholar] [CrossRef]
  42. Deshmukh, S.S. Geology of the area around Taljhari and Berhait, Rajmahal hills, Santhal Paraganas, Bihar. Rep. 22nd Intern. Geol. Cong. 1964, 7, 61–84. [Google Scholar]
  43. Ghose, N.C. Pyroclastic rocks of India in space and time. In Proceedings of the Indian Geological Congress (IGC); Indian Geological Congress: Roorkee, India, 2000. [Google Scholar]
  44. Le Bas, M.J.; Le Maitre, R.W.; Streckeisen, A.; Zanettin, B. A chemical classification of volcanic rocks based on the total alkali-silica diagram. J. Petrol. 1986, 27, 745–750. [Google Scholar] [CrossRef] [Green Version]
  45. McDonough, W.F.; Sun, S.-S. The composition of the Earth. Chem. Geol. 1995, 120, 223–253. [Google Scholar] [CrossRef]
  46. Chatterjee, N. An assembly of the Indian Shield at c. 1.0 Ga and shearing at c. 876-784 Ma in Eastern India: Insights from contrasting P-T paths, and burial and exhumation rates of metapelitic granulites. Precamb. Res. 2018, 317, 117–136. [Google Scholar] [CrossRef]
  47. Chatterjee, N.; Nicolaysen, K. An intercontinental correlation of the mid-Neoproterozoic Eastern Indian Tectonic Zone: Evidence from the gneissic clasts in Elan Bank conglomerate, Kerguelen Plateau. Contrib. Mineral. Petrol. 2012, 163, 789–806. [Google Scholar] [CrossRef]
  48. Chatterjee, N.; Banerjee, M.; Bhattacharya, A.; Maji, A.K. Monazite chronology, metamorphism–anatexis and tectonic relevance of the mid-Neoproterozoic Eastern Indian Tectonic Zone. Precamb. Res. 2010, 179, 99–120. [Google Scholar] [CrossRef]
  49. Biswas, S.K. Mesozoic volcanism in the east coast basins of India. Ind. J. Geol. 1996, 68, 237–254. [Google Scholar]
  50. Desa, M.A.; Ramana, M.V.; Ramprasad, T.; Anuradha, M.; Lall, M.V.; Kumar, B.J.P. Geophysical signatures over and around the northern segment of the 85oE Ridge, Mahanadi offshore, eastern continental margin of India and their tectonic implications. J. Asian Earth Sci. 2013, 73, 460–472. [Google Scholar] [CrossRef]
  51. Ratheesh-Kumar, R.T.; Windley, B.F.; Sajeev, K. Tectonic inheritance of the Indian Shield: New insights from its elastic thickness structure. Tectonophysics 2014, 615–616, 40–52. [Google Scholar] [CrossRef]
  52. Armstrong, J.T. CITZAF—A package for correction programs for the quantitative electron microbeam X-ray analysis of thick polished materials, thin-films and particles. Microbeam Anal. 1995, 4, 177–200. [Google Scholar]
  53. Morimoto, N.; Fabries, J.; Ferguson, A.K.; Ginzburg, I.V.; Ross, M.; Seifert, F.A.; Zussman, L.; Aoki, K.; Gottardi, G. Nomenclature of pyroxenes. Mineral. Mag. 1988, 52, 535–550. [Google Scholar] [CrossRef]
  54. Putirka, K.D. Thermometers and barometers for volcanic systems. Rev. Mineral. Geochem. Mineral. Soc. Amer. 2008, 69, 61–120. [Google Scholar] [CrossRef]
  55. Higgins, O.; Sheldrake, T.; Caricchi, L. Machine learning thermobarometry and chemometry using amphibole and clinopyroxene: A window into the roots of an arc volcano (Mount Liamuiga, Saint Kitts). Contrib. Mineral. Petrol. 2022, 177, 10. [Google Scholar] [CrossRef]
  56. Jorgenson, C.; Higgins, O.; Petrelli, M.; Bégué, F.; Caricchi, L. A machine learning-based approach to clinopyroxene thermobarometry: Model optimization and distribution for use in Earth sciences. J. Geophys. Res. Solid Earth 2022, 127, e2021JB022904. [Google Scholar] [CrossRef] [PubMed]
  57. Putirka, K.D.; Johnson, M.; Kinzler, R.; Walker, D. Thermobarometry of mafic igneous rocks based on clinopyroxene-liquid equilibria, 0–30 kbar. Contrib. Mineral. Petrol. 1996, 123, 92–108. [Google Scholar] [CrossRef]
  58. Blundy, J.; Melekhova, E.; Ziberna, L.; Humphreys, M.C.S.; Cerantola, V.; Brooker, R.A.; McCammon, C.A.; Pichavant, M.; Ulmer, P. Effect of redox on Fe-Mg-Mn exchange between olivine and melt and an oxybarometer for basalts. Contrib. Mineral. Petrol. 2020, 175, 103. [Google Scholar] [CrossRef]
  59. Andersen, D.J.; Lindsley, D.H. Internally consistent solution models for Fe-Mg-Mn-Ti spinels: Fe-Ti oxides. Amer. Mineral. 1988, 73, 714–726. [Google Scholar]
  60. Putirka, K.D. Rates and styles of planetary cooling on Earth, Moon, Mars, and Vesta, using new models for oxygen fugacity, ferric-ferrous ratios, olivine-liquid Fe-Mg exchange, and mantle potential temperature. Amer. Mineral. 2016, 101, 819–840. [Google Scholar] [CrossRef]
  61. Kinzler, R.J.; Grove, T.L. Primary magmas of mid-ocean ridge basalts.1. Experiments and methods. J. Geophys. Res. 1992, 97, 6885–6906. [Google Scholar] [CrossRef]
  62. Yang, H.J.; Kinzler, R.J.; Grove, T.L. Experiments and models of anhydrous, basaltic olivine-plagioclase-augite saturated melts from 0.001 to 10 kbar. Contrib. Mineral. Petrol. 1996, 124, 1–18. [Google Scholar] [CrossRef]
  63. Tormey, D.R.; Grove, T.L.; Bryan, W.B. Experimental petrology of normal MORB near the Kane Fracture Zone: 22–25° N, mid-Atlantic Ridge. Contrib. Mineral. Petrol. 1987, 96, 121–139. [Google Scholar] [CrossRef]
  64. Grove, T.L. Corrections to expressions for calculating mineral components in “Origin of calc-alkaline series lavas at Medicine Lake volcano by fractionation, assimilation and mixing” and “Experimental petrology of normal MORB near the Kane Fracture Zone: 22°–25° N, mid-Atlantic ridge”. Contrib. Mineral. Petrol. 1993, 114, 422–424. [Google Scholar]
  65. Grove, T.L.; Kinzler, R.; Bryan, W. Fractionation of mid-ocean ridge basalt (MORB). In Mantle Flow and Melt Generation at Mid-Ocean Ridges; Phipps Morgan, J., Blackman, D., Sinton, J., Eds.; Geophysical Monograph 71; American Geophysical Union: Washington, DC, USA, 1992; pp. 281–310. [Google Scholar]
  66. Till, C.B.; Grove, T.L.; Krawczynski, M.J. A melting model for variably metasomatized plagioclase and spinel lherzolite. J. Geophys. Res. 2012, 117, B06206. [Google Scholar] [CrossRef]
  67. Chatterjee, N.; Sheth, H. Origin of the Powai ankaramite, and the composition, P-T conditions of equilibration and evolution of the primary magmas of the Deccan tholeiites. Contrib. Mineral. Petrol. 2015, 169, 32. [Google Scholar] [CrossRef]
  68. Till, C.B. A review and update of mantle thermobarometry for primitive arc magmas. Amer. Mineral. 2017, 102, 931–947. [Google Scholar]
  69. Krein, S.B.; Molitor, Z.J.; Grove, T.L. ReversePetrogen: A Multiphase dry reverse fractional crystallization-mantle melting thermobarometer applied to 13,589 mid-ocean ridge basalt glasses. J. Geophys. Res. Solid Earth 2021, 126, e2020JB021292. [Google Scholar]
  70. Roeder, P.L.; Emslie, R.F. Olivine liquid equilibrium. Contrib. Mineral. Petrol. 1970, 29, 275–289. [Google Scholar] [CrossRef]
  71. Thompson, R.N.; Gibson, S.A. Subcontinental mantle plumes, hotspots and preexisting thinspots. J. Geol. Soc. Lond. 1991, 148, 973–977. [Google Scholar] [CrossRef]
  72. Ebinger, C.J.; Sleep, N.H. Cenozoic magmatism throughout East Africa resulting from impact of a single plume. Nature 1998, 395, 788–791. [Google Scholar] [CrossRef]
  73. Duggen, S.; Hoernle, K.A.; Hauff, F.; Kluegel, A.; Bouabdellah, M.; Thirlwall, M.F. Flow of Canary mantle plume material through a subcontinental lithospheric corridor beneath Africa to the Mediterranean. Geology 2009, 37, 283–286. [Google Scholar] [CrossRef]
  74. Begg, G.C.; Hronsky, J.A.M.; Arndt, N.T.; Griffin, W.L.; O’Reilly, S.Y.; Hayward, N. Lithospheric, cratonic, and geodynamic setting of Ni–Cu–PGE sulfide deposits. Econ. Geol. 2010, 105, 1057–1070. [Google Scholar] [CrossRef]
  75. Begg, G.C.; Hronsky, J.M.A.; Griffin, W.L.; O’Reilly, S.Y. Global- to deposit-scale controls on orthomagmatic Ni-Cu(-PGE) and PGE reef ore formation. In Processes and Ore Deposits of Ultramafic-Mafic Magmas through Space and Time; Mondal, S., Griffin, W., Eds.; Elsevier: Amsterdam, The Netherlands, 2017. [Google Scholar]
  76. Royer, J.-Y.; Coffin, M.F. Jurassic to Eocene plate tectonic reconstructions in the Kerguelen Plateau region. Proc. Ocean Drill. Prog. Sci. Res. 1992, 120, 917–928. [Google Scholar]
  77. Marshall, J.F.; Lee, C.S. Basin framework and resource potential of the Abrolhos sub-basin. In Geophysical Yearbook 1988–1989; Wolf, K.H., Paine, A.G.L., Eds.; Australian Bureau of Mineral Resources: Canberra, Australia, 1989; pp. 63–67. [Google Scholar]
Figure 1. Location of the different parts of Kerguelen LIP. The continental landmasses are shown in black and the submarine plateaus and ridges are shown in grey.
Figure 1. Location of the different parts of Kerguelen LIP. The continental landmasses are shown in black and the submarine plateaus and ridges are shown in grey.
Minerals 13 00426 g001
Figure 2. (A) Location of the Rajmahal Traps, Sylhet Traps, and Raniganj–Giridih–Koderma (RGK) dikes on the Indian Shield. CGC—Chotanagpur Gneiss Complex, EGB—Eastern Ghats belt, EITZ—Eastern Indian Tectonic Zone. (B) Detailed geological map of the Rajmahal Traps (after [32]). The Bouguer gravity contours (−60 to +25 mGal) are after [36].
Figure 2. (A) Location of the Rajmahal Traps, Sylhet Traps, and Raniganj–Giridih–Koderma (RGK) dikes on the Indian Shield. CGC—Chotanagpur Gneiss Complex, EGB—Eastern Ghats belt, EITZ—Eastern Indian Tectonic Zone. (B) Detailed geological map of the Rajmahal Traps (after [32]). The Bouguer gravity contours (−60 to +25 mGal) are after [36].
Minerals 13 00426 g002
Figure 3. Bulk chemical compositions of the Rajmahal basalts and basaltic andesites (data from [10,21]). (a) Total alkali versus SiO2 after [44]. (b) Ti/Zr versus Zr/Y. (c) Primitive mantle-normalized [45] incompatible elements.
Figure 3. Bulk chemical compositions of the Rajmahal basalts and basaltic andesites (data from [10,21]). (a) Total alkali versus SiO2 after [44]. (b) Ti/Zr versus Zr/Y. (c) Primitive mantle-normalized [45] incompatible elements.
Minerals 13 00426 g003
Figure 4. Cross-polarized light images of (a) Maharajpur basalt showing plagioclase (Pl) and clinopyroxene (Cpx) phenocrysts, and (b) Tinpahar basalt showing Pl and olivine (Ol) phenocrysts, and backscattered electron images of (c) Taljhari basalt showing Pl and Cpx phenocrysts, and (d) Kherwa basalt showing Pl and Cpx microphenocrysts surrounded by a fine-grained groundmass comprising Cpx, Pl, ilmenite (Ilm), and hematite (Hem).
Figure 4. Cross-polarized light images of (a) Maharajpur basalt showing plagioclase (Pl) and clinopyroxene (Cpx) phenocrysts, and (b) Tinpahar basalt showing Pl and olivine (Ol) phenocrysts, and backscattered electron images of (c) Taljhari basalt showing Pl and Cpx phenocrysts, and (d) Kherwa basalt showing Pl and Cpx microphenocrysts surrounded by a fine-grained groundmass comprising Cpx, Pl, ilmenite (Ilm), and hematite (Hem).
Minerals 13 00426 g004
Figure 5. Chemical composition of minerals in the Rajmahal basalts from the northeast, northwest, and central sectors (see Figures S1–S3 for further details). (a,b) Average augite (En—enstatite, Di—diopside, Hd—hedenbergite, Wo—wollastonite, Quad—quadrilateral components, Jd—jadeite, and Ae—aegirine) plotted according to [53]. (c) Plagioclase (Ab—albite, Olg—oligoclase, Ane—andesine, Lab—labradorite, Byt—bytownite, An—anorthite, Or—orthoclase). (d) Oxides.
Figure 5. Chemical composition of minerals in the Rajmahal basalts from the northeast, northwest, and central sectors (see Figures S1–S3 for further details). (a,b) Average augite (En—enstatite, Di—diopside, Hd—hedenbergite, Wo—wollastonite, Quad—quadrilateral components, Jd—jadeite, and Ae—aegirine) plotted according to [53]. (c) Plagioclase (Ab—albite, Olg—oligoclase, Ane—andesine, Lab—labradorite, Byt—bytownite, An—anorthite, Or—orthoclase). (d) Oxides.
Minerals 13 00426 g005
Figure 6. Portions of the pseudoternary projections of the basalt tetrahedron: (a) Ol-Pl-Cpx from Qz, and (b) Ol-Cpx-Qz from Pl (after [63,64]) showing the compositions of the Rajmahal Group I basalts and their corresponding melts (determined by subtracting plagioclase). For clarity, only the multiple saturation points for plagioclase lherzolite (PL-MSP at 1 bar and 2–10 kbar with 2 kbar intervals, after [59]) of the average basalt and spinel lherzolite (SL-MSP at 9–15 kbar with 3 kbar intervals, after [63]) of the average primary magma are shown. The solid lines with arrowheads represent examples of model reverse fractionation paths along the Ol + Pl + Cpx and Ol + Pl cotectic boundaries, and the dashed lines with arrowheads in (b) represent part of the inferred Ol + Pl + Cpx cotectic boundary at 1 bar and 2 kbar.
Figure 6. Portions of the pseudoternary projections of the basalt tetrahedron: (a) Ol-Pl-Cpx from Qz, and (b) Ol-Cpx-Qz from Pl (after [63,64]) showing the compositions of the Rajmahal Group I basalts and their corresponding melts (determined by subtracting plagioclase). For clarity, only the multiple saturation points for plagioclase lherzolite (PL-MSP at 1 bar and 2–10 kbar with 2 kbar intervals, after [59]) of the average basalt and spinel lherzolite (SL-MSP at 9–15 kbar with 3 kbar intervals, after [63]) of the average primary magma are shown. The solid lines with arrowheads represent examples of model reverse fractionation paths along the Ol + Pl + Cpx and Ol + Pl cotectic boundaries, and the dashed lines with arrowheads in (b) represent part of the inferred Ol + Pl + Cpx cotectic boundary at 1 bar and 2 kbar.
Minerals 13 00426 g006
Figure 7. Bivariate plots showing the variation of the major oxides (ae) and Ni with MgO (f) for the Rajmahal Group I basalts (corrected melt compositions) and their primary magmas. The compositions of melt in which clinopyroxene first crystallizes (end of stage 1) are also plotted. The direction of reverse fractionation with increasing MgO is shown by arrows at the bottom.
Figure 7. Bivariate plots showing the variation of the major oxides (ae) and Ni with MgO (f) for the Rajmahal Group I basalts (corrected melt compositions) and their primary magmas. The compositions of melt in which clinopyroxene first crystallizes (end of stage 1) are also plotted. The direction of reverse fractionation with increasing MgO is shown by arrows at the bottom.
Minerals 13 00426 g007
Table 1. Mineral assemblages in the Rajmahal basalts.
Table 1. Mineral assemblages in the Rajmahal basalts.
RegionLocation AugPgtPlHemIlmOther
Northeast
SahebganjRB88-12mphmphphen gm
Moti JharnaRB88-19phen phen gm
MaharajpurRB88-16phenmphphengmgm
RB88-17phen phengmgm
RB88-18phen phengmgm
TaljhariEB89-154phen phen
RB88-24phen phengmgm
TinpaharRB88-31phen phen gmMgh
RB88-32phen phengmgm
RB88-39mph phen Ol
RB88-35mph phengmgmCumm
KherwaEB89-112mphmphphengm Zeol
EB89-117mphmphphengmgmZeol
EB89-118mphmphphengm
BerhaitEB89-121phenmphphengmgm
EB89-123mph phen gm
EB89-127mph mphgm
EB89-128mphmphphengmgm
Northwest
DariyachakDAR-2GBphen phen gmMgh
DAR-5YLphenmphphengmgm
DAR10-90phen phengmgm
Central
Litipara-Amrapara Rd.LA1phen phengmgm
LA2phen phengmgm
LA4phen phengmgm
LA5phen phengmgm
LA6phen phengmgm
Litipara-Simlong Rd.TT1phen phengmgm
TT2phen phengm
TT3phenmphphengmgm
TT4phen phengmgm
Litipara-Hiranpur Rd.LH1phen phengmgm
AmraparaBL1phen phengmgmKfs, Zeol
BL3phen phengmgmKfs, Zeol
KundpaharKP1phen phengmgm
KP2mph phengmgmKfs, Zeol
KP3mph phengmgm
KP4mph phen
KP5phen phengmgm
KP6phen phengmgm
KP7phen phengmgm
Aug—augite, Pgt—pigeonite, Pl—plagioclase, Hem—hematite, Ilm—ilmenite, Mgh—maghemite, Ol—olivine, Cumm—cummingtonite, Kfs—K-feldspar, Zeol—zeolite, phen—phenocryst, mph—microphenocryst, gm—groundmass.
Table 2. Composition range of minerals in the Rajmahal basalts.
Table 2. Composition range of minerals in the Rajmahal basalts.
A. Augite and Plagioclase
SectorNortheastNorthwestCentral
corerimcorerimcorerim
Augite
En45–5649–5248–5449–5047–5449–54
Fs8–2110–159–1211–157–157–19
Wo29–3934–3936–4037–4036–4131–41
Mg#68–8776–8380–8577–8276–8973–88
Plagioclase
An57–8258–7561–7145–6660–8151–77
Ab18–4225–4128–3834–5419–3922–48
Or0–40–30–10–10–10–1
B. Other minerals
All sectors
IlmeniteHematite Pigeonite Olivine
Hem0–743–92En42–74Fo71–76
Ilm85–976–53Fs18–46Fa24–30
Pph1–10–12Wo4–12
Gk1–90–6Mg#48–81
En—enstatite, Fs—ferrosilite, Wo—wollastonite, Mg# = 100.molar Mg/(Mg + Fe), An—anorthite, Ab—albite, Or—orthoclase, Hem—hematite, Ilm—ilmenite, Pph—pyrophanite, Gk—geikelite, Fo—forsterite, Fa—fayalite.
Table 3. Clinopyroxene thermobarometry of the Rajmahal basalts.
Table 3. Clinopyroxene thermobarometry of the Rajmahal basalts.
LocationSamplePTPTPTPT
kbar°Ckbar°Ckbar°Ckbar°C
AvgAvg
Equations (32a) and (32d) [54][55][56]
Northeast sector
SahibganjRB88-120.311560.311562.0 ± 1.01138 ± 440.0 ± 0.01139 ± 21
Moti JharnaRB88-193.211803.211804.6 ± 2.01151 ± 190.0 ± 2.31157 ± 25
MaharajpurRB88-160.00111191.611452.0 ± 0.71133 ± 460.0 ± 0.11125 ± 23
RB88-183.21172
TaljhariEB89-1540.00111470.00111472.0 ± 0.51137 ± 520.0 ± 0.21141 ± 23
RB88-240.0011147
TinpaharRB88-314.612095.411954.5 ± 2.21127 ± 431.8 ± 4.51147 ± 33
RB88-324.41196
RB88-397.51218
RB88-355.01156
KherwaEB89-1120.00111111.311363.7 ± 2.11099 ± 620.0 ± 0.71119 ± 43
EB89-1172.71160
BerhaitEB89-1216.612131.811603.3 ± 2.61148 ± 320.4 ± 1.01134 ± 34
EB89-1230.0011127
EB89-1270.11147
EB89-1280.51152
Northwest sector
DariachakDAR-2GB2.011700.711342.0 ± 0.81093 ± 420.1 ± 0.31137 ± 28
DAR-5YL0.0011100
DAR10-900.0011131
Central sector
Litipara-AmraparaLA12.211601.511602.0 ± 0.61137 ± 500.0 ± 0.61139 ± 29
LA23.01177
LA40.81155
LA50.61156
LA60.81154
Litipara-SimlongTT12.811851.111633.1 ± 1.91141 ± 320.0 ± 0.51149 ± 25
TT20.041162
TT31.71166
TT40.0011138
Litipara-HiranpurLH10.111550.111554.1 ± 2.61150 ± 270.0 ± 1.01149 ± 28
AmraparaBL12.311881.911823.8 ± 2.41096 ± 510.0 ± 1.21152 ± 48
BL31.61175
KundpaharKP13.611921.911622.3 ± 1.11123 ± 520.0 ± 0.51136 ± 27
KP23.31190
KP32.01157
KP40.71156
KP52.51161
KP60.0011125
KP71.31153
Table 4. Clinopyroxene-bulk and whole rock thermobarometry of the Rajmahal basalts.
Table 4. Clinopyroxene-bulk and whole rock thermobarometry of the Rajmahal basalts.
LocationSampleP (kbar)T (°C)KD aPl (%) b
Equations (P1) and (T1) [57]:
Central sector
Litipara-AmraparaLA24.611850.29
Litipara-SimlongTT31.611660.32
Litipara-HiranpurLH10.00111370.31
[61], and Equation (16) [54]:
Northeast sector
Moti Jharna88-2141188 4
Tinpahar88-3021168 7
88-423.51183 4
Northwest sector
S of Dariachak cRJ1-25-111168 6
S of Dariachak dRJ1-26-70.0011147 4
Central sector
Litipara-SimlongTT321169 5
Litipara-HiranpurLH10.0011158 10.5
W of Pakur eRJ1-30-30.0011158 3.5
RJ1-30-40.0011156 3.5
KundpaharKP60.0011157 10
a Fe2+-Mg exchange coefficient for Cpx-bulk, b amount of plagioclase subtracted from bulk to determine melt composition, c Lalmatia, d Bejam Pahar, e Dhanbad village.
Table 5. P-T of crystallization and last equilibration of primary magma with mantle for the Rajmahal Group I basalts.
Table 5. P-T of crystallization and last equilibration of primary magma with mantle for the Rajmahal Group I basalts.
CrystallizationStage 1 aStage 2 aPrimary Magma b
P (kbar)T (°C)F%OlPlCpxF%OlPlP (kbar)T (°C)
Northeast sector
88-214118860.713.752.433.912.228.671.491286
88-302116860.313.052.634.419.830.269.891284
88-423.5118357.813.053.433.614.830.269.891284
Northwest sector
RJ1-25-11116846.812.253.134.722.930.269.891280
RJ1-26-70.001114757.810.952.636.430.330.269.891280
Central sector
TT32116953.812.553.633.919.829.470.681268
LH10.001115854.213.051.935.226.028.671.48.51275
RJ1-30-30.001115855.612.150.237.729.630.269.891281
RJ1-30-40.001115656.510.952.636.429.630.269.891281
KP60.001115754.211.851.636.628.130.269.891282
a Percent fractionation and proportions of phases added. b Primary magmas (Mg# 73–76) are equilibrated with olivine Fo90–91.5, Cpx Mg# 91.5–93, and Opx Mg# 91–92.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chatterjee, N.; Ghose, N.C. Thermobarometry of the Rajmahal Continental Flood Basalts and Their Primary Magmas: Implications for the Magmatic Plumbing System. Minerals 2023, 13, 426. https://doi.org/10.3390/min13030426

AMA Style

Chatterjee N, Ghose NC. Thermobarometry of the Rajmahal Continental Flood Basalts and Their Primary Magmas: Implications for the Magmatic Plumbing System. Minerals. 2023; 13(3):426. https://doi.org/10.3390/min13030426

Chicago/Turabian Style

Chatterjee, Nilanjan, and Naresh C. Ghose. 2023. "Thermobarometry of the Rajmahal Continental Flood Basalts and Their Primary Magmas: Implications for the Magmatic Plumbing System" Minerals 13, no. 3: 426. https://doi.org/10.3390/min13030426

APA Style

Chatterjee, N., & Ghose, N. C. (2023). Thermobarometry of the Rajmahal Continental Flood Basalts and Their Primary Magmas: Implications for the Magmatic Plumbing System. Minerals, 13(3), 426. https://doi.org/10.3390/min13030426

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop