Petrography, Mineralogy, and Geochemistry of Thermally Altered Coal in the Tashan Coal Mine, Datong Coalﬁeld, China

: A suite of coal samples near a diabase dike were collected to investigate the petrographic, mineralogical, and geochemical characteristics of thermally altered coal in Proximate analysis, vitrinite reﬂectance measurement, and petrographic analysis were applied to identify and characterize the alteration halo; optical microscope observation, qualitative X-ray diffractometry, and SEM-EDS were applied to study the phases, occurrence, and composition of minerals; XRF, ICP-MS, and AFS were applied to determine concentrations of major and trace elements; and the occurrence modes of elements were studied by correlation and hierarchical cluster analysis as well as SEM-EDS. The results demonstrated that the 3.6 m dike has caused an alteration halo of approximately 2 m in diameter. In addition, the thermally altered coals were characterized by high vitrinite reﬂectance, low volatile matter, and the occurrence of thermally altered organic particles. Dolomite and ankerite in the thermally altered coal may be derived from hydrothermal ﬂuids, while muscovite and tobelite may be transformed from a kaolinite precursor. The average concentration of Sr in the Tashan thermally altered coal reached 1714 µ g/g, which is over 12 times that of the Chinese coal; the phosphate minerals and Sr-bearing kaolinite account for this signiﬁcant enrichment. The cluster analysis classiﬁed elements with geochemical associations into four groups: group 1 and 2 were associated with aluminosilicates, clays, and carbonates and exhibited enrichment in the coal/rock contact zone, indicating that the dike may be the source of the elements; group 3 included P 2 O 5 , Sr, Ba, and Be, which ﬂuctuate in coals, suggesting that their concentrations were inﬂuenced by multiple-factors; group 4 did not manifest obvious variations in coals, implying that the coal itself was the source. Contributions: Conceptualization, methodology, software, resources, curation, writing—original preparation, X.S. and H.M.; writing—review and editing, visualization, H.M.; supervision, X.S.; project administration,

The heating and expansion of igneous intrusions and the cooling and contraction of the surrounding rocks result in massive joints and fractures in coals [22], inducing enhanced permeability and declined stability of coal seams, which may lead to gas outbursts and roof falls [23]. Igneous intrusions reduce the economic viability of coal mines because of the difficult mining conditions they create. In addition, thermally altered coal, which is not marketable, may be abandoned in the practical mining process since it generally has high ash yields, low calorific values, and difficulties with combustion.
Coal is a type of sedimentary rock that is sensitive to changing temperatures and stress that may cause physical and chemical variations. Heating from igneous intrusions petrographically alters the coal. Thermally altered coal is macroscopically dull and dense [10,14];

Coal-Bearing Strata
The coal-bearing strata of the Datong coalfield include the Jurassic Datong, Lower Permian Shanxi, and Upper Carboniferous Taiyuan formations (Figure 2). At present, the Jurassic coal seam has been exhausted, and exploitation has turned to the underlying Carboniferous-Permian coal seams [21]. The Shanxi formation conformably overlies the Taiyuan formation. Its thickness ranges up to 96.1 m, with an average of 67.2 m. The coal seams include Shan 1 , Shan 2 , Shan 3 , and Shan 4 , of which only the Shan 4 coal seam is mineable, with an average thickness of 2.63 m.

Igneous Intrusion
Datong Coalfield was intruded by igneous intrusions in both the Indosinian (Late Permian to the end of the Triassic) and the Yanshan (Late Jurassic to the Early Cretaceous) periods [47][48][49][50][51]. Intrusions dating back to the Indosinian period are dominated by mafic to ultramafic lamprophyre and carbonatite [49,51], whereas rocks that intruded in the Yanshan period are mainly diabase and basaltic andesite [52]. The detailed lithological characteristics of the igneous intrusions in the Datong coalfield were described in Chen et al. [21]. In the practical mining process, sills and dikes were encountered in the Shan 4 , 2, 3 −1 , 3 −2 , 5 (3)(4)(5), and 8 coal seams in the Tashan Coal Mine. The large-scale igneous intrusions in coal seam affect the safety of coal mines by causing roof instability. For example, two severe roof falling incidents induced by igneous intrusion emplacement were recorded in the Tashan Coal Mine in 2005 and 2020. Moreover, they play a deleterious role in economic efficiency by reducing the mineable thickness and altering coal into cinder or natural coke, which are not marketable or combustible.

Sampling
The 8222 and 5222 tunnels of the Tashan Coal Mine parallel to each other and encountered the same diabase dike. The detailed characteristics of the dike were described in Song et al. [53]. Although we did not obtain photographs of the dike in the 5222 tunnel because of poor lighting conditions, the overlying dike in the 8222 tunnel exhibited the same characteristics ( Figure 3). The 5222 dike is 3.6 m in thickness, comprising a 2.6 m dark green core and two pale chill zones on each side. A suite of 16 coal samples was collected along the 5222 tunnel following the Chinese Standard Method GB/T 482-2008 [54]. The coal samples were collected at essentially the same height along the transect at closely spaced intervals in the vicinity of the dike and larger intervals farther away from the dike, up to a distance of approximately 20 m. All collected samples were immediately sealed in plastic bags to minimize contamination and oxidation.

Methods
The coal samples were subjected to proximate analysis according to ASTM Standard for moisture (ASTM D3173/D3173M-17a [55]), ash yield (ASTM D3174-12(2018)e1 [56]), and volatile matter (ASTM D3175-20 [57]). Particulate blocks for petrographic analysis were prepared according to the standard procedures described in ISO 7404-2 [58]. The mean maximum and random reflectance of vitrinite were measured using a Leica DM4500P microscope equipped with a DETA V4000 Spectrometer according to ISO 7404-5 [59]. Point counts were conducted under incident white light, immersion oil, at a 500× magnification using 500 counts per block and two blocks per sample to quantify the contents of the different macerals and evaluate the extent of thermal alteration. The identification of macerals followed the ICCP system 1994 proposed by the International Committee for Coal and Organic Petrology for vitrinite [60], inertinite [61], and liptinite macerals [62]. The classification of organic particles in thermally altered coals was performed using the ICCP classification system [63].
Optical microscopic observation, powder X-ray diffractometry (XRD), and scanning electron microscope analysis (SEM) were applied to study the mineralogical characteristics of the Tashan coal. Qualitative XRD analysis was applied to confirm the main mineralogical phases. Powdered coal samples were examined on a Bruker D8 advance powder X-ray Cu-Kα radiation diffractometer. The XRD pattern was recorded over a 2θ interval of 5-70 • with a step size of 0.02 • . The resulting diffractogram peaks were qualitatively analyzed using HighScore Plus software. Scanning electron microscopes (Hitachi SU8010 and Bruker Quanta FEG 450, 1kV and 20 kV accelerating voltage) in conjunction with an energydispersive X-ray spectrometer (SEM-EDS) were used to study the morphology, distribution, and composition of the minerals in the coal.
Samples were crushed and ground to pass 200 mesh for geochemical analysis. Concentrations of major and trace elements in the Tashan coal were determined at the Analytical Laboratory Beijing Research Institute of Uranium Geology following the methods described in Chinese Standards GB/T 14506-2010 [64]. The contents of the major element oxides, including SiO 2 , Al 2 O 3 , Fe 2 O 3 , MgO, CaO, Na 2 O, K 2 O, MnO, TiO 2 , and P 2 O 5 , were determined by X-ray fluorescence spectrometry (PANalytical Axios-max); high resolution inductively coupled plasma mass spectrometry (Thermo Fisher ELEMENT XR) was used to determine the concentrations of elements, such as Li, Be, Sc, V, Cr, Co, Ni, Cu, Zn, Ga, Rb, Sr, Mo, Cd, In, Sb, Cs, Ba, W, Tl, Pb, Bi, Th, U, Nb, Ta, Zr, and Hf; the concentrations of Se and As in the Tashan coal were determined using liquid chromatography-atomic fluorescence spectrometry (LC-AFS 6500).

Identification and Characterization of the Alteration Halo
An alteration halo was induced in the ambient area of the dike because of thermal effects, in which the samples exhibited high vitrinite reflectance, abnormal coal composition (ash and volatile matter), and changes in petrographic characteristics. Identification and characterization of the alteration halo were performed following the standard coal classification methods described in ISO 11760-2018 [65] and ASTM D 388-19a [66].

Vitrinite Reflectance
Vitrinite reflectance increased sharply in the immediate area of the intrusive body (Table 1, Figure 4), indicating an upgrading coal rank approaching the dike [3,67]. According to ISO 11760-2018 [65], T1-T3 are high-rank coals (or anthracite), whereas T1 and T2 are high-rank B coal (or anthracite B coal with a mean vitrinite random reflectance of ≥3.0% and <4.0 %), T3 is high-rank C coal (or anthracite C coal with a mean vitrinite random reflectance of ≥2.0% and <3.0%); T4 and T5 are medium rank B coal (or bituminous B coal with a mean vitrinite random reflectance of ≥1.0% and <1.4%); T6-T16 are medium rank C coals (or bituminous C coal with a mean vitrinite random reflectance of ≥0.6% and <1.0%).
The vitrinite maximum reflectance (R max ) is also indicative of coal rank upgrading. The coal within 1.6 m of the dike (T1-T5) had a R max higher than 1%, whereas the R max fluctuated between 0.67% and 0.94% in the more distant areas.
The volatile matter yield (dry, ash-free basis) of the Tashan coal ranged from 7.36% to 49.92% and increased with increasing distance from the intrusive body (Table 1, Figure 4). As suggested by the American standard ASTM D 388-19a [66], samples T1 and T3 are semianthracite, T2 is anthracite, T4 and T5 are low volatile bituminous coal, and T6-T16 are medium to high volatile bituminous coal.

Petrographic Characteristics
Intruding magma caused morphological and optical changes in the organic particles of thermally altered coals; for example, liptinite becomes difficult to distinguish [28,29,[33][34][35], and porous and massive altered organic particles as well as newly formed components occur [63]. These altered organic particles indicate that coals experienced plastic phases at high temperatures during the intrusive event [36], differentiating thermally altered coals from unaltered coals. Macerals of the Tashan thermally altered coal exhibited various melting and plasticity characteristics ( Figure 5); for instance, vitrinite was converted into porous coke in the adjacent area of the dike; the edges of the macerals became rounded, which may have been caused by melting; and fissures and cracks became prevalent in the coal, inducing a 'fish bone' structure in some isolated carbonized vitrinite [68]. Quantitative analysis was carried out to evaluate the extent of thermal alteration ( Figure 6, Table 2). For example, the contents of non-altered organic particles ranged from 32.4% to 43.4% in T1-T4 and reached a peak at 75.4% in T5; the contents of altered organic particles, however, exhibited a reverse trend compared with non-altered particles, ranging between 56.1% and 67.2% in T1-T4, followed by a sharp decrease to 23.8% in T5; and newly formed organic particles only accounted for a small proportion of particles, and no obvious trend was observed in T1-T5. The quantification of the organic particles in thermally altered coal demonstrated that thermal alteration dissipated at a further distance; fewer altered organic particles in T5 could be a sign of the heat pulse running out.  Macerals in T6-T16 showed prominent isotropy, indicating little or no thermal alteration. Table 3 summarizes the point counting analysis data of T6-T16, and the results indicate that vitrinite, inertinite, and liptinite in the unaltered coal accounted for an average of 65.7%, 29%, and 5.3%, respectively. The vitrinite group mainly comprised collodetrinite, telinite, and collotellinite, with average contents of 33.1%, 17.6%, and 10.8%, respectively; the inertinite group was mainly composed of semifusinite and fusinite, with average contents of 12.1% and 9.1%, respectively, whereas macrinite and micrinite only accounted for a small proportion of the contents; and liptinite was represented by sporinite. According to the classification by petrographic composition described in ISO 11760-2018 [65], most of the Tashan coals are of moderately high vitrinite, except T8, T11, and T12, which are medium vitrinite coals.

Characterization of the Alteration Halo
Both the ISO and ASTM classification methods indicated an elevated coal rank in the vicinity of the dike; however, the results were not very consistent. It was suggested that the T1-T3 samples were in the anthracite stage on the basis of ISO 11760-2018 [65], while ASTM D388-19a [66] classified T1 and T3 as semianthracite. This is likely because T1 and T3 have considerable quantities of carbonate minerals that can decompose during the volatile matter determination process, affecting the results to some extent [69,70]. In the present study, two methods as well as thermally altered coal petrographic analysis were combined to identify and characterize the alteration halo (Table 4). Accordingly, thermally altered coals had an R max of >1%; coals adjacent to the intrusion (T1-T3, <0.6 m) were altered to anthracite or semianthracite with low volatile matter and moderately high ash; T4 and T5, which were slightly farther away from the dike, were altered to bituminous B coal with low volatile matter; and the unaltered coal was bituminous C coal with medium-high volatile matter, low-medium ash, and an R max lower than 1%. Thermally altered organic particles were only observed in T1-T5, whereas macerals in T6-T16 did not manifest thermal alteration. Overall, proximate analysis, vitrinite reflectance, and petrographic characteristics consistently suggested that the edge of the alteration halo is located between T5 and T6.    Table 5 shows the main crystalline phases identified from the diffractograms of the Tashan coal. Minerals found in the unaltered coal were dominated by kaolinite, quartz, calcite, and dolomite and were less abundant in ankerite, analcime, diopside, and boehmite. Minerals identified in the thermally altered coal included kaolinite, dolomite, ankerite, quartz, muscovite, tobelite (NH 4 -illite), and anhydrite. Clay minerals in the Tashan coal are dominated by kaolinite, which was identified in all samples. Authigenic kaolinite occurring as cell fillings or replacing the organic constituent was found in both thermally altered and unaltered coals, for example, in the form of detrital pellets as well as book-like and vermicular aggregates (Figure 7a-c). This occurrence may indicate that authigenic kaolinite formed during peat accumulation or during early diagenetic processes [71,72]. However, kaolinite filling in the cleats and fractures as well as in devolatilization vacuoles was also observed in thermally altered coals (Figure 7d-e), indicating its epigenetic origin.
Kaolinite can be transformed into chlorite, muscovite, and illite with elevated coal ranks [38,[73][74][75][76]. According to the XRD results, muscovite occurred in the thermally altered coals T1-T4, which is consistent with the conclusion drawn by Golab and Carr [44] that muscovite is indicative of thermal alteration. In the present study, muscovite was intimately intergrown with kaolinite (Figure 7f), and the identification of muscovite was supported by SEM-EDS observations. The occurrence modes of muscovite indicate its epigenetic formation by the transformation of a kaolinite precursor.
The XRD pattern indicates the occurrence of tobelite (NH 4 -Illite) in T4 and T5, which was further demonstrated by SEM analysis. NH 4 -Illite has been found in some high-rank coals [76][77][78]. It was suggested by Dai et al. [77] that ammonian illite may form because of the interaction between kaolinite and nitrogen released from the organic matter at high temperatures during an intrusive event. In the present study, cleat-filling tobelite with flake morphology was observed through SEM analysis (Figure 7g). The identification of tobelite was supported by the SEM-EDS spectra in Figure 7h, and similar results were reported by Dai et al. [79]. The occurrence modes of tobelite differentiate it from other clays and may suggest an epigenetic origin.

Quartz
According to the XRD results, quartz occurred in most samples and was lacking in some samples (e.g., T4, T5, T8, T12, and T15), which may be because the quartz content was below the detection limit. However, quartz was rarely observed by optical analysis except in the T1 sampled directly adjacent to the intrusion, where, for example, epigenetic quartz occurring as cell-and fracture-fillings was observed (Figure 8). Chen et al. [37] reported quartz veins in thermally altered coals intruded by syenite porphyry and lamprophyre sills from the Huainan Coalfield, where the epigenetic quartz was derived from Si-rich hydrothermal fluids. In the present study, the restricted epigenetic mineralization of quartz reflects that the mafic diabase dike was unlikely to provide massive silicon-rich hydrothermal fluids during the intrusive event.

3.
Carbonate minerals Dolomite and ankerite were more prevalent than calcite in the thermally altered coal, while carbonates in the unaltered coal were dominated by calcite ( Table 5). Because of a lack of epigenetic mineralization, carbonates in unaltered coal were mainly of syn-depositional origins, such as the calcite filling in the fusinite (Figure 9). However, most dolomite and ankerite in thermally altered coal formed epigenetically, and this view has been supported by many studies [2,16,22,33,76,[80][81][82]. It is generally suggested that epigenetic carbonate minerals form through the reaction of magma or hydrothermal solution with CO and CO 2 released from coal at high temperatures during an intrusive event [2,37,71,72]. Dai et al. [77] suggested that during coal metamorphism, calcium released from the organic matter may interact with CO 2 in the pore water, forming carbonates. In the Tashan coals, epigenetic dolomite and ankerite included considerable quantities of Fe, Mg, and Ca, indicating that these minerals may be deposited by hydrothermal fluids. Epigenetic carbonate filling in the devolatilization vacuoles and cleats of macerals was observed in the Tashan coals, providing the following insights into the interactions between coal and rocks during the emplacement event. (i) The occurrence of thermal alteration microtextures demonstrates temperatures of more than 400 • C [83,84]. Volatile matter was released at this temperature (see Table 1, reduced V daf value of thermally altered coal), and devolatilization vacuoles formed and finally served as accommodation for deposited minerals (Figure 10a,b). (ii) Dikes intruded vertically and may serve as relatively longlived magma conduits [26], and heat expansion caused by the supplement of magmas followed by the shrinkage of host rocks eventually caused fractures in coal, for example, the development of prismatic fracturing perpendicular to the heating surface [24]. Joints and fractures in the coal acted as paths for magmas and hydrothermal fluids and were finally filled by carbonates (Figure 10c,d). (iii) Organic inclusions in carbonate minerals may reflect the alteration of magma on coal seams during emplacement (Figure 10e,f), that is, when a large amount of the magma invades quickly, organic debris may be carried and transported and eventually mixed in epigenetic minerals.

Major Elements in the Tashan Coal
The contents of major element oxides in the Tashan coal are shown in Table 6. Compared with unaltered coals, the thermally altered coals had enhanced average contents of major elements because of the devolatilization and interaction with magma ( Figure 11). SiO 2 and Al 2 O 3 were major components of the Tashan coal and averaged 7.26% and 6.47% in the thermally altered coals, respectively, slightly higher than the averages of 4.7% and 4.65% in the unaltered coals ( Table 6). The average contents of MgO and CaO in the thermally altered coals were 0.75% and 1.93%, respectively, 6.22 and 2.04 times the averages for the unaltered coals (Table 6). P 2 O 5 accounted for 0.31% on average in the thermally altered coal, slightly higher than the average of 0.18% in the unaltered coals. Na 2 O and K 2 O were enriched in the samples immediately adjacent to the intrusion, thus inducing increased average values in thermally altered coal ( Table 6). The contents of Fe 2 O 3 and MnO exhibited similar trends at different distances from the dike ( Figure 11); although they tended to be enriched in the coal/rock contact zone, their high values in samples T10 and T13 led to virtually equal average contents in the thermally altered and unaltered coal ( Table 6). The average contents of TiO 2 did not exhibit obvious differences in the thermally altered and unaltered coals (Table 6). The average contents of MgO, CaO, and P 2 O 5 in the thermally altered coal were significantly higher than those in average Chinese coals [45,85]-3.41, 1.57, and 3.40 times higher, respectively-while the Al 2 O 3 content was 1.08 times that of Chinese coals. However, the average contents of SiO 2 , Fe 2 O 3 , Na 2 O, K 2 O, MnO, and TiO 2 were lower than those of Chinese coals. The average P 2 O 5 content in the unaltered coal was 1.92 times higher than that in the Chinese coals, whereas the contents of other elements were lower. Table 7 summarizes the concentrations of trace elements in the Tashan coal. Lithium, V, Cr, Rb, Sr, Cs, Ba, W, Tl, Th, Nb, Ta, and As were more enriched in the thermally altered than in the unaltered coal, whereas Co, Ni, Ga, Mo, Cd, In, and Sb exhibited depletions in the thermally altered coal. The average concentrations of Be, Cu, Zn, Pb, Bi, U, Zr, Hf, Se, and Sc were virtually equal in both coals.

Trace Elements in the Tashan Coal
In the current study, the average element concentrations in the North China Carboniferous-Permian coal [85,86], Chinese coal [45], and average world coal [87] were collected to calculate concentration coefficients (CC) of elements in the Tashan coal; for example, CC i and CC i ' were used to denote the corresponding concentration coefficients in the thermally altered and unaltered coal, respectively (Table 8, Figure 12). The enrichment of elements in the Tashan coal was classified according to the method proposed by Dai et al. [88].   Figure 12), the Tashan thermally altered coal was significantly enriched in Sr (CC 1 = 11.1) and slightly enriched in Li (CC 1 = 3.3). In the unaltered coal, Sr (CC 1 ' = 6.0) was enriched, and Li (CC 1 ' = 2.2) was slightly enriched.
Compared with average Chinese coals, strontium was significantly enriched in both thermally altered (CC 2 = 12.2) and unaltered coals (CC 2 ' = 6.5). Ga was slightly enriched in both coals, with concentration coefficients of 2.5 and 2.9 in the thermally altered and unaltered coals, respectively.
Compared with average world coals, the Tashan thermally altered coal was significantly enriched in Sr (CC 3   Nd, no data; a data with * are from Tang and Huang [86], the others are from Ren et al. [85]; b data are from Dai et al. [45]; c data are from Ketris and Yudovich [87]; d average concentration of trace elements in thermally altered coal; e average concentration of trace elements in unaltered coal; CC i and CC i ' respectively denote the concentration coefficient of the corresponding element in thermally altered and unaltered coal calculated on the basis of the average element concentration of the thermally altered and unaltered coal samples versus the respective value established for the C-P coals in North China, Chinese coal, and world coal.

Modes of Occurrence of Elements
The modes of occurrence of elements in coal are important, as this parameter dictates the behavior of the element during coal combustion, beneficiation, conversion, weathering, leaching, or any other chemical reaction that the coal undergoes [89]. Studies on the modes of occurrence of elements in thermally altered coal may provide information on the sources of elements and elemental migration during the intrusive event. In this study, correlation and hierarchical cluster analysis, as well as SEM-EDS, were used to study the affinities and geochemical associations of elements in the thermally altered coals.

Correlation and hierarchical cluster analysis
Correlation between element concentrations and ash yields may provide preliminary insights into the elemental inorganic-organic affinities [77,90,91]. Hierarchical cluster analysis based on Pearson's correlation, which categorized elements in the Tashan thermally altered coal into four groups, was useful to analyze the geochemical associations and potential sources of the elements in the coal (Figure 13). Pearson's correlation coefficients between elements and ash yields, as well as inter-elements, are summarized in Table 9.   Table 9. Pearson's correlations between elements and ash yields and selected elements.  Figure 13). In this group, CaO, Na 2 O, Cr, and Ni were strongly correlated with the ash yield (r ash ≥ 0.6), suggesting their inorganic affinity (Table 9); MgO, SiO 2 , Al 2 O 3 , MnO, Fe 2 O 3 , and As had lower but still positive correlation coefficients with the ash yield (0.6 > r ash ≥ 0.3), suggesting that these elements are mainly associated with minerals in coal; TiO 2 , V, and Zn were weakly correlated with the ash yield (0.3 > r ash ≥ 0), suggesting that the modes of occurrence of these elements are complex. Cu, In, W, and Tl had negative correlation coefficients with the ash yield; however, these elements are not necessarily associated with organics. All of the group 1 elements had significant correlation coefficients with the aluminosilicates, and most of them were strongly correlated with carbonates, except V, Cu, W, and Tl. Therefore, it is deduced that the group 1 elements are mainly associated with aluminosilicates and carbonates in coal.

Correlation with Ash Yield
Group 2 included K 2 O, Rb, Co, and Cs ( Figure 13). K 2 O and Rb were strongly correlated with the ash yield (rash ≥ 0.6, Table 9), while Cs and Co had lower but still positive correlation coefficients (0.6 > r ash ≥ 0.3). The close correlation between K 2 O and Rb (r = 0.977) suggests that they both from illite or muscovite in the coal, which is consistent with Finkelman [92], who stated that substantial amounts of rubidium may be present in illite and mixed-layer clays, as rubidium can readily substitute for potassium. Finkelman et al. [93] suggested that Cs is geochemically associated with K and Rb. In the present study, the association of Cs and Co (r = 0.972) suggests that they are both from the same potassium-bearing minerals, which may be muscovite and illite.
Group 3 included P 2 O 5 , Sr, Ba, and Be ( Figure 13). Elements in this group were strongly correlated with the ash yield (r ash = 0.867-0.967, Table 9), except Be (r ash = 0.200). The high correlation coefficients between P 2 O 5 , Sr, and Ba suggest that these elements may be associated with phosphates in coal, such as goyazite, gorceixite, crandallite-group minerals, and apatite [39,93]. As for Be, it only exhibited significant correlations with Co (r = 0.616) and was weakly correlated with the ash yield. Previous studies have suggested organic and aluminosilicate associations of Be in coal [93,94]. It is thus inferred that beryllium in the Tashan coal may be associated with clays and organics; however, further studies are needed.
Group 4 included Ga, Sc, Cd, Nb, Ta, Mo, Pb, Hf, Sb, Se, Zr, U, Li, Bi, and Th ( Figure 13), which were negatively correlated with the ash yield (ranging from −0.394 to −1, Table 9). However, these elements are not necessarily associated with organic matter. Lead, Bi, Se, and Cd may be associated with sulfides and selenides; for example, Pb and Se can exist together as clausthalite [95]. These minerals may occur as fine-grained particles within or shielded by organics [93] and only account for a small proportion, which is below the XRD detection limit. Zirconium, Ta, Th, Hf, and U may be associated with accessory zircon. Lithium, however, may be associated with clays and accumulated at the edge of the alteration halo; therefore, it exhibits a reverse trend with the ash yields.

SEM-EDS Analysis
While statistical analysis provides useful information on the modes of occurrence in coal, it should be combined with other methods for reliable interpretations [39,45,91,96,97]. In the present study, the direct determination of some mineral compositions in the Tashan coals was accomplished by SEM-EDS.
The SEM-EDS analysis confirmed the modes of occurrence of some elements in the Tashan coals (Table 10, Figure 14). For example, Si and Al were mainly associated with clays, while quartz was another source of Si in coal. Iron, Ca, and Mg were mainly associated with ubiquitous carbonates in the thermally altered coals; however, a substantial amount of Fe was also detected in pyrite. In addition, Ca and Mg were also detected in anhydrite, and some Ca occurred in apatite; P was detected in goyazite and apatite; and Zr, Nb, Hf, and Ta were associated with zircon.

Strontium in the Tashan Coal
Concentrations of strontium in the Tashan coal were extremely high (peak at 3067 µg/g in sample T4) when compared with other coals. Despite the relatively low concentration in T5, the average concentration of Sr in the thermally altered coal was 1714 µg/g, which is more than 12 times that of Chinese coal (140 µg/g on average). Finkelman [92] suggested that the modes of occurrence of strontium can be both organic and inorganic; inorganic Sr may be associated with phosphates such as crandallite-group minerals and goyazite, barite, celestite, carbonates, clay minerals, and zeolites [17,92,96,98]. The quantification of the modes of occurrence of strontium by the sequential leaching procedure demonstrated that in bituminous coals, 50% of the Sr was associated with phosphates and carbonates, 25% was associated with silicates, and 25% was associated with organics [93]. In the present study, the high correlation coefficients between the elements in group 3 reflect that Sr is associated with phosphates, which was confirmed by the presence of goyazite by SEM ( Figure 7g). Moreover, SEM-EDS analysis demonstrated that Sr is also associated with clays, such as kaolinite ( Figure 15). A total of 22 clay minerals were analyzed by SEM-EDS, and Sr was detected in 13 of them, with a maximum concentration of 23.46 wt.% in kaolinite (Table 10). Although SEM-EDS analysis only provides semi-quantitative elemental compositions of the minerals in coal, it indicates that the abundant Sr-bearing clays resulted in high strontium values in the Tashan coal.

Element Sources
As discussed above, the interaction between rocks and coals has caused changes in the concentrations of elements, which may enrich, deplete, or have no impact in the contact zone. The enriched elements in the thermally altered coal have two potential sources: the elements either migrated out from the intrusion and enriched in coals and/or were enriched through residual accumulation during pyrolysis of the coal [99], and the contribution of each source differs depending on the mineral affinity. Elements that do not exhibit significant variations in the altered and unaltered zone may primarily have a coal origin; otherwise, they would be derived from the intrusion or driven off from the alteration halo [43]. Hierarchical cluster analysis categorizes elements that exhibit similar trends into the same group; these elements may have the same source ( Figure 16). The group 1 elements fluctuated in the altered and unaltered zone and increased sharply in the contact zone (Figure 16a), which indicates that intrusion is the source of these elements. It is generally suggested that minerals are the major hosts of the vast majority of elements present in coal [39], and thus the input of these elements from the intrusive body can be represented by epigenetic mineralization derived from hydrothermal fluids.
The intrusive rock itself may be altered by coal [37,44,53]. It is evident that the dike in the area adjacent to the coal became pale, whereas the color was dark green in the area far away from coal (Figure 3a). This phenomenon may provide more evidence for element migration. Song et al. [53] reported that mafic and ferrous minerals, such as pyroxene ((Ca,Na,Li) (Mg,Fe 2+ ,Fe 3+ ,Mn,CrAl) Si 2 O 6 ), olivine ((Mg,Fe) 2 SiO 4 ), and magnetite (FeFe 2 O 4 ), interact with CO 2 during the intrusive event. It is assumed that alteration of these minerals and subsequent elemental migration may cause the pale margins of the dike and the enrichment of group 1 elements in the thermally altered coal.
Concentrations of the group 2 elements remained relatively stable in the unaltered coal but increased in the contact zone (Figure 16 b), which may be due to an input of these elements from the intrusion.
The group 3 elements fluctuated greatly in the altered and unaltered coal. However, these elements were enriched in the contact zone (Figure 16c). The alteration of minerals in the intrusion, such as apatite (Ca 5 (PO 4 ) 3 (F,Cl,OH)), and the migration of these elements may play a part in the elevated group 3 concentrations; however, other factors may have also affected the concentrations of these elements. For example, it is possible that some of the group 3 elements contained in the roof rock or the intrusion were leached out by groundwater and introduced into the adjacent coal [37]. The trend of group 4 element exhibits a different pattern from the above elements; since these elements were not enriched or depleted with distance (Figure 16d), the concentration pattern may indicate that the coal itself is the source of these elements and that they responded neutrally to the intrusive event.

Conclusions
Petrographic, mineralogical, and geochemical compositions of the Tashan coal were investigated in this paper, and the following conclusions were drawn: (1) The 3.6 m diabase dike has caused an alteration halo of approximately 2 m. Within the area, coals exhibit an elevated vitrinite reflectance and ash yield; volatile matter is driven off; and thermally altered organic particles occur.
(2) More crystalline phases were identified in the thermally altered coal than in the unaltered coal. Epigenetic carbonates such as dolomite and ankerite may have been deposited by hydrothermal fluids. Muscovite and tobelite may have been transformed from a kaolinite precursor.
(3) Significantly enriched strontium in the Tashan coal is associated with phosphate minerals such as goyazite as well as the Sr-bearing clays.
(4) Elements associated with aluminosilicates, clays, and carbonates are enriched in the coal/rock contact zone; the enrichment was caused by epigenetic mineralization and element migration. On the other hand, P 2 O 5 , Sr, Ba, and Be may be controlled by multiple factors, such as hydrothermal fluids and groundwater, while elements exhibiting no variations in thermally altered and unaltered coals indicate that the coal itself is the source. Funding: No funding was received to assist with the preparation of this manuscript. Data Availability Statement: All data generated or analyzed during this study are included in this published article.