1. Introduction
Shape fabric analysis is an important aspect to unravel the complex geological structures and determine the state of strain in rocks (Woodcock, Reference Woodcock1977; Cobbold & Gapais, Reference Cobbold and Gapais1979). Various methods have been used to investigate shape fabrics, such as stereographic projection, eigenvalue calculations, and anisotropy of magnetic susceptibility (AMS) measurements (e.g. Woodcock, Reference Woodcock1977; Cobbold & Gapais, Reference Cobbold and Gapais1979; Launeau et al. Reference Launeau, Bouchez and Benn1990; Simpson & De Paor, Reference Simpson and De Paor1993; Tikoff & Fossen, Reference Tikoff and Fossen1995; Wallis, Reference Wallis1995; Launeau & Robin, Reference Launeau and Robin1996; Olivier et al. Reference Olivier, Saint Blanquat, Gleizes, Leblanc, Bouchez, Hutton and Stephen1997; Grasemann et al. Reference Grasemann, Fritz and Vannay1999; Gomez-Rivas et al. Reference Gomez-Rivas, Bons, Griera, Carreras, Druguet and Evans2007). Amongst these methods, the eigenvalue method is a powerful tool which provides detailed information about the shape as well as the strength of the fabrics. It also helps to quantify the degree of randomness, both arithmetically and graphically (Williams & Chapman, Reference Williams and Chapman1979; Woodcock & Naylor, Reference Woodcock and Naylor1983).
AMS is often measured to quantify the fabric of low-anisotropy rocks, i.e. rocks which do not show any mesoscopic field foliation such as granite, quartzite etc. (Owens & Bamford, Reference Owens and Bamford1976; Tarling & Hrouda, Reference Tarling and Hrouda1993; Borradaile & Henry, Reference Borradaile and Henry1997; Bouchez Reference Bouchez, Bouchez, Hutton and Stephen1997; Borradaile & Jackson, Reference Borradaile, Jackson and Martín-Hernández2004; Mondal & Mamtani, Reference Mondal and Mamtani2013). AMS helps to understand the deformation and strain variation in a rock and establish the tectonic history of a region (Mamtani & Greiling, Reference Mamtani and Greiling2005; Almqvist et al. Reference Almqvist, Henry, Jackson, Werner and Lagroix2014; Ferré et al. Reference Ferré, Gébelin, Till, Sassier and Burmeister2014; Mondal & Mamtani, Reference Mondal and Mamtani2014 and references therein). AMS involves inducing a magnetic field in a sample in different directions and measurement of the induced magnetization in each direction (Tarling & Hrouda, Reference Tarling and Hrouda1993). The orientation and magnitude of the three principal axes of the AMS ellipsoid allow the magnetic planar and/or linear fabric to be determined.
Although AMS provides information about the magnetic fabrics, it is challenging to relate it to deformation in case of poly-deformed rocks. In superposed deformation, the early formed magnetic foliation in rocks may be reoriented and become very weak. Thus, it is crucial to quantify these fabrics in light of regional deformation. The eigenvalue method helps to analyse such fabrics using few parameters. Collectively, most of the previous studies performed on the geological application of AMS methods show a lack of clear assessment of the correlation between eigenvector analysis of AMS data and its connection with the mechanistic issues related to regional tectonics. Of a few studies that have addressed such topic, Parés et al. (Reference Parés, van der Pluijm and Dinarès-Turell1999) used eigenvector analysis of AMS data to link magnetic fabric development with the progressive deformation in mud rocks. In this regard, Pueyo et al. (Reference Pueyo, Román-Berdiel, Bouchez, Casas and Larrasoaña2004) also discussed the utilization of normalized eigenvalues in order to analyse magnetic fabric statistics of paramagnetic granites.
In the present paper, we use published AMS data from five different granitoids of India and analyse them using the eigenvalue method to quantify their degree of relative randomness. The present investigation uses the following steps: (1) calculation of eigenvalues and corresponding eigenvectors from the ‘orientation tensor matrix’; and (2) determination of shape and strength parameters from the magnetic foliation data. The study graphically quantifies the degree of randomness of poles to magnetic foliation.
2. Shape fabric analysis using eigenvalue method – the background
In geology, the characterization of directional data along with their relative attitude, randomness and shape are crucial. For a large-array dataset, overlap between data is common and their representation becomes challenging. Therefore, an approach consisting in specifying directional features by a limited number of representative parameters should avoid this large-array problem. Shape fabric analysis is very useful method in representing such a dataset. Attributing these representative parameters can be solved using a procedure that involves a tensor matrix determined from the dataset, and from which eigenvectors and eigenvalues may be calculated (Watson, Reference Watson1966; Mark, Reference Mark1973, Reference Mark1974; Owens, Reference Owens1973; Mark & Andrews, Reference Mark and Andrews1975; Woodcock & Naylor, Reference Woodcock and Naylor1983).
Before diving into main purpose of the current study, we would like to briefly refer to some well-known previously proposed methods of shape fabric analysis. In order to depict the preferred orientation of tectonically strained objects, Flinn (Reference Flinn1965) established a diagram, having the ‘maximum/intermediate’ and ‘intermediate/minimum’ as its vertical and horizontal axes, respectively, which has been of immense interest to structural geologists till today. Subsequently, Woodcock (Reference Woodcock1977) introduced another spherical data analysis method, termed as the ‘Flinn–Woodcock plot’, by integrating the shape of ellipsoid with the distribution type (i.e. cluster/great-circle girdle) of the orientational data. In addition to these, other interesting methods of analysing the fabric shapes were proposed by Jelínek (Reference Jelínek1981), which not only described the type of magnetic anisotropy, but also addressed some important issues associated with the anisotropic concentration of orientational data. Borradaile (Reference Borradaile2003) provided an in-depth review of several methods of spherical-orientation data analysis, relating to sampled AMS data distribution in time, space and orientation.
However, it may be noted that these studies were all restricted to the fabric shapes originating from the single deformation event and provided scope to further study the more complex, multi-modal data distributions, which are the product of some specific regional tectonics.
In this present study, we perform the eigenvalue analysis method and analyse magnetic fabrics of different syntectonic granitic plutons from different parts of India to comment on the variation of magnetic shape fabric from a coaxially deformed region to a non-coaxially deformed one. As Scheidegger (Reference Scheidegger1965) and Woodcock (Reference Woodcock1977) provide a detailed description of the theoretical basis of this eigenvalue method, we present here only the salient aspects of the principle involved in shape fabric analysis. The eigenvalue method assumes that each orientation is represented by a unit vector. Therefore, a 3 × 3 matrix, called an ‘orientation tensor matrix’, is formed by taking the summation of the cross-products of the direction cosines of unit vectors (Scheidegger, Reference Scheidegger1965).
$$a = \left\lfloor {\matrix{
{\sum {{l^2}} } \hfill & {\sum {lm} } \hfill & {\sum {ln} } \hfill\cr
{\sum {ml} } \hfill & {\sum {{m^2}} } \hfill & {\sum {mn} } \hfill\cr
{\sum {nl} } \hfill & {\sum {nm} } \hfill & {\sum {{n^2}} } \hfill\cr} } \right\rfloor $$
A normalized form of this matrix (A) is given below:
$$A = a/N.$$
where, N is the number of data; l, m and n are the direction cosines of a particular unit vector. Further, the eigenvectors (v1, v2 and v3) and their corresponding eigenvalues (λ1, λ2 and λ3) are determined from the above tensor matrix. The maximum eigenvector v1 represents the direction along the minimum ‘moment of inertia’ of the distribution, while v3 indicates maximum ‘moment of inertia’ (Watson, Reference Watson1966). The eigenvalues are used mainly in their normalized form:
$${{\rm{S}}_1} + {{\rm{S}}_2} + {{\rm{S}}_3} = 1,\,{\rm{where}}\,{{\rm{S}}_i} = {\lambda _i}/{\rm{N}}.$$
These normalized eigenvalues are directly related to the shape of the fabric, as the clusters and girdles in equal-area projection of a distribution consist of S1 > S2 ∼ S3 and S1 ∼ S2 > S3, respectively (Watson, Reference Watson1966).
The last step of this analysis involves the calculation of two parameters from the graphical presentations of S1, S2 and S3: (i) The shape parameter (K), which can be obtained by plotting the ln(S2/S3) vs ln(S1/S2), similar to the Flinn diagram (Flinn, Reference Flinn1962); and (ii) the strength parameter (C) defined by ln(S1/S3) which quantifies the randomness of the distribution.
In the present paper, we use the above two parameters, and position of the eigenvectors on the equal-area projection of the poles to the magnetic foliations corresponding to the primary dataset coming from different granitoids of India. The randomness is quantified, and this is followed by an assessment of the shape of the distribution data in the light of their regional tectonics. Figure 1 shows the stepwise procedure of the above-described calculations applied to the poles to magnetic foliations (K3) represented by their plunge (P) and trend (T) coming from the Chakradharpur granite (India).
Stepwise procedure for eigenvector analysis (Woodcock & Naylor, Reference Woodcock and Naylor1983) from pole to magnetic foliation (K3) data (Mamtani et al. Reference Mamtani, Pal and Greiling2013).
3. Available data
As already mentioned, we have used the AMS data from five granite bodies, namely Chakradharpur, Godhra, Malanjkhand, J. N. Kote and Chitradurga, located in various parts of India (see Fig. 2). Extensive AMS studies have been carried out in the above regions. All the granites have been inferred to be syntectonic with regional deformation and associated with regional-scale shear zones in their respective vicinities. Below we provide available data which are only salient for the present analysis.
Location map of the various younger granites from Indian continent. (a–d) Regional maps of Chakradharpur (Mamtani et al. Reference Mamtani, Pal and Greiling2013) (a), Malanjkhand (Majumder & Mamtani, Reference Majumder and Mamtani2009) (b), Godhra (after Sen & Mamtani, Reference Sen and Mamtani2006) (c) and J. N. Kote (Mondal & Mamtani, Reference Mondal and Mamtani2014) (d) granitoids. (e, f) Map of Chitradurga granites (Mondal, Reference Mondal2018). Dashed boxes in (e) and (f) are the northern and southern part of the granite, respectively. Inset in (a) shows the locations of these granitoids in the map of India. The lower-hemisphere equal-area projection in each panel shows the pole (K3) to magnetic foliation and their corresponding symmetry arguments. Planes of symmetry (m) are shown as dashed red lines in each of the stereonets. Software StereoNet (Allmendinger et al. Reference Allmendinger, Cardozo and Fisher2013; Cardozo and Allmendinger, Reference Cardozo and Allmendinger2013) was used for all lower-hemisphere equal-area projection and contouring (http://www.geo.cornell.edu/geology/faculty/RWA/programs/stereonet.html). Dashed black lines in (a), (b), (d), (e) and (f) show the trend of regional shear zone.
3.a. Chakradharpur granitoid
The magnetic fabric in this Precambrian granitoid has developed syntectonically along with the evolution of the Singhbhum Shear Zone (SSZ) that lies to its south (Mamtani et al. Reference Mamtani, Pal and Greiling2013). The Chakradharpur granitoid lacks well-developed magmatic fabrics, and the mean magnetic susceptibility (Km) varies from 67.6 to 659 µSI (Table 1). It has been inferred that biotite is the main paramagnetic phase contributing to its AMS (Mamtani et al. Reference Mamtani, Pal and Greiling2013). The mean magnetic foliation plane (K1K2) has an orientation of N54°E (strike) with a steep to vertical dip (Fig. 2a). Mamtani (Reference Mamtani2014) has performed a 2D vorticity analysis and concluded that the fabric of this granite body is dominated by pure shear (kinematic vorticity number Wk = 0.58) and associated with the evolution of nearby shear zone.
Characterization of every granitic body
Km = mean susceptibility; σK = standard deviation of magnetic susceptibility; K1 = magnetic lineation; K3 = pole to magnetic foliation – shaded regions for each studied granite depict % of sample sites, having ferromagnetic character. References for available data: Chakradharpur (Mamtani et al. Reference Mamtani, Pal and Greiling2013); Malanjkhand (Majumder & Mamtani, Reference Majumder and Mamtani2009); Godhra (Sen & Mamtani, Reference Sen and Mamtani2006); J. N. Kote (Mondal & Mamtani, Reference Mondal and Mamtani2014); Chitradurga (South) and Chitradurga (North) (Mondal, Reference Mondal2018).
3.b. Malanjkhand granitoid
AMS measurements have been performed in the ∼2.48 Ga Malanjkhand body (Fig. 2b) in order to evaluate the time relationship between fabric development in the granite and the regional tectonics (Majumder & Mamtani, Reference Majumder and Mamtani2009). The Central Indian Suture (CIS) that defines the southern margin of the Central Indian Tectonic Zone (CITZ) demarcates the NW boundary of this granite body. The magnetic fabric trajectory that can be drawn in the Malanjkhand granite is interpreted to be related to synmagmatic deformation (Majumder & Mamtani, Reference Majumder and Mamtani2009). Mamtani (Reference Mamtani2014), has previously identified two prominent sectors (domain-I and domain-II) in Malanjkhand granite based on their spatial disposition with respect to the adjacent CITZ. The Km shows a variation from 190.8 to 3490.5 µSI in domain-I. However, for domain-II, this range is found from 268.47 to 5789.62 µSI (Table 1). It is also interpreted that domain-I and domain-II are pure- (Wk = 0.98) and simple-shear (Wk = 0.34) dominated, respectively.
3.c. Godhra granitoid
The Godhra body (955 ± 20 Ma) is known to have developed its fabric synchronously with Grenvillian-age tectonic rejuvenation of the Central Indian Tectonic Zone (CITZ) that lies to its south (Fig. 2c). AMS studies have revealed that the paramagnetic minerals, namely biotite and in some samples hornblende, and the ferromagnetic (sensu lato) magnetite are the important phases that contribute to the AMS (Mamtani & Greiling, Reference Mamtani and Greiling2005). The magnetic foliation in the Godhra granite body is dominantly ENE–WSW in orientation (Fig. 2b), i.e. parallel to the CITZ (Sen & Mamtani, Reference Sen and Mamtani2006). Km values vary widely from 169.07 to 5756 µSI. Recently, Mamtani (Reference Mamtani2014) concluded that the late-stage fabric in this magmatic body was dominated by pure shear (Wk = 0).
3.d. J. N. Kote granitoid
This Archaean-age granite body (Fig. 2d) is located in the vicinity of the Chitradurga granite (Fig. 2e) and was emplaced syntectonically during the evolution of the adjacent Chitradurga Shear Zone (CSZ). AMS was performed on this granite (Mondal, Reference Mondal2018) and the magnetic fabric is found to be parallel to the CSZ. Their susceptibility (Km), being ∼69.85 µSI, indicates a dominant contribution of paramagnetic phases. The magnetic fabric of the J. N. Kote body is interpreted to be dominantly due to simple shearing (Wk = 0.80) that occurred during its emplacement along the granite–TTG contact (Mondal, Reference Mondal2018).
3.e. Chitradurga granitoid
The fabric in granite (∼2.6 Ga) from the Chitradurga region (Western Dharwar Craton, south India; Fig. 2e, f) is analysed using AMS study (Mondal, Reference Mondal2018; Mondal et al. Reference Mondal, Bhowmick, Das and Patsa2020). The microstructural investigation on the granite shows a progressive textural overprint from magmatic, through high-T to low-T solid-state deformation textures. The mean magnetic foliation in the rocks of the region is dominantly NW–SE-striking. Km values range from 138.42 to 1289.06 µSI in the southern region, while in the northern region they vary from 116.69 to 585.14 µSI. The vorticity analysis from magnetic fabric in the southern region of the Chitradurga granite reveals that the NW–SE-oriented fabric formed under pure shear condition (Wk = 0.06; Mondal, Reference Mondal2018; dashed box in Fig. 2e). However, the northern region of the granite is closed to the adjacent CSZ and is inferred to be controlled by simple shearing (Wk = 88; dashed box in Fig. 2f).
Based on the available magnetic dataset, each granitic body is characterized in Table 1 (see table caption for references) by its emplacement as well as deformation age(s), number of AMS stations, mean bulk susceptibility (Km and the information on standard deviation), and declination and inclination of mean K1 and K3. This characterization, inturn, embodies the basis for two initial important assumptions made for this study which are discussed in Section 4.
We also use the magnetic foliation data to quantify the degree of randomness of the fabric under this study. The study also helps us to understand the mode of shearing (simple/pure) responsible for the development of the magnetic fabrics in the corresponding areas. Details of this application of magnetic fabric using the eigenvalue method are presented in Section 5.
4. Initial assumptions made in the current analysis
Prior to discussing the application of the present study, we epitomize the twofold aspects of our assumptions in the following subsections.
4.a. Behaviour of granitic bodies under respective regional deformation
The granitoids analysed in this study are replete with structural fabric which has developed during their syntectonic emplacement on account of corresponding regional tectonics. These granitic plutons show superimposition of low-temperature deformation textures over high-temperature and that developed uniformly when the granites cooled from high temperature to low temperature.
In the case of Malanjkhand (Majumder & Mamtani, Reference Majumder and Mamtani2009) and Godhra granite (Mamtani & Greiling, Reference Mamtani and Greiling2005), the magnetic foliations (K3) not only fit well with the field fabric orientation of the two granites, but also show similarities to that of the older gneiss adjacent to their margin. The magnetic lineation (K1) is found to be consistent and uniform even in a large domain (kilometre-scale). Furthermore, the presence (having overall consistency in stretching lineation at km scale) and absence of mylonite, respectively, in domain-I and domain-II of the Malanjkhand granite (shown in Fig. 2b) support the uniform and unique behaviour of each, while responding to the syntectonic cooling that succeeded their emplacement. A discussion on the kinematics of Chakradharpur granite (Mamtani et al. Reference Mamtani, Pal and Greiling2013) reveals that the granite preserves an oblique relationship between its magnetic and field foliation, interpreting the development of the former in the crystallization stage. Besides, the field foliation of the granite has been found parallel to the adjacent shear zone (SSZ) orientation located in its southern margin, indicating that the magnetic and field fabrics of this granite ‘froze’ as succeeding events under the progressive deformation linked to its syntectonic crystallization. According to Mondal (Reference Mondal2018), the other three studied granites i.e. Chitradurga (South), Chitradurga (North) and J. N. Kote, also preserve micro- and mesoscale features that are similar to the Malanjkhand and Godhra granite. This study confirms the existence of parallelism between magnetic foliations in these and a regional trend of the adjacent shear zone (CSZ) as well as field foliations measured from adjacent metasedimentary rocks. Moreover, shape preferred orientation (SPO) analysis, performed by the same author on recrystallized sub-grains of stretched quartz, demonstrates that their long-axis orientation exactly traces the magnetic foliations, which are again found to be in good agreement with the regional stress states associated with these three granites (Mondal, Reference Mondal2018).
Bearing in mind the above-discussed information related to all of our studied granites (or part, as in Malanjkhand pluton), we assume that they behaved in a uniform and homogeneous way in order to record an overall shear type (pure/simple) during their corresponding regional deformation. Indeed, the pluton-scale estimation of different kinematic vorticity numbers (see their values in Section 3) from the granitoids or any domain of it is also found to be in line with our assumption.
In this connection, considering the main characteristics of K3 data on the stereonet, such as well-defined girdles or strong maxima (the latter in our case) at the pluton scale, we observe that there exist three mutually orthogonal planes of symmetry for all the analysed granites (Fig. 2 insets). As suggested by Paterson and Weiss (Reference Paterson and Weiss1961), this feature of the stereonet pattern can be used to interpret the overall orthorhombic symmetry (Bingham, Reference Bingham1974) of structural fabric at larger scale (pluton scale, here), even if any deviation (monoclinic/triclinic/axial, etc.) may exist at the metric or decametric site scale as a result of natural variability, including local response to deformation. Likewise, extrapolation of any small-scale symmetry pattern seems not to be always rational when assessing the structural dataset over a comparatively large scale.
Since no field observations have marked occurrences of significant local stress-induced perturbations in any of the studied plutons, we trust in the mutual orthogonality of three symmetry planes to infer the overall orthorhombicity of the magnetic foliation dataset at the pluton scale for each granite.
4.b. AMS ellipsoid and its relation to the rock fabric
The relationship between AMS ellipsoid and rock fabric (or mineral preferred orientations) has long been an absorbing subject of discussion, and provides the final basis of our assumptions in the present study. Several research (Rochette, Reference Rochette1987; Bouchez, Reference Bouchez, Bouchez, Hutton and Stephen1997), that thematically targeted granitic bodies with regard to this particular issue, clearly revealed that this relationship can be granted for case where paramagnetic minerals make a dominant contribution to the bulk magnetic susceptibility (K) of a granitic body. Although the interpretation of ferromagnetic and paramagnetic data may differ substantially (Trindade et al. Reference Trindade, Raposo, Ernesto and Siqueira1999; Terrinha et al. Reference Terrinha, Pueyo, Aranguren, Kullberg, Kullberg, Casas-Sainz and Azevedo2017). Especially when the contribution of iron-phyllosilicates (such as biotite, amphibole, etc.) defines the K of a granite, their magneto-crystalline anisotropy gives rise to parallelism between their magnetic axes and crystallographic axes (and therefore their shape axis) (Martín-Hernández and Hirt, Reference Martín-Hernández and Hirt2003). However, in comparison with such granites where AMS can directly be utilized as a strain indicator, drawing similar correlations is somewhat complicated for ferromagnetic-dominated (magnetite-bearing) granites. This is due to the functioning of other factors such as a strong control of magnetite content on K (Benn et al. Reference Benn, Rochette, Bouchez and Hattori1993) and the ‘interaction anisotropy’ (Grégoire et al. Reference Grégoire, de Saint-Blanquat, Nédélec and Bouchez1995). In the present study, information on bulk mean-magnetic susceptibility (Km and its standard deviation) for different granitoids is subcategorized in Table 1 according to the percentage of their AMS sites whose samples show either a paramagnetic or ferromagnetic character. Bouchez (Reference Bouchez, Bouchez, Hutton and Stephen1997) stated that granites with Km < 500 µSI can be interpreted as paramagnetic granites whereas the others (Km > 500 µSI) are referred to as ferromagnetic. Based on this, Chakradharpur, Malanjkhand (domain-I), J. N. Kote and Chitradurga (North) granites (column 1 of Table 1) show that samples from ∼90% of these sites are paramagnetic, whereas very few (∼8–12 %) appear to be ferromagnetic, so they behave as overall paramagnetic-dominated granites. However, for the rest of the granites (such as Malanjkhand (domain-II), Godhra and Chitradurga granite (South)), ferromagnetic samples are observed in a substantial percentage (∼58 %) of sites, with paramagnetic ones in <42–44 % of sites. It has also been previously established for all granites that the iron-phyllosilicate such as biotite (and hornblende in a few samples) was the main paramagnetic phase in contributing to the development of magnetic susceptibility, whilst magnetite (sensu lato) acted as an important ferromagnetic contributor to the AMS (see individual references for each granite in Table 1 caption). Taking account of the probable influence of magnetite content on the magnetic fabric, we observe that magnetic anisotropy (Km vs P’; after Jelínek, Reference Jelínek1981) scatterings for ferromagnetic samples do not provide any one-to-one relationships for Malanjkhand (domain-II), Godhra and Chitradurga (South) granite (Supplementary figure 1, available online at https://doi.org/10.1017/S0016756822000747). Furthermore, lower-hemisphere equal-area projections of K3 orientations for all ferromagnetic samples confirm their resemblance with respect to that in paramagnetic samples (Supplementary figure 2, available online athttps://doi.org/10.1017/S0016756822000747). This information therefore allows us to conclude that so far as the orientations of magnetic foliation are considered, the magnetite content of these granites does not remarkably create noise in the contribution of paramagnetic phases to the AMS.
In addition, ore petrography of the paramagnetic samples clearly revealed that the opaque phases present in all the studied granites are goethite, pyrite and hematite, excluding magnetite. These studies also concluded that magnetite (for all the granites) displays ‘Verwey Transition’ (Tarling & Hrouda, Reference Tarling and Hrouda1993) in the cooling experiments, indicating their Multidomain (MD) characteristics. This finally helps us to eliminate the possibility of inverse fabric development in the studied granitoids. Based on the above discussions, we have preferably assumed that the magnetic fabric dataset used in the current study records mineral (or shape) preferred orientation and hence provides the strain information directly for each of the studied plutons.
5. Results
The published magnetic data (poles to magnetic foliation) from Chakradharpur, Malanjkhand (domain-I and -II), Godhra, J. N. Kote, northern and southern Chitradurga granite are analysed using the eigenvalue method. In Fig. 2 (insets), the positional individuality of characteristic eigenvectors, calculated (arithmetically) from orientation tensor matrices (see Fig. 1, third and fourth steps) using the available dataset, is depicted on the contoured lower-hemisphere equal-area projection of poles to magnetic foliation for each granitoid, separately.
It may be noted that all the granitoids show orthorhombic fabric symmetry classes, where the position of the maximum eigenvector (V1), for each granitoid exhibits their tendency to overall trace the densest area of data distribution on stereonets (Fig. 2, inset).
In Fig. 3, the attitude of the plane containing the intermediate and minimum eigenvector (i.e. the V2–V3 plane) is separately portrayed to determine its relative disposition with respect to mean magnetic foliation (MF) and shear zone, for individual granitoid. It is clearly evident from Fig. 3b1 and e that the strike of the V2–V3 planes for the pure-shear-dominated regions, such as Malanjkhand domain-II and southern Chitradurga, exactly mimics the strike of MF planes. Conversely, in the simple-shear-dominated regions such as Malanjkhand domain-I, J. N. Kote and northern Chiradurga, the relationship between V2–V3 and MF does not hold with the aforementioned observation, as a distinct shift of the strike of the V2–V3 planes from that of the MFs may be noted in Figure. 3b2, c and f, respectively. The angular difference between the strikes of V2–V3 and MF in these three granitoids is found to be 10°, 5° and 9°, respectively. Interestingly, however, the other two studied granitoids i.e. Chakradharpur and Godhra (the shaded region in the Fig. 3), exhibit an eccentric relationship between their V2–V3 and MF directions, contrary to our observations both in Figure 3b1 and e and in Figure 3b2, c and f. Despite being a pure-shear-dominated region, the V2–V3 plane of Godhra granite shows a marked angular variation (6°) from its MF plane’s strike. Also, Chakradharpur, which is a general-shear-dominated region, reflects greater angular variation (∼11°) between the strikes of its V2–V3 and MF than what we previously observed in other simple-shear-dominated regions. Figure 3 (inset), where the angular relationships (between V2–V3 and MF) of all granitoids are plotted against their vorticity numbers, also clearly reflects such eccentricities observed in the two granitoids (see the shaded elliptical region in Fig. 3 in-set). For the other granitoids, the angular shift of the V2–V3 plane from MF’s orientation linearly increases with a general increase in their vorticity number, i.e. the stronger the component of simple shear under which a particular region deforms, the greater is the angular shift of the V2–V3 plane from its MF plane’s strike.
Superimposition of V2–V3 planes and best-fit great circles (see legends) on lower-hemisphere equal-area projections of pole to magnetic foliation (K3) for each individual granitoid, depicting their relations with the corresponding mean magnetic foliation planes (MF) and the flow apophyses directions of extension (Ae = shear zone). The positional uniqueness of eigenvectors amid the K3 data distributions is also exhibited in each panel. (a), (b1), (b2) (c), (d), (e) and (f) denote Chakradharpur, Malanjkhand (domain-II), Malanjkhand (domain-I), J.N.Kote, Godhra, Chitradurga southern and northern granite, respectively. The eigen parameters are achieved by writing algorithms in MATLAB interface.
Seeking an explanation of this eccentricity observed in the above-mentioned two granitoids (Chakradharpur and Godhra), we analyse the modality of their data (K3) distributions in Fig. 4. The histogram plots in Fig. 4 reveal that the declinations of K3 appear to have been bimodal for Chakradharpur and Godhra (Fig. 4a, d), while an overall unimodality is prevalent for the other granitoids. In the next section, we elaborate on the role of the varying modality of data distributions to explain the connection between the orientations of V2–V3 and MF for the studied granitoids in light of their different deformation regimes.
(a), (b1), (b2), (c), (d), (e) and (f) are the histograms of declinations of magnetic foliation of Chakradharpur, Malanjkhand domain-II and -I, J. N. Kote, Godhra, Chitradurga southern and northern granite, respectively. (a) and (d) clearly show the bimodal distributions, while the rest remain overall unimodal.
Figure 5 depicts the Cartesian plot (in ln space) of S2/S3 vs S1/S2. The shape parameter (K) is equal to ln(S1/S2)/ln(S2/S3). By analogy to the Flinn diagram, each plot in Figure 5 is divided into two regions by the line K = 1. The region K > 1 depicts the cluster shape of the distribution, while K < 1 denotes the girdle shape. It is envisaged that the K values of all the granitoids show the cluster shape of the distribution, except for the J. N. Kote granite. The K values for Chakradharpur and Chitrdurga are found to be comparable and lie close to a uniaxial clustering. The K value for the Godhra granite tends toward 1 and suggests a feeble clustering of the K3. Figure 6a presents a ternary diagram based on three indices, P (Point or Cluster), G (Girdle) and R (Random), that represent the three end members of fabric distribution patterns obtained from the normalized eigenvalues (S1, S2 and S3). To measure their relative values, we use the following relationships (Vollmer, Reference Vollmer1990):
$${\rm{P}} = \left( {{{\rm{S}}_1} - {{\rm{S}}_2}} \right)$$
$${\rm{G}} = 2\left( {{{\rm{S}}_2} - {{\rm{S}}_3}} \right)$$
$${\rm{R}} = 3{{\rm{S}}_3}$$
The eigenvalue ratio graphs show the shape of magnetic foliation (K3) data for all the granitoids. K denotes the shape parameter. The graph reads exactly same conventions of legends for different granites as is shown in Fig. 6.
(a) Ternary diagram, based on the indexes P (Point or Cluster), G (Girdle), R (Random) for all the granites. OriginLab (a data analysis and graphing software) was used to construct the ternary diagram. (b) The ln(S1/S2) vs ln(S2/S3) plot. C denotes the strength parameter. Note that Godhra and J. N. Kote granite show lowest and highest C values, respectively.
The above indices range from 0 to 1 with P + G + R = 1. All the studied granitoids show a K3 distribution pattern of less girdle (<50 %). The granitoid from Chitradurga North shows the most random K3 distribution, while the Chakradharpur granitoid displays the least random distribution. The strength parameter (C) of K3 is presented in Figure 6b. Any point in each solid line in this figure represents an equal strength value that increases away from the origin. A dataset having a perfectly uniform distribution will lie on the origin. Figure 6b shows that the J. N. Kote and Godhra bodies have the highest and lowest C values, respectively, thus indicating a maximum and minimum strength of their K3 distribution. The detailed discussion and conclusions from the above results are presented in Sections 6 and 7, respectively
6. Discussion
6.a. Estimating eigenvectors of magnetic fabric: a geometric approach coupled with mechanistic issues
Considering an individual magnetic foliation as a unit vector and following Scheidegger (Reference Scheidegger1965), we primarily construct the 3 × 3 ‘orientation tensor matrix’ (see also the discussion in Section 2) to compute three eigenvectors (V1, V2, V3) and the normalized eigenvalues (S1, S2, S3) associated with them. The above method has been applied for all the study areas, and therefore we obtain a different set of characteristic eigenvectors for each granitoid. Accordingly, assuming the individual datum of each distribution as a single point of unit mass within a sphere, an ellipsoid can be constructed. The three axes of the ellipsoid thus represent the direction of maximum, intermediate and minimum ‘moment of inertia’ of the assumed dissemination of masses. In such a scenario, Watson (Reference Watson1966) proposed that the direction of the maximum (V1) and minimum (V3) eigenvectors derived from an ‘orientation tensor matrix’ should follow the directions at which the ‘moment of inertia’ of the scattered masses is minimized and maximized, respectively.
For the magnetic fabric distributions in Malanjkhand (domain-II) and Chitradurga (south), the computed maximum eigenvectors (V1) show consistency in tracing the direction of the mean pole to K3, which leads to the coincidence between the strikes of V2–V3 and MF in these regions. Besides, the minimum eigenvectors (V3) define the poles to their best-fit great circles (Fig. 3). Hence, the characteristics of the fabric eigenvectors obtained from these two granites conform very well to Watson’s (Reference Watson1966) findings. However, in the rest of the cases, strikes of the V2–V3 planes show distinct angular variations with respect to the MF, thus not only contradicting our previous observations as well as Watson’s (Reference Watson1966) propositions, but also instantly drawing additional attention to the mechanistic issues associated with the origin of their magnetic fabrics. Earlier AMS studies suggested that orientation of the mean magnetic foliation (MF) can be treated as an Instantaneous Stretching Axis (ISAmax) for syntectonic granitoids, that deformed under a steady-state material flow system (Mamtani et al. Reference Mamtani, Pal and Greiling2013; Mamtani, Reference Mamtani2014; Mondal, Reference Mondal2018). Considering the analytical determination of the 2D vorticity of flow (Xypolias, Reference Xypolias2010), these investigations also estimated vorticity numbers in the studied granites, using (1) the angle between strikes of different planes (magnetic and field fabrics) on a horizontal plane, and (2) the parallelism between the direction of the shear zone and the extensional apophysis (Ae) of flow. It has been recorded that material flow in the Malanjkhand (domain-II), Chitradurga (south) and Godhra granitoid was close to the pure shear deformation, while simple shear was the dominant mechanism in the Malanjkhand (domain-I), J. N. Kote and Chitradurga (North) regions. Chakradharpur was the only granitoid that experienced general shear.
When correlating the above information with the obtained results, it can be physically accepted that the orientation of V2–V3, having been an estimation of the direction at which ‘moment of inertia’ of the distribution is maximized, should mimic the ISAmax of a particular steady-state flow under pure shear regime and retain its orientation throughout the all stages of progressive deformation. This is why the strikes of V2–V3 planes exactly coincide with that of MFs in Figure 3b1 and e. Contrasting this scenario, when the non-coaxial progressive deformations are considered, the rotational component of shear seems plausible to rotate the maximum eigenvector (V1) from the direction of minimized ‘moment of inertia’. Therefore, the rotational characteristics of eigenvectors with respect to a fixed coordinate on the stereonet physically explains the shift of V2–V3 orientation from MF (=ISAmax) in the simple-shear-dominated granitoids (e.g. Fig. 3b2, c and f). However, the results obtained from the Chakradharpur and Godhra granitoids (shaded region in Fig. 3) do not fit well with the above explanation as their V2–V3 planes also depict recognizable angular variations from ISAmax orientations in spite of their flows’ being under general and pure shear regime, respectively. We would like to state that such discrepancies are attributed to the unique bimodal scattering (see Fig. 4a, d) pattern of their magnetic fabric distributions in comparison with rest of the granitoids. It may be noted that both the bimodality and the substantial amount of rotational strain contributed to shifting the V2–V3 from MF (or ISAmax) for Chakradharpur granitoid, while the former was considered to be the sole reason behind the noted strike mismatch of V2–V3 and MF in Godhra. We elaborate on these deficiencies of using the eigenvalue approach, which are related to the scattering pattern and symmetry classes of fabric data distribution, further in the next subsection (6.3).
It also seems important to mention here that the position of the V2–V3 plane may conveniently be used in the 2D vorticity analysis of any flow under pure shear deformation, as it stays the same along with the direction of ISAmax and does not rotate with respect to flow apophyses of extension (Ae = shear zone) throughout the entire history of progressive deformation. Nevertheless, the rotational tendency of the V2–V3 plane with respect to ISAmax as well as Ae during the progressive stages of deformation under any sub-simple/simple shear confirms their inadequacy when performing vorticity analysis from them. Therefore, we conclude that estimation of fabric eigenvectors as well as determination of their connection with respect to the direction of ISAmax (=MF plane), at least after checking the mode and symmetry of any fabric distribution over the stereonet, can be utilized as a powerful tool to determine the type of shear mechanism (i.e. pure or simple) associated with the steady-state material flow of any particular region.
6.b. Normalized eigenvalue ratios: an indicator of the strength of a deformation event
The strength parameter diagram (Fig. 6b) shows that none of the C values of the granitoids lie at the origin, which confirms that no granitoid has magnetic foliation with a perfect uniform distribution. The arithmetical procedure used in the present paper to evaluate the strength parameters certainly shows its limitation when approaching with a structural dataset from poly-deformed regions. Since all the studied granitoids are known to have experienced more than one deformational event, we have consciously tried to incorporate the magnetic foliation data, as they manifest the late-stage deformation event only. It is evident from the figure that CJ.N.Kote > CChitradurga(north) > CChitradurga(south) > CChakradharpur > CMalanjkhand > CGodhra. This represents the sequence of the strength of K3 data of all granitoid bodies.
Although there exists a prominent variation of sample numbers (n) in K3 datasets for different granitoids, we assert that this would not anyway influence the reliability of obtained ‘C’ values in determining the strength parameters of magnetic fabrics. In accordance with Section 2 (‘methodology’ part), it may be noted that ‘C’ values are not computed directly from the eigenvalues (i.e. λ i ) but from their ‘normalized’ values (i.e. S i ), where S i = −λ i /n. Therefore, taking the summation over ‘n’ numbers of normalized direction cosines’ product (i.e. components of the ‘orientation tensor matrix’) ultimately gives rise to S1 + S2 + S3 = 1. In the current study, calculated values of S i for all the granitoids truly satisfy this relation. This is why we infer that the variation of sample number (n) does not at all affect the calculated C values [=(ln(S1/S3)] as well as our interpretation regarding the comparison of randomness for different granitoids.
6.c. Deficiencies of using eigenvalue approach in magnetic fabric analysis
The inconsistencies (associated with Fig. 3a, d) discussed in Subsection 6.1 clearly demand an overall realization of the limitations of using the eigenvalue method when analysing any structural fabric data that are relevant to coaxial or non-coaxial flows. In the present paper, all the analysis so far deals with the fabrics, which exhibit clusters with mainly orthorhombic symmetry. That the possibility of differences in modality may exist with individual clusters further leads us to plot the histograms of declinations of K3 in Fig. 4, separately for each granitoid. It may be noted that notwithstanding the rest of the granitoids, which more-or-less show unimodality, Chakradharpur and Godhra exhibit (Fig. 4a, d) a clear bimodal data distribution. This observation is found to be compatible with the observation made by Woodcock (Reference Woodcock1977) that for the bimodal or multimodal distributions, maximum eigenvector may not tend to exactly match the direction of the minimized ‘moment of inertia’. We consider this to be one of the reasons behind the mismatch of strike orientation between the V2–V3 plane and ISAmax in these two granitoids (Fig. 3a, d). Apart from the modality of the data distributions, we would like to add that the other important factor which can induce such a mismatch is the symmetry of the fabric distributions. As stated before, all the granitoids exhibit overall orthorhombic symmetry, which is another essential condition (in addition to unimodality and coaxiality of deformation) behind the observed coincidence of V2–V3 with the direction of maximized ‘moment of inertia’ (=ISAmax/MF) in pure-shear-dominated regions (Fig. 3b1, e). Considering other fabric symmetries such as axial, spherical, monoclinic and triclinic (Turner & Weiss, Reference Turner and Weiss1963, pp. 43–4), the eigenvalue approach would be useful in identifying the ISAmax in only the first two cases as previous investigations confirmed the coaxiality between V1 and the minimum ‘moment of inertia’ in such symmetries (Woodcock & Naylor, Reference Woodcock and Naylor1983). However, since the eigenvector analysis imposes an orthorhombic symmetry on the analysed data even if the distribution has a different kind of symmetry, extra attention should be paid when approaching with the monoclinic and triclinic data distribution. Therefore, although the examples of atypical granitoids, i.e. Chakradharpur and Godhra, serve the purpose of studying significant deficiencies of using the eigenvalue method in any structural fabric data analysis, we insist that, for any unimodal distribution of orthorhombic symmetry, estimation of fabric eigenvectors can be used to directly assess the style of deformation mechanism (pure/general/simple), whilst their corresponding eigenvalues provide some important aspects to decipher the fabric shape as well as strength of that deformation. In addition, considering the inadequacies of the eigenvalue approach in the above-mentioned complexities, we must not deny the necessity of future research in this realm, which would quantify the shift of the maximum eigenvector from the minimized ‘moment of inertia’ and their relationship with the mechanistic issues involved in such cases.
7. Conclusions
In the current study, we have used the eigenvalue method to determine the shape and strength of the magnetic foliations in various granitoids from India. The study also highlights the relationship between the positions of the eigenvectors and the mode of shearing that are responsible for the development of the magnetic fabric.
Below we summarize the major findings of the present study:
-
1) The positional uniqueness of characteristic eigenvectors (V1, V2 and V3) of magnetic fabrics can be concluded to be an important geometrical aspect, which directly indicates the mode of shearing (coaxial/general/non-coaxial)-related mechanistic issues associated with any particular deformed region. In the case of an orthorhombic unimodal fabric distribution of a syntectonic granitoid deformed under pure shear regime, the orientation of the V2–V3 plane exactly mimics that of the instantaneous stretching axis (ISAmax) of the steady-state material flow system associated with that deformation.
-
2) In contrast to the pure-shear-dominated granitoids, the V2––V3 plane shows its rotating tendency with respect to ISAmax and flow apophysis of extension (=shear zone) in the case of simple-shear-dominated granitoids.
-
3) The P, G and R ternary diagram is a representation which can simultaneously quantify and compare the randomness of the magnetic fabric distributions of different syntectonic granitoids.
-
4) K and C are two parameters which can be used to classify several datasets according to their shape and strength. These are also two representatives of the uniformity of any dataset.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S0016756822000747
Acknowledgements
The present research is funded by Indian Statistical Institute. Suggestions provided by Professor H N Woodcock and discussions with Professor Manish A Mamtani helped to improve the manuscript considerably. Professor Jean-Luc Bouchez and other, anonymous, reviewers are gratefully acknowledged for providing critical reviews on the earlier version of the manuscript. Rajdeep Mondal is thanked for helping with drawings and discussions. Detailed reviews by Dr. Emilio L Pueyo and an anonymous reviewer helped to improve the paper considerably. Editorial handling by Professor Olivier Lacombe is greatly appreciated. Logistic supports provided by Gourav Das, Sirshendu Kumar Biswas, Amlan Ghosh, Sutanu Das and Swarnasree Mondal are acknowledged.
Conflicts of interest
None.
Open access
