Indication of Sr Isotopes on Weathering Process of Carbonate Rocks in Karst Area of Southwest China

: Based on the determination of the major and trace element content, and the Sr isotope composition of soils from limestone and dolomite proﬁles in the karst area of Southwest (SW) this study discussed the distribution and migration characteristics of the elements and the inﬂuencing factors of Sr isotope fractionation, in order to better understand the chemical weathering and pedogenesis process, as well as to explore the material source of laterite. The chemical weathering analysis results (Chemical Index of Alteration, mass balance coefﬁcient, and A-CN-K ternary) indicate that the weathering intensity of dolomite proﬁle PB is higher than that of limestone proﬁles CZ and CY. The correlation analysis between Sr isotopes and P 2 O 5 suggests that in addition to illite, apatite in the dolomite proﬁle also affects the 87 Sr/ 86 Sr composition. Sr content and 87 Sr/ 86 Sr demonstrate two stages in the weathering process of the carbonate, carbonate dissolution stage, and residual silicate weathering stage. As the carbonate minerals dissolve at the beginning of weathering, the Sr content decreases and 87 Sr/ 86 Sr increases slightly. After the decomposition of carbonate, the Sr content remains unchanged and 87 Sr/ 86 Sr increases. Finally, the study suggests that the material source of laterite is more likely to be underlying bedrock rather than the aeolian source. and the accumulation of K-bearing minerals is more obvious. In a A-CN-K ternary diagram the soil samples of dolomite proﬁle PB are located above the illite, corresponding to a higher CIA value > 90%, indicating that these samples have experienced strong weathering. Most of the soil samples of limestone weathering proﬁles are located below or close to illite (not beyond the scope of illite), and the corresponding CIA is 75–85%, indicating that the weathering of limestone soil samples has not yet reached the strong lateralization stage. (PB-3-PB-8),


Introduction
Chemical weathering is an important soil forming process, which will lead to the decomposition of rocks and minerals to form a weathering crust [1][2][3]. It involves complex chemical, mechanical, and biological processes, which can control the geochemical cycle of various elements, provide nutrients in the soil, and affect the global carbon cycle and longterm climate change [4][5][6][7][8]. Studying geochemical characteristics is helpful to understand the laws and control mechanisms of element release, migration, and distribution during the chemical weathering process [9,10].
Major and trace elements are powerful tools in the study of rock weathering and pedogenesis. Many chemical indicators composed of major elements are commonly used to evaluate the weathering degree [11][12][13][14]. Recently, proxies 4Si-M + − R 2+ , ∆4Si vs. R 3+ / (R 3+ + R 2+ + M + ) were proposed by Meunier et al. and were considered to be more accurate in distinguishing chemical weathering trends [11,15]. The immobile elements were utilized to address the material sources of sedimentary rocks, terra rossa, loess, and deep-sea sediments [16][17][18]. Due to the fact that rare earth elements (REE) have special chemical properties and unique system variations, they have been used to clarify the chemical weathering process of soil and sediments, such as mineral dissolution, and redox fluctuations [13,16,19].
Sr is an important trace component of rock-forming minerals [20][21][22][23]. Sr isotopic composition hardly undergoes any significant fractionation during geochemical processes, and its changes mainly reflect changes in material sources [24][25][26]; therefore, Sr isotopes have been used to study material sources of soils and terrestrial deposits, nutrient cycling in ecosystems, and water/rock reactions in surface or groundwater systems [27]. Julie et al. used Sr isotopes to constrain the contribution of soil mineral development, weathering fluxes, atmospheric deposition, and to elucidate the weathering processes of individual minerals in a small tropical granitoid watershed in the Luquillo Mountains of Puerto Rico [28]. Sr isotopes were also used by Santoni et al. to identify the water-rock interactions and mixing rates within a complex granite-carbonate coastal aquifer under high touristic pressure [23]. Each mineral has an individual 87 Sr/ 86 Sr ratio, and different weathering rates of the mineral cause variations in Sr isotopic composition during the weathering process; therefore, Sr isotopes have been used as tracers to address a wealth of geoscience issues in the weathering process of rock [29][30][31][32].
Numerous studies on rock weathering with Sr isotopes have been widely carried out [22,33]; however, most of them are focused on igneous rocks, and few have considered carbonate rocks [34,35]. The distribution area of red soil developed on carbonate rocks in the world is very extensive. The karst region of SW China covers approximately 42.6 × 10 4 km 2 and is one of the largest contiguous exposed carbonate distribution areas in the world [12,14]. The exploration of the pedogenesis process on the weathered carbonate rock in karst areas is the basis for revealing the geochemical cycle of elements in the karst environment and research on supergene resources; however, the material source and formation process of the red soil covering the bedrock is still controversial, and it has attracted wide and strong interest from soil scientists over the years [18,34].
In the southwest-central Guizhou Province, the widely distributed carbonate rocks are covered by thick red soil due to the humid and mild monsoon climate; therefore, the natural setting facilitates the study of carbonate rocks and lateritic soils. In this study, we investigated the Sr isotopic composition and element geochemistry of soil profiles to explore Sr isotopic characteristics and influencing factors in the soil profiles overlying dolomite and limestone, and to explain the material source of soil in this area.

Study Area and Sampling
The study karst area is located in the central part of the Guizhou Province, Southwest (SW) China ( Figure 1). The climate in SW China is humid, and the karst landform is widely developed in the Guizhou and Yunnan provinces. The karst terrain in this area is one of the largest continuously distributed karst landforms in the world. The study area has a subtropical monsoon climate, controlled by the East Asian summer monsoon, the Indian Ocean monsoon, and the Tibetan Plateau monsoon. The annual precipitation is 1100-1300 mm, and the rainy season is from May to October, accounting for 70-85% of the annual precipitation. The weather is relatively mild, and the annual average temperature is above 15 • C [36]. The PB profile (26 • [18,36]. The sampling points and regional lithology are shown in Figure 1. The vertical soil profiles were established due to road construction by an excavator. The surface soil of the profiles is loose and covered with lush vegetation. The samples were collected in order from the surface soil to the bedrock. A total of 9 soil samples and 4 bedrock samples were collected from PB profile. PB-1, PB-2, and PB-3 were collected at 5 cm, 10 cm, and 20 cm away from the ground, then soil samples were taken at every 40 cm and recorded as PB-4-PB-9. Four bedrock samples were collected at 290 cm, 390 cm, 590 cm, and 1090 cm, respectively, and were recorded as PB-10 to PB-13. CY-1 and CZ-1 were collected at 40 cm and 30 cm away from the ground, then 4 soil samples (CY-2 to CY-5) and 5 soil samples (CZ-2 to CZ-6) were collected at intervals of 40 cm and 30 cm, The vertical soil profiles were established due to road construction by an excavator. The surface soil of the profiles is loose and covered with lush vegetation. The samples were collected in order from the surface soil to the bedrock. A total of 9 soil samples and 4 bedrock samples were collected from PB profile. PB-1, PB-2, and PB-3 were collected at 5 cm, 10 cm, and 20 cm away from the ground, then soil samples were taken at every 40 cm and recorded as PB-4-PB-9. Four bedrock samples were collected at 290 cm, 390 cm, 590 cm, and 1090 cm, respectively, and were recorded as PB-10 to PB-13. CY-1 and CZ-1 were collected at 40 cm and 30 cm away from the ground, then 4 soil samples (CY-2 to CY-5) and 5 soil samples (CZ-2 to CZ-6) were collected at intervals of 40 cm and 30 cm, respectively. Subsequently, bedrock samples were collected at 250 cm and 210 cm away from the surface and recorded as CY-6 and CZ-7.

Analytical Methods
A total of 26 samples from three profiles were collected, including 20 soil samples and 6 bedrock samples. The soil samples were air-dried to remove impurities such as plant roots, debris, and conglomerates, and placed in a sealed polyethylene bag for storage, pre-

Analytical Methods
A total of 26 samples from three profiles were collected, including 20 soil samples and 6 bedrock samples. The soil samples were air-dried to remove impurities such as plant roots, debris, and conglomerates, and placed in a sealed polyethylene bag for storage, pretreatment, and analysis. Air-dried soil samples for whole-rock geochemical analysis were ground in an agate mortar and sieved through a 200-mesh (74 µm) light sieve.
The major elements were determined via X-ray fluorescence spectrometry (XRF) (PW2400 X-ray fluorescence spectrometer, PANalytical Company, Westborough, MA, USA) in accordance with the GB/T14506.28-2010 silicate rock chemical analytical procedure [37]. The trace elements and REE contents were analyzed using high-resolution inductively coupled plasma mass spectroscopy (HR-ICP-MS) (Element I, Thermo Fisher, Waltham, WA, USA) according to the GB/T14506.30-2010 [38]. An ICP-MS procedure was conducted at a temperature of 20 • C and a humidity of 30%. The analytical precision values for major and trace elements were less than 5% and 10%, respectively. The major and trace element data of all samples were provided by the Beijing Research Institute of Uranium Geology. Sr is determined by the isotope dilution method at the State Key Laboratory of Environmental Geochemistry, Institute of Geochemistry, Chinese Academy of Sciences. The powder samples (sieved through 200 mesh) were placed into a low-pressure digestion tank and the analytical grade HF + HNO 3 + HCl mixed acid was added to digest the samples. The dissolved sample is evaporated so that it becomes dry and converted to perchlorate, and 1.8-2.0 mol/L hydrochloric acid was used as elution acid. Sr isotopes were measured by multicollector Thermal Ionization Mass Spectrometry (TIMS) (IsoProbe-T, GV Company, U.K.). TIMS typically provides a value of the 87 Sr/ 86 Sr ratio of 0.710235 ± 0.000014 (2σ, n = 26) for the NBS987 Sr standard, and the measured isotope data are precise to 0.003% or better.
The Chemical Index of Alteration is widely used to measure the chemical weathering degree [39,40]. It can demonstrate the degree of feldspar decomposition and clay mineral variation in soils. The calculation method is: CIA = Al 2 O 3 /(Al 2 O 3 + CaO* + Na 2 O + K 2 O). Al 2 O 3 , CaO, K 2 O, and Na 2 O are all molar numbers of oxides, in which CaO* is the molar content of silicate minerals, excluding the content of CaO in carbonates and phosphates [41]. Usually when the molar number of CaO is greater than Na 2 O, n(CaO*) = n(Na 2 O) can be considered, whereas n(CaO*) = n(CaO) can be considered when the molar number of CaO is less than Na 2 O [41].
Mass balance coefficient (τ) is one of the best methods to indicate the depletion or enrichment of elements in soil samples, and it can provide information about REE mobility during geochemical processes [13]. The coefficient τ was calculated by the equation: where C j,w , C j,p and C i,w , C i,p represent the concentrations of element j and the relatively immobile element i in the soil samples and the bedrock, respectively. A positive τ i,j , value indicates that element j is enriched in soils relative to that in bedrock, whereas a negative value indicates the depletion. And zero of τ i,j indicates the immobility of the element j in soil. We used Th as the immobile reference element in this study.

Results
The vertical changes in the concentration of several selected elements and other geochemical parameters are shown in Figure 2. Multiple elements have a very low content in the bedrock samples (PB-10, CY-6, CZ-7) of three profiles, but Sr content is still relatively high. The concentration of Al 2 O 3 , Fe 2 O 3 , and P 2 O 5 in the dolomite profile PB is higher than in the limestone profiles CY and CZ, whereas SiO 2 , K 2 O, and Sr are lower. CIA is a common proxy for determining the weathering intensity. Figure (Table 3). At the same time, the PB profile has a higher degree of fractionation between light and heavy rare earth elements with (La/Yb) N = 10.49-13.10. In contrast, the content of REE in limestone profiles CY and CZ is lower than in the PB profile, as well as the fractionation degree between light rare earth elements (LREE, La, Ce, Pr, Nd, Pm, Sm, and Eu) and heavy rare earth elements (HREE, Gd, Tb, Dy, Ho, Er, Tm, Yb, and Lu). The content of REE in the CY profile ranges from 111. 18  CY profile, and 10.07-10.64 in CZ profile) shows that compared with HREE, the enrichment of LREE in the three profiles is higher. The chondrite-normalized REE pattern displays a right-leaning trend, which is also indicative of the enrichment of LREE (Figure 3). In the REE pattern, the PB profile shows a significant Ce anomaly (δCe = 0.85-1.7), whereas the Ce anomaly is not obviously in the limestone profile CZ and CY.

Major and Rare Earth Elements Indicate Chemical Weathering of Dolomite and Limestone
Chemical weathering is a process in which primary minerals are dissolved or transformed into other minerals and plays an important role in the process of soil formation. As a commonly used weathering indicator, CIA can directly reflect the degree of soil weathering [12,14]. It is generally believed that when 50 < CIA < 65 is reflected, the primary chemical weathering occurs under cold and dry climate conditions; when 65 < CIA < 85, it indicated intermediate weathering under warm and humid conditions; when CIA > 85, it is expressed as strong weathering under hot and rainy climate conditions [15]. The CIA value of the dolomite profile PB is higher than the limestone profiles. The CIA values of the soil A-CN-K (Al 2 O 3 -CaO* + Na 2 O-K 2 O) ternary diagram is usually used to study the weathering process, reflecting the changes of major elements and minerals, and indicating the trend of chemical weathering [42,43]. In the A-CN-K diagram (Figure 4), the bedrock is close to the CN apex and soil samples are close to the A-K line, and they gradually close to the A apex from the lower part to the upper part along with the profiles. CaO and Na 2 O in the bedrock are lost rapidly after the beginning of dissolution, which indicates the dissolution of feldspar minerals (plagioclase and potash feldspar) and the formation of clay minerals (illite). In the soil layer, illite further develops into kaolinite and K 2 O is gradually lost and Al 2 O 3 is enriched [43]; therefore, the weathering trend reflected in the A-CN-K diagram can be divided into two stages, stage I: the dissolution of carbonate minerals and leaching of soluble elements, and stage II: the formation and alteration process of clay minerals. Although dolomite and limestone profiles both have two stages of weathering, they show different weathering trends. In the first stage of weathering, dolomite profile PB shows that the weathering trend from bedrock to rock-soil interface is closer to the A-CN line, whereas the weathering trend from rock-soil interface to soil samples is closer to the A-K line and intense weathering products such as kaolinite, chlorite, and gibbsite [44,45]. The first stage in limestone profiles CY and CZ deviates more from the A-CN line than the PB profile, and the accumulation of K-bearing minerals is more obvious. In a A-CN-K ternary diagram the soil samples of dolomite profile PB are located above the illite, corresponding to a higher CIA value > 90%, indicating that these samples have experienced strong weathering. Most of the soil samples of limestone weathering profiles are located below or close to illite (not beyond the scope of illite), and the corresponding CIA is 75-85%, indicating that the weathering of limestone soil samples has not yet reached the strong lateralization stage. mary chemical weathering occurs under cold and dry climate conditions; when 65 < CIA < 85, it indicated intermediate weathering under warm and humid conditions; when CIA > 85, it is expressed as strong weathering under hot and rainy climate conditions [15]. The CIA value of the dolomite profile PB is higher than the limestone profiles. The CIA values of the soil samples in CY and CZ profiles are between 75-85, indicating that they have experienced relatively strong weathering, whereas the CIA values of almost all samples in the dolomite profile PB are > 90, indicating that they have experienced strong chemical weathering and leaching processes. These processes lead to large consumption of CaO, Na2O, K2O, and other soluble components in the study area.
A-CN-K (Al2O3-CaO* + Na2O-K2O) ternary diagram is usually used to study the weathering process, reflecting the changes of major elements and minerals, and indicating the trend of chemical weathering [42,43]. In the A-CN-K diagram (Figure 4), the bedrock is close to the CN apex and soil samples are close to the A-K line, and they gradually close to the A apex from the lower part to the upper part along with the profiles. CaO and Na2O in the bedrock are lost rapidly after the beginning of dissolution, which indicates the dissolution of feldspar minerals (plagioclase and potash feldspar) and the formation of clay minerals (illite). In the soil layer, illite further develops into kaolinite and K2O is gradually lost and Al2O3 is enriched [43]; therefore, the weathering trend reflected in the A-CN-K diagram can be divided into two stages, stage I: the dissolution of carbonate minerals and leaching of soluble elements, and stage II: the formation and alteration process of clay minerals. Although dolomite and limestone profiles both have two stages of weathering, they show different weathering trends. In the first stage of weathering, dolomite profile PB shows that the weathering trend from bedrock to rock-soil interface is closer to the A-CN line, whereas the weathering trend from rocksoil interface to soil samples is closer to the A-K line and intense weathering products such as kaolinite, chlorite, and gibbsite [44,45]. The first stage in limestone profiles CY and CZ deviates more from the A-CN line than the PB profile, and the accumulation of K-bearing minerals is more obvious. In a A-CN-K ternary diagram the soil samples of dolomite profile PB are located above the illite, corresponding to a higher CIA value >9 0%, indicating that these samples have experienced strong weathering. Most of the soil samples of limestone weathering profiles are located below or close to illite (not beyond the scope of illite), and the corresponding CIA is 75-85%, indicating that the weathering of limestone soil samples has not yet reached the strong lateralization stage. During chemical weathering, REE-containing minerals are dissolved, and REE is released into the soil and transferred through the soil solution [9,16]. The content of REE usually decreases with the increase in weathering degree [13,19]. According to the REE mass balance coefficient figure (Figure 5), both LREE and HREE in the soil samples of the PB profile have been greatly lost, and the τ is close to −1, which is consistent with the Sustainability 2022, 14, 4822 9 of 16 strong degree of weathering (CIA > 90). The τ of LREE and HREE in the limestone profile CY are concentrated at −0.5-0.8, indicating that the loss of REE is not as strong as that in the PB profile. Meanwhile, the REE in the CZ profile shows a slight depletion with a τ range from −0.5 to 0. All these indicate that the weathering degree of the dolomite profile is higher than that of the limestone. The behavior of the Ce element in the PB profile is different from other elements, so we calculated the Ce anomaly and showed the results in Figure 5a. There is a positive Ce anomaly in the upper part of the PB profile, δCe = 1.5-1.8 (PB-3-PB-8), the strong positive anomaly does not appear in the limestone profiles, and the Ce anomaly in the CY and CZ profiles is not obvious. leased into the soil and transferred through the soil solution [9,16]. The content of R usually decreases with the increase in weathering degree [13,19]. According to the R mass balance coefficient figure (Figure 5), both LREE and HREE in the soil samples of PB profile have been greatly lost , and the τ is close to −1, which is consistent with strong degree of weathering (CIA > 90). The τ of LREE and HREE in the limestone pro CY are concentrated at −0.5-0.8, indicating that the loss of REE is not as strong as that the PB profile. Meanwhile, the REE in the CZ profile shows a slight depletion with range from −0.5 to 0. All these indicate that the weathering degree of the dolomite pro is higher than that of the limestone. The behavior of the Ce element in the PB profile different from other elements, so we calculated the Ce anomaly and showed the results Figure 5a. There is a positive Ce anomaly in the upper part of the PB profile, δCe = 1.5-(PB-3-PB-8), the strong positive anomaly does not appear in the limestone profiles, a the Ce anomaly in the CY and CZ profiles is not obvious. Therefore, dolomite may have a faster weathering rate than limestone. The soil fr the dolomite profile has experienced more intense weathering and is in the intense later ization stage. In contrast, the limestone profiles have a lower degree of chemical weath ing and lateralization. Therefore, dolomite may have a faster weathering rate than limestone. The soil from the dolomite profile has experienced more intense weathering and is in the intense lateralization stage. In contrast, the limestone profiles have a lower degree of chemical weathering and lateralization.

Influencing Factors of 87 Sr/ 86 Sr during Soil Chemical Weathering Process
Sr isotopes are rarely affected by near-surface chemical, physical, and biological processes [44,45]; however, the weathering and decomposition of Sr-containing minerals may cause changes in the Sr isotopic composition [46][47][48][49]. The soil has three soluble reservoirs of Sr: carbonate minerals, clay minerals that adsorbed exchangeable strontium, and phosphorus minerals formed during diagenesis [50][51][52]. In the hot and humid subtropical climate of the Guizhou Province, the weathering profile contains few carbonate minerals, so the latter two are the main reservoirs of Sr in the study profiles. The 87 Sr/ 86 Sr ratios of the three profiles increase sharply from the rock-soil interface to the soil, then decrease upwards in the soil layer with an increasing weathering degree (Figure 2). This decrease is particularly obvious in the CY profile. The ratios of Rb/Sr and 87 Sr/ 86 Sr both gradually decrease upwards, which is contrary to the increase of Rb/Sr and 87 Sr/ 86 Sr in the granitic weathering profiles. This may be due to the existence of a large number of clay minerals that can absorb exchangeable Sr [32,53].
In the A-CN-K ternary diagram (Figure 4), the soil sample of the limestone profile falls in the area of illite, indicating that the main clay minerals in soil are illite and kaolinite formations. Illite is often regarded as a factor that affects the change of 87 Sr/ 86 Sr [31,35]. As a soluble element, the content of K 2 O in the laterite profile is closely related to the decomposition and formation of K-containing minerals such as K-feldspar and illite [23,53]. The K 2 O element content in the limestone profiles CY and CZ decreases with depth, which is caused by the gradual decomposition of illite due to the increase in weathering intensity. Due to the ionic radius of Rb being similar to that of K and the electricity valence being equal, K-containing minerals are usually rich in Rb. 87 Sr is derived from 87 Rb through β decay, so K-containing minerals usually have a higher 87 Sr/ 86 Sr [48]. In the soil layer, illite decomposes, resulting in 87 Sr/ 86 Sr becoming smaller. The good positive correlation between the K element and 87 Sr/ 86 Sr also indicates weathering, and the decomposition of illite controls the 87 Sr/ 86 Sr composition of the profile (Figure 6a). Sr isotopes are rarely affected by near-surface chemical, physical, and biological p cesses [44,45]; however, the weathering and decomposition of Sr-containing minerals m cause changes in the Sr isotopic composition [46][47][48][49]. The soil has three soluble reservo of Sr: carbonate minerals, clay minerals that adsorbed exchangeable strontium, and ph phorus minerals formed during diagenesis [50][51][52]. In the hot and humid subtropical mate of the Guizhou Province, the weathering profile contains few carbonate minerals the latter two are the main reservoirs of Sr in the study profiles. The 87 Sr/ 86 Sr ratios of three profiles increase sharply from the rock-soil interface to the soil, then decrease wards in the soil layer with an increasing weathering degree (Figure 2). This decreas particularly obvious in the CY profile. The ratios of Rb/Sr and 87 Sr/ 86 Sr both gradually crease upwards, which is contrary to the increase of Rb/Sr and 87 Sr/ 86 Sr in the gran weathering profiles. This may be due to the existence of a large number of clay miner that can absorb exchangeable Sr [32,53].
In the A-CN-K ternary diagram (Figure 4), the soil sample of the limestone pro falls in the area of illite, indicating that the main clay minerals in soil are illite and kaolin formations. Illite is often regarded as a factor that affects the change of 87 Sr/ 86 Sr [31,35]. a soluble element, the content of K2O in the laterite profile is closely related to the deco position and formation of K-containing minerals such as K-feldspar and illite [23,53]. T K2O element content in the limestone profiles CY and CZ decreases with depth, which caused by the gradual decomposition of illite due to the increase in weathering intens Due to the ionic radius of Rb being similar to that of K and the electricity valence be equal, K-containing minerals are usually rich in Rb. 87 Sr is derived from 87 Rb throug decay, so K-containing minerals usually have a higher 87 Sr/ 86 Sr [48]. In the soil layer, il decomposes, resulting in 87 Sr/ 86 Sr becoming smaller. The good positive correlation tween the K element and 87 Sr/ 86 Sr also indicates weathering, and the decomposition illite controls the 87 Sr/ 86 Sr composition of the profile (Figure 6a). The content of K in the dolomite profile PB is very low, and soil samples that fal the A-CN-K ternary diagram are above the illite, indicating that the main clay minera the soil is kaolinite, so there may be other factors affecting 87 Sr/ 86 Sr.
The PB profile is a typical carbonate weathering profile. It can be seen that CaO a good correlation with 87 Sr/ 86 Sr ( Figure 6b); however, in the strong weathering and lea ing environment, the carbonate rocks are leached and decomposed, and the CaO-conta ing mineral may not be carbonate rock but apatite [31]. The consistency of the Ce anom The content of K in the dolomite profile PB is very low, and soil samples that fall in the A-CN-K ternary diagram are above the illite, indicating that the main clay mineral in the soil is kaolinite, so there may be other factors affecting 87 Sr/ 86 Sr.
The PB profile is a typical carbonate weathering profile. It can be seen that CaO has a good correlation with 87 Sr/ 86 Sr ( Figure 6b); however, in the strong weathering and leaching environment, the carbonate rocks are leached and decomposed, and the CaOcontaining mineral may not be carbonate rock but apatite [31]. The consistency of the Ce anomaly and P 2 O 5 content change ( Figure 5), and the high P 2 O 5 content in the PB profile, indicates the possibility of apatite, and previous studies on Guangxi red soil also proved that the presence of apatite may be affected 87 Sr/ 86 Sr. We analyzed the correlation between P 2 O 5 and 87 Sr/ 86 Sr and found that only P 2 O 5 and 87 Sr/ 86 Sr in the PB profile have a good correlation, which proves that apatite is also a factor that affects the Sr isotope in the PB profile ( Figure 6).

Significance of 87 Sr/ 86 Sr for Evaluating the Soil Chemical Weathering Process
The Sr isotopic composition of the bedrock samples and the soil samples are different. The bedrock contains Sr in the main minerals, dolomite/calcite, which are dissolved during the weathering process, resulting in Sr leaching, which is the significant difference between soil and bedrock [26,54]. The Sr content and 87 Sr/ 86 Sr present two weathering steps during the weathering process of carbonate rocks (Figure 7). From the bedrock to the soil-rock interface, the Sr content decreases significantly and the 87 Sr/ 86 Sr ratios increase significantly; from the rock-soil interface to the soil samples, the Sr content changes a little, and 87 Sr/ 86 Sr increases moderately. It shows that the weathering process of carbonate rocks can be divided into two stages: Stage I represents the leaching process of carbonate minerals, and Stage II represents the weathering process of saprolite (Figure 7a). This may be related to two different sources of Sr in carbonate rocks. One is carbonate minerals with high Sr content and low 87 Sr/ 86 Sr, and the other is silicate minerals with low Sr content and high 87 Sr/ 86 Sr [34]. Carbonate rocks are more easily decomposed during the weathering process, so 87 Sr/ 86 Sr from the bedrock to the soil will suddenly increase.
Sustainability 2022, 14, x FOR PEER REVIEW 12 of and P2O5 content change ( Figure 5), and the high P2O5 content in the PB profile, indica the possibility of apatite, and previous studies on Guangxi red soil also proved that the pr ence of apatite may be affected 87 Sr/ 86 Sr. We analyzed the correlation between P2O5 and 87 Sr/ 8 and found that only P2O5 and 87 Sr/ 86 Sr in the PB profile have a good correlation, which prov that apatite is also a factor that affects the Sr isotope in the PB profile ( Figure 6).

Significance of 87 Sr/ 86 Sr for Evaluating the Soil Chemical Weathering Process
The Sr isotopic composition of the bedrock samples and the soil samples are differe The bedrock contains Sr in the main minerals, dolomite/calcite, which are dissolved du ing the weathering process, resulting in Sr leaching, which is the significant differen between soil and bedrock [26,54]. The Sr content and 87 Sr/ 86 Sr present two weatheri steps during the weathering process of carbonate rocks (Figure 7). From the bedrock the soil-rock interface, the Sr content decreases significantly and the 87 Sr/ 86 Sr ratios crease significantly; from the rock-soil interface to the soil samples, the Sr content chang a little, and 87 Sr/ 86 Sr increases moderately. It shows that the weathering process of c bonate rocks can be divided into two stages: Stage I represents the leaching process carbonate minerals, and Stage II represents the weathering process of saprolite (Figu 7a). This may be related to two different sources of Sr in carbonate rocks. One is carbona minerals with high Sr content and low 87 Sr/ 86 Sr, and the other is silicate minerals with lo Sr content and high 87 Sr/ 86 Sr [34]. Carbonate rocks are more easily decomposed during t weathering process, so 87 Sr/ 86 Sr from the bedrock to the soil will suddenly increase. Niobium (Nb) is a stable element and has almost no migration during weathering [5 57]; therefore, the concentration of Nb in the soil solution is very low, and the Nb/Sr ratio close to zero. As the degree of weathering increases, the Nb concentration and Nb/Sr ratio the soil are expected to increase [34]. For this reason, 87 Sr/ 86 Sr vs. Nb/Sr diagrams are co monly used to assess Sr endmembers and reveal the evolution of weathering profiles. The is a strong correlation between Nb/Sr and 87 Sr/ 86 Sr in the limestone profile CY and CZ (R 0.89,0.92), showing a trend that Nb/Sr and 87 Sr/ 86 Sr increase with the weathering (Figure 7 Nb/Sr and 87 Sr/ 86 Sr increase rapidly from bedrock to rock-soil interface, but the changes to t soil layer are not obvious. The dolomite profile PB has obvious two-stage weathering tren Fresh carbonate bedrock has a relatively low 87 Sr/ 86 Sr ratio and high Sr content. It is prefere tially weathered under warm and humid conditions [35]; therefore, clastic silicate miner end members with a higher 87 Sr/ 86 Sr ratio to control the Sr isotopic composition of saprolite. the first stage from the bedrock to the rock-soil interface, the 87 Sr/ 86 Sr ratio increased and t Niobium (Nb) is a stable element and has almost no migration during weathering [55][56][57]; therefore, the concentration of Nb in the soil solution is very low, and the Nb/Sr ratio is close to zero. As the degree of weathering increases, the Nb concentration and Nb/Sr ratio in the soil are expected to increase [34]. For this reason, 87 Sr/ 86 Sr vs. Nb/Sr diagrams are commonly used to assess Sr endmembers and reveal the evolution of weathering profiles. There is a strong correlation between Nb/Sr and 87 Sr/ 86 Sr in the limestone profile CY and CZ (R 2 = 0.89, 0.92), showing a trend that Nb/Sr and 87 Sr/ 86 Sr increase with the weathering (Figure 7b). Nb/Sr and 87 Sr/ 86 Sr increase rapidly from bedrock to rock-soil interface, but the changes to the soil layer are not obvious. The dolomite profile PB has obvious two-stage weathering trends. Fresh carbonate bedrock has a relatively low 87 Sr/ 86 Sr ratio and high Sr content. It is preferentially weathered under warm and humid conditions [35]; therefore, clastic silicate minerals end members with a higher 87 Sr/ 86 Sr ratio to control the Sr isotopic composition of saprolite. In the first stage from the bedrock to the rock-soil interface, the 87 Sr/ 86 Sr ratio increased and the Sr content decreased (Nb/Sr increased) (Figure 7b). As dolomite contains more clastic impurities (silicates), the degree of Sr accumulation in the dolomite rock-soil interface after the bedrock dissolves is higher than that of the limestone profile. In the second stage, carbonate is depleted (CaO concentration is usually <0.5 wt.%), and silicate components constitute the main body of the weathering layer. In this process, silicate minerals are affected by long-term strong weathering and transform into more stable minerals, such as kaolinite and gibbsite [44,55]. The exchangeable Sr existing in the silicate minerals is partially leached; therefore, the Sr and 87 Sr/ 86 Sr in the soil are reduced compared with the bedrock. Dolomite contains more silicate minerals and has experienced more intense weathering. Sr leaching is stronger than that in limestone soil samples, resulting in lower Sr content in the dolomite profile.

87 Sr/ 86 Sr Indication of Material Source of Soil Overlying Carbonate Rock
For a long time, the material source of laterite overlying carbonate rock has been controversial. They may have multiple origins, such as aeolian, and they may develop from clastic rock or in-situ weathering of carbonate rock. In SW China, some studies have also proposed the possibility that the overlying red soil comes from Chinese loess or Emeishan basalt; however, in our research, it is difficult to find evidence to support this hypothesis.
In the chondrite-normalized REE pattern, except for the strong Ce anomaly in the PB profile, the patterns of REE in the soil samples of the three profiles are very similar to the patterns in the respective bedrock samples, showing a trend of tilting to the right. This means that LREE are relatively enriched in HREE. This similarity indicates that the soil samples have a good inheritance from the bedrock and it supports the hypothesis of the homology relationship between the carbonate rock and the overlying weathered soil.
In the open superficial environment, Rb-containing minerals (mainly mica and potassium feldspar) are more resistant to decomposition than Ca (and Sr)-containing minerals (mainly feldspar and calcite) [24,58]. Due to different weathering rates of different minerals, the weathering process usually leads to changes in Rb and Sr concentrations and leaching rates of radioactive and non-radioactive Sr [25,29,59,60]. Generally, minerals with a high Rb/Sr ratio have good weather resistance and a high 87 Sr/ 86 Sr ratio; therefore, the ratios of Rb/Sr and 87 Sr/ 86 Sr of soil profile samples increase with the increase of weathering intensity. The Rb/Sr and 87Sr/86Sr ratios of the study profiles are significantly positively correlated (Figure 8a). The strong linear relationship indicates that the radiogenic 87 Sr is mainly derived from 87 Rb decay, and the weathering decomposition of Rb-containing minerals leads to the loss of Rb and radiogenic 87 Sr/ 86 Sr from the weathering system, whereas the contribution of the heterogeneous is very limited. Compared with profiles PB and CZ, the ratios of Rb/Sr and 87 Sr/ 86 Sr in profile CY are higher, which may be due to the weaker weathering degree and residual minerals with a high Rb/Sr ratio.
The study profile was compared with some similar lateritic weathering profiles in the Guizhou and Guangxi Provinces (Figure 8b). All of the data comes from published reports [18,34,35]. Some previous studies and our PB profile are in the same area (a distance of more than 10 km), with the same climate and similar stratigraphic age. Although these profiles and our profiles have similar ranges of Rb/Sr and 87 Sr/ 86 Sr ratios, their slopes and intercepts are not the same, indicating that their material sources are not a uniform aeolian material source, but should be evolved from the underlying bedrock.
In Figure 8b, we found that soil samples from Triassic and Permian bedrocks have different distribution ranges. The soils overlying the Permian bedrock have lower Rb/Sr and 87 Sr/ 86 Sr ratios while those overlying the Triassic bedrock have higher Rb/Sr and 87 Sr/ 86 Sr ratios, bounded by 87 Sr/ 86 Sr = 0.716. This difference may be due to the change of continental debris source in the South China Sedimentary Basin during the Permian to Triassic period [35]. The changes in the sources of these continents may be related to plate tectonics; therefore, the Sr isotope records in red soil may provide evidence for the movement of the South China block in the Late Paleozoic. mainly derived from 87 Rb decay, and the weathering decomposition of Rb-contain minerals leads to the loss of Rb and radiogenic 87 Sr/ 86 Sr from the weathering sys whereas the contribution of the heterogeneous is very limited. Compared with profile and CZ, the ratios of Rb/Sr and 87 Sr/ 86 Sr in profile CY are higher, which may be due to weaker weathering degree and residual minerals with a high Rb/Sr ratio.  [18]; data of profiles T-QZ1, T-QZ2, and P-QZ3 are from Liu et al. [34]) and the Guilin province profile data (P-PG and P-BY1 are from Wei et al. [35]).

Conclusions
This study assessed the Sr isotope composition and weathering processes of the carbonate rock (dolomite and limestone) weathering crust. The weathering degree of the dolomite weathering profile PB was higher than that of the limestone weathering profiles CZ and CY. A-CN-K ternary suggested the main clay mineral in the CZ and CY profiles is illite, and the clay minerals in the PB profile are changing to gibbsite and kaolin. The mineral composition had a significant impact on the Sr isotopic characteristics. The decomposition of illite and apatite in the carbonate rock profiles could affect the Sr isotopic composition in soils. Two stages existed in the weathering process of the carbonate rock profiles. In the first stage, the Sr content decreased and 87 Sr/ 86 Sr slightly increased due to the carbonate minerals rapidly dissolving from the bedrock to the rock-soil interface. In the second stage, the main body of the weathering process from the rock-soil interface to the soil layer is composed of clastic silicate minerals, which resulted in the Sr content slightly decreasing and 87 Sr/ 86 Sr increasing significantly. The ratios of Rb/Sr and 87 Sr/ 86 Sr divided the soil profiles into different groups, indicating that their parent materials are different and the possibility of exogenous soil origin (aeolian) of the region is very low; therefore, the red soil overlying the carbonate rock is more likely to be developed from the underlying bedrock.