4.1 Distribution of rare earth elements and provenance tracing in the deep phaeozems
The relative abundance of rare earth elements in the soil varies during the soil formation due to the combined influence of soil-forming parent material, climatic conditions, humus, and clay components[48]. Chondrite is among the initial materials on earth. Therefore, it neither contains light nor heavy rare earth elements and hence can be used as a reference for determining the degree of fractionation in the samples. The data on the rare earth elements in phaeozems was normalized using chondrite for accounting for either odd or even numbers of elements influencing the relative abundance of rare earth elements[49], as shown in Table 3. The distribution of rare earth elements in the phaeozems of Arongqi is represented in Fig. 5. The distribution pattern indicates that all the sample points show a right-sloping pattern, and this pattern is consistent with that previously observed in the Chinese soils[38], North American shales (Haskin and Paster, 1979), and the Earth's crust, implying that the source of deep phaeozems is terrestrial. The slope of the curve from La to Eu is steeper, whereas the slope of the curve from Eu to Lu is gentle. Ce has a distinct valley shape, while Eu has a weak valley shape, both of which show negative anomalies. The distribution of rare earth elements in the deep phaeozems shows that LREE is relatively abundant and exhibit typical terrestrial rare earth characteristics[50–51].
(La/Yb)N method is used to characterize the degree of differentiation between light and heavy rare earths. A ratio greater than one implies enrichment of LREE, whereas a ratio less than one suggests depletion of HREE[52]. As observed from Table 3, the highest value of (La/Yb)N is 10.41, the minimum value is 5.92, and the average value is 7.58, all of which are greater than 1, implying that the LREE and HREE are well differentiated in the deep phaeozems of Arongqi. (La/Sm)N is a coefficient reflecting the degree of fractionation of LREE[53], and (La/Yb)N is a coefficient reflecting the degree of fractionation of HREE. It can be observed from Table 3 that the (La/Sm)N values of the deep phaeozems range from 2.54 to 4.89 with a mean value of 3.54, and there is a huge variation in the fractionation of LREE. The (La/Sm)N value ranges between 0.90 and 1.63, with a mean of 1.14, indicating that the fractionation of HREE in the deep phaeozems is not significantly different. However, the fractionation of LREE is more pronounced than that of HREE in the deep phaeozems samples.
The analysis of rare earth elements is an important way of deposit tracing and can reliably maintain the geochemical information about the originating place[54]. The rare earth element standards are frequently used to characterize the provenance quality. Numerous studies have demonstrated that the deposits with the same provenance frequently exhibit comparable rare earth element distribution curves[53]. A low similarity in the rare earth model curves is observed with a wide fluctuation in the distribution of LREE. Moreover, sometimes, dissimilar model curves are observed, and characteristics of multiple provenances become evident. Studies have demonstrated that when the parent rock is a medium-acid felsic rock such as granite, ∑REE of deposits is relatively high, the Eu negative anomaly is evident, and the LREE/HREE ratio is high. When the parent rock is basic, ∑REE is relatively low, the negative Eu anomaly is not evident, and the HREE are relatively abundant[55–56]. Based on this characteristic, it is hypothesized that the pedogenic parent rock in the deep phaeozems contains granite components. This can be verified by analyzing the provenance properties of deep phaeozems using the La/Yb-ΣREE map and La/Yb-ΣREE, as shown in Fig. 6. The deep phaeozem samples were concentrated in granite and alkaline basalt areas, proving the previous theory that the deep phaeozems may have granite components. However, granite and alkaline basalt rocks were not discovered in the bottom bedrock drilling and surrounding mountains[29], indicating that the bottom and surrounding bedrock were not the origins of phaeozems. The rare earth elements and types of volcanic rocks from the Longjiang Formation (K11) and Guanghua Formation (K1gn) do not conform to the provenance properties[29]. Figure 6 demonstrates that the material source for deep phaeozems is a mixture of rocks, and the findings agree with the REE distribution pattern for the deep phaeozems. Therefore, it can be stated that the parent material in the deep phaeozems is transportable. This is consistent with the inference that the accumulation-type parent material of phaeozems gets deposited on the surface in the study area[29].
4.2 δCe, δEu anomaly, and anomaly index (Ceanom) in the deep phaeozems
Deep phaeozems differ from the surface-layer phaeozems in some aspects, such as minimal human interference during the formation of deep phaeozems and the rigid environment surrounding the deep phaeozems. The rare earth elements, also well-known as environmental response factors, react to the environmental changes throughout the soil formation process[40]. The variable rare earth elements, Eu and Ce, are susceptible to redox reactions in specific natural environmental conditions and thus, exhibit geochemical properties distinct from those of other rare earth elements. δCe and δEu are often used to reflect upon the soil-forming environment[43]. The δCe and δEu values represent the anomalies of Ce and Eu, respectively. δCe and δEu < 1 indicate negative anomalies, while δCe and δEu > 1 indicate positive anomalies[42,57]. Under oxidizing conditions, Ce3+ is easily oxidized to more stable Ce4+. Thus, Ce is separated from other trivalent rare earth elements, showing a positive δCe anomaly. On the contrary, under reducing conditions, δCe shows a negative anomaly as Eu3+ is reduced to Eu2+. Eu element mainly exists in the form of Eu3+; therefore, the reducing environment results in the loss of Eu, leading to δEu negative anomaly[57]. As shown in Table 3, both δCe and δEu in the deep phaeozems show negative anomalies, implying a depletion of Ce and Eu in the deep phaeozems under reducing conditions in the presence of water.
The anomaly index (Ceanom) method can also be deployed to understand the environment surrounding the deep phaeozems[32]. Several research and verification studies have highlighted that this index is a good indicator of paleoenvironment restoration[57]. Ceanom >-0.1 indicates an anoxic reducing environment, whereas Ceanom < -0.1 indicates an oxidizing environment. As shown in Table 3, the Ceanom values for the deep phaeozem samples range from 3.49 to 4.20, with a mean of 3.90. This value is much greater than − 0.1, implying anoxic reducing conditions in the deep layer, and is consistent with the inferences drawn based on δCe and δEu data described above.
Table 3
Different parameters related to rare earth elements in the deep phaeozem samples from Arongqi
| REE | LREE | HREE | L/H | La/Yb | La/Sm | Gd/Yb | δCe | δEu | Ceanom |
Min | 283.32 | 207.17 | 66.35 | 2.71 | 5.92 | 2.54 | 0.9 | 0.36 | 0.51 | 3.49 |
Max | 535.62 | 405.63 | 129.99 | 3.87 | 10.41 | 4.89 | 1.63 | 0.46 | 1.03 | 4.2 |
Average | 391 | 298.71 | 92.29 | 3.24 | 7.58 | 3.54 | 1.14 | 0.41 | 0.86 | 3.9 |
4.3 Salinity of water bodies formed by the deep phaeozems
Sr/Ba and Th/U were analyzed to characterize the salinity of water bodies formed by the deep phaeozems. Both Sr and Ba exist in water in the form of soluble bicarbonate, sulfate, and chloride, and Sr has a stronger migration ability than Ba. When the water salinity starts increasing, Ba precipitates, and further, when the salinity reaches a particular level, Sr starts precipitating, Sr/Ba increases, and Sr/Ba is positively correlated with the water salinity. According to one study, Sr/Ba in the freshwater environment is less than 0.6, Sr/Ba in brackish water is between 0.6 and 1, and Sr/Ba in salty water is greater than 1. Th/U can also be used as a major means for the measurement of paleosalinity. Th is easily adsorbed by the clay minerals, and U is easily leached or oxidized. Th/U is less than 2 in saltwater habitats, whereas it is higher in freshwater environments[58–60]. The Th/U and Sr/Ba ratios are shown in Fig. 7. All the samples have a Th/U value greater than 2, with a mean value of 2.87. Most of the Sr/Ba values are less than 0.6, except Sr/Ba values of the deep phaeozems at HT1-6 and HT0-2, which are between 0.6 and 1 with a mean value of 0.38. Both the Th/U and Sr/Ba values indicate that the freshwater water body is formed by the deep phaeozems. A comprehensive analysis in conjunction with the investigation of types of borehole deposits further suggests that the water body formed by the deep phaeozems is Arun riverine system.
4.4 Paleoclimate records of the deep phaeozems
The concentration and composition of trace elements in the deposits can provide valuable information about paleoclimate. Therefore, several characteristic elements are frequently used for the reconstruction of the paleoclimate[59–62]. At present, the widely used element indices are Rb/Sr[63], climate index (CI), and Sr/Cu[61]. Rb has a relatively strong geochemical connection under weathering conditions. Sr is easily oxidized and leached out. Strong weathering leads to the leaching of Sr under warm and humid climatic conditions. However, Sr leaching is minimal because of weak weathering under arid climatic conditions. A lower Rb/Sr value implies a dry climate, and a higher Rb/Sr value indicates a warmer and humid climate. Sr is a typical dry-prone element, and a high Sr content indicates a dry climate. Cu is a typical humidity-prone element, and a high Cu content indicates a humid climate. Sr/Cu < 10 represents a humid climate, and Sr/Cu > 10 represents a dry and hot climate[61,63]. High concentrations of elements in the deposits, especially Ni, Zn, A1, Fe, and Co, indicate a humid and warm climate. Under dry conditions, a significant proportion of Sr, Ca, Mg, and Na undergo precipitation because of the loss of water by evaporation and increased alkalinity. Thus, higher Sr, Ca, Mg, and Na indicate an arid climate. CI is used to describe the changes in the climate system. Lower CI implies a dry climate, while higher CI indicates a humid one[61]. CI provides an accurate description of both modern and paleoclimates[39,64]and can be calculated by Eq. 7[61]. The climate patterns based on the values of CI are shown in Table 4.
CI = (Fe + Mn + Cr + Ni + V + Co)/(Ca + Mg + Sr + Ba + K + Na) (7)
Table 4
Climate patterns based on the climate index
Climate index | < 0.2 | 0.2–0.4 | 0.4–0.6 | 0.6–0.8 | 0.8–1.0 |
Type of Paleoclimate | dry heat | dry heat | Semi dry heat-Semi humid | Semi humid | Humid |
Rb/Sr values in the deep phaeozems of Arongqi are shown in Table 5 and Fig. 8. The maximum and minimum values of Rb/Sr are 0.53 and 0.04, respectively, with a mean value of 0.23. The statistics reveal that the Rb/Sr ratio is 0.44 when CI is 0.2, the ratio is 0.47 when CI lies in the range of 0.2–0.4, and the ratio is 0.09 when CI lies in the range of 0.4–0.6. When CI is in the range of 0–0.4, the Rb/Sr ratio is 0.91, and it indicates an arid climate, as shown in Fig. 8. Sr/Cu ratio in the deep phaeozems of Arongqi is much greater than 10. The minimum and maximum values of Sr/Cu are 10.52 and 84.36, respectively, with a mean value of 31.02. Sr/Cu value indicates that the recorded climate in the deep phaeozems of Arongqi was arid. Table 5 and Fig. 8 show that the CI of the deep phaeozems ranges from 0.58 to 0.05, with a mean value of 0.22. The statistics show that the proportion of CI less than 0.2 is 0.56, the proportion of CI in the range 0.2–0.4 is 0.31, and the proportion of CI in the range 0.4–0.6 is 0.13. The proportion of CI in the range 0–0.4 is 0.87, with a mean value of 0.29, implying a dry-hot-semi-dry-hot climate. Rb/Sr, Sr/Cu, and CI all indicate that the climate recorded in the deep phaeozems of Arongqi was dry-hot-semi-dry-hot.
Table 5
Rb/Sr and climate index (CI) in the deep phaeozems of Arongqi
Proportion | 0–0.2 | 0.2–0.4 | 0.4–0.6 | 0.6–0.8 | 0.8–1.0 |
Rb/Sr | 0.44 | 0.47 | 0.09 | 0 | 0 |
CI | 0.56 | 0.31 | 0.13 | 0 | 0 |
4.5 Formation of the phaeozems from various parent materials, including fluvial deposits
The salinity of water bodies was recorded in the deep phaeozems based on the contact boundary between river deposits and deep phaeozems in quaternary boreholes. Moreover, the source of materials, climate information, and the formation of deep phaeozems were also similarly discussed. The study concluded that the deep phaeozems were formed from an early riverine marsh ecosystem under an arid climate system.
The formation of phaeozems has been discussed for a long time, and several hypotheses have been put forward, including the anthropogenic origin of phaeozems from the loess parent material at Neolithic sites[65–67], natural origin of surface phaeozems from the loess parent material[68–70], the origin of phaeozems from the zeolite basalt residues[71], and origin of phaeozems from the volcano deposits[72]. However, there is no evidence that the phaeozems of anthropogenic origin exist in Neolithic sites of Northeast China; therefore, only phaeozems of natural origin were considered in this study. The phaeozems are formed as a result of the interactions of several elements, and the underlying mechanism is complex. As described previously, other elements also play a greater part in the formation of phaeozems from river parent material deposits under arid conditions. The environment and mechanism of phaeozem formation are determined by the combined effect of many factors, with each factor playing a different role. Therefore, this study considers the phaeozems of natural origin and explores multiple formation mechanisms of phaeozems derived from various parent materials. As shown in Table 6, the phaeozems formed from different parent materials are compared to understand the formation and evolution of phaeozems and the major controlling factors determining the formation of phaeozems.
Table 6
Formation of phaeozems from different parent materials
Phaeozem formation | Loess parent material[69] | Fluvial deposit parent material | Basalt Residual parent material [71] | Andesitic residual parent materials [72] |
Based on the type of parent material | Loess | Fluvial deposit | Basalt (weathered basalt) | Andesitic materials |
Modern distributed vegetation | Modern forest grassland landscape | Modern grassland landscape | Modern forest landscape | Modern forest landscape |
Elevation (m) | 220–300 | 320–350 | 850–1150 | Around 500 |
Physiographic or Topographic position | Hill, mountain root | river terraces | Hill slope | Above the island |
Climate zone | Temperate | Temperate | Tropical | Tropical |
Phaeozem records vegetation characteristics | Xerophytic, mesophytic, hygrophytic, aquatic | Xerophytic, hygrophytic, aquatic | | |
Phaeozem records climatic characteristics | Warm and wet | Arid - semi-arid | Hot and humid | Hot and humid |
Characteristics of soil formation | Warm climate and humid | Arid climate | Presence of zeolite minerals | Coral organic matter |
Factors controlling soil formation in this study | Coniferous forests, broad-leaved forests (deciduous), shrubs and terrestrial herbs, ferns | Coniferous forests, shrubs and terrestrial herbs, ferns | forest | Modern agricultural landscape |
4.5.1 Zonal mechanism
Parent materials, namely, fluvial deposits and loess[69], are concentrated in the warm and moist regions of the horizontal distribution zone. A previous study[2]demonstrates that the phaeozems are extensively dispersed across latitudes between 30° and 54°, with the highest distribution at 46°. According to that study, the temperate continental monsoon climate, characterized by long, cold, and dry winters as well as short, hot, and rainy summers, is conducive to the production and accumulation of organic matter in the environment. The current conditions are a result of the climatic conditions in the region. However, climate does not play a major role in the preservation and distribution of phaeozems[69]. The phaeozems formed from basalt deposits and volcanic deposits as parent materials are concentrated in the tropics[69,73−74], because of the strong chemical weathering and the presence of special minerals in the tropics. It is important to note that climate does not play a crucial part in this process[71–72]. The vertical zone mechanism is demonstrated by the distribution of phaeozems from different parent materials and the distribution geomorphology of phaeozems. All the phaeozems are situated on the mountain slopes or intermountain valleys with relatively elevated topography. The topography, climate, landform, and hydrology are found to be distinct at different altitudes. The phaeozems from the loess parent material are dispersed between 220 and 300 m, whereas those from the basalt residual parent material are distributed between 850 and 1150 m in the mountains. Based on this topography, the climate is vertically differentiated, and the vertical climate prevails. The topographical features such as inter-mountain depressions at the base of slopes are conducive to the storage of water from precipitation. Climate, geomorphology, and hydrology determine different environmental conditions at different altitudes. A combination of all these factors defines the occurrence, development, and growth of phaeozems at various altitudes.
4.5.2 Commonality in formation mechanism
The factors controlling the formation of phaeozems from diverse parent materials are extraordinarily complicated. Generally, the formation is the combined result of several distinct mechanisms, of which one or two are prominent.
According to the previous studies[44,69], the climate is warm and humid, and the degree of chemical weathering is low during the formation of the phaeozems from the loess parent material. If the conditions are warm and humid, grassland vegetation ensures an adequate source of organic matter and chemical weathering. The low degree of weathering causes mineral and organic matter in the soil to decompose slowly, thus providing a stable environment for the formation of phaeozems. Although the phaeozems from zeolite basalt parent material are formed in a tropical climate, the high degree of chemical weathering in a tropical climate is unfavorable for the formation of phaeozems. The presence of zeolite not only prevents the conversion of montmorillonite to kaolin but also maintains the alkali saturation level in the environment. This allows the preservation and accumulation of organic matter, thus leading to the formation of phaeozems in tropical climates[69,71]. The formation of phaeozems from volcanic deposits takes place because the clay minerals in volcanic deposits are unfavorable to soil infiltration. This helps in the storage of soil physical and chemical indicators and ensures the presence of organic matter. Thus, phaeozems were formed due to the combined effects of multiple factors[72]. In this study, the formation of phaeozems from river parent material is due to the swampy environment that ensures the source of organic matter and to the arid climate, which is the major contributing factor. In comparison of this study with the previous studies, it can be stated that the formation of phaeozems from various parent materials is a result of adequate sources of organic matter with common properties and a favorable environment for the preservation of organic matter. A comparative study also helps to verify the finding that the phaeozems are formed in warm and humid regions[44,69]; however, this may not be entirely accurate. The factors controlling the formation of phaeozems from various parent materials vary, and several mechanisms play essential roles in regulating the organic matter content.