Isotopic and Geochemical Signatures of Dolostones and Their Implications for Carbonate Incipient Weathering Processes in the Datangpo Region, Guizhou, China

1. Introduction

Subtle variations in the global carbon cycle can lead to significant fluctuations in atmospheric CO₂ concentrations, thereby influencing the stability of the global climate[1,2,3]. To mitigate the greenhouse effect and combat global warming, it is essential to have a comprehensive understanding of the mechanisms that govern the global carbon cycle, including its sources, sinks, and their spatio-temporal distributions[1,2]. Previous research has extensively focused on terrestrial carbon sinks, particularly in relation to forest carbon sequestration and agricultural soil erosion sinks[4,5]. However, even after accounting for net forest carbon sequestration within terrestrial carbon sink systems, the global carbon cycle remains unbalanced, highlighting the well-documented phenomenon known as the ‘missing carbon sink’ [6,7]. Recently, carbonate weathering has emerged as a crucial component of the global carbon sink system, significantly influencing the regulation of global climate change [1,6,8]. Through biogeochemical processes such as the global water cycle and aquatic photosynthesis, the dissolution of carbonate rocks contributes to carbon sequestration over both short and long time scales, potentially addressing the gaps in the global carbon sink framework[2,6,8]. During the weathering-induced dissolution process, carbonate rocks absorb atmospheric CO₂. Compared to other minerals, carbonate minerals exhibit a higher susceptibility to dissolution, thereby providing substantial capacity for carbon sequestration[1]. The reaction equation of carbonate weathering is shown as follows: CaCO₃(s)+CO₂(g)+H₂O(l)⇌Ca²⁺(aq)+2HCO₃⁻(aq). This is a reversible reaction process. Under laboratory conditions, the reverse reaction occurs easily, releasing CO₂ again. However, in natural conditions, the addition of exogenic waters such as atmospheric precipitation and surface runoff dilutes HCO₃^- and Ca²⁺, preventing the reaction from reaching equilibrium[9]. During the weathering process, CO₂ is converted into the form of HCO₃^-, which is transported to the ocean, rivers, aquatic plants, and caves where it is sequestered, with only a small portion being re-released[10,11]. Carbonate weathering contributes approximately 0.5–0.9 Pg C yr⁻¹ to the global carbon cycle, which is about 0.45–0.82 times that of forest ecosystem carbon sinks[12,13]. Therefore, carbonate weathering is considered an important missing carbon source. As carbonate weathering serves as a vital pathway for reducing atmospheric CO₂ and maintaining equilibrium between carbon sources and sinks, understanding its entire processes is critical for research on global climate change and achieving carbon neutrality goals.

Extensive research has already been conducted on carbonate weathering, leading to significant advancements in understanding weathering rates, soil formation, and mineralization dynamics[8,11]. Traditional studies have relied on weathered soil profiles to elucidate elemental migration patterns during the weathering process[14,15]. Utilizing weathered soil profiles to reveal carbonate weathering is difficult because soil composition is complex. In addition to in situ weathering, soil may also be influenced by dust deposition, the transport of external water, and soil filtration. These factors can significantly hinder our understanding of carbonate weathering[14,15,16]. In addition, the intense weathering of carbonate rocks often makes it challenging to preserve the incipient weathering layer within soil profiles, thereby hindering researchers' ability to understand the incipient weathering process[14,15]. To our knowledge, there is a research gap in understanding the carbonate incipient weathering processes, as well as the associated elemental and isotopic changes. Previous researchers have not yet provided a clear definition of 'carbonate incipient weathering'. Given that carbonate incipient weathering is a crucial part of the overall carbonate weathering process, studying it is of great importance for understanding laterite formation and the role of carbonate weathering in carbon dioxide sequestration.

Guizhou Province, a paradigm of karst topography replete with carbonate rocks, is of particular interest. In the Datangpo region, mining and road construction have generated numerous artificial carbonate rock exposures. Notably, the dolomite lens of the Liangjiehe Formation, interbedded within sandstone, is relatively unaffected by other factors. Having been exposed to the air for merely a few years, it has formed a carbonate spheroidal weathering crust but has not yet weathered into soil. This makes it an ideal subject for studying the carbonate incipient weathering processes. The mineralogical and geochemical characteristics of spheroidal-weathered dolostones from the Liangjiehe Formation may be affected by carbonate incipient weathering, thereby documenting fluid sources and environmental alterations throughout these processes. With advancements in in-situ analytical techniques, we are now able to investigate the detailed mineralogy, geochemistry, and isotopic geochemistry of spheroidal-weathered dolostones, from core to rim, to elucidate the processes and mechanisms involved in carbonate incipient weathering.

Consequently, this study collected lenticular, spheroidal-weathered dolostones from the Liangjiehe formation in Zhaiying Town, Songtao County, Guizhou Province, China. Comprehensive investigations were conducted into the mineralogy, geochemistry, and isotopic geochemistry to explore the processes associated with carbonate incipient weathering.

2. Geological setting

The study area is located in the Datangpo region of Songtao County, Guizhou Province, China. Geologically, it lies at the intersection of the southeastern margin of the Yangtze Block, within the Nanhua Basin, and the Jiangnan Orogenic Belt (Figure 1a & Figure 1b). This extensional basin was formed as a result of the global breakup of the Rodinia supercontinent[17,18]. The basin contains a series of carbonate rocks, including rhodochrosite and dolomite[19]. The stratigraphic sequence in the study area is part of the Nanhua System (Figure 1b), which includes the following formations, listed from oldest to youngest: Liangjiehe Formation, Tiesiaoao Formation, Datangpo Formation, and Nantuo Formation[17]. An angular unconformity is observed between the Liangjiehe Formation of the Nanhua System and the underlying Banxi Group of the Qingbaikou System[17,20]. The Nantuo Formation of the Nanhua System is in para-unconformity contact with the Doushantuo Formation[21]. The low-temperature sedimentary sequence of the Nanhua System records two glacial events. The mixed rocks and sandstones of the Liangjiehe and Tiesiaoao Formations are indicative of the Sturtian glaciation (~720-663 Ma), which is characterized by poor roundness[22]. In contrast, the glaciomarine mixed rocks, siltstones, and sandstones of the Nantuo Formation document the Marinoan glaciation (~654-635 Ma), with thicknesses ranging from 60 to 200 meters[19,21]. The Datangpo Formation consists of interglacial sediments deposited between the Sturtian and Marinoan glaciations. It spans approximately 10 million years (~663-654 Ma), is conformably in contact with both the underlying Tiesiaoao Formation and the overlying Nantuo Formation, and has a thickness that varies between 30 and 700 meters[19].

3. Material and Method

The samples used in this study were primarily obtained from the Liangjiehe (Figure 1c). These samples are fragments from spheroidal-weathered dolostones. Within the sandstone of the Liangjiehe Formation, there is a lenticular spheroidal-weathered dolostone measuring approximately 126 cm in length and 32 cm in width (Figure 1f). Its surface displays a reddish-brown weathering crust with a laminated texture (Figure 2a & Figure 2b). The lenticular dolostone was located beside a mine-pit road and have been exposed due to the recent construction of the mine-pit road. The lenticular dolostone shows typical spheroidal weathering crust but has not yet weathered into soil. When struck with a hammer, it reveals three distinct layers from the inside out (Figure 2a): a bluish-grey layer (P zone), a yellowish-brown layer (S1 zone), and a reddish-brown layer (S2 zone).

3.1. Petrographic Observation

To observe mineral-phase transformations, samples include three-layer of spheroidally weathered dolostone were made into thin sections and analyzed using a microscope for mineral identification. Those sections were prepared using the Struers Labotom-15 cutting machine, selecting appropriate cross for sectioning. These were then flattened and polished with the Struers Tegramin polisher. Crystal Bond was applied to adhere the polished setions to glass slides, followed by cutting 200 μm sub-sections with the Struers Accutom-50 and final polishing to a thickness of approximately 40-60 μm.

After preparation, these sections were photographed using a stereo microscope (Leica EZ4W) to identify minerals species. Polarizing microscopy (Zeiss Scope.A1) was employed to observe and document the characteristics and positions.

3.2. Scanning Electron Microscopy (SEM)

To elucidate mineral phase variations in different spheroidal weathered crust layer, SEM-EDS analysis was performed. SEM experiments were conducted at the Guangdong Marine Resources and Coastal Engineering Laboratory, where gold-sputtered samples were mounted on conductive tape and placed in the evacuated chamber of a VEGA3 TESCAN field-emission SEM. Images were captured in backscatter mode at 20 kV and 1 nA, with a scan speed of approximately 90 s per image. EDS spectra, acquired for approximately 20 s each, aided in mineral identification.

3.3. X-ray Diffraction (XRD) Mineral Identification

XRD of powder samples determined mineral content changes in different spheroidal weathered crust layer. For powder samples XRD, about 50 mg of powder is drilled with a dental drill in-situ of spheroidal weathered crust layer. The powder is ground to below 200 mesh with an agate mortar, placed on a slide, flattened, and then tested by XRD. The powder samples XRD is performed on a PANalytical X'Pert Pro MDP with a Cu target, at 60 kV, 55 mA, covering 6- 50° (2θ), with a step size of 0.02°/s. XRD data is collected via X'pert Data Collector and analyzed using HighScore Plus for mineral phase identification.

Micro-XRD analysis was conducted to identify carbonate species changes in different spheroidal weathered crust layer. A fragment, containing P zone、S1 zone and S2 zone, was selected and sliced to a thickness about 100 mm. The slices were attached to glass fiber tubes. The rotation range of the platform was adjusted. Under the lens, the dolomite matrix in different spheroidal weathered crust layer was located in situ for XRD diffraction experiments.The D/MAX RAPID II XRD was used for micro-XRD, which equipped with a Mo radiation, operated at 50 kV and 30 mA, scanning diffraction angles from 2° to 44° at a step size of 0.02°/s. Diffraction patterns were converted to angle-intensity series using 2DP software, and mineral identification was performed with the PDXL software provided by Rigaku Corporation[23,24].

3.4. Trace and REEs Analysis

To elucidate the migration of elements in different spheroidal weathered crust layer, this study employed LA-ICP-MS to quantify the elemental contents, the points of LA-ICP-MS are shown in Figure 3b. Given that carbonate rocks primarily consist of Ca(Mg)CO₃, the major elements of focus were Ca and Mg. Additionally, considering the established correlation between carbonate rocks and manganese deposits in the Datangpo area, where carbonate lens occurrences are often accompanied by manganese ore layers[17,18], the study also investigated the contents of Mn and Fe. To track the changes in redox conditions and REEs content during the weathering process, the contents of Ni, Co, and REEs were also tested.

In-situ LA-ICP-MS analysis was performed. Samples were analyzed using the GeoLas Pro 193-nm ArF excimer laser ablation system coupled with an Agilent 7900 ICP-MS. Helium was used as the carrier gas, with a spot size of 32 μm, energy density of 5 J/cm², and repetition rate of 5 Hz. NIST610 and NIST612 served as external standards for trace element calibration every ten points, and data were processed using Iolite software.

3.5. Carbon and Oxygen Isotope Analysis

To facilitate a more comprehensive comparison and analysis of carbon and oxygen isotope, this study additionally collected unaltered dolomite (TJT1) in the study strata and carbonaceous mudstones (TZYY) from the Datangpo Formation strata for carbon and oxygen isotope testing (Figure 1c). The unaltered dolomite were collected to serve as a reference for primary mineral and geochemistry characters. The carbonaceous mudstones were sampled to investigate potential links between weathering and manganese mineralization (Figure 1d & Figure 2g & Figure 3a).

The positions of C-O isotope analysis are shown in Figure 3a-3c. Each samples, approximately 200 μg, were obtained by using a dental drill in situ drilled from section. These powders were loaded into glass tubes, evacuated with high-purity helium, and reacted with 100% pure phosphoric acid for 24 hours to convert it to CO₂[25,26,27,28]. The GasBench II system was used to introduce the gas into a Thermo Scientific 253 Plus isotope ratio mass spectrometer for analysis. Calibration was performed using the international standard NBS-18 (δ¹³C = -5.014‰, δ¹⁸O = -23.2‰), and results were reported in δ units (‰) relative to the Vienna Pee Dee Belemnite (V-PDB) standard.

3.6. Calculate Method

Migration Coefficient

To represent elemental accumulation and depletion, migration coefficient were calculated by using the formula :

RE = (X - Y) / Y

(1)

Where X represents the average elemental content in the current stage, and Y represents the average elemental content in pristine dolomite. Given the heterogeneity of carbonates and potential measurement errors associated with instrumentation, elements with an absolute RE value greater than 0.5 were defined as significantly depleted/enriched.

Cerium (Ce) Anomaly

The Ce anomaly, a commonly used indicator in REEs studies, serves to distinguish redox environments[29,30,31]. In this research, the Ce anomaly is calculated using the formula:

Ce/Ce* = 2Ce_N/(La_N + Pr_N)

(2)

Where N represents normalization values based on the Post-Archean Australian Shale (PAAS)[32].

4. Result

4.1. SEM and Microscope

Under SEM, P zone exhibits a denser texture with a dolomite matrix and occasional quartz grains (Figure 4a). When examined under optical microscope, pristine dolomite appears bluish-gray, with a fine-grained and compact microcrystalline dolomite matrix, intercalated with occasional quartz grains (Figure 4b).

SEM exhibits that dolomite in S1 zone has a relatively loose matrix composed of a mixture of dolomite and weathering residue, with occasional fractures filled with calcite (Figure 4c). Microscopically, dolomite in S1 zone is predominantly yellowish-brown, with a matrix comprising creamy fine-grained dolomite and yellowish-brown fine-grained weathering residues (Figure 4d). Occasional fractures are filled with creamy calcite, displaying yellowish-brown hues.

SEM analysis of dolomite in S2 zone indicates an even looser matrix, dominated by weathering residues with lesser dolomite and more frequent quartz grains (Figure 4e). Fractures are wider and primarily filled with weathering residues. Under microscope, dolomite in S2 zone appears reddish-brown, with quartz grains more prevalent, a reddish-brown weathering residues matrix, and occasional creamy dolomite remnants (Figure 4f). Fractures are predominantly filled with pale yellow material.

4.2. XRD Pattern

The powder XRD test results are presented in Figure 5a. By comparing the XRD data with standard mineral cards (PDF-2024), three minerals were identified: dolomite, quartz, and feldspar. The prominent peaks of dolomite (1 1 $\overline{4}$), (1 1 3), and (1 1 0) were all detected in the XRD spectra, along with the strong quartz peaks (1 0 1), (1 0 0), and (1 1 2). The intense peak of feldspar (0 0 2) was also present, confirming the presence of dolomite, quartz, and feldspar as the dominant minerals in the rock. Notably, the quartz peak intensities follow the trend P zone < S1 zone < S2 zone (Figure 5a), consistent with SEM and microscopic observations, indicating an increase in quartz content from P zone to S2 zone. As the degree of weathering intensifies increase, the content and peak intensities of feldspar increase.

In-situ XRD identified that the peak position of dolomite (1 1 $\overline{4}$) at different weathering stages didn't change significantly, which indicates that the dolostone matrix's mineral phase remained stable, with only dissolution occurring(Figure 5b).

4.3. LA-ICP-MS Data

In S1 zone, notable depletion was observed for Mg, Ca, Co, and Ni, while significant enrichment was evident for Mo (Figure 6a & Figure 7). In S2 zone, notable enrichment was observed for Mn, Mo, and Zn, while significant depletion was evident for Mg and Ca (Figure 6b & Figure 7). The REEs patterns in P zone, S1 zone, and S2 zone exhibit similar trends (Figure 8), with the ratio of light REEs (LREEs) vs. heavy REEs (HREE) varying minimally (P zone =3.74 ppm, S1 zone = 2.86 ppm, S2 zone = 3.21 ppm) (Table 1). This indicates that significant fractionation of REEs did not occur during the weathering process. The total REEs content exhibits a trend of initial decrease followed by an increase, as evident from the REEs distribution map (Figure 9). Specifically, higher REEs concentrations are observed in the P zone, which decrease in S1 zone, and then slightly increase again in S2 zone (Table 2). The overall Ce anomaly ranges from 0.76 to 1.99. Notably, a significant increase in the Ce/Ce* ratio is observed from P zone to S2 zone (Figure 7f).

4.4. Carbon and Oxygen Isotopes

The δ¹³C values of unaltered dolostone (TJT1) ranged from -6.07 to -2.76‰, with an average of -3.99 ‰, while the δ¹⁸O values ranged from -14.73 to -12.40‰, averaging -13.32‰. For the carbonaceous mudstones (TZYY), the δ¹³C values ranged from -9.65 to -7.99‰, averaging -8.55‰, and the δ¹⁸O values ranged from -13.63 to -10.35‰, averaging -11.68 ‰ (Table.3).

P zone exhibited δ¹³C values ranging from -5.26 to -3.35‰, averaging -4.18‰, and δ¹⁸O values ranging from -13.79 to -12.83‰, averaging -13.38 ‰, consistent with the isotopic signatures of the unaltered dolomite (TJT1)(Figure 10). S1 zone showed δ¹³C values ranging from -7.61 to -2.52‰, averaging -4.87‰, and δ¹⁸O values ranging from -12.22 to -8.06‰, averaging -10.49‰. S2 zone displayed δ¹³C values ranging from -16.56 to -12.43‰, averaging -14.17 ‰, and δ¹⁸O values ranging from -8.46 to -7.03‰, averaging -7.93 ‰ (Figure 11).

5. Discussion

5.1. Petrograph and Geochemistry Features Reveal Two Stages of Carbonate Incipient Weathering

The spherically weathered dolostone exhibits three distinct layers in hand specimens (Figure 1f & Figure 2a & Figure 3b). The innermost layer is bluish-grey and represents pristine dolomite (P zone). The middle layer is yellowish-brown (S1 zone), indicating dolostone that has just begun to weather. The outermost layer is reddish-brown, corresponding to dolostone that has been subjected to prolonged weathering (S2 zone). Besides color differences, the S1 and S2 zones also show significant differences in petrography, redox conditions, and acidity-alkalinity conditions.

Under the microscope, the S1 zone is relatively dense with few fractures, while the S2 zone is looser and more fractured (Figure 4). In terms of redox conditions, the Ce/Ce* ratio increases from the S1 zone to the S2 zone (Figure 7f). Ce is a redox-sensitive element predominantly existing as Ce³⁺ and Ce⁴⁺. In oxidizing aquatic environments, Ce³⁺ is oxidized to Ce⁴⁺, which is insoluble and tends to precipitate[33]. During carbonate weathering with fluid involvement, more Ce³⁺ precipitates from the fluid in more oxidizing areas, raising the local Ce content and thus the Ce/Ce* ratio. Thus, the increasing Ce/Ce* from the S1 to the S2 zone indicates a more oxidizing environment. From an acidity-alkalinity perspective, the S1 zone shows substantial Ni and Co loss (Figure 6a). The geochemical properties of Co and Ni indicate their solubility in acidic oxidizing environments and precipitation under alkaline conditions[34].This indicates that the S1 zone has an acidic environment. In the S2 zone, Ni and Co no longer show significant loss, but Mn shows significant enrichment (Figure 6b). Mn(II) is more readily oxidized to insoluble Mn(III) or Mn(IV) oxides under neutral or near-neutral conditions[35]. Thus, the S2 zone has a near-neutral environment.

Given the significant differences in rock properties, porosity, redox, and acidity-alkalinity conditions between the S1 and S2 zones, this study suggests taht they record two distinct stages of carbonate incipient weathering. In Stage I, carbonates interact with atmospheric CO₂ and water to produce weakly acidic carbonic acid (H₂CO₃). This reaction initiates the leaching of Ca²⁺ and Mg²⁺ from the carbonates via interstitial water, resulting in dolomite dissolution and structural disintegration. During this stage, elements such as Ni and Co are solubilized in the weakly acidic interstitial water and subsequently removed, leading to their depletion. Stage II is characterized by an increased production of HCO₃^⁻ as the reaction between carbonates, water, and CO₂ continues. This process gradually elevates the pH of the interstitial water as HCO₃^⁻ accumulation, which will hinder Ni and Co depletion. Concurrently, as the pH increases, which Mn(II) in the interstitial water reacts with oxygen to form manganese oxides or hydroxides.

5.2. Carbonate Incipient Weathering Indicates Lateritic Soil Formation Model

Carbonate weathering is often associated with laterite formation, yet the academic community still disputes the origin of laterite. Currently, four main hypotheses exist regarding laterite formation[36,37,38]:(1) In-situ weathering, where laterite forms as carbonate weathers in place, leaching carbonate minerals and leaving residual materials; (2) Dissolution-metasomatism, where laterite forms because carbonate reacting with external fluids rich in Al, Si, Fe, etc., forming clay minerals via metasomatism; (3) Clastic sedimentation, proposing that laterite is deposited by eolian or fluvial processes; and (4) Polygenetic theory, attributing laterite formation to a combination of the above factors.

This study, examining the petrology and geochemistry of spheroidal weathered dolostone, suggests that in the study area, laterite most likely originates from in-situ weathering of carbonate rocks. Supporting evidence includes: (1) Under microscopes and SEM, the rock matrix was observed to loosen and residual materials to accumulate with intensifying weathering, characteristic of in-situ weathering hypothesis (Figure 4). (2) Powder XRD analysis of P zone revealed over 50% dolomite, alongside significant quartz, feldspar, etc., indicating sufficient material in dolomite for laterite formation (Figure 5a). (3) No remarkable Al or Fe enrichment was detected in the S1 and S2 zones, implying limited metasomatic interaction between fluids and carbonate minerals (Figure 6). (4) The REE distribution patterns of the S1 and S2 zones resemble those of P zone, suggesting minimal external-source sediment input (Figure 8).

Thus, the overlying laterite on carbonate rocks in the study area is most plausibly formed by in-situ weathering of carbonate rocks.

5.3. C-O Isotopes Variation During Carbonate Incipient Weathering

In studies of stalagmites, oxygen isotope values of carbonate are believed to correlate with atmospheric precipitation, serving as indicators of monsoon rainfall and the effects of monsoon climate[39,40,41]. During the formation of carbonate rocks, oxygen isotopes predominantly inherit the isotopic signal from atmospheric precipitation[39]. Notably, variations in the oxygen isotopes of atmospheric precipitation display a significant latitude effect, with higher latitudes resulting in greater isotopic fractionation and consequently lighter oxygen isotope values[39].

The dolostone of this study were formed during the Neoproterozoic Era, a period characterized by the presence of numerous glacial deposits in the Yangtze Block. This suggests that the latitude of the Yangtze Block during the Neoproterozoic was higher than it is today[21]. Consequently, due to the latitude effect on atmospheric precipitation, carbonates formed during the Neoproterozoic are expected to exhibit lighter oxygen isotope values compared to modern carbonates. The range of oxygen isotope values in the P-zone of this study, range from -14.73 to -12.40‰, is indeed lighter than those of modern carbonates formed in Zhijin Cave (-10 to -8‰)[42], Guiizhou province. It is aligning with the isotopic characteristics of carbonates formed during the Neoproterozoic. In the S2 zone, the range of oxygen isotope values is -8.46 to -7.04‰, with an average of -7.83‰ (Figure 10 & Figure 11), which is consistent with the isotopic characteristics of modern carbonates from Zhijin Cave . These modern carbonates are in isotopic equilibrium with atmospheric precipitation, inheriting the isotopic signature of current atmospheric precipitation[42]. The S1 zone displays oxygen isotope values that are intermediate between those of P zone and S2 zone, indicating a gradual increase in weathering intensity from the S1 zone to the S2 zone, accompanied by a progressive alignment with the oxygen isotope values of modern carbonates. This suggests that variations in carbonate oxygen isotope values during weathering are influenced by atmospheric precipitation (Figure 11).

The carbon isotope composition of carbonates can provide insights into the composition and origin of the fluids involved in their formation[26,27,43]. During the formation of carbonates, the carbon isotopic signature is closely linked to the isotopic composition of dissolved inorganic carbon (DIC) in the fluids utilized, with minimal isotopic fractionation occurring during the precipitation of carbonates from DIC[44,45]. The carbon isotopic signature of DIC is primarily determined by its source. There are three main sources of DIC: (1) dissolution from surrounding rocks; (2) the reaction between atmospheric CO₂ and water; and (3) the reaction between CO₂ produced by plant and root respiration and water[1,46,47]. The pristine dolostone in this study is embedded within the sandstone strata of the Liangjiehe Formation (Figure 1c). Given that the surrounding rocks are sandstone, the influence of DIC from these rocks is considered minimal. Moreover, the dolostone was encased by sandstone of the Liangjiehe Formation, which suggests limited exposure to the atmosphere during their formation. This indicates that the DIC during the dolostone formation was primarily derived from the reaction between CO₂ produced by plant and root respiration and water. The CO₂ generated by vegetation and root respiration in the overlying strata reacted with water to form DIC, which was subsequently transported by pore water to precipitate as the dolostone. The similarity in REEs partitioning patterns between the pristine dolomite and pore water further supports this hypothesis (Figure 8).

The influence of vegetation on DIC is regulated by two factors: vegetation cover and the ratio of C₃ to C₄ plants. Increased vegetation cover results in a higher proportion of organic matter in the soil, which enhances respiration and organic matter decomposition[1,14]. This decomposition typically produces CO₂ with lighter δ¹³C values, leading to lighter δ¹³C values for DIC. Additionally, the C₃/C₄ plant ratio affects DIC values, as different plant types (C₃ and C₄) employ distinct photosynthetic pathways, resulting in the release of CO₂ with varying δ¹³C values[48]. C₃ plants (e.g., trees and most shrubs) exhibit lower δ¹³C values (-32 to -25‰), whereas C₄ plants (e.g., corn, sorghum, and prairie grasses) have higher δ¹³C values (-14 to -10‰). C₃ plants thrive in moist and cool climates, while C₄ plants are adapted to hot and arid environments. An increase in the C₃/C₄ plant ratio would consequently lead to lighter δ¹³C values for DIC[49,50].

During the Neoproterozoic Period, the Yangtze Block was characterized as a relatively cold region, exhibiting lower vegetation cover compared to contemporary conditions [51]. Consequently, carbonates formed during the Neoproterozoic should display heavier carbon isotope values than those of modern carbonates. The carbon isotope values of P zone (-6.07 to -2.76‰) are heavier than those of stalagmites from Zhijin Cave (-11 to -6‰), aligning with the isotopic characteristics expected for carbonates formed during the Neoproterozoic Period. In the S2 zone, carbon isotope values fluctuate between -16.56 and -12.43‰, with an average of -14.17‰ (Figure 10 & Figure 11). Atmospheric CO₂ typically has carbon isotope values ranging from -12 to -8‰. During rapid carbonate precipitation, carbon isotope values can be approximately 2‰ lighter than those of DIC[52]. Therefore, the carbon isotope signature of S2 zone’s carbonates is consistent with those expected from rapid precipitation influenced by atmospheric CO₂. S1 zone exhibits carbon isotope values that are intermediate between those of P zone and S2 zone (Figure 10 & Figure 11), indicating a gradual increase in weathering intensity from the S1 to the S2 zone, along with a progressive approach toward the carbon isotope values of atmospheric CO₂. This implies that variations in carbon isotope values during weathering are influenced by the carbon isotopic composition of atmospheric CO₂.

In summary, this study demonstrates that the oxygen and carbon isotope compositions of Neoproterozoic pristine dolomite lenses were governed by the oxygen isotope composition of atmospheric precipitation and the carbon isotope composition of DIC, which was influenced by the vegetation cover decrease during the Neoproterozoic Period. Due to the cold climate and high latitudes prevalent during the Neoproterozoic Period, carbon isotopes in carbonates were heavier than those formed in modern times, while oxygen isotopes were lighter. During incipient weathering processes, the oxygen isotope composition of carbonates was influenced by modern atmospheric precipitation, whereas the carbon isotope composition was affected by modern atmospheric CO₂.

5.4. Carbonate Re-Precipitation During Carbonate Incipient Weathering

During carbonate incipient weathering, a significant change in C-O isotopic values was observed in this study. The isotopes were influenced by atmospheric precipitation oxygen and atmospheric CO₂ carbon isotopes. If only carbonate dissolution occurred during carbonate incipient weathering, C-O isotopic values would remain stable. However, these changes were observed in this study, which indicates re-precipitation of carbonates during carbonate incipient weathering pprocesses. Atmospheric precipitation and CO₂ reacted with carbonates to form HCO₃^-, which then re-precipitated as carbonate minerals. This re-precipitation process may release CO₂, which hasn't been considered in the carbon-sink calculations of carbonate weathering. This could lead to an overestimation of the carbon-sequestration capacity of carbonate weathering and affect global carbon-reservoir calculations.

The proportion of re-precipitated carbonates in the initial weathering zone is neither too high nor too low. Powder XRD data shows that carbonates are barely detectable in the weathering zone (Figure 5a). This might be due to low carbonate content (< 5%) or the amorphous state of newly-formed carbonates. Geochemical data shows that Ca ions were lost by about 85% in the S2 zone (Figure 6), so newly-formed carbonates would account for no more than 15%. However, the detectable carbonate C-O isotopic signals in the MAT-253 Plus indicate that carbonate content isn't too low. Further research is needed to quantitatively describe the proportion of re-precipitated carbonates in the carbonate incipient weathering process.

6. Conclusion

In recent years, carbonate weathering has emerged as a significant global carbon sink, highlighting the critical need to investigate the mechanisms underlying the weathering of carbonate rocks. Research on these mechanisms primarily relies on weathered soil profiles; however, intense carbonate weathering often leads to the absence of the incipient weathering horizon within these profiles, which hinders a comprehensive understanding of incipient weathering processes and their regulatory mechanisms. In this study, we utilize surface-weathered dolostone to elucidate the processes of incipient weathering in carbonates. The following conclusions were drawn:

1. The δ¹³C and δ¹⁸O isotope ratios of pristine dolomite, ranging from -5.26 to -3.35‰ and -13.79 to -12.83‰, respectively, suggest that dolomite formation occurred under high-latitude, cold climatic conditions during the Nanhua Period.

2. The study identifies two distinct stages of incipient weathering in dolomite. In Stage I, there is a significant decrease of Ni and Co indicate the presence of a weakly acidic microenvironment. The fluctuating δ¹³C values (-7.61 to -2.52‰) and δ¹⁸O values (-12.22 to -8.06‰) , which indicate carbonate is re-precipitated. In Stage II, there is an enrichment of Mn which reflects a transition to a neutral environment. The δ¹³C values ranging from -16.56 to -12.43‰ and δ¹⁸O values from -8.46 to -7.03‰, which are influenced by atmospheric CO₂ and atmospheric precipitation.

3. The mineralogical, geochemical, and carbon-oxygen isotopic characteristics of the surface-weathered carbonate rock provide a systematic documentation of the incipient weathering processes of carbonate rocks, positioning it as a potential proxy for studying these processes.

This study focuses on the spheroidal weathered dolostone in Guizhou, China. However, the observed phenomena are not region-specific and are likely to be widespread. We believe that similar phenomena occur elsewhere on Earth. In this study, we observed carbonate re-precipitation during the incipient weathering of carbonate. But the proportion of re-precipitation carbonates wasn't accurately quantified. Future research will further explore this to better assess carbonate weathering contributions to the global carbon cycle. We also look forward to more studies on the spheroidal weathering of carbonate rocks from other regions to uncover the patterns of carbonate incipient weathering.