Using the high-temperature phase transition of iron sulfide minerals as an indicator of fault slip temperature

The transformation of pyrite into pyrrhotite above 500 °C was observed in the Chelungpu fault zone, which formed as a result of the 1999 Chi-Chi earthquake in Taiwan. Similarly, pyrite transformation to pyrrhotite at approximately 640 °C was observed during the Tohoku earthquake in Japan. In this study, we investigated the high-temperature phase-transition of iron sulfide minerals (greigite) under anaerobic conditions. We simulated mineral phase transformations during fault movement with the aim of determining the temperature of fault slip. The techniques used in this study included thermogravimetry and differential thermal analysis (TG/DTA) and in situ X-ray diffraction (XRD). We found diversification between 520 °C and 630 °C in the TG/DTA curves that signifies the transformation of pyrite into pyrrhotite. Furthermore, the in situ XRD results confirmed the sequence in which greigite underwent phase transitions to gradually transform into pyrite and pyrrhotite at approximately 320 °C. Greigite completely changed into pyrite and pyrrhotite at 450 °C. Finally, pyrite was completely transformed into pyrrhotite at 580 °C. Our results reveal the temperature and sequence in which the phase transitions of greigite occur, and indicate that this may be used to constrain the temperature of fault-slip. This conclusion is supported by field observations made following the Tohoku and Chi-Chi earthquakes.

www.nature.com/scientificreports www.nature.com/scientificreports/ material 19,20 , the characterization of magnetic minerals 21,22 , high-speed rotary shear experiments 23,24 , among others. Greigite and pyrite are very stable under pressures of up to 3.4 GPa and 40 GPa 25,26 , respectively, which correspond to depths of ~125 km and ~1450 km beneath the earth's surface, respectively. The epicenters of high magnitude earthquakes that have occurred in the last 30 years (e.g., Tohoku earthquake in Japan, Wenchuan earthquake in China, Chi-Chi earthquake in Taiwan, and Hanshin earthquake in Japan) were all located at depths less than 30 km, and are classified as shallow focus earthquakes. Offshore drilling research found that pyrite was converted into pyrrhotite as a result of high-temperature friction in the fault-slip zone in the offshore sediments of northeastern Japan (Tohoku earthquake) 27 . A study of the 1999 Chi-Chi earthquake in Taiwan also observed that high temperatures generated by the fault slip caused pyrite to transform into pyrrhotite 28 . Therefore, the high-temperature phase transition of iron sulfide minerals such as greigite (Fe 3 S 4 ) may be another method for measuring fault slip temperature.
We therefore investigated the high-temperature phase transition of pyrite in nature, caused by fault slip, by designing a series of experiments in which greigite undergoes a high-temperature phase transition in an anaerobic environment. Our motivation for using greigite (rather than pyrite) as the starting material is that greigite is considered to be an intermediary in the conversion of mackinawite (FeS) to pyrite (FeS 2 ) [29][30][31][32] . Although greigite is a metastable iron sulfide mineral, it has been found in a deep core from Qinghai Lake in China, where it was preserved in Pliocene-age sediments 2 . A study of sedimentary rocks from southwestern Taiwan also found a stratum containing magnetic minerals such as greigite and pyrrhotite 33 . Greigite is an interesting iron sulfide mineral, very few studies that discuss the greigite-pyrite-pyrrhotite transition under high-temperature conditions have been published. Therefore, our intent is to show that the sequence in which the phase transitions occur corresponds to our prediction, namely that greigite first transforms to pyrite, which is then converted to pyrrhotite. The second transformation in this sequence resembles the phase transition of pyrite as it occurs in nature when it is caused by friction due to fault slip.

Methods
In this study, greigite was first synthesized via a hydrothermal method 34 . The phase transition temperature of greigite was then obtained by carrying out thermogravimetric and differential thermal analyses (TG-DTA). Moreover, the phase diagram of iron and sulfur was calculated by using Fact Sage software. The high-temperature phase transformation and the transformation pathway of greigite under anaerobic conditions was observed using synchrotron radiation in situ X-ray diffraction (XRD) measurements. The related details are as follows.
Synthesis. Synthetic greigite was fabricated by hydrothermal synthesis according to the following procedure.
The solution was deployed, filtered, and then stored in an anaerobic glove box. Iron (II) sulfate heptahydrate (FeSO 4 •7H 2 O) and thioacetamide (C 2 H 5 NS) were formulated into 0.2 and 0.3 M solutions, respectively, such that the volume of each was 50 mL. The solutions were then mixed and loaded into an autoclave. The anaerobic conditions in the autoclave were maintained by injecting nitrogen gas. The mixture was heated to 115 °C and held at this temperature for 2 h to fabricate greigite. After cooling, the solution was filtered with deionized water and an acetic acid filter and then dried in a freeze dryer for more than 12 h to obtain greigite powders.
Analysis. The fabricated greigite was confirmed to consist of the pure phase by XRD (Bruker D8 Advance with a wavelength of Cu Kα, a scanning rate of 1°/min, and a scanning step of 0.02°; the working voltage was 30 kV and the current was 20 mA). Next, we observed the weight change and exothermic heat release reaction of the sample at high temperatures by TG-DTA (Netzsch STA 449 F3 at a heating rate of 10 °C/min, from 25 to 900 °C). We used 5N ultra-high-purity nitrogen and an oxygen trapping system to reduce the oxygen content in the chamber. Finally, the sequence of greigite phase changes was examined by in situ XRD with synchrotron radiation. The in situ XRD diffraction patterns were recorded in the BL01C2 and BL17A1 beamlines of the National Synchrotron Radiation Research Center (NSRRC) in Taiwan. The energy of the storage ring was 1.5 GeV and the current was 360 mA. The wavelength of the beamline was 1.03 Å and 1.32 Å for BL01C2 and BL17A1, respectively. The exposure time was 72/192s. Both beamlines used a MAR345 detector, which is a two-dimensional area detector. The original 2D data were integrated into a one-dimensional diffraction pattern by GSASII, and the wavelength was converted to Cu Kα by Winplot software for a clear comparison.

XRD. The diffraction peaks of the synthesized samples were compared to the values on the International
Centre for Diffraction Data (ICDD) card 01-089-1998 (Fig. 1); however, the diffraction pattern was shifted to a slightly higher angle relative to that on the card. Greigite crystallizes in the cubic crystal system and the measured lattice parameter of the crystallographic axes is 9.845 Å (9.876 Å on the ICDD card). The observed differences in the pattern may have been due to defects resulting from the mismatch in the Fe 2+ /Fe 3+ ratio in the greigite sample. Figure 2 shows the TG-DTA results for greigite (Fe 3 S 4 ). The DTA curve does not change significantly with temperature. The change in the TG curve below 100 °C is the result of desorption of water from the sample. The TG curve shows a small decrease at 300 °C and a large change at approximately 590 °C, indicating weight loss during heating. Assuming that the weight loss was the loss of sulfur from the iron sulfide mineral, the variation in the ratio between iron and sulfur can be calculated (as marked in Fig. 3 with green arrows). Figure 3 presents the phase diagram of iron and sulfur calculated by Fact Sage software (the software algorithm 35,36 is based on the well-known Gibbs-Helmholtz relation: G = H − TS, where H is enthalpy, S is entropy, and G is free energy calculated from this formula), with an ordinate temperature and abscissa of the molar ratio of sulfur in the iron sulfide mineral (sulfur accounts for 57.1% of greigite). Greigite can be replaced with a stable mineral in the phase diagram because it is a metastable mineral. The iron/sulfur ratio as a function of temperature was calculated according to the weight change observed in the TG-DTA experiment, which is www.nature.com/scientificreports www.nature.com/scientificreports/ marked on the figure with green arrows. The calculation result showed that greigite undergoes conversion to pyrite at 317 °C and is converted to pyrrhotite between 600 and 700 °C. The most stable mineral phase would exist at the lowest free energy based on this software-simulation concept. Therefore, the phase transformation occurred in the sequence of greigite → greigite + pyrrhotite + pyrite → pyrrhotite + pyrite → pyrrhotite, which is consistent with the TG-DTA results.

Phase diagram.
In situ XRD. Figure 4 presents the XRD pattern of greigite and its phase transition at different temperatures.
The pure phase of greigite was stable from room temperature to 300 °C. According to the figure, pyrite and pyrrhotite (hexagonal structure) began to appear at 320 °C. The (311), (400), and (511) diffraction peaks of greigite disappeared at approximately 450 °C, and the (200) and (023) diffraction peaks of pyrite disappeared at approximately 580 °C. Only the pyrrhotite phase remained when the temperature reached 900 °C. After returning to room temperature, the final product of the heating process, pyrrhotite (monoclinic structure), was present 37 . This result suggests that the high-temperature phase transition is irreversible, and that pyrrhotite can exist as a stable phase under anaerobic conditions. Figure 5 illustrates the formation of iron sulfide and the sequence of high-temperature phase transitions of greigite. Mackinawite is well known to be the initial precipitate of iron sulfide minerals during the reaction of Fe 2+ with S 2− solutions. In a slightly oxidizing  www.nature.com/scientificreports www.nature.com/scientificreports/ or hydrogen sulfide rich environment, mackinawite changes into greigite and then into pyrite (i.e., pyritization). There are a lot of locations where greigite was found, for example in Qinghai Lake (China) 2 , in mudstones (Taiwan) 19 , among others. The deposition rate in these areas 2,19 was extremely rapid, thus the proportion of organic carbon and hydrogen sulfate in the sediments was diluted; this resulted in stagnant pyritization, causing the greigite to be preserved. A study of the 1999 Chi-Chi earthquake in Taiwan indicated that the high temperatures generated by the fault slip caused pyrite to transform into pyrrhotite at a transition temperature exceeding 500 °C 28 . Core drilling off the Pacific coast near the site of the Tohoku earthquake (2011) also determined that pyrite changed into pyrrhotite near the sliding surface and its corresponding transition temperature was greater than 640 °C 27 . Because pyrite or greigite was very stable up to 3.4 GPa, but the pressure of these fault-slips was less than 0.8 GPa and the duration for the fault-slip was very short (a few seconds or minutes, meaning the time effect is slight). Moreover, it is reported that the reaction of pyrite is primarily controlled by the temperature during the fault movement (pressure effect is minor) 27 . In this study, we observed via in situ XRD, that pyrite was completely transformed into pyrrhotite at temperatures above 580 °C. The TG-DTA results showed that pyrite transformed into pyrrhotite between 520 and 630 °C. Our phase transition temperature results were therefore similar to those  www.nature.com/scientificreports www.nature.com/scientificreports/ that were found in nature (i.e., the fault-slip temperatures in the Chi-Chi and Tohoku earthquakes). Our study suggests that the phase transition temperature at which greigite undergoes phase transsition into pyrrhotite can provide an indicator of the fault-slip temperature.

Conclusion
The temperature of a fault slip is difficult to determine. In this study, we synthesized greigite to investigate its high-temperature phase transitions and to simulate the temperature of fault movement. In the high-temperature phase transition process in an anaerobic environment, greigite started to transform into pyrite and pyrrhotite at 300 °C. At 600 °C, only the pyrrhotite phase was observed to remain. Our experimental results are similar to the observations made in the field, indicating that greigite can be used as an indicator of the fault-slip temperature. Figure 5. Diagram comparing the phase transition sequence and temperature of greigite, including pyritization, the preservation of greigite in Qinghai Lake, phase-change sequence during the Chi-Chi and Tohoku earthquakes, and the phase transformation of greigite determined in this study.