Diagenesis of Magnetic Minerals in Active/Relict Methane Seep: Constraints From Rock Magnetism and Mineralogical Records From Bay of Bengal

In this study, we conducted a comprehensive investigation of rock magnetic, mineralogical, and sedimentological records of sediment cores supplemented by a high resolution seismic data to elucidate the controls of structural and diagenetic (early vs. late) processes on the sediment magnetism in active and relict cold seep sites in the Bay of Bengal. Two distinct sediment magnetic zones (Z-I and Z-II) are defined based on the down-core variations in rock magnetic properties. The sediment magnetism is carried by complex magnetic mineral assemblages of detrital (titanomagnetite, titanohematite) and authigenic (fine-grained greigite) minerals. Overall, the magnetic susceptibility varies over one order of magnitude with highest values found in relict core. Uppermost sediment magnetic zone (Z-I) is characterized by higher concentration of magnetite as seen through elevated values of magnetic susceptibility (χlf) and saturation isothermal remanent magnetization (SIRM). A systematic gradual decrease of χlf and IRM1T in Z-I is attributed to the progressive diagenetic dissolution of iron oxides and subsequent precipitation of iron sulfides. Magnetic grain size diagnostic (ARM/IRM1T) parameter decreases initially due to the preferential dissolution of fine-grained magnetite in the sulfidic zone (Z-I), and increases later in response to the authigenic formation of magnetite and greigite in methanic zone (Z-II). Distinct low S-ratio and χlf values in methanic zone of relict core is due to increased relative contribution from highly preserved coercive magnetic (titanohematite) grains of detrital origin which survived in the diagenetic processes. A strong linkage between occurrence of authigenic carbonates and greigite formation is observed. Two plausible mechanisms are proposed to explain the formation and preservation of greigite in Z-I and Z-II: 1) decline in methane flux due to massive hydrate accumulation within the active fault system and formation of authigenic carbonate crust in the sub-surface sediments hindered the supply of upward migrating fluid/gas; thereby limiting the sulfide production which preferentially enhanced greigite formation in Z-I and 2) restricted supply of downward diffusing sulfide by the carbonate layers in the uppermost sediments created a sulfide deficient zone which inhibited the pyritization and favoured the formation of greigite in the methanic zone (Z-II).

Methane migrates from deep-seated gas reservoirs and interacts with downward diffusing seawater sulfate to trigger multiple zones of anaerobic oxidation of methane (AOM) Knittel and Boetius, 2009;Lin et al., 2018) driven by sulfate reduction (Mazumdar et al., 2019). At seep sites, non-steady and multiple diagenetic fronts would yield a complex pattern of magnetic mineral assemblages (Roberts, 2015;Reilly et al., 2020). Recent study by Reilly et al. (2020) provides new insights on magnetic mineral diagenesis through a shallow sulfate-methane transition zone (SMTZ) in a sediment core from Perseverance Drift, Northwestern Weddell Sea, Antartica. They reported that non-steady state disturbances including fluctuations in carbon flux and sedimentation can shift the pore-water geochemical profiles resulting in preservation of authigenic ferrimagnetic iron sulfides. In contrast, steady-state diagenetic processes at relict sites would lead to dissolution of detrital magnetic minerals followed by subsequent diagenetic formation of iron sulfide (pyrite) which can significantly reduce magnetic susceptibility of the bulk sediments (Canfield and Berner, 1987;Roberts and Weaver, 2005;Mohamed et al., 2011;Kars and Kodama, 2015;Roberts, 2015;Amiel et al., 2020;Lin et al., 2020;Reilly et al., 2020).
Magnetic minerals are potential recorders of the geological and methane-related geochemical processes and provide vital information on the methane seepage dynamics and diagenetic processes in marine gas hydrate sedimentary system (Musgrave et al., 2006;Larrasoaña et al., 2007;Kars and Kodama, 2015;Roberts, 2015;Badesab et al., 2019). Previous rock magnetic studies on sediment cores overlying shallow and deep-seated gas hydrate deposits in the K-G basin mainly aimed at delineating the magnetic signatures of detrital and diagenetic processes associated with evolution of gas-hydrate system (Dewangan et al., 2013;Badesab et al., 2017). For example, a magnetic based proxy for deciphering paleo-methane seepage events and sediment dispersal patterns, tracking of the rapid sedimentation events and its control on the magnetic mineral diagenesis and evolution of gas hydrate deposits was developed (Usapkar et al., 2014;Badesab et al., 2017;Badesab et al., 2019;Badesab et al., 2020b). Furthermore, the mechanism of formation and preservation of greigite in deep methanic sediments were unravelled by Badesab et al. (2020c). The occurrence of silicate-hosted greigite inclusion were reported by Badesab et al. (2020a). They further highlighted its importance in understanding magnetic mineral diagenesis in the gas hydrate system. More recently, Badesab et al. (2020c) conducted a comprehensive rock magnetic and mineralogical investigation on a seep-impacted sediment core from a newly discovered site and successfully delineated the processes governing the methane seepage dynamics and evolution of shallow gas hydrate (2-3 mbsf) system in the K-G basin. However, a dedicated rock magnetic study focussing on understanding the complex diagenetic (early and late) processes controlled by variability in methane fluxes and underlying fault/fractures at active and relict seep dominated sedimentary systems is still lacking.
Sediment cores retrieved during gas hydrate discovery cruise in 2018 in the K-G basin provides an excellent opportunity to establish a magnetic mineral inventory and investigate the diagenetic processes at active/relict and paleo-seep sites associated with shallow and deep-seated marine gas hydrates respectively. The active seep sites (SPC-01 to SPC-06) are mainly characterized by gas flares as observed in the water-column images, higher methane flux and presence of live chemosynthesis dependent organisms predominantly Bivalvia, Gastropods, and Polychaete (Mazumdar et al., 2019;Dewangan et al., 2020). Relict seep site (Stn-11/GC-02) is characterized by the absence of gas flares, presence of non-living shell fragment, Calyptogena Sp. and occurrence light brown colored authigenic carbonates at 48 cm below sea floor at this site. Contrastingly, paleo-seep site (MD161/Stn-8) is marked by the presence of methane-derived carbonates, non-living with chemosynthetic clams (Calyptogena sp.) and authigenic carbonates within 17-23 mbsf in this core which provided evidence of seepage of methane and sulfide-bearing fluids to the seafloor in the geological past (Mazumdar et al., 2009). In this study, we conducted a comprehensive investigation of rock magnetic, mineralogical, and sedimentological records of sediment cores supplemented by high resolution seismic data to elucidate the controls of structural and diagenetic (early vs. late) processes on the sediment magnetism in an active (SSD-45/Stn-4/SPC-01 to SPC-06) and relict (SSD-45/Stn-11/GC-02) cold seep sites in the Bay of Bengal.

GEOLOGY OF KRISHNA-GODAVARI BASIN
Krishna-Godavari (K-G) basin is situated along the eastern continental margin of India ( Figure 1). The K-G basin covers onshore area of 28,000 km 2 and offshore area about 145,000 km 2 and extends from Ongole in South to Vishakhapatnam in north (Ojha and Dubey, 2006). Average sediment thickness is about 3-5 km in onshore region, which increases to about 8 km in the offshore region (Prabhakar and Zutshi, 1993;Bastia, 2007). The Krishna and Godavari rivers and their tributaries deliver major sediment load to the K-G basin. Annual sediment load delivered by the Krishna and Godavari rivers is estimated to be around 67.7 × 10 6 and 170 × 10 6 metric ton respectively (Biksham and Subramanian, 1988;Ramesh and Subramanian, 1988). Both the river systems flows through the Deccan trap basalt and Precambrian metamorphic rocks and supplies large amount of detrital magnetite-rich sediments to the K-G basin (Ramesh and Subramanian 1988;Sangode et al., 2007). Sediments in the K-G basin are primarily composed of montmorillonite clay with traces of illite and kaolinite (Rao, 1991). Neo-tectonic activities in the K-G basin resulted in the formation of geomorphic structures including shale diapirs, mounds, and faults. Migration of methane through the fault system favoured the accumulation of gas hydrates in the K-G basin (Dewangan et al., 2010).

Sampling and Measurements
During gas hydrate exploration cruise (SSD-45) of Council of Scientific and Industrial Research -National Institute of Oceanography (CSIR-NIO) onboard R/V Sindhu Sadhana, six short spade cores (SPC-01 to SPC-06; Table. 1) were retrieved from active site (SSD-45/Stn-4) and one gravity core (Stn-11/GC-02) from relict seep site (as evident through the presence of non-living shell fragment, Calyptogena Sp. and authigenic carbonates at 48 cmbsf in this core) in the K-G basin ( Figure 1). Upon recovery, sediment cores from active sites smelled strongly of hydrogen sulfide (H 2 S) suggesting a shallow SMTZ that motivated this study. Magnetic samples were taken onboard immediately after collection of pore water samples to minimize sediment alteration due to warmer temperatures and oxygen exposure. The cores were subsampled in the presence of high purity nitrogen flushing to minimize the exposure of samples to atmospheric conditions, which could oxidize hydrogen sulfide or iron mono-sulfide. Samples for magnetic analysis were collected by pushing nonmagnetic plastic bottles into the split halves at ∼2 cm intervals. After sub-sampling, magnetic samples were transferred to a cold room (4°C) for storage. Magnetic measurements were carried out on wet sediment samples (packed in 25 mm plastic sample bottled) at the Paleomagnetic laboratory of CSIR-NIO and specialized rock magnetic measurements on dry samples at Center for Advanced Marine Core Research (CMCR), Kochi University, Japan.

Magnetic Measurements
Low-field magnetic susceptibility measurements were carried out using a Bartington Instruments MS2B dual frequency susceptibility meter. The susceptibility was measured at low (lf) and high (hf) frequencies χ lf 0.47 kHz and χ hf 4.7 kHz. The frequency-dependent susceptibility was calculated as χ fd % (χ lf − χ hf )/χ lf × 100%. All remanent magnetizations were measured using an AGICO JR-6A spinner magnetometer. An anhysteretic remanent magnetization (ARM) was imparted using a 100 mT alternating magnetic field superimposed on a fixed direct current (DC) bias field of 50 μT. An isothermal remanent magnetization (IRM) was imparted using an inducing field of +1T in the forward direction by a MMPM10 pulse magnetizer. S-ratio is calculated as the ratio of IRM at −300 mT and SIRM (IRM −300 mT /IRM 1T ) (Thompson and Oldfield, 1986). Thermomagnetic measurements were conducted on eight selected samples using a Natsuhara Giken (Model NMB-89) magnetic balance at CMCR with a heating rate of 10°C/min in a 0.3-T field.
Hysteresis loops, first-order reversal curves (FORC) (Pike et al., 1999), and back-field curves were also measured for 14 selected samples with a maximum field of 1T (number of loops per FORC: 91, averaging time: 200 ms, field increment: 4 mT, and a slew rate limit: 1T/s) at CMCR. Processing of FORC diagrams were done using the FORCinel software (Harrison and Feinberg, 2008). FORC distributions were processed using a smoothing factor (SF) of 5, but SF 6 were used where possible. Low-temperature magnetic measurements were conducted on selected samples using a Quantum Design Magnetic Properties Measurement System (MPMS) at CMCR. A room temperature SIRM (RT-SIRM) was imparted at room temperature (300 K) in 2.5 T. Samples were then cooled to 5 K and warmed back to 300 K in a zero magnetic field. A low temperature SIRM (LT-SIRM) was then imparted at 5 K in 2.5 T. Samples were warmed up to 300 K in a zero magnetic field (termed "ZFC" for zero fieldcooled). Samples were then cooled to 5 K in the presence of a 2.5 T magnetic field. A LT-SIRM was again imparted at 5 K, and samples were warmed to 300 K in a zero magnetic field (termed "FC" for field-cooled).

Bulk Sediment Grain Size Analyses
Sediment grain size measurements were carried out using a Malvern Mastersizer 2000 Laser Particle Size Analyzer at CSIR-NIO. Before analyses, desalination of sediment samples were carried out by washing repeatedly with deionized water followed by decarbonisation using dilute HCl (1N) acid. For removal of organic carbon, samples were treated with 10% H 2 O 2 and later, Sodium hexametaphosphate was added to sediments for better dispersal of the sediment fractions. Grain size values are reported as volume%.

Mineralogical Analysis
Magnetic particles were extracted from bulk sediments following the separation method of Petersen et al. (1986). Bulk sediment samples were disintegrated into 500 ml of distilled water with addition of 200 mg of a dispersing agent (sodium hexametaphosphate) and suspensions were treated with ultrasonics. Further, sediment suspension was allowed to continuously flow through the magnetic extraction set-up. The extraction method is based on the principle that high magnetic field gradients are created by inserting magnets on the teflon coated magnetic finger. During continuous circulation, magnetic particles within the bulk sediments get attracted towards the magnetic finger. These magnetic fractions are removed and washed several times in a small volume of distilled water to remove clay fractions and other contaminants. Magnetic particles were mounted on a carbon tape and were carbon coated for imaging. Images of magnetic grains were captured in secondary electron (SE) imaging mode at energy levels between 15 and 20 keV using a scanning electron microscope (SEM; JEOL JSM-5800 LV) at CSIR-NIO. Energy dispersive X-ray spectroscopy (EDS) probe attached to the microscope was used to determine the composition of magnetic minerals. Magnetic mineralogy was determined on the extracted magnetic particles using a Rigaku X-Ray Diffractometer (Ultima IV). The samples were allowed to run from 15°to 65°of 2θ at 1°/min scan speed using Cu Kα radiation (λ 1.5414 Å).

Down-Core Changes in Rock Magnetic Parameters
The magnetic susceptibility (χ lf ) profile of all spade cores (SPC-01 to SPC-06) from active seep and a gravity core from relict (SSD-45/Stn-11/GC-02) sites along with IRM 1T , ARM/IRM 1T , IRM 1T /χ lf and S-ratio are shown in Figure 2. Two distinct sediment magnetic zones (Z-I and Z-II) are defined based on down-core variations of χ lf in all cores ( Thompson, 1999;Peters and Dekkers, 2003). A general trend of downcore decrease suggesting coarsening in magnetic grain size is observed till the bottom of Z-I in all cores except SPC-02 ( Figure 2I). In SPC-01, SPC-03, SPC-04, SPC-06, and GC-02 ARM/IRM 1T first decreases probably due to preferential dissolution of fine-grained magnetite in Z-I, and rises downcore in response to the authigenic formation of greigite ( Figures Thompson, 1999;Peters and Dekkers, 2003;Snowball, 1991;Snowball and Thompson, 1990). A trend of initial rise in IRM 1T /χ lf values are noticed in Z-I of cores SPC-01, SPC-02, SPC-03, SPC-05, and GC-02 suggest the presence of fine grained magnetic particles (Figures 2D,J,P,A2). After initial rise, a gradual decline in IRM 1T /χ lf values till end of Z-II of the same cores is noticed (Figures 2D,J,P,A2,A14). Overall, S-ratios vary between 0.82 and 0.99 for majority of the samples suggesting that bulk magnetic mineralogy is characterized by both soft as well as hard coercive minerals in the cores from active seep sites (Figures 2E,K,Q,W,A3,A9; Frank and Nowaczyk, 2008). A systematic downcore decrease in χ fd (%) is seen in all cores from active sites ( Figure 2F,L,R,X,A4,A10). High values in Z-I indicate higher concentration of fine-grained magnetic particles, while relatively lower values in Z-II suggest the relative dominance of coarse magnetic particles. Interestingly, we noticed an opposite trend in χ fd (%) profile of relict site. Lower χ fd (%) values were found in Z-I and vice-versa ( Figure 2A16). Magnetic record of relict site (SSD-45/Stn-11/GC-02) showed similar magnetic zonation (Z-I, Z-II), but exhibit different downcore trends in ARM/IRM 1T , IRM 1T /χ lf and S-ratio ( Figure 2A11-A16). Three sediment intervals (11 cmbsf, 13 cmbsf, 15 cmbsf) within Z-I showed higher χ lf , IRM 1T and IRM 1T /χ lf suggesting the presence of fine-grained ferrimagnetic iron sulfides in these samples ( Figures 2A11,A12,A14). An anomalous increase in ARM/IRM 1T between 39 and 57 cmbsf suggests the dominance of fine-grained magnetic particles in this interval ( Figure 2A13). A major shift in S-ratio at 29 cmbsf followed by a minor rise in IRM 1T downcore (29-33 cmbsf) indicates that the relict site has much higher concentration of high-coercive (coarse-grained) minerals ( Figures 2A12,A15) as compared to other sites.

MAGNETIC MINERALOGY AND GRAIN SIZE DIAGNOSTIC PARAMETERS Thermomagnetic Measurements
Figures 3A-H show the thermomagnetic curves of the representative sediment samples from Z-I to Z-II of cores SPC-01/02/03/05. A steady decrease in χ between 552 and 633°C in majority of the samples imply that the bulk magnetic mineralogy is governed by magnetite (Dunlop et al., 1997). Two samples from SPC-01 (21 cmbsf) ( Figure 3B) and SPC-02 (41 cmbsf) ( Figure 3D) showed a major drop in χ at 695 and 673°C respectively. A minor increase in χ between 357 and 473°C for most of the samples could be either due to the presence of titanomagnetite with a wide range of Ti-contents (Lattard et al., 2006) or because of conversion of paramagnetic minerals into magnetite during heating process (Hirt and Gehring, 1991;Pan et al., 2000;Passier et al., 2001;Figures 3A-H).

FORC Diagrams
FORC diagram provided additional information on the type of magnetic minerals and their domain states. Representative FORC diagrams for Z-I and Z-II samples of cores SPC-01/ 02/03/06 are shown in Figure 4A-N. In majority of the samples, FORC distributions are characterized by closed contours with a peak at Bc ∼10-12 mT which suggest that magnetic particles exhibit vortex state to multidomain (MD) type behavior (Figures 4A-N;Roberts et al., 2017;Lascu et al., 2018;Roberts et al., 2020).

Low-Temperature Measurements
Results of the low-temperature magnetic measurements displaying characteristic magnetic signatures with respect to the sediment magnetic zones (Z-I, Z-II) ( Figure 5). A typical Verwey transition (Tv) temperature indicative of presence of magnetite ( Figures 5A,B,D-F) is noticed in RT-SIRM, ZFC-FC, and in the first derivative of magnetization curves (Verwey, 1939;Verwey, 1947;Özdemir et al., 2002;Chang et al., 2016a). Samples from magnetically reduced zone (Z-II) do not exhibit any Tv. The displayed LT magnetization curve ( Figure 5C) at 29 cmbsf of core SPC-01 is indicative of maghemite (Özdemir and Dunlop, 2010). Approximately 62% loss of RT-SIRM induced at 5 K is seen during warming between 5 and 30 k ( Figure 5C) indicating the presence of SP size magnetic particles. We do not observe any Besnus transition (typical for pyrrhotite; Besnus and Meyer, 1964;Dekkers, 1989;Larrasoaña et al., 2007) in the studied samples.

Grain Size Distributions of Sediments, Chemosynthetic Communities and Authigenic Carbonates
Bulk sediment grain size in the studied cores mainly composed of silt and clay size fractions ( Figure 9). Sand content is relatively low and

Correlation Between Magnetic Parameters of Active and Relict Methane-Seep Sites
A clear distinction in terms of magnetic mineral concentration, grain size and mineralogy is seen between active and relict seep sites Frontiers in Earth Science | www.frontiersin.org March 2021 | Volume 9 | Article 638594 11 ( Figure 10). A positive correlation between χ lf and IRM 1T is observed in all the samples, but showed different slopes which can be attributed to the presence of different magnetic mineralogies ( Figure 10A). Contrastingly, anomalous lower values (0.11-0.55) of S-ratio and wide range of χ lf displayed by samples from relict site indicates the predominance of highly coercive minerals ( Figure 10C). Coupling between magnetic grain size (ARM/IRM 1T ) and physical grain size is noticed in majority of samples from active and relict seep sites ( Figure 10D). We observed two populations (A and B). The samples from group-A showed uniformity in physical grain size, but exhibit wide range of magnetic grain size, while group-B samples follow equivalent trend and showed slight variation in magnetic and large scattering in physical grain size ( Figure 10D).

DISCUSSION
Room, low and high temperature based rock magnetic data, FORC diagrams, XRD and SEM-EDS observations provide vital information about complex magnetic mineral assemblages in the studied cores. Data revealed that the sediment cores from active seep sites have undergone significant diagenetic dissolution in Z-I, while antigenically formed fine-grained magnetite and greigite is predominant in Z-II. Contrastingly, sulfitic (Z-I) and methanic (Z-II) zones of relict core are devoid of authigenic magnetite and mainly dominated by large amount of coarsegrained titanomagnetite and highly coercive (titanohematite) grains of detrital origin which most likely survived the diagenetic dissolution and remained preserved in these zones ( Figures 8I-P). Here, we analyse these information to elucidate the controls of structural and diagenetic (early vs. late) processes on the sediment magnetic signals in an active/relict cold seep sites in the Bay of Bengal.

Magnetic Mineral Inventory of Sediments at Active/Relict Cold Seep Sites in the K-G Basin
Magnetic mineral assemblages in the sediment cores from active and relict seep sites from K-G basin, primarily consists of detrital (titanomagnetite, titanohematite), diagenetic (greigite/pyrite) and authigenic (  Frontiers in Earth Science | www.frontiersin.org March 2021 | Volume 9 | Article 638594 8I-P). Higher values of grain size indicative (ARM/IRM 1T ) ratio (Maher, 1988) in Z-II of cores SPC-01/03/04/06, and GC-02 reveals high concentrations of fine SD/PSD magnetite ( Figures  2C,O,U,A7,A13). ARM/IRM 1T first decreases because of preferential dissolution of fine-grained magnetite in the upper part of the reduced zone (Z-I), and increases later in response to the authigenic formation of fine-grained authigenic magnetite (Rodelli et al., 2019;Lin et al., 2020) and greigite in deeper positions downcore (Z-II). Diagenetically formed magnetic iron sulfide (greigite) are also found to co-exist with detrital titanomagnetite grains and are spatially distributed in both sediment magnetic zones (Z-I and Z-II) of cores from active and relict seep sites (Figures 6, 7, 8I-P, 11). These can be explained based on the fact that few titanomagnetite grains might have survived the early diagenetic dissolution and got rapidly buried and remain preserved due to increased sedimentation. Verwey transition (T v ) temperatures has been successfully utilized as a diagnostic for discriminating detrital magnetite in sulfitic marine sedimentary systems. For example in Oman continental margin (Chang et al., 2016a), offshore of Western Australia , and Ria de Vigo, Spain (Mohamed et al., 2011). Low-temperature magnetic data of SPC-01 and SPC-02 revealed the presence of magnetite ( Figures 5A,B,D-F; Moskowitz et al., 1993;Özdemir et al., 1993;Verwey, 1947;Muxworthy and McClelland, 2000;Prozorov et al., 2007;Chang et al., 2013;Chang et al., 2014;Chang et al., 2016a;Chang et al., 2016b). At relict seep site (GC-02), magnetic mineralogy in the high χ lf zone (Z-I) comprised of abundant detrital coarse-grained titanomagnetite and titanohematite particles ( Figures 8I,J). Within sulfidic zone, minor peak of high χ lf , χ fd %, and IRM 1T /χ lf between 11 and 13 cmbsf reflect the diagenetic formation of greigite in this interval ( Figures  2A11,A14,A16, 8K-L). Reduced methane flux due to either closure of underlying fractures/fault or cut-off in the supply of downward diffusing sulfate concentrations most likely interrupted the reductive diagenetic process and rather favored the formation and preservation of greigite by limiting H 2 S production at this site Neretin et al., 2004). Overgrowth of sub-micron sized particles on the larger titanomagnetite grain in Z-II is suggestive of authigenic growth of greigite probably due to weakening of anaerobic oxidation of methane (AOM) process in GC-02 ( Figure 8O).
Variability in rock magnetic parameters in each core reflects the different degree of diagenetic disturbances created by methane seepage within the studied sites in the K-G basin. Similar observations were reported from Vestnesa Ridge, Svalbard margin (Sztybor and Rasmussen, 2017) and Oregon continental margin (Karlin, 1990). Variations in S-ratio can be attributed to the changes in magnetic mineralogy ( Figures  2E,K In spade cores (SPC-01 to SPC-06), we did not observe any signatures of higher coercivity (hematite) in SEM-EDS and XRD data (Figures 6, 7). In fact, XRD data confirmed the presence of greigite in samples exhibiting low S-ratio in both Z-I and Z-II (Figures 6B,C,E,F,K). The S-ratio drop in Z-II of relict core can be well-explained by the occurrence of numerous titanohematite grains ( Figure 2A15, 8M,O,P). Preservation of titanomagnetite grains ( Figures 8M-P) in Z-II could be either due to weakened AOM (Reidinger et al., 2005) or higher resistance against dissolution (Nowaczyk, 2011;Rodelli et al., 2019) which might have hindered the early diagenetic processes and led to the preservation of detrital magnetic particles. Significant drop in S-ratio ( Figure 2A15) followed by high χ fd % ( Figure 2A16) in Z-II of relict core (SSD-45/Stn-11/GC-02) can be attributed to the presence of mixture of detrital titanohematite ( Figures 8M,O,P), authigenic greigite (Dillon and Bleil, 2006;Larrasoaña et al., 2007;Rodelli et al., 2019) and hexagonal 3C pyrrhotite (Horng, 2018;Horng and Roberts, 2018). Similar observations were recently reported from Perseverance Drift, North-western Weddell Sea (Reilly et al., 2020).
Authigenically formed ferrimagnetic iron sulfides (greigite and pyrrhotite) are well-reported in methane-rich sedimentary system, for instance in Nankai Trough (Horng, 2018;Kars and Kodama, 2015), Hydrate Ridge (Larrasoaña et al., 2007), and offshore Taiwan (Horng and Roberts, 2018). FORC diagrams and low-temperature magnetic measurements provided additional constraints on the magnetic mineralogy and domain states. FORC distributions showed the mixture of vortex to multidomain magnetic particles ( Figures 4A-N). In natural system, greigite can commonly occur as magnetostatically interacting framboids. A micromagnetic modelling study involving hysteresis and FORC simulations demonstrated that magnetostatic interactions could remarkably alter the FORC response of greigite framboids and can form patterns similar to vortex and MD (Rodelli et al., 2018;Valdez-Grijalva et al., 2020). It is quite possible that greigite clusters observed in SPC-02 ( Figures 7E,F) and SPC-04 ( Figure 7J) might be contributing significantly towards the observed FORC (vortex to multi-domain) distributions. Low temperature magnetic data indicated the presence of magnetite in Z-I and Z-II ( Figures 5A,B,D-F). Similar to our findings, Korff et al. (2016) reported cyclic magnetite dissolution and preservation of detrital magnetite in methanic sediments of the abyssal northwest Pacific Ocean. Sample from the zone of reduced magnetic susceptibility (Z-II) of SPC-01 did not show any indication of a Verwey transition in RT-SIRM, ZFC-FC, and the first derivative of magnetization curves ( Figure 5C). We observed that nearly 62% of the LT-SIRM applied at 5 K is lost during warming between 5 and 30 K ( Figure 5C). This observation indicates the presence of superparamagnetic (SP) magnetic particles in these samples (Tarduno, 1995;Passier and Dekkers, 2002). Similar to our observation, Kars and Kodama (2015) reported a related pattern of sudden decrease in LT-SIRM from 5 to ∼30 K and about ∼50-60% loss of the applied remanence in the sediment core samples from Nankai Trough, offshore Japan. They linked the loss of remanence to the presence of SP particles which get unblocked very rapidly. The SP magnetic particles in this sample ( Figure 5C) is most likely iron sulfides like greigite, pyrrhotite, titanohematite (Horng, 2018;Rodelli et al., 2019;Roberts et al., 2020). This interpretation can be explained by the fact that fine-grained magnetic particles in this sample would get readily dissolved in the sulfitic environment (as seen through lowest χ lf in Z-II). Contrastingly, the SP greigite and pyrrhotite nanoparticles are thermodynamically more stable compared to the SP magnetite nanoparticles and would rather remain preserved Roberts et al., 2018;Rodelli et al., 2018). Lack of LT transitions in the magnetization derivation curves in this sample ( Figure 5C) is persistent and further confirms the presence of greigite (through high χ fd %, Figure 2F) and pyrrhotite (Chang et al., 2009;Roberts et al., 2011).

Early Versus Late Diagenesis of Magnetic Minerals in Active and Relict Methane Seep Sites in the Bay of Bengal: Constraints From Rock Magnetism and Mineralogical Proxies
In active cold seep sedimentary system, temporal movement of diagenetic fronts due to variability in methane fluxes and sulfate reduction rates favours non-steady state diagenesis. Such process could generate complex magnetic signatures in sediments involving magnetite reduction and authigenic formation of magnetite (Amiel et al., 2020;Lin et al., 2020;Riedinger et al., 2014) and iron sulfides (Canfield and Berner, 1987;Larrasoaña et al., 2003;Liu et al., 2012;Riedinger et al., 2014;Roberts, 2015). Rock magnetic data showed progressive downcore dissolution of magnetite in sediment cores from active (Figures 2A-A10) and relict ( Figures 2A11-A16; Karlin and Levi, 1983) sites in K-G basin. A gradual downcore decrease in ARM/IRM 1T and χ lf within Z-I in cores SPC-01/03/04/05/06 suggest substantial dissolution of fine-grained magnetic particles due to early diagenesis ( Figures 2C,O,U,A1,A7; Rodelli et al., 2019). A slight increase in ARM/IRM 1T in Z-II of SPC-01/03/04/06 ( Figures 2C,O,U,A7) can be attributed to the presence of finegrained authigenic magnetite (Karlin et al., 1987;Roberts, 2015;Roberts et al., 2018;Rodelli et al., 2019;Lin et al., 2020). A similar steady increase of ARM/IRM 1T in sulfidic zone of a sediment cores from deep methanic zone of Southern Eastern Mediterranean continental shelf (Amiel et al., 2020) and Niger deep-sea fan was attributed to the gradual authigenic growth of Frontiers in Earth Science | www.frontiersin.org March 2021 | Volume 9 | Article 638594 SD magnetite (Dillon and Bleil, 2006). A downcore increase in ARM/IRM 1T in Z-II of active seep sites (SPC-01/03/04/06) could also be due to presence of authigenic magnetite formed as a result of microbial iron reduction process in methanic sediments which might have permitted the newly formed magnetite to survive the sulfidic dissolution and favoured preservation (Rodelli et al., 2018;Rodelli et al., 2019;Lin et al., 2020;Figures 2C,O,U,A7). SEM-EDS and XRD data of cores showed the presence of greigite in methanic zone (Z-II) from active seep sites ( Figures 6C,E,F,K, 7E,F,K,L). Following the concept put forward by Lin et al. (2020), we propose that following mechanism for authigenic greigite formation in Z-II of active seep sites. Increased supply of reactive iron (secondary iron oxyhydroxides) produced during iron-sulfide mineral oxidation due to decline in methane seepage intensities (Lin et al., 2016;Lin et al., 2020) and limitation of sulfide in methanic zone (below SMTZ) probably favoured the greigite formation in this zone. During early diagenesis, sulfate gets readily depleted through decomposition of organic matter (Karlin, 1990a). We put forward the plausible mechanism for the preservation of fine-grained authigenic magnetite in Z-II as seen through a trend of high ARM/IRM 1T in SPC-01/03/04/06 ( Figure 2C,O,U,A7). Finegrained magnetite (SD-type) could also occur as magnetic inclusions in silicates which most likely protected the magnetite from diagenetic dissolution (Hounslow and Morton 2004;Feinberg et al., 2005;Tarduno et al., 2006;Muxworthy and Evans, 2013;Chang et al., 2016a;Chang et al., 2016b) and could explain the preservation of fine-grained magnetite in Z-II. Another possibility is that authigenic magnetite formation coupled to microbial iron reduction (Lin et al., 2020) within methanic sediments might have allowed the freshly formed finegrained magnetite to survive sulfidic dissolution fronts and rather favoured preservation as evident through high ARM/IRM 1T values in Z-II (methanic) of cores SPC-01/03/04/06 from active seep sites. However, we do not rule out the possibility of presence of fine-grained titanomagnetite in Z-II. Similar findings were reported from rock magnetic studies on sediment cores from continental margin of Oman northwestern Arabian Sea (Chang et al., 2016b) and northeast part of Japan Sea (Yamazaki et al., 2003). At relict site, an anomalous (39-57 cmbsf) interval within (Z-II) showed increased ARM/ IRM 1T values and significant drop in S-ratio (29-63 cmbsf) indicating the preservation of numerous highly coercive magnetic particles as confirmed through SEM observations ( Figures 2A13,A15, 8M−P). We proposed that overgrowth of sub-micron sized authigenically formed SP greigite particles on the larger titanomagnetite grain (at 43 cmbsf, Figure 8O) might have contributed towards the observed increase in IRM 1T /χ lf ( Figure 2A14). The presence of SP greigite is further confirmed by the observed increase in χ fd % in Z-II of GC-02 ( Figure 2A16).
Reduced magnetization at gas hydrate vents has been established in Cascadia accretionary margin off Vancouver Island (Novosel et al., 2015). At active seep sites, methane is being transported from deep-seated gas reservoir up to the seafloor and the intensities of methane flux governs the SMTZ location. Presence of abundant live chemosynthetic communities (Bathymodiolus sp. and Acharax sp.) provide evidences of active methane seepage at the studied site ( Figures 8A-F). Elevated methane flux will shift the SMTZ fronts upward and vice-versa (Borowski et al., 1996). Production of hydrogen sulfide via AOM at SMTZ causes the dissolution of iron oxides, and precipitation of iron sulfides (Reidinger et al., 2014). Diagenetic dissolution of magnetite followed by the subsequent transformation to pyrite will create a distinct χ lf minima compared to ferrimagnetic greigite which would significantly increase χ lf (Dewangan et al., 2013;Kars and Kodama, 2015;Badesab et al., 2019). Several researchers established that variability in the SMTZ depth can be used to quantify methane flux (Housen and Musgrave, 1996;Neretin et al., 2004;Novosel et al., 2005;Reidinger et al., 2005;Musgrave et al., 2006;Larrasoaña et al., 2007;Roberts, 2015). A close linkage between the depth of SMTZ and minima in χ lf has been wellestablished and increasingly utilized as a proxy for detecting SMTZ fronts in methanic sediments for example in Cascadia margin (Housen and Musgrave, 1996;Larrasoaña et al., 2007), offshore of Argentina and Uruguay (Riedinger et al., 2014), Nankai trough, Japan (Kars and Kodama, 2015), and offshore New Zealand (Rowan and Roberts, 2006). In the analyzed sediment cores, the distinct χ lf minima most probably reflect recent SMTZ front, which varies between 5 cmbsf and 16 cmbsf (Figure 2). The shallowest χ lf minima is noticed in active seep core SPC-05 ( Figure 2Y) and deepest at SPC-03 ( Figure 2M) and relict site GC-02 ( Figure 2A11), respectively. Such variability in SMTZ depth provides evidences on the differential rate of diagenesis constrained by fluctuations in methane fluxes. In such nonsteady state diagenetic scenario, for instance at active seep sites (SPC-01 to SPC-06), different degree of diagenetic disturbances created by changing methane flux could variably modify and preserve the primary magnetic mineral assemblages in sediments. SEM observations indicate the presence of abundant coarse-grained titanomagnetite and titanohematite in Z-I ( Figures  7A,D,G, 8I-P) and Z-II ( Figures 7C,I, 8M-P). More recently Amiel et al. (2020) conducted a detailed magnetic investigation on the sediment cores from SE Mediterranean Continental Shelf to examine the influence of early diagenesis in methanic environment. They reported that rapid sedimentation and nonsteady state conditions triggered by variability in methane flux favoured the survival of abundant iron-oxide minerals in the sulfidic and deep methanic zone (Riedinger et al., 2014). In the present study, survival and preservation of these minerals in Z-I and Z-II can be explained by the fact that titanohematite and titanomagnetite are more stable and offers strong resistance to reductive dissolution induced by early and late diagenetic processes (Channell and Hawthorne, 1990;Poulton et al., 2004;Nowaczyk, 2011;Rodelli et al., 2019) or due to increased sedimentation (Riedinger et al., 2014;Amiel et al., 2020). Stability against diagenetic alteration increases together with elevated Ti content and, in turn, results in a higher magnetic stability. Dissolution of titanohematite involves a noticeably slower reaction compared to magnetite (Cornell and Schwertmann, 1996;Yamazaki et al., 2003;Liu et al., 2004;Rodelli et al., 2019). Previous studies have demonstrated that increasing titanium content stabilizes titanomagnetite compounds because Ti 4+ substitution correspondingly reduces the number of Fe 3+, the acting electron acceptor under anaerobic conditions (Karlin, 1990;Garming et al., 2005). During non-steady diagenesis, formation of secondary minerals takes places whenever pore-water chemistry changes due to upward or downward diffusion of dissolved hydrogen sulfide  and could lead to the formation and preservation of authigenic magnetite and greigite. Microbial iron reduction within methanic sediments is the most likely mechanism for authigenic magnetite formation (Rodelli et al., 2018;Lin et al., 2020).
More recently Reilly et al. (2020) and Rodelli et al. (2019) revealed that high coercivity magnetic minerals respond slowly to the diagenetic dissolution through entire SMTZ and can be a potential source for ferric iron via microbial iron reduction process coupled to AOM (Beal et al., 2009;Riedinger et al., 2014;Amiel et al., 2020). They further concluded that increased concentration of high coercivity minerals in methanic zone (below SMTZ) can produce large amount of ferric (Fe 3+ ) via iron reduction where sulfide is limited and favours the authigenic precipitation of greigite. In the relict core (SSD-45/Stn-11/GC-02) greigite and titanohematite are found to occur throughout the sulfidic (Z-I) and methanic (Z-II) ( Figures 2A11,A16, 8I-P, 11A-D). It is highly possible that diagenetic dissolution of titaohematites during iron reduction generated enough (Fe 3+ ) and provided conducive environment for greigite formation in the intervals where sulfide got exhausted in core SSD-45/Stn-11/GC-02. Our SEM ( Figures 8K,L,O) and XRD ( Figures 11B,C) observations followed by high values of χ fd % ( Figure 2A16) support the findings of Reilly et al. (2020). Bivariate magnetic plots ( Figure 10) and FORC diagrams ( Figure 12) also helped to evaluate the influence of steady and non-steady state diagenetic processes on the magnetic record of sediment cores from active (SPC-01 to SPC-06), relict (Stn-11/GC-02) and paleo (MD161/ Stn-8; Mazumdar et al., 2009) seep sites in the K-G basin. Scattered data plot between IRM 1T /χ lf and χ lf reflects a broad range of SP greigite ( Figures 2A16, 8K,L,O) particles and a larger variations in their concentrations ( Figure 10B). A scattered plot between magnetic (ARM/IRM 1T ) and physical (mean) grain size showed a systematic pattern of fining in magnetic grain size with increase in diagenetic dissolution of detrital magnetic minerals ( Figure 10D). Non-steady state diagenetic processes can potentially affect the sediment magnetic parameters. For example, abrupt change in the sedimentation rate can cause migration of the SMTZ (e.g., Fu et al., 2008;Riedinger et al., 2014), temporal changes in the organic matter load into the sediment can affect the depths of the different respiration processes (Larrasoaña et al., 2003;Rowan et al., 2009), changes in the sediment composition can lead to the accumulation of magnetic minerals in the sediment and change in the magnetic parameters. These mechanisms require time-variations in the flux of the sediments from Krishna and Godavari rivers, which are main source of detrital magnetic particles to the studied area. Multiple layers of mass transport deposits (MTD's) generated as a result of frequently occurring high sedimentation events have been identified in the shallow and deep offshore regions of the K-G basin (Ramana et al., 2007;Ramprasad et al., 2011;Yamamoto et al., 2018). Recently Badesab et al. (2019) reported a high resolution sediment magnetic record of high sedimentation events in the K-G basin. Samples from relict site (SSD-45/ Stn-11/GC-02) showed remarkably high χ lf and IRM 1T ( Figures 2A11,A12) compared to active seep sites (SPC-01 to SPC-06) probably due to non-steady state diagenesis driven by high sedimentation events. It is most likely that coarse titanomagnetite and titanohematite grains survived diagenesis due to 1) its inherent property to resist the dissolution (Nowaczyk, 2011;Rodelli et al., 2019), and 2) rapid burial due to increased sedimentation (Badesab et al., 2019). Detailed rock magnetic and microscope analyses conducted on the sediment cores from continental margin of Oman and northern California by Rowan et al. (2009) reported that newly formed greigite particles are small initially and would display thermally unstable SP behaviour unless it grows fully through the stable SD blocking volume. They attributed the magnetic enhancement to the late diagenetic growth of magnetically stable greigite. Greigite can display both SD-type behaviour (Roberts, 1995) as well as ultrafine-grained SP (Rowan and Roberts, 2006;Rowan et al., 2009). Differences in the growth of greigite particles from active (SPC-01 to SPC-06) and paleo (MD161/Stn-8; Badesab et al., 2020a) samples is clearly reflected in the FORC diagrams ( Figures 12A-D). FORC diagram of samples from active seep site (SPC-01) showed the magnetic particles exhibit a mixed (vortex state to multidomain) type behavior ( Figures 12A,B; Roberts et al., 2017;Roberts et al., 2020). In contrast, FORC diagrams of paleo-seep samples showed broad distribution with concentric and elongated contour indicative of presence of magnetically stable greigite probably formed as a result of late diagenetic growth ( Figure 12C; Rowan et al., 2009) and biogenic mineralization ( Figure 12D; Rodelli et al., 2018;Badesab et al., 2020a) at site MD161/Stn-8.

Structural Control on the Formation and Preservation of Greigite
Methane venting episodes linked to active and paleo-cold seeps are well documented in the K-G basin (Mazumdar et al., 2009;Dewangan et al., 2020). High-resolution seismic data across active seep site (SSD-45/Stn-4; present study) showed the presence of well-developed seismic chimney and faults, which facilitated the migration of deep-seated gas through the gas hydrate stability zone (GHSZ) up to the sea floor and favored gas hydrate accumulation in the K-G basin (Dewangan et al., 2020; Figure 1B). At the same site, shallow (2-3 mbsf) gas hydrates were recovered during a gas hydrate exploration cruise in 2018 (Mazumdar et al., 2019). Vertical migration of methane through the sedimentary column produces multiple AOM's resulting in the growth of diagenetic minerals including pyrite, greigite and pyrrhotite (Roberts, 2015). Interplay between the rates downward diffusing sulfate and upward migrating methane concentrations control the amount of hydrogen sulfide which is a key factor constraining the greigite formation (Kasten et al., 1998;Jørgensen et al., 2004;Neretin et al., 2004). For instance, enhanced methane flux would generate large amount of hydrogen sulfide which could drive pyritization process to complete resulting in the precipitation of pyrite, while limitation of sulfide due to decline in methane would favor greigite/pyrrhotite formation (Liu et al., 2004). At active site, presence of greigite in both sediment magnetic zones of cores SPC-01/02/04 has been confirmed through SEM-EDS and XRD data (Figures 2, 6, 7). We try to investigate the structural control on the variability in methane fluxes and subsequent diagenetic formation of greigite in the studied cores. Control of opening and closing dynamics of the underlying fracture/fault system on methane migration in the K-G basin has been well-established Dewangan et al., 2020). Spatial distribution of authigenic carbonates in cores SPC-01, SPC-04, SPC-06 ( Figure 13) provide evidence on several AOM's and is indicative of past fluctuations in methane fluxes. We put forward following plausible mechanisms to explain the formation and preservation of greigite in Z-I and Z-II: 1) decline in methane flux either due to massive hydrate accumulation (Mazumdar et al., 2019;Dewangan et al., 2020) or hindering of upward migrating fluid/gas by the carbonate layers in the subsurface sediments most likely resulted in low sulfide production and preferentially favoured formation and preservation of greigite in the Z-I. Similar observations were reported in fault-controlled cold seep-hydrate Woolsey Mound in the Northern Gulf of Mexico (Simonetti et al., 2013), gas hydrate field offshore Vancouver Island (Riedel et al., 2002), Vestnesa Ridge, W-Svalbard Margin (Singhroha et al., 2020), and in the K-G basin (Badesab et al., 2017;Dewangan et al., 2011). They explained that plugging of an active fault system due to massive hydrate formation can also cause significant drop in methane flux. 2) Occurrence of greigite below the carbonate layers in cores SPC-01, SPC-04, SPC-06 ( Figure 13) in Z-II suggest that supply of downward diffusing sulfide is most likely hindered and subsequently a sulfide deficient zone is formed below the carbonate layer which inhibited the pyritization and rather favoured the formation and preservation of greigite. A detailed rock magnetic analysis of a 30 m long sediment core from paleoseep site in the Bay of Bengal reported the presence of both diagenetic and biogenic greigite (within 17-23 mbsf) i.e., beneath a thick layer of authigenic carbonate formed as a result of intense methane seepage event (Mazumdar et al., 2009;Dewangan et al., 2013;Badesab et al., 2020a). More recent study has demonstrated that elevated supply of Fe 2+ via iron reduction produced by high coercivity minerals in sulfide limited zone (Z-II) can also favour the authigenic growth of greigite (Rodelli et al., 2019;Reilly et al., 2020). Greigite presence in low S-ratio intervals of Z-II in relict core (SSD-45/Stn-11/GC-02) can also be linked with such process (Figures 2A15, 8O, 11C,D). The presence of greigite is further confirmed by the observed increase in χ fd % ( Figure 2A16). However additional geochemical and microbiological data are warranted to resolve the mechanism.
A conceptual model explaining the evolution of magnetic mineralogies at active (Stn-4/SPC-01 to SPC-06) and relict (Stn-11/GC-02) seep sites is developed (Figure 14). Shallow gas hydrates (2-3 mbsf) and authigenic carbonates (throughout the cores) were recovered at active seep sites (Figure 13; Mazumdar et al., 2019). In sulfidic (Z-I), detrital magnetic particles minerals supplied through Krishna and Godavari river systems reacted with the H 2 S produced by bacterial activity via organic matter decomposition and AOMcoupled sulfate reduction which led to the diagenetic dissolution of magnetic particles and subsequent precipitation of iron sulfide (pyrite) with intermediate precursors such as greigite and pyrrhotite (Berner, 1984;Karlin, 1990a). SEM-EDS and XRD data coupled with χ fd % confirmed the presence of greigite in both sulfidic (Z-I) and methanic (Z-II) (Figures 2, 6, 7). Most likely mechanism which can explain the formation and preservation of greigite in Z-I and Z-II is as follows. Decline in methane flux due to massive hydrate accumulation along the active fault system and formation of authigenic carbonate crust in the sub-surface sediments might have hindered the supply of upward migrating fluid/gas and limited the sulfide production which preferentially enhanced greigite formation in Z-I. Secondly, restricted supply of downward diffusing sulfide by the carbonate layers in the uppermost sediments might be a responsible process for creating a sulfide deficient zone beneath. These might have inhibited the pyritization process and rather favoured the formation of greigite in methanic zone (Z-II). Samples from relict site (SSD-45/Stn-11/GC-02) showed notably high χ lf and IRM 1T ( Figure 2A11,A12) compared to active seep sites which is most probably due to non-steady state diagenesis driven by high sedimentation events in the K-G basin. Coarse titanomagnetite and titanohematite grains might have survived diagenesis either due to its inherent property to resist the dissolution (Nowaczyk, 2011;Rodelli et al., 2019), or because of rapid burial due to increased sedimentation (Badesab et al., 2019).

CONCLUSION
A set of rock magnetic, mineralogical and sedimentological records of the sediment cores were used to characterize the changes in the magnetic concentration, grain size, and mineralogy in a newly discovered active (Stn-4/SPC-01 to SPC-06) and relict (Stn-11/GC-02) cold seep sites in the K-G basin, Bay of Bengal. Sediment magnetism is mainly carried by complex magnetic mineral assemblages of detrital (titanomagnetite, titanohematite) and authigenic (magnetite, greigite) minerals. Topmost sulfidic sediment magnetic zone (Z-I) is marked by a higher concentration of magnetic minerals as evident in rock magnetic parameters. A systematic downcore decrease in χ lf and IRM 1T in this zone is due to the progressive diagenetic dissolution of iron oxide and subsequent precipitation of iron sulfides. Magnetogranulometric proxy (ARM/IRM 1T ) provided valuable insights on the preferential dissolution of fine-grained magnetite in the sulfidic zone (Z-I), and authigenic formation of magnetite in methanic zone (Z-II) of all studied cores. At relict site, elevated concentration of highly coercive magnetic (titanohematite) grains relative to soft minerals (as seen through lower S-ratio and χ lf ) which survived diagenesis are found in the methanic zone. A strong linkage between occurrence of authigenic carbonates and greigite formation and preservation is noticed. We summarized our observations into a conceptual model which provides insights on the evolution of different magnetic mineralogies constrained by variable diagenetic processes at active and relict seep sites in the K-G basin, Bay of Bengal.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: CSIR-NIO's data repository (https:// publication-data.nio.org/s/pb7ZQHybKm5DqDj).

AUTHOR CONTRIBUTIONS
VG, FB, PD, and MK developed the concept for the manuscript. PD and VG have contributed towards sediment core collection and sampling during the research expedition (SSD-

ACKNOWLEDGMENTS
Director, CSIR-NIO, Secretary MOES, and National Gas Hydrate Program (NGHP) are thanked for supporting this study. We greatly acknowledge the support of Dr Aninda Mazumdar, participants of scientific cruise SSD-045 onboard research vessel R/V Sindhu Sadhana of CSIR-NIO, ship cell team and crews onboard for providing assistance in collection of sediment cores and sampling. Director of the Center for Advanced Marine Core Research (CMCR), Kochi University, Japan is thanked for allowing access to the Paleomagnetic laboratory during a research visit of FB as a part of INSA-JSPS fellowship. Myriam Kars is thanked for providing support to FB during the fellowship period at CMCR. High-end rock magnetic analyses at CMCR were supported by DST-JSPS under the scheme JSPS Postdoctoral Fellowship for Overseas Researchers (Strategic Program; JSPS grant GR 17202) and Establishment of Young Researcher Fellowship Program 2017-2018 (IA/Indo-Japanese/F-1/2018) scheme awarded to F. B. We thank Chisa Nishimori for providing laboratory and technical support in Kochi. Thanks to Areef Sardar, Girish Prabhu, Kaveri Salunke and Prem Piyush Panda for their support in carrying out mineralogical analyses at CSIR-NIO. We thank associate editor (Luigi Jovane) and all three reviewers for constructive comments that improved the paper. The data used in the present study are available on CSIR-NIO's data repository (https://publication-data. nio.org/s/pb7ZQHybKm5DqDj). This is CSIR-NIO publication no. 6671.