Abstract

A significant linear correlation between Ba and Eu contents has been observed in a set of mudrock formations in the lower Cambrian Hetang Formation in South China. Combined with petrographic features and reinvestigation of published data obtained from the same laboratory, the analytical interference of Ba can be excluded as a main factor of this correlation. Electron microprobe analyses (EMPA) have indicated that Ba-rich minerals (hyalophane and cymrite) precipitated from hydrothermal fluids account for the total Ba content. It is evident that significant contribution of both Ba and Eu from hydrothermal fluids is the cause of this linear relationship, while this coexistence is ultimately controlled by their similar chemical properties and abundances in rock sequences along the fluid path. This interpretation was supported by similar linear Ba-Eu correlations retrieved from barite-bearing chert in another area of the Yangtze platform during the same geological period. Therefore, linear Ba-Eu correlations may be more common than previously thought in this area, where hydrothermal activities were active during the specific time interval.

1. Introduction

The rare earth elements (REEs) include lanthanum (La) and the f-block elements (cerium through lutetium) in the periodic table, whose atomic number ranges from 57 to 71. Scandium and yttrium are also included in this group as they have ionic radii similar to the lighter f-block elements and are found coexisting in the same ores (Atwood, 2012). The REEs are insoluble in most geological settings and resistant to remobilization beyond the mineralogical scale during weathering, diagenetic, and metamorphic processes. That makes them important tracers for characterizing and understanding a variety of geochemical “reservoirs” ([1, 2]; Atwood, 2012).

The most stable oxidation state of all REEs is the REE (III) form, except for europium (Eu) and cerium (Ce) that are also present in the form of Eu2+ and Ce4+, respectively [3]. The distinctive redox chemistries of Eu and Ce are of considerable importance to understanding geochemical conditions with unique insights into magmatic, aqueous, and sedimentary processes ([1, 2, 46]; Murray et al., 1990; Atwood, 2012). The ionic radius of Eu2+ is about 17% larger than that of Eu3+, which leads to different substitution behaviors of Eu2+ than other trivalent REEs (e.g., Sm and Gd) and subsequent anomalous REE patterns (i.e., Eu anomalies). Europium occurs in the form of Eu2+ under highly reducing conditions. It exists only within certain magmatic and/or hydrothermal environments that are rarely found at the surface of the Earth. An important geological example is that Eu can be highly concentrated in plagioclase feldspar by substituting into the calcium site. The anomalous Eu abundance in magmatic rocks indicates relatively shallow igneous partial melting or fractional crystallization processes (Atwood, 2012). Positive Eu anomalies have also been widely reported in modern deep-sea hydrothermal systems [79] and in massive deposits of submarine hydrothermal exhalative origin [1012]. Consequently, it has been used as an effective tracer for tracking and characterizing hydrothermal activities in sedimentary rocks and ore deposits [1319]. In addition, in spite of no direct evidence showing europium reduction under surface conditions, positive Eu anomalies might be present at highly reducing conditions during diagenesis. For example, the trace and rare earth element abundances in phosphate nodules from the lower Cambrian strata in South China have shown positive Eu anomalies, possibly implying the replacement of Ca2+ by Eu2+ under reducing conditions [20].

Accurate measurement of Eu abundances in geological samples is a prerequisite to effectively obtaining important information of geochemical reactions and environmental conditions in natural processes recorded by elements. The determination of Eu abundances by inductively coupled plasma-mass spectrometry (ICP-MS), however, is interfered by the coexisting element of barium (Ba). This is due to similar mass/charge ratios of Ba oxide (BaO) and hydroxide (BaOH+) species that partially overlap those of Eu isotopes (151Eu and 153Eu), with 151Eu being less affected by BaO and BaOH+ than 153Eu [2123].

A routine assessment of the interference level of Ba is to determine the yield of Ba oxide and hydroxide species using a single-element solution of Ba (e.g., 500 ng/ml Ba in BaCl2 solution) and compare the yield of potential interfering species (e.g., 135Ba16O) with the Eu isotope that has the same mass number (e.g., 151Eu) in a single-element solution of Eu [21, 24]. The previous results indicated that 1000 ng/g of Ba could cause an apparent Eu concentration of 0.22 ng/g in 151Eu [21]. An alternative way to evaluate the extent of the interference is analyzing the correlation between Ba and Eu in the samples. Plots of Ba and Eu have been widely used in previous studies, and linear Ba-Eu correlation has been commonly interpreted as an indicator of apparent interference of Ba with Eu, with less attention on other possible causes ([5, 20, 25]; Ling et al., 2013).

To better understand the ultimate controlling factors, if any, of Ba-Eu correlation other than analytical errors embedded with analytical techniques, we investigated the trace and rare earth elemental contents of a drilling profile in the lower Cambrian Hetang Formation, South China. Hydrothermal activities and Ba enrichment have been reported in previous studies for this formation and its stratigraphic equivalents in other areas [16, 2528]. In this study, a series of petrological and mineralogical analytical techniques are used to quantify the abundances of Eu and Ba and elucidate their relationship with hydrothermal activities.

2. Geological Setting and Samples

During the transition between the Ediacaran period and the Cambrian period, the Yangtze platform in South China gradually evolved from a rift to a passive continental margin basin in response to extensional tectonism [29, 30]. With sedimentological features, the sedimentary sequences can be classified into three distinct facies: shallow water, transitional, and slope to deep basin facies (Figure 1). The shallow water facies is represented by thick carbonate strata, while the slope to deep basin facies consists of a series of interbedded black chert and shale. The transitional facies is composed of both carbonate and black shale [20, 29]. The occurrence of ore deposit beds containing Ni-Mo-PGE-Au sulfides has been reported in this black shale-chert sequence along the NE transitional belt at the margins of the Yangtze platform in lower Cambrian [15].

The samples in this study were taken from a drilling section (Well Wanning 2) located at Ningguo, Anhui province (Figure 1), covering lithological units of the Xijianshan and Hetang formations. The sedimentary strata of the Hetang Formation include mudrock, chert, and carbonates. The major mineralogy of mudrock varies, including carbonaceous materials, silicates, or Ca-containing minerals. The Hetang Formation is underlain by the Xijianshan Formation which is mainly composed of dolostone and chert. Carbonate and quartz veins are common in this set of strata, indicating enhanced fluid migration during diagenesis. The samples cover the whole depth of the Hetang Formation and are extended into the Xijianshan Formation (Figure 2). To ensure accurate measurements, two batches of samples with similar mineralogy and petrographic features at selected depths were collected for analyses (see below).

3. Methods

The trace element analysis was carried out on a Finnigan Element II ICP-MS at the State Key Laboratory for Mineral Deposits Research at Nanjing University. For carbonate rocks, 100 mg of powdered samples was weighted in Teflon beakers and dissolved by 1 M HAc under ultrasonic bathing for 2 hours. For mudrock and chert, 50 mg of powdered sample was dissolved by a mixed HF-HNO3 solution in high-pressure Teflon bombs. For calcareous or dolomitic mudrock, the carbonate fractions were dissolved by 1 M HAc and the solid residues were then dissolved by the same mixed acid solution in high-pressure Teflon bombs. After dissolution, the solutions were centrifuged and the supernatants were transferred to clean Teflon beakers for evaporation. The solutions were added with 1 ml of HNO3 and dried three times to remove HAc or HF. The residues were then dissolved in 5% HNO3 solution for analysis. A solution containing Rh in 10 ppb was used as an internal standard to monitor signal drifting during measurements. The detailed analytical procedure was described in the work of Gao et al. [31], and parameters for the instrument were presented in Table 1. Analytical uncertainties are estimated to be less than 10%. The value of Eu/Eu is calculated as , where EuPAAS, SmPAAS, and GdPAAS represent the abundances of Eu, Sm, and Gd, respectively, normalized by the post-Archean Australian Shale (PAAS) standard [32].

The electron microprobe analysis (EMPA) was carried out on polished, carbon-coated thin sections to study the petrological and mineralogical features in the samples. It was conducted using a JEOL JXA-8100 electron microprobe with a wavelength-dispersive system at the State Key Laboratory for Mineral Deposits Research at Nanjing University. The instrument was operated with three crystal spectrometers, with the accelerating voltage and the specimen current of 15 kV and 20 mA, respectively. The diameter of the beam spot used for quantitative element analysis was 1 μm. The standards used were natural minerals supplied by the American National Standards Institute (ANSI). The detection limit of the microprobe analysis was approximately 0.002% for each element.

4. Results

4.1. Petrological and Mineralogical Features

The locations of samples analyzed by EMPA within the profile are labelled in Figure 2. Minerals containing Ba and K, such as hyalophane and cymrite, were found in all of the mudrock and chert samples with backscattered electron (BSE) observation and EMPA analyses. For each sample analyzed by EMPA, the occurrence of Ba-rich minerals, in hydrothermal veins and/or in matrix, is also specified (Figure 2).

Hyalophane is present as the matrix in all the samples analyzed, while hydrothermal veins that contain hyalophane are observed in most samples. In matrices, the occurrence of hyalophane includes fillings of pore spaces in chert samples (Figure 3(a)) and thin rims around clastic K-feldspar grains (Figures 3(b)3(d)). In hydrothermal veins, hyalophane is present as irregular or amorphous aggregates. The close association of hyalophane in veins with the whole matrix is very common (Figure 3). In particular, the texture of hyalophane-bearing veins penetrating into the matrix is observed, with abundant hyalophane rims formed around clastic K-feldspar grains (Figures 3(b) and 3(c)). Cymrite grains show euhedral tabular forms, primarily coexisting with hyalophane in the matrix in a number of chert samples (Figure 3(a)).

4.2. Trace and Rare Earth Elements

Samples over the whole profile, Batch 1, were measured by ICP-MS for abundances of trace and rare earth elements, including carbonates, mudrock, and chert (Table 2). To limit analytical errors from instruments, a second batch of samples, i.e., Batch 2, with similar mineral compositions and petrological features were also used for analyses (Table 3).

There are 7 carbonate samples in Batch 1, including dolostone and limestone. The concentration of Ba and Eu in carbonates ranges from 76.67 to 326.92 ppm (Figure 4(a)) and 0.19 to 0.45 ppm (Figure 4(b)), respectively. Slight negative Eu anomalies are observed, with Eu/Eu values from 0.71 to 0.99 (Figure 4(c)). There is a covariation of Sm and Gd concentrations in carbonate samples, with both being in the range of 0.50 to 3.00 ppm (Figure 4(d)).

The abundances of Ba and Eu in mudrock and chert are much higher than in carbonate rocks. In Batch 1, the concentration of Ba ranges from 431.47 to 20386.00 ppm, approximately 40 times higher than that in carbonates on average (Figure 4(a)). The concentration of Eu is in the range of 0.27 to 6.22 ppm, which is about 9 times higher than that in carbonates (Figure 4(b)). Values of Eu/Eu range from 1.12 to 14.17, suggesting pronounced positive Eu anomalies (Figure 4(c)). In Batch 2, the concentration of Ba and Eu ranges from 1818.70 ppm to 17316.74 ppm (Figure 4(a)) and 0.86 ppm to 3.71 ppm (Figure 4(b)), respectively. The average concentrations of Ba and Eu are 30 and 7 times higher than those in carbonates, respectively, with Eu/Eu values ranging from 1.08 to 11.20 (Figure 4(c)). The covariation of Sm and Gd can also be observed in mudrock and chert samples in both batches. The range becomes broader for both batches, varying from 0.78 to 16.15 ppm (Figure 4(d)).

5. Discussion

5.1. Ba-Eu Correlation

In previous studies, to evaluate the interference of Ba with Eu concentrations, correlations of Ba and Eu were verified with different elemental ratios, including the relationship between Ba/Nd and Eu/Eu ([5]; Ling et al., 2013), Eu/Eu and Ba/Sm [20], and Ba and Eu/Eu[25]. Among them, instead of the absolute concentration of Eu, the ratio of Eu/Eu was the one that is frequently used to evaluate the Eu abundance. Due to intrinsic chemical properties and complex redox conditions, Sm and Gd behave differently than Eu in natural environments. Changes in their abundances may not result from the same geological event(s) that impacted the Eu speciation and its distribution in different phases. Therefore, the involvement of Sm and Gd in evaluating Ba interference with the ratio of Eu/Eu may obscure the original Eu distribution and authentic Ba-Eu correlation. The abundances of Ba and Eu are evaluated directly in this study.

The concentrations of Ba and Eu in mudrock and chert samples in Batch 1 show a significant linear correlation, with the coefficient of determination () of the linear regression line being 0.88 (Figure 5(a)). In carbonate samples within the same batch, however, no linear Ba-Eu correlation is observed (Figure 5(b)). Concentrations of Ba and Eu in mudrock and chert samples in Batch 2 have also shown linearity between them, with an of 0.86 (Figure 5(c)). By combining data from both Batch 1 and Batch 2, a linear Ba-Eu correlation with an overall of 0.85 suggests a consistent relation between Ba and Eu in mudrock and chert samples (Figure 5(d)), with the Ba/Eu ratios ranging from 797.32 to 4806.16.

The abundance of Eu (i.e., EuT), reported as the result collected by ICP-MS, has two potential sources: the authentic Eu content in samples (EuA) and the signal error caused by Ba interference (EuF). The value of EuA is the direct result of involved geological processes, while EuF is linearly correlated with the Ba concentration. If EuF is much higher than EuA, a linear Ba-EuT correlation most likely indicates major interferences from Ba. This is how linear Ba-EuT correlations were interpreted in previous studies ([5, 20, 25]; Ling et al., 2013). However, if Ba and EuA are enriched by the same geological processes and the abundance of EuA is dominant in the total Eu, a Ba-Eu linear correlation is a genuine record of geochemical processes involving both elements. Both possibilities are discussed in the following sections to unravel the cause of linear Ba-EuT correlations in this study.

5.2. Evaluation of Analytical Interferences

The interference of Ba in the determination of Eu was a systematic error caused by the overlap of peaks from Ba oxide and hydroxide species with Eu isotopes during ICP-MS analysis [21]. Due to the difference in system errors of individual instruments and internal standards used, the ratio of the concentration of Ba over EuF could vary significantly between different instruments. For example, the Ba/EuF ratios of 1000 : 0.22 and 1000 : 0.10 from ICP-MS analyses were reported by Dulski [21] and Smirnova et al. [24], respectively. Despite this discrepancy, for the same instrument, the Ba/EuF ratio is believed to be constant under similar working conditions [21].

As reported above in “Introduction,” to validate the extent of Ba interference on Eu, the most effective method is to determine the yield of Eu caused by single-element solutions of Ba oxide or hydroxide with varying concentrations [21, 24]. Such experiment has been conducted on a Finnigan Element XR ICP-MS in the same laboratory, with identical pretreatment procedures and similar analytical conditions to the Finnigan Element II ICP-MS used in this study [33]. The results revealed a Ba/EuF ratio of 400000 : 1 to 1400000 : 1 for the low-resolution (300) mode analysis and much higher Ba/EuF ratios for the mid-resolution (4000) and high-resolution (10000) mode analyses. In this case, the Ba/EuF ratio was supposed to be at least 100 times higher than the Ba/EuA ratio of analyzed samples in this study, demonstrating that the extent of Ba interference on Eu determination was likely minimal.

In addition, previously published data of Ba and EuT concentrations obtained from the same instrument under similar working conditions are also collected to evaluate the Ba/EuF ratio. In a study of black shale in the Niutitang Formation (lower Cambrian) in South China, the concentrations of Ba and EuT were reported to be in the ranges of 659 to 13686 ppm and 0.33 to 3.97 ppm, respectively [15] (Figure 6(a)). The Ba/EuT ratio was calculated to be varied from 391.54 to 8146.43. Both the Ba and Eu concentrations in black shale, along with their ratios (Ba/EuT), are in a similar range as in samples of this study (Tables 2 and 3). However, no linear correlation between Ba and EuT can be extracted in spite of a weak covariation trend, suggesting that EuF is not a major component of the reported EuT values.

Another set of data comparable to this study was reported on phosphate nodules of the Mufushan Formation, which is also in the period of lower Cambrian in South China [20]. The Ba and EuT concentrations are in a narrower range than other reported data, from 817 to 6144 ppm and 0.05 to 11.42 ppm, respectively (Figure 6(b)). The Ba/EuT ratio varies significantly from 95.53 to 59040.00. There is no linear correlation between Ba and EuT (Figure 6(b)). Therefore, considering the similar range of Ba concentrations and Ba/EuT ratios obtained in this study, systematic analytical errors can be excluded as a main factor attributing to the significant linear correlation between Ba and Eu.

Furthermore, the concentrations of Sm and Eu from most samples in this study have also shown a positive correlation (Figure 7(a)). A similar relationship can be observed between Eu and Gd as well (Figure 7(b)). Unlike Eu, the interference of Ba on Sm and Gd is believed to be negligible because of differences in mass over charge ratio [21]. Therefore, the linear Ba-Eu correlation, coupled with Eu-Sm and Eu-Gd linearities, suggests the presence of an intrinsic relationship between Ba and Eu which is mainly caused by geological processes. The positive Eu-Sm and Eu-Gd correlations in this study may result from the impact of the geological processes and the original elemental feature of the sediments.

5.3. Evidence of Ba-Eu Coexistence

Europium is commonly enriched in plagioclase as a result of substitution into the calcium site. Similar to Ca2+, Ba2+ can also be substituted by Eu2+ due to the same electronic charge and similar cation radius [3436]. Therefore, Eu may be enriched in minerals or geological fluids that have high abundances of Ba. A linear Ba-Eu correlation can be present as a result of their substitutional relationship, particularly when Ba2+ is the dominant cation in minerals (or geological fluids) as shown above (Figures 4(a), 4(c) and 4(d)).

The linear relationship between Ba and Eu suggests that the majority of Ba and Eu in this set of rock formations are closely related during geological processes. The BSE observation and elemental analysis by EMPA provide additional evidence showing the highly concentrated areas of Ba in the samples. The Ba-rich minerals, hyalophane in particular, are pervasive in the rock sequence and present in the matrix of all the samples analyzed by EMPA (Figure 2). Hyalophane is believed to be formed in limited geological environments, and its origin is mostly from hydrothermal fluids (McSwiggen et al., 1994; Moro et al., 2001).

In the rock matrix, hyalophane mainly occurs as rims of clastic K-feldspars (Figure 3(c)). Both hyalophane and cymrite are also present in intergranular pores (Figure 3(a)). The hydrothermal veins containing hyalophane are common textures, shown penetrating into the matrix (Figures 3(b) and 3(c)). Overall, the occurrence of such Ba-rich minerals suggests that they were formed by infiltration of late-stage hydrothermal fluids into the rocks. Hyalophane of hydrothermal origin with similar occurrences (e.g., hyalophane rims around K-feldspar) was also found in black shale in the Devonian period from the Selwyn Basin, Canada [37].

In this study, both barium and europium elements are concentrated by late-stage hydrothermal fluids passing through surrounding sedimentary sequences. They were then accumulated in hyalophane and cymrite during diagenesis. This scenario agrees well with the consistently high concentration of Ba in the chert and mudrock samples analyzed, with an average value of 6431.04 ppm (Tables 2 and 3). It is also confirmed by the fact that, as hyalophane and cymrite are not included in the trace element and rare earth element analyses of carbonate rocks because they are insoluble in HAC, lower Ba concentrations and nonlinear Ba-Eu correlation are obtained (Figures 4(a) and 5(b)).

Enrichment of Eu has been reported in modern deep-sea hydrothermal systems [79] and in massive deposits of submarine hydrothermal exhalative origin [1012]. High abundances of both Ba and Eu in hydrothermal fluids are present in a variety of hydrothermal environments [7, 3840].

Wang et al. [25] investigated the occurrences and textures of chert on the marginal zone of the Yangtze platform in the western Hunan area, during the Ediacaran-Cambrian transition. Four lithotypes of chert are identified in the study, including mounded, vein, brecciated, and bedded chert. Further geochemical investigations have shown prevalent positive Eu anomalies in chert and homogenization temperature of fluid inclusions in the range of 120-180°C. Consequently, the chert was suggested to be formed by hydrothermal fluids. The variation of occurrences and textures was considered to be a reflection of different formation conditions. The mounded chert was formed at or near the vent fields, while the formation of brecciated chert occurred in large-scale fracture systems. The chert in veins was precipitated from silica-rich hydrothermal fluids ascending along fractures, and the bedded chert was formed in a quiet water environment away from vent fields. This wide range of occurrences of hydrothermal chert was ascribed to the extensional tectonism in the Yangtze Block during the Ediacaran-Cambrian transition [25, 41].

Although the positive Eu anomaly was observed in most chert samples, barite was only found in mounded and brecciated chert. Only in those two lithotypes of chert can a prominent linear Ba-Eu correlation be extracted (Figure 8(a)). No significant Ba-Eu correlation can be recognized in the other two types of chert (Figure 8(b)), although the concentration of Ba and the Ba/Eu ratio of bedded chert are in a similar range as those of mounded and brecciated chert.

The results above provide another line of evidence supporting that the distinct Ba-Eu correlations in different lithotypes of chert are caused by the source of elements and environmental conditions. In areas close to vent fields or having large-scale fracture systems, where the mounded and brecciated chert were formed, chemical components in hydrothermal fluids can be preserved. However, bedded chert was suggested to form far away from vent fields and in a quiet water environment, where the original correlation between Ba and Eu contents within hydrothermal fluids could be erased by secondary processes. Therefore, we suggest that the linear Ba-Eu correlation in mounded and brecciated chert most likely is being caused by significant contribution from Ba-rich hydrothermal fluids.

6. Conclusions

Trace elemental analysis of the lower Cambrian Hetang Formation in Anhui Province has shown a linear relationship between Ba and Eu contents in mudrock and chert formations. Combined with petrological and mineralogical observation, significant contribution of Ba and Eu from Ba-rich hydrothermal fluids is proposed to account for that correlation. Positive Eu anomaly in the rocks is suggested to be an effective record of the influence of hydrothermal activities.

The ratio of Ba/Eu provides important information for evaluating the interference of Ba on determination of Eu concentrations during analysis. Relatively low Ba/Eu ratios and significant linear Ba-Eu correlations may suggest that Ba and Eu were originated from the same hydrothermal source. In addition, petrological and mineralogical investigations are essential in identifying the occurrences of Ba and Eu in rocks and facilitating understanding of Ba-Eu correlations.

At the western boundary of the Yangtze platform, pronounced linear Ba-Eu correlations were also present in barite-bearing chert (mounded and brecciated chert) formed by Ba-rich hydrothermal fluids. We suggest that such linear Ba-Eu correlations may be common in geological environments, which is ultimately constrained by the availability of both elements in rock sequences and the timing of hydrothermal activities.

Data Availability

The trace and rare earth element data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was financially supported by the National Key Research and Development Program (Grant no. 2017YFC0603101), National Natural Science Foundation of China (Grant nos. 41802026, 41830425, 41890845, and 41621003), China Postdoctoral Science Foundation (Grant no. 2018M633553), and 111 Project (D17013). Research support to QF from the NSF CAREER Program under award OCE-1652481 and the American Chemical Society Petroleum Research Fund 54474-DNI2 is also acknowledged. Special thanks are due to Dr. Wenlan Zhang and Qian Liu for their help with analytical work.