Abstract

The mineral chemistry of twenty chlorite samples from the United States Geological Survey (USGS) spectral library and two other regions, having a wide range of Fe and Mg contents and relatively constant Al and Si contents, was studied via infrared (IR) spectroscopy, near-infrared (NIR) spectroscopy, and X-ray fluorescence (XRF) analysis. Five absorption features of the twenty samples near 4525, 4440, 4361, 4270, and 4182 cm−1 were observed, and two diagnostic features at 4440 and 4280 cm−1 were recognized. Assignments of the two diagnostic features were made for two combination bands ( and ) by regression with IR fundamental absorptions. Furthermore, the determinant factors of the NIR band position were found by comparing the band positions with relative components. The results showed that Fe/(Fe + Mg) values are negatively correlated with the two NIR combination bands. The findings provide an interpretation of the NIR band formation and demonstrate a simple way to use NIR spectroscopy to discriminate between chlorites with different components. More importantly, spectroscopic detection of mineral chemical variations in chlorites provides geologists with a tool with which to collect information on hydrothermal alteration zones from hyperspectral-resolution remote sensing data.

1. Introduction

Chlorites are ubiquitous ferromagnesian phyllosilicates most commonly found in epimetamorphic rocks as hydrothermal alteration products and after erosion in sediments together with various clay minerals [1, 2]. The crystal structure of the chlorite group minerals can be described as a 2 : 1-type hydrous aluminosilicate (talc-like layer) with the octahedral sheet “sandwiched” between two opposite tetrahedral sheets and linked by an extra octahedral sheet (brucite-like layer). The simplified structural formula of chlorite minerals can be described as Mg6Mg6(Si8O20)(OH)4(OH)16, in which sheets of talc (Mg6(Si8O20)(OH)4) and brucite (Mg6(OH)12) are included (Figure 1). The brucite-like and talc-like sheets are bonded to one another by long hydrogen bonds between the oxygen atoms from the siloxane sheet of the talc-like layer and the hydroxy groups of the brucite-like layer. Cation substitution is very common in chlorites and leads to a wide range of chemical compositions. In some cases, partial substitution of the Mg2+ by Fe3+ occurs in the brucite-like layers accompanied by a coupled charge compensation, whereas substitution of Si4+ by Al3+ occurs in the tetrahedral sheet of the talc-like layer. These substitutions produce a variety of clinochlorites with a structural formula of Mg4Al2Mg6(Si6Al2O20)(OH)4(OH)12. The presence of Al in both the octahedral and tetrahedral sheets is necessary to ensure similar crystal cell parameters for the formation of a stable structure. Mg2+ can also be replaced by Fe2+, leading to a type of Fe-rich chlorite known as ripidolite with a theoretical formula of (Fe2+,Mg)3Al3(Fe2+,Mg)6(Si5Al3O20)(OH)4(OH)16 [3].

Most of the available vibrational spectral studies have focused on the infrared, Raman, and other vibrational spectra of chlorites, and some correlations between spectral features and chemical composition have been published. Tuddenham and Lyon [4] observed the relationship between the substitution amount of Al for Si in the tetrahedral positions and the wavenumber of the Si-O stretching. In addition, Stubican and Roy [5] stated that the substitution of Al for Si was closely related to the position of the strongest Si-O band in the 665 to 685 cm−1 region. Hayashi and Oinuma [6] reported a band shift from 540 to 560 cm−1 with an increasing octahedral Al content but an opposite shift with increasing Mg and Fe contents. Moreover, Hayashi and Oinuma observed a shift toward lower frequencies of the band between 620 and 692 cm−1 with increasing octahedral Mg and Fe and decreasing octahedral Al content. Additionally, they observed another shift toward lower frequency in the OH stretching bands from 3400 to 3436 cm−1 and from 3560 to 3586 cm−1 with increasing Fe content. These shifts in OH stretching bands were interpreted as being due to a shorter -OH interlayer distance in Fe-rich chlorites. The two broad bands near 3560 and 3420 cm−1 were ascribed to the interlayer OH with a weak shoulder band around 3620 cm−1 associated with the inner OH of the 2 : 1 sheet [7]. The two relatively strong bands were generally assigned to the (SiAl)O-OH and (AlAl)O-OH vibrations [8], and the intensities were mainly determined by the tetrahedral sheet composition, while the exact band positions were also related to the composition of the octahedral interlayer hydroxide sheet.

The near-infrared (NIR) spectra of chlorites have also been studied by some researchers, but they mostly focused on the application of NIR spectra to discriminate chlorites from other minerals. Post and Crawford [9] observed a shift in the overtone band from 7102 cm−1 to 7205 cm−1 and suggested that this shift might make the identification of chlorite by remote sensing more difficult. Yang et al. [10] and Laakso et al. [11] observed a hydroxyl absorption shift (4424 cm−1 to 4456 cm−1) toward higher frequency with increasing Mg and toward lower frequency with increasing Fe, but they did not discuss the hydroxyl absorption band near 4273 cm−1. Petit et al. [12] assigned the NIR features of octahedral smectites (2 : 1 layer) near 4550 cm−1 and 4370 cm−1 to the combination of hydroxyl stretching and bending modes. As with the layer minerals of 2 : 1, the NIR features of chlorites may also be related to the combination of hydroxyl stretching and bending in the octahedral sheet.

IR studies and NIR studies are historically distinct from one another. The former have mainly focused on functional groups and spectral interpretation, whereas the latter have focused on statistical analysis. The use of spectral units of cm−1 (wavenumber) in IR and nm (wavelength) in NIR is symptomatic of the historical divide separating these research communities.

The aim of this paper is to compare NIR data with IR data with respect to the Fe/(Fe + Mg) values and identify changes that occur in the NIR spectra as a function of the IR features and cation substitutions. The NIR reflectance technique, which is particularly fast and efficient for identifying both minerals bearing hydroxyl moiety and carbonate minerals, would contribute to detecting the gossans and hydrothermal products on the earth’s surface when mapping the land using multi- or hyperspectral remote sensing [13].

2. Materials and Methods

2.1. Samples

Mineral component and IR and NIR data for eight samples were collected from the USGS spectral library. The descriptions and conditions of these chlorites have been published (see references in Table 1), and only the prominent features will be summarized and discussed here. The IR features were used to assign the combination bands of chlorite in the NIR region. Another twelve samples were collected from the two regions (Nachtai and Arjin regions) in Northwestern China (Table 2). Samples were ground, and chlorite minerals were purified. NIR and XRF were used to collect spectra and determine the mineral components. The mineral component data of these twenty samples were used to elucidate the relationship between spectral features and certain components by comparing the IR and NIR band positions.

2.2. NIR Spectroscopy

All NIR spectra were obtained using a PANalytical ASD FieldSpec Pro® 3 spectrometer (hereafter referred as ASD) that records spectra from the 350 to 2500 nm wavelength (4000 to 28,571 cm−1) region with a spectral resolution of 10 nm and a sampling interval of 1 nm in the shortwave infrared (1300–2500 nm) region. The spectrometer was connected to a contact probe with an internal halogen bulb, which ensures stable illumination conditions during data collection. The raw values of at-sensor radiance were converted to surface reflectance values using a Spectralon™ reflectance panel (i.e., the “white reference,” SRT-99-100, Labsphere Inc., North Sutton, New Hampshire), which is a commercially available plate made of polytetrafluoroethylene [16]. Finally, these relative reflectance values were converted to absolute reflectance values by multiplying the relative reflectance value for each wavelength with the reflectance factor obtained from the calibration certificate of the Spectralon panel, in accordance with the procedure of Clark et al. [17]. The wavelengths of the spectra were converted to wavenumbers to facilitate the comparison with the IR data.

2.3. XRF Analysis

Twelve samples from the two regions were ground and concentrated and then analyzed with a PANalytical Axios Wavelength-Dispersive X-ray Fluorescence Spectrometer (WDXRF) to determine the oxide component. This method can determine the contents of Al2O3, SiO, FeO, and MgO (Table 2). The XRF spectroscopy instrument is equipped with an array of six analyzing crystals and fitted with a 4 kW rhodium target. A vacuum value between 5 and 100 Pa was used as the medium for the analyses to avoid interactions of X-rays with air particles. Ground samples with Li2B4O7 were placed into yellow platinum crucibles. Then, the samples were melted for 10 min in the smelter and refrigerated for XRF measurement.

2.4. Spectral Component Analysis

Spectral manipulation, which included baseline adjustment, smoothing, and normalization, was processed using the program written by VBA in Microsoft Excel software. The absorption features of chlorites located between 2250 and 2340 nm (4440 and 4270 cm−1) can be explicitly separated into their constituent absorption bands using the Gaussian-Lorentz model [18], and the mathematical shape description of the individual absorptions is accurate [1924]. The Gaussian-Lorentz model was developed and validated by empirical studies of isolated vibration absorption bands in both transmission and reflectance spectra of autunite, nontronite, and smectite [25, 26]. The Gaussian-Lorentz model fitting was undertaken until reproducible results were obtained with an R2 of greater than 0.995.

3. Results

3.1. MIR Bands

The MIR spectra collected from the USGS spectral library are dominated by Si-O absorption bands at approximately 1070 and 1030 cm−1, and these features are not related to NIR spectra [27]; they are therefore not mentioned in the following discussion.

3.1.1. Bands in the 3664–3406 cm−1 Spectral Range

The IR spectra ranging from 3664 to 3406 cm−1 are shown in Figure 2, and Table 3 shows comparisons with some chlorite minerals reported in the literature in order to assist with the assignment of the observed bands for chlorite spectra. A relatively weak shoulder band was observed in the IR spectra near 3662 cm−1, which also occurred on the IR spectra of talc and was attributed to the stretching of OH groups from the 2 : 1 layer [28, 29]. Two additional wide absorptions which occurred from approximately 3554 to 3423 cm−1 were very strong and overlapping in large widths at half height. IR spectra showed an increase in the frequency of these bands for the Fe-Mg substitution at constant Al/(Al + Si) and a stable frequency for Al-Si substitution at constant Fe/(Fe + Mg). The bands that occurred near 3554 cm−1 were assigned to the hydroxyl stretch in (AlAl)O-OH (), and the bands that occurred near 3423 cm−1 were assigned to (SiAl)O-OH stretching () [7, 30]. Furthermore, Prieto et al. [30] observed that an increase of Fe content induced a shift of these two bands toward lower frequency. This relationship was also validated using USGS samples, as shown in Figure 3.

3.1.2. Bands in the 1033–428 cm−1 Spectral Range

The IR spectra of the eight samples collected from the USGS spectral library displayed three intense bands at 1060, 1030, and 956 cm−1 (Figure 4 and Table 3). These bands were assigned to Si-O stretching in 2 : 1 layers [13, 30] and were also observed in other phyllosilicates with 2 : 1 layers [31, 32]. The second type of absorption in this region was located at approximately 850, 757, and 649 cm−1. These bands were assigned to (AlAl)O-OH, (SiAl)O-OH, and (SiSi)O-OH bending absorptions (, , and ), respectively [13, 33]. The frequency of 850 cm−1 bands shifted toward higher frequency with an increasing Fe content. In the 400–550 cm−1 region, the IR spectra displayed three bands near 555, 472, and 432 cm−1. These bands were assigned to (Fe,Mg)-O-Si or Al-O-Si bending [30].

3.2. NIR Bands

For many natural phyllosilicates, the combination bands occurring in the NIR region are broad and overlapping (Figure 5). In the spectral component analysis, four or five independent bands were observed at approximately 4525, 4440, 4361, 4270, and 4182 cm−1 (2210, 2252, 2293, 2341, and 2391 nm), which are the average band positions of all twenty samples. The two typical samples (HL-01 and HL-13), including 4 and 5 independent absorption bands, are shown in Figure 6. The weak band near 4525 cm−1 was assigned to the combination of OH stretching and deformation vibrations of Al-OH groups in 2 : 1 layers, and this band was also observed at a similar frequency in kaolinites [34]. Another weak band near 4361 cm−1 was tentatively attributed to an Mg-OH combination band [10]. The absorptions near 4440 and 4270 cm−1 are intense and have been successfully used to detect chlorites in using remotely sensed data as a diagnostic feature [35, 36]. The two absorption features occurred from 4424 to 4451 cm−1 and from 4240 to 4298 cm−1. The band at 4440 cm−1 was closely related to the IR bands at 3554 and 850 cm−1, and positive linear correlations were established between the NIR band at 4440 cm−1 and IR bands at 3554 and 850 cm−1 following the regressions shown in Figures 7(a) and 7(b). These regressions indicated that the NIR band (hydroxyl combination band) near 4440 cm−1 could be assigned to the combination of (AlAl)O-OH stretching () and (AlAl)O-OH bending () in the IR region. The absorption of the band at 4280 cm−1 was likely related to the IR bands at 3427 and 750 cm−1. Another two positive correlations existed between the NIR band at 4270 cm−1 and IR bands at 3423 and 757 cm−1, as shown in Figures 7(c) and 7(d). These regressions illustrated that the frequency of the NIR band near 4270 cm−1 could be induced by the combination of (SiAl)O-OH stretching () and (SiAl)O-OH bending modes ().

The energy required for the combination band is the sum of the stretching and bending bands in the case of evenly spaced energy levels. Although the wavenumber is proportional to the energy, the combination band should occur at the wavenumber of the sum of the stretching and bending fundamentals. However, due to the anharmonic feature of vibrations, the combination bands appear at a wavenumber that is higher than the sum of the fundamental bands [37]. The observed NIR band assignments are commonly determined by analogy with the IR and bands [26, 38]. Nevertheless, to interpret spectra, researchers often need to establish a relationship that allows them to estimate the wavenumbers of the combination bands in the NIR from the wavenumbers of the fundamentals ( and ) in the IR and vice versa [38]. The anharmonicity constant X, which is determined by the nature of the oscillator and may be affected by the hydroxyl environment (formation of hydrogen bonds, cations, intermolecular interactions, electronegativity, and so forth), is calculated as follows [39]: where is the wavenumber of the combination band in the NIR region and and are the wavenumbers of the stretching and bending bands, respectively, in the IR region. The values of the anharmonicity constant X are calculated from the frequencies (wavenumbers) of the fundamental vibration and the combination vibration for the same OH group (Table 4). The anharmonicity constant values of the 4440 cm−1 band range from 16 to 70 cm −1, with a mean value of 36 cm−1. The anharmonicity constant values of the 4270 cm−1 band range from 85 to 109 cm−1, with a mean value of 94 cm−1.

From (1), the relations between the combination band () and the fundamental bands ( and ) are as follows:

All along the chemical series, the combination bands shifted progressively from 4423 cm−1 to 4450 cm−1 and from 4240 cm−1 to 4298 cm−1. Figure 8 shows two direct linear correlations between band positions and Fe/(Fe + Mg) values for the twenty chlorite samples collected from the USGS spectral library and the two regions, and this relationship indicated that both absorptions observed in the NIR region shifted toward higher frequency as the Fe/(Fe + Mg) value decreased.

These correlations revealed that the Mg-Fe substitution would be the dominant factor in NIR band positions. In other words, the band positions of the two diagnostic combination bands could be used to discriminate among different chlorites with a variety of Fe-Mg contents. As Curtis et al. reported [40], chlorites with higher Fe/(Fe + Mg) values (>0.5) would be chamosite and those with lower values (<0.5) would be clinochlore. Assuming that the absorption coefficients for the two hydroxy combination bands are similar, it is possible to distinguish different chlorites using their Fe/(Fe + Mg) values. The frequency of chamosite absorption should be lower than 4440 cm−1, and the frequency of clinochlore absorption should be higher than 4440 cm−1. In addition, the frequency of chamosite absorption should be lower than 4270 cm−1, and the frequency of clinochlore absorption should be higher than 4270 cm−1.

4. Discussion

Chlorite samples with various Fe-Mg substitutions were studied to establish a correlation between the fundamental vibrations (, , , and ) and combination bands ( and ) of structural OH groups in chlorite minerals. Direct linear relations were found between the wavenumbers of the fundamental absorptions and the combination bands, making it possible to identify which stretching and bending bands in the IR region form the combination bands in the NIR region in chlorites. It has been shown that these relationships can be used for any other minerals, especially for the OH-bearing minerals and carbonate minerals. The collected and experimental data are well fitted with (2) and (3) based on the anharmonic vibration theory. This theory was first used by Petit et al. [39] in studying first overtones in talcs, whereas we used this method to determine the relationship between combination bands and fundamental bands in chlorites for the first time.

From Table 4, the IR OH stretching absorption bands of the samples collected from USGS spectral library range from 3565.44 to 3543.46 cm−1 (covering 22 cm−1) and from 3452.80 to 3406.57 cm−1 (covering 46 cm−1). The IR OH bending absorption features range from 763.87 to 748.73 cm−1 (covering 15 cm−1) and from 653.91 to 642.25 cm−1 (covering 12 cm−1). The results of IR band positions indicated that the NIR combination band positions might be dominated by IR OH stretching bands.

The Fe-Mg substitution with Fe/(Fe + Mg) values between 0.03 and 0.83 at relatively constant Al/(Al + Si) values ranging from 0.13 to 0.44 was accompanied by both a shift toward lower frequency and a widening of the IR OH stretching bands. These results confirm the hypothesis of Petit et al. [34] that IR OH stretching band frequency is dependent on the amount of Fe-Mg substitution. According to Shirozu [8], Fe substitution for Mg in chlorites indirectly weakens the surplus negative charge of the surface oxygens and results in an increase of the O-OH distance and further decreases the frequency of the OH IR stretching band. As the combination band is induced by OH stretching and bending, the NIR OH combination band positions are similarly dependent on the amount of Fe-Mg substitution.

The correlation established above between OH combination bands and fundamental bands provided a technique for determining the cationic environment of the OH group and finding the corresponding OH fundamental bands in the IR region for the interpretation of the OH combination in the NIR region. Compared to IR spectroscopy and Raman spectroscopy, NIR spectroscopy is almost nondestructive, fast, and easy to perform (no preparation of samples) and has a high sensitivity to the hydroxyl group environment. The results clearly show the analytical efficiency of the NIR spectra technique for clay minerals. Nevertheless, the whole interpretation of the observed absorption features needs a further quantitative approach and complementary studies using a rigorous mathematical model. The Raman spectra of the fundamental vibrations of the OH group could also be added, and the methods developed by Koga et al. could be applied [41]. However, this method would require considerable knowledge about the components of the samples. Thus, the relationship given above can be generally applied to better distinguish the presence of certain types of chlorite minerals.

5. Conclusions

The assignment of NIR bands in chlorite spectra and the dominant impact factor of chlorite spectra were achieved by comparing NIR absorption features with IR absorption features and XRF component results. The application of NIR spectroscopy for the study of chlorites shows great potential for understanding the interactions between the fundamental absorptions and combinational absorptions. A number of conclusions can be drawn based on the NIR spectra: (a) chlorites are characterized by two NIR absorptions; the high wavenumber band near 4440 cm−1 could be attributed to the combination of (AlAl)O-OH stretching and (AlAl)O-OH bending in the IR region, and the lower wavenumber band near 4280 cm−1 could be induced by the combination of (SiAl)O-OH stretching and (SiAl)O-OH bending in the IR region. (b) The positions of the two diagnostic NIR absorptions both have negative correlations with Fe/(Fe + Mg) values. The NIR absorption features of Fe-rich chlorites may occur at lower frequencies than 4440 cm−1 and 4270 cm−1, and the NIR absorption features of Mg-rich chlorites may occur at positions higher than these two frequencies.

The interpretation of the NIR spectra of chlorites is also important for mapping hydroxyl-bearing minerals using remotely sensed hyperspectral images and geological field surveys. Spectral analysis with spectral images or field spectrometry is generally used to discriminate compositional variations within altered mineral series, which is important, as mineral composition may systematically vary in an alteration system as a function of the temperature and composition of the hydrothermal fluids and with proximity to zones of mineralization. Mapping these alterations can allow remote sensing researchers and geologists to locate samples within an alteration system to distinguish a certain mineral from other minerals and to define alteration relationships [11, 42, 43]. Accurate interpretations of spectra could contribute to more accurate quantitative analyses of phyllosilicates and to a better understanding of their mineralization relationships.

Conflicts of Interest

The authors declare no conflict of interest.

Acknowledgments

This research was conducted as part of the Fundamental Research Grant Scheme of the National Natural Science Foundation of China (vote no. 41502312), the Chinese Ministry of Science and Technology (vote no. 2015BAB05B03-01), and the China Geological Survey Foundation (vote no. DD20160002). The authors are thankful to the USGS spectral library (https://speclab.cr.usgs.gov/spectral-lib.html), the Key Laboratory for Geohazards in Loess Areas, the Chinese Ministry of Land and Resources, and the Xi’an Center of China Geological Survey for their contributions to this investigation.