Introduction

The so-called “geochemical effect of lanthanides” phenomenon was established by the authors 10 years ago. This effect is a geochemical reflection of lanthanide contraction in the periodic table of elements, i.e., any quantifiable geochemical behaviors of lanthanides (for example: the normalized values of chondrite for lanthanides as well as the parameters characterizing the migration, distribution, and occurrence modes of lanthanides in geological bodies can all be considered as the results of quantifiable geochemical behaviors) usually have a strong linear relationship with certain atomic structural parameters (such as the radii of their trivalent ions)1,2,3,4,5,6,7,8. To investigate this relationship, lanthanides can be divided either into two parts (LREE and HREE), or take the 14 elements of lanthanides can be evaluated comprehensively. This relationship is described either by a linear regression function or through other functions, such as quadratic polynomial function, log function, exponential function, or power function. Usually, there are significant differences in the manifestations of the geochemical effect of lanthanides among geological bodies of different genetic types5,6. These differences likely result from the bodies’ unique geochemical states and processes, and so useful geological information may be acquired from them.

Within the currently available literature, there is no corresponding background information about this phenomenon. A few textbooks have mentioned that the complexation potential of the lanthanide elements enhances with an increase in their atomic numbers. However, these publications do not thoroughly probe into the atomic structural parameters of the lanthanide elements, such as the radii of their trivalent ions9.

Since the discovery of the geochemical effect of lanthanides, this research topic has received considerable attention from the authors, who have focused on the scientific significance and applications of this effect. In a recent article, we proposed a geochemical effect equation for lanthanides and explored the geochemical significance of certain parameters in all types of function regression equations8. However, being limited by the conditions at the time, we were unsuccessful in recognizing the geological significance of goodness of fit, which is a highly essential parameter within these regression equations. Concurrently, the explanation for the scientific significance of certain evaluated parameters was not sufficiently thorough. In this context, we took the late Paleozoic C2 coal seam of Handan Coalfield in Hebei Province of China as an ideal example and further evaluated the geological and geochemical significance of all parameters in the geochemical effect regression equation for lanthanides.

Geological setting and previous research

The Fengfeng Mine of Handan Coalfield (Hebei Province, China) was selected as the study area (Fig. 1), which contains a series of late Paleozoic coal-bearing strata (Fig. 2). Owing to the influence of Yanshanian magmatism on the C2 coal seam of the Shanxi Formation in the Early Permian, a series of thermal metamorphisms occurred almost continuously from south to north in the same coal seam; as a result, the coal gradually changed from low-rank bituminous to high-rank meta-anthracite (Fig. 1). To understand the influence of magmatism and coal's thermal metamorphism on the migration and distribution of elements in coal, systematic sampling and analyses have been conducted on coals with different ranks of metamorphism (Rmax from 0.89% to 7.41%) in the C2 coal seam. The important parameters, including the concentration of lanthanides are shown in Table 1. Sampling and analysis methods were exhibited in the literature10. The previous work in this field has established a useful foundation for our study.

Figure 1
figure 1

Location of the C2 coal seam and sampling ___location in the Fengfeng mine of the Handan Coalfield, Hebei Province (Dai and Ren, 200710).

Full size image
Figure 2
figure 2

Profile of the sedimentary sequence of coal-bearing strata in the Fengfeng mine of the Handan Coalfield, Hebei Province (Dai and Ren, 200710).

Full size image
Table 1 Metamorphism, lanthanide abundance and related parameters of Late Paleozoic C2 coal seam in Fengfeng mine, China (Dai and Ren, 200710).
Full size table

Research ideas and methods applied in this study

After numbering 12 coal samples according to the ranks of the C2 coal seam from low to high as No. 1–No. 12, their lanthanides were normalized with chondrite, and the results of normalization and the radii of their trivalent ions were subjected to a regression analysis with a linear function, quadratic polynomial function, log function, exponential function, and power function. The key obtained parameters are shown in Table 2. On this basis, the relationships between these parameters and other known geochemical parameters were utilized in investigating the geological and geochemical significance of these parameters.

Table 2 Important parameters obtained by regression of normalized lanthanide chondrite and the radius of their trivalent ion using various functions in the coals from Late Paleozoic C2 coal seam in Fengfeng mine, China.
Full size table

Results and discussion

Correlation of parameters in the regression equation and parameters of REE

We took the relationship of lanthanides normalized with chondrite and the radii of their trivalent ions as an example and confirmed the previous findings. Specifically, the regression coefficients of each function (i.e., slope in the case of linear regression) from this study had significantly positive correlations with ∑REE and ∑LREE, and the power exponent of the power function has a significantly linear positive correlation with LREE/HREE8. In addition, this study also identified significant positive correlations between the intercept values obtained from log function regression and ∑REE, with a goodness of fit value of above 0.99 (Fig. 3a) and a value of 1 for some coal samples. The regression coefficients of the log function also showed significant positive correlations with the LREE concentration (with a goodness of fit of 0.9228, Fig. 3b). On this basis, a new regression equation was proposed:

$$F = k_{1} \sum LREE \, \ln \, r \, + \, k_{2} \sum REE,$$
(1)

where F is one of the geochemical behaviors of lanthanides, r is the radius of trivalent ion for lanthanides, and k1 and k2 are both proportionality coefficients. In the case of normalization of lanthanides with chondrite (the normalization process with chondrite was also taken as a geochemical process of the lanthanides, and more specifically their differentiation results relative to chondrite), the value of k2 can be approximated to 1.

Figure 3
figure 3

(a) Correlation between the intercept values obtained from log function regression for the geochemical effect of lanthanides in the C2 coal seam in Fengfeng Mine of Handan Coalfield in Hebei Province and the total concentration of REE, and (b) the correlation between the regression coefficients and concentration of LREE. The original data of REE refer to Table 1 (Dai Shifeng, 200710), the regression coefficients refer to Table 2.

Full size image

In regression functions, few regression parameters have a correlation with ∑HREE. In this study, the regression coefficients of the power function for lanthanides in C2 coal seam of Handan and the radii of their trivalent ions had a significant positive correlation with ∑HREE (with a goodness of fit of 0.7033). Therefore, revising the equation as the following, produces a formula that is better than its counterpart in the literature8 (F = − k1∑REE rk2 f).

$$F \, = \, - K(k_{0}^{\prime } \sum LREE \, + \, k_{1}^{\prime } \sum HREE)r^{k2\; \, f} ,$$
(2)

where r represents the radii of trivalent ions of lanthanides, f is the differentiation degree of light (LREE) and heavy rare earth (HREE) elements, i.e., LREE/HREE, k0, k1, and k2 are all proportionality coefficients, and \(k_{0}^{\prime }\) and \(k_{1}^{\prime }\) are proportionality coefficients produced after extracting the common factor K.

Three parameter types in the geochemical effect regression equation for lanthanides

Certain parameters listed in the above formulas, including regression coefficients, intercept values, and power exponents in the regression equation, are referred to as “parameters in common” in this paper, as the relationships of these parameters with ∑REE, ∑LREE, ∑HREE, and LREE/HREE are universal. Their geochemical significance conforms to the geochemical effect of lanthanides in coal and fits within geological bodies other than coal, such as the geochemical effect of lanthanides in magmatic rock (to be discussed later). Regardless of whether these samples were obtained from the same coal seam (or cognate magmatic rock, metamorphic rock or sedimentary rock, and so on) or different coal seams (or different sources and different kinds of magmatic rock, metamorphic rock, or sedimentary rock, and so on), the geochemical parameter characteristics of lanthanides (i.e., ∑REE, LREE/HREE, and so on) can always be described with the above parameters. Thus, the rule is “universally applicable”. This suggests that any geochemical behaviors or processes of lanthanides are limited by grand heat-tectonic events in geology and are strictly restrained by the internal microstructures of their ions in statistics. Moreover, these behaviors or processes are always closely related to the parameters of lanthanides, such as ∑REE, ∑LREE, ∑HREE, LREE/HREE, and so on. The geochemical behavior of every single lanthanide has a strict functional relationship with the radius of its trivalent ion, the total quantities of ∑REE and ∑LREE, and the differentiation degree of heavy and light rare earth elements. In other words, these are not purely individual behaviors. However, it is entirely possible that the values and ± of the proportionality coefficients in the above formulas (such as K, \(k_{0}^{\prime }\), \(k_{1}^{\prime }\), k1, k2, and so on) might vary with the different geochemical processes experienced by lanthanides, their different occurrence modes, their different geochemical environments, and so forth. Thus, they can be regarded as the “individual parameters” or “state parameters” in the above regression equations. For example, the complexation ability of lanthanides increases constantly with the reduction of the ionic radius. In this process, k1 in formula (1) should be negative. In contrast, in the fractional crystallization process of magma, the relationships of the LREE and HREE after normalization with chondrite and the radii of their trivalent ions usually presented a significant double-pitch-up pattern, i.e., the K values of the regression coefficients in these linear relationships should all be positive6,9. The latter research will show that the “deviation value” of lanthanides is also an “individual parameter”.

In this study, the linear regression coefficients in the relationship of the LREE after normalization with chondrite and the radius of its trivalent ion (hereinafter called “regression coefficients of geochemical effect equation for chondrite-linear lanthanides”) had significant or even highly significant positive correlations with concentration of Al, Si, La, Ce, and Pr (Table 3, with goodness of fit values of 0.4919, 0.5062, 0.9838, 0.8617, and 0.5559, respectively). Yb/La implies that the presence of multistage hydrothermal activity may also lead to the reallocation of trace elements, including rare earth element of coal-bearing rocks between coal seams and partings. Compared to partings, the high elemental ratios in coal seams (Yb/La, Nb/Ta, and Zr/Hf) are mainly caused by the re-precipitation of these elements (such as Yb, Nb, and Zr) in coal seams. In partings, these elements (Yb, Nb, and Zr) are more active in the leaching process and are easily leached, adsorbed by the underlying organic matter, and then precipitated in the underlying coal seam11. In this study, the regression coefficients in the exponential function of REE have a strong logarithmic relationship with Yb/La (Fig. 4a), whereas the coefficient of independent variable X in the exponential function shows a significant negative correlation with Yb/La (Fig. 4b). The power exponent obtained from the power function regression presents a similar relationship as shown in Fig. 4b. In addition, it indicates that the “parameters in common” such as certain functional regression coefficients, the coefficient of independent variable X (note that X refers to the radius of trivalent ions of lanthanides, which is regarded as a variable according to the regression formula) in the exponential function and power exponent of the geochemical effect of lanthanides, can also reflect geochemical processes under certain conditions; they also display certain characteristics of “state parameters”. This is because the ratios of certain elements inside lanthanides themselves have an environmental significance. Therefore, parameters such as the regression coefficients of exponential function for REE share the characteristics of both “parameters in common” and “individual parameters”, and these kinds of parameters can thus be called “two-sided parameters”.

Table 3 Ash yield (%, db), major and minor (wt %, db), and trace elements (μg/g, db) concentrations in the coals from the Fengfeng–Handan Coalfield (Dai and Ren, 200710) and related ratios.
Full size table
Figure 4
figure 4

(a) Correlation between regression coefficients in the exponential function of geochemical effect for chondrite–lanthanides and Yb/La, and (b) the correlation between the coefficient of independent variable X (i.e., the coefficient of the radius of trivalent ion for lanthanides) and Yb/La in the C2 coal seam. The original data of REE refer to Table 1 (Dai Shifeng, 200710), the regression coefficients refer to Table 2.

Full size image

The goodness of fit in the above regression equation, which is also an “individual parameter”, likely depends upon the sources, geological origin, geological environment, geochemical process, occurrence mode, and other factors related to lanthanides. For example, the goodness of fit of a geochemical effect equation for chondrite-linear lanthanides of REE in the C2 coal seam shows a negative correlation with its LREE/HREE (with a goodness of fit of 0.7191). However, in other cases, this relationship may not hold true; one such instance was for the 616 basalt samples obtained from a sea bed12 (the detailed description is giving in another paper). This relationship cannot be established if lanthanides are normalized with a depleted mantle or enriched mantle, which represents a different geochemical process. However, in the case of No. 8 coal seam of Taiyuan Formation in Xishan of Taiyuan, Shanxi, China in the literature13, there is indeed such a relationship, but it is positive. Regardless of whether the 12 coal samples in C2 coal seam of Handan, or the five samples in No. 8 coal seam in Xishan of Taiyuan, Shanxi in the literature are used 13, all originate from the same coal seam of different metamorphic degrees. This likely implies that the linear goodness of fit is closely related to the differentiation degree of heavy and light rare earth elements only when the lanthanides are from same source and have undergone different evolutionary or metamorphic phases.

This study also found that the goodness of fit of the geochemical effect equation for chondrite-linear lanthanides of LREE also showed negative correlations with the concentrations of Al, Si, and La (with goodness of fit values of 0.4098, 0.5963, and 0.5328, respectively); that of REE presented similar correlations with the concentrations of Al, Si, and La (with goodness of fit values of 0.4084, 0.634, and 0.5328, respectively) and was positively correlated with the concentration of Zr (with a goodness of fit of 0.4147, which can reach 0.6858 after removing sample No. 5). The goodness of fit of the geochemical effect equation for chondrite-quadratic polynomial lanthanides of REE also showed negative correlations with the concentrations of Al, Si, and La (with goodness of fit values of 0.4296, 0.6369, and 0.5087, respectively), and exhibited a positive correlation with the concentration of U when the power function was used for fitting (with a goodness of fit of 0.6876). In general, in many cases, the goodness of fit of function regression usually shows a positive or negative correlation with several elements of lanthanides. This suggests that, in this example, some geochemical states (geochemical processes, occurrence modes, or origin evolution types) of lanthanides are related to the activity rules of some lanthanide elements as well as those of non-lanthanide elements especially Al and Si.

“Deviation value” of lanthanides and its geochemical significance

Both previous studies and this present study indicate that, in some cases, geochemical behaviors of lanthanides in certain geological bodies should be more thoroughly described using a linear geochemical effect equation for lanthanides. In other cases, non-linear functions such as the log function, exponential function, power function, and quadratic polynomials are more appropriate (namely, the goodness of fit of the regression function is higher). At present, it is still difficult to offer a detailed and accurate explanation on this phenomenon. However, we believe that it might have been caused by different “geological or geochemical states”. In theory, the higher the goodness of fit for a functional relationship of the geochemical effect of lanthanides, the more singular the geochemical process, or “geological or geochemical states” such as source, origin, occurrence mode of the corresponding lanthanides. Certainly, in this study, it is only a scientific speculation or inference. Regarding the specific correspondence between the functional regression relationship and the geochemical state of lanthanides, studies and statistical analysis are still needed for further clarification. Based on the above speculation and inference and taking the geochemical effect equation for chondrite-linear lanthanides as an example, the investigation was conducted in two parts, i.e., LREE and HREE. The most ideal circumstance for the geochemical processes or occurrence modes characterized is a goodness of fit of 1. In general, this goodness of fit cannot be 1, because the geochemical behaviors or states of lanthanides, regardless of their distribution, migration process, or occurrence modes, are not purely singular. Therefore, the geochemical behavioral parameters of some lanthanides obtained from a regression formula with a goodness of fit 1 are only the “ideal value”. In most circumstances, the actual values of parameters for some geochemical behaviors of lanthanides do not match ideal values. This “deviation value” can be obtained from the difference between the parameter value obtained from the actual goodness of fit and the ideal value obtained by assuming the goodness of fit of 1. Table 4 shows the “deviation values” calculated from the geochemical effect equation for chondrite-linear lanthanides in 12 samples of the C2 coal seam in Handan. Tables 1 and 4 show that the deviation values of some elements for LREE are related to the anomaly of Ce. For example, the deviation value of Ce shows a negative correlation with δCe (goodness of fit = 0.6327), and positive correlations with ∑REE and ∑LREE (goodness of fit = 0.6571 and 0.6724, respectively). The same is true for the deviation value of Sm (the goodness of fit values for ∑REE and ∑LREE = 0.7915 and 0.737, respectively), which also has a positive correlation with ∑HREE (goodness of fit = 0.5898). The deviation value of Ho shows a negative correlation with ∑HREE (goodness of fit b = 0.44), while the deviation value of Yb exhibits a positive correlation with it (the goodness of fit being 0.4909). Thus, the deviation values of some elements for LREE can reflect their geological environments (such as the oxidation–reduction environments) and that the magnitudes of such deviation values are restrained overall by ∑REE, ∑LREE, or even ∑HREE. The deviation values of La and Eu also show positive correlations with LREE/HREE (goodness of fit = 0.614 and 0.8127, respectively), whereas those of Pr and Nd show negative correlations with LREE/HREE (goodness of fit = 0.7703 and 0.4152, respectively). This suggests that the differentiation degree of heavy and light rare earth is also one of the reasons causing the “deviation” of lanthanides.

Table 4 “Deviation value” of lanthanide elements and related parameters in the coals from Late Paleozoic C2 coal seam in Fengfeng mine, China (μg/g).
Full size table

In general, the ratios of many elements can indicate sedimentary environments and geological origins in coal geology. For instance, a smaller ratio of Th/U in coal generally reflects a stronger hydrothermal influence on a coal seam14. High ratios of Yb/La, Nb/Ta, and Zr/Hf are mainly caused by the re-precipitation of Yb, Nb and Zr in the coal seam. In partings, these elements are active in the leaching process and can be easily leached, adsorbed by the underlying organic matter, and precipitated in the underlying coal seam11. The ratio of Nb/Ta can be used as the indicator to identify different evolutionary processes of magma9,15. The ratio of Sr/Ba increases with the distance away from the coast and can qualitatively reflect paleosalinity; specifically, the ratio of Sr/Ba < 1 in freshwater sediments, > 1 in marine sediments, and between 0.6 and 1 in blackish water facies16. The ash index (SiO2 + Al2O3)/(MgO + CaO + Fe2O3) mainly reflects the degree of reduction for a peat bog. When it is high, the reduction of media is weak in peat precipitation, and the aqueous nature of freshwater is evident; when it is low, the reduction of media is strong in peat precipitation17.

Based on the covariant figure of “deviation values” of lanthanides and some of the above elements' ratios, the 12 continuously metamorphous coal samples in the C2 coal seam of Handan can be divided into two different series. For instance, based on the covariant figure of the deviation value of the La-ash index in Fig. 5, these continuously metamorphous coal samples can be divided into series A (No. 1–No. 7), and series B (No. 8–No. 12). The literature10 indicates that samples No. 1–No. 7, from the coal seam are uninfluenced or slightly influenced by magmatic-hydrothermal fluids in C2, whereas samples No. 8–No. 12, from the coal seam are strongly influenced by magmatic-hydrothermal fluids in C2. The deviation values of La or Lu vs. the ash index can also constitute a similar covariant figure. Moreover, based on the covariant figure of the goodness of fit of chondrite-linear lanthanides and the ash index, these continuously metamorphous coal samples can also be divided into series A and series B. This suggests that the “deviation value” and goodness of fit are both useful “individual parameters” or “states parameters” that can indicate indicating their geochemical environments. Some “parameters in common”, such as the mentioned slope value (regression coefficient) of chondrite-linear lanthanides, power value after power function regression, and exponential value after exponent function regression, can all be used to distinguish the two series of A and B in the C2 coal seam in the covariant figure plotted from the “parameters in common” and ash index. Therefore, these “parameters in common” in fact fall within the scope of “two-sided parameters” described above. Other “parameters in common” observed by the predecessors, such as δEu, can also constitute a good covariant figure capable of identifying geological environments (Fig. 6). In Fig. 6, the samples of A series in the C2 coal seam, being uninfluenced or slightly influenced by magmatic-hydrothermal fluids are basically above the regression line, whereas those of B series in the C2 coal seam strongly influenced by magmatic-hydrothermal fluids are all below the regression line. Even the covariant figure constituted by ∑REE, ∑LREE, LREE/HREE, and the ash index can be used to distinguish A and B in the C2 coal seam. This result can also be inferred from the formulas (1) and (2).

Figure 5
figure 5

Covariant figure of deviation value of La–ash index in the C2 coal seam. The deviation value of La refers to Table 4, and the ash index refers to Table 3.

Full size image
Figure 6
figure 6

Covariant figure of δEu–ash index in the C2 coal seam. The value of δEu refers to Table 1, and the ash index refers to Table 3.

Full size image

The sixteen elements of B, F, Cl, Br, Hg, As, Co, Cu, Ni, Pb, Sr, Mg, Ca, Mn, Zn, and U in Table 3 and the ratios with certain geological or geochemical environmental significance (such as Sr/Ba, Th/U, and Th/Co), as well as other element ratios (such as Nb/Ta), can constitute a covariant figure capable of indicating a geological environment or origin through the “deviation values” of lanthanides or ash indexes. The above elements can be found within the hydrothermal fluid of magma10. For instance, the covariant figure of Ca%—ash index can reflect a strong logarithmic negative correlation between the indexes but also clearly distinguish between A and B of the C2 coal seam (Fig. 7). Certain covariant figures can also reflect slight differences in environmental influence. For example, although the covariant figure of the Lu deviation value-Co can clearly distinguish between A and B, it fails to classify sample No. 8, which is also strongly influenced by the hydrothermal fluids of magma, into Series B. According to a previous study (Figures 10 and 11, see reference)10, although some transition elements such as Co are found in the hydrothermal fluids of magma, in the phase of sample No. 8, the hydrothermalism is not strong enough. Herein, compared to the samples prior to No. 7, which are uninfluenced or slightly influenced by hydrothermal fluids of magma, there are few changes in elements, such as Co. Thus, sample No. 8 is not separated from Series A (Fig. 8).

Figure 7
figure 7

Covariant figure of Ca–ash index in the C2 coal seam. The concentration of Ca and the ash index refer to Table 3.

Full size image
Figure 8
figure 8

Covariant figure of the deviation value of Lu–Co in the C2 coal seam. The deviation value of Lu refers to Table 4, and the concentration of Co refers to Table 3.

Full size image

Environmental indication significance of parameter Rmax

Rmax, an important parameter used to characterize the metamorphic degrees of coal in the Handan C2 coal seam in Hebei, has no evident correlation with any parameter of the geochemical effect equation for lanthanides. However, the covariant figure constituted by Rmax with other parameters still has certain environmental significance. For example, the covariant figure of the Rmax—ash index reflects the significant negative correlation between the Rmax and ash index while clearly distinguishing between A and B for the C2 coal seam (Fig. 9). In turn, this completely indicates several super-metamorphic coal forms in the C2 coal seam with the superposition of the magmatic-hydrothermal process in Yanshanian based on plutonic metamorphism. For series A of samples No. 1–No. 7, the ash index reflects the changes in the reduction degree of peat bog upon the formation of C2 coal seam17. However, series B of samples No. 8–No. 12, it is more likely to reflect the result of magmatic-hydrothermal influence on the coal seam. Specifically, the large-scale injection of Ca and Mg in the hydrotherm causes the ash index of the coal seam to present an overall declining trend. The same applies to the covariant figures constituted by the above elements or elements' ratios related to magma hydrothermal fluid with the ash index or some parameters of the geochemical effect of lanthanides. Similarly, the covariant figure of Rmax-Th/U and Rmax-Th/Co can reflect the power function negative correlation between them and clearly distinguish between A and B for the C2 coal seam. Moreover, although Rmax has no correlation with Sr/Ba, it can clearly distinguish between these two series as well.

Figure 9
figure 9

Covariant figure of Rmax–ash index in the C2 coal seam. The value of Rmax refers to Table 1, and the ash index refers to Table 3.

Full size image

Conclusions

  1. 1.

    The geochemical significance of the parameters in the geochemical effect regression equation for lanthanides was systematically explained, and new descriptive formulas characterizing the geochemical behaviors of lanthanides were proposed. We speculate that different geochemical processes of lanthanides can be described with different regression functions; the higher the goodness of fit, the purer and more singular the different geochemical behaviors or processes of lanthanides will be.

  2. 2.

    Based on the above speculation, the concept of lanthanide deviation was proposed, and the parameters in the geochemical effect regression equation for lanthanides were roughly divided into three types: The first type is parameters acting in common, and its geochemical significance fits within any geochemical process experienced by any geological body other than coal. The second type is individual or state parameters, which are probably only applicable to certain geological bodies or geochemical environment that have stronger geochemical, environmental, and genetic significances. However, the difference between individual parameters and parameters in common is not “either this or that”. In this study, some parameters acted in common, such as the regression coefficients of some functions. Moreover, they had environmental indication significance under certain conditions. Therefore, the concept of “two-sided parameters” was established.

  3. 3.

    In this study, lanthanides had certain interactions with the elements from magmatic fluid or certain non-magmatic-hydrothermal-sourced elements such as Al and Si. Although the sources and changes of ∑REE and ∑LREE are not related with the magmatic-hydrothermal effect, they have certain interactions with the elements from the magmatic hydrothermal fluid. This manifests that the covariant figure with these elements can also be used to indicate different metamorphic series in the C2 coal seam. Moreover, this covariant figure conforms to the scientific logic inside regression formulas proposed in this study.

  4. 4.

    Rmax in the C2 coal seam of Handan reflected the metamorphic degree of coal and indicated the influence of post-modification activities to coal through the combination with other indexes. For example, in the covariant figure of Rmax and Sr/Ba or the ash index, the degree of the influence of later magmatic-hydrothermal fluid on the C2 coal seam was clearly reflected.