Trace Metal and Cd Isotope Systematics of the Basal Datangpo Formation, Yangtze Platform (South China) Indicate Restrained (Bio)Geochemical Metal Cycling in Cryogenian Seawater

: The behaviour of bioavailable trace metals and their stable isotopes in the modern oceans is controlled by uptake into phototrophic organisms and adsorption on and incorporation into marine authigenic minerals. Among other bioessential metals, Cd and its stable isotopes have recently been used in carbonate lithologies as novel tracer for changes in the paleo primary productivity and (bio)geochemical cycling. However, many marine sediments that were deposited during geologically highly relevant episodes and which, thus, urgently require study for a better understanding of the paleo environment are rather composed of a mixture of organic matter (OM), and detrital and authigenic minerals. In this study, we present Cd concentrations and their isotopic compositions as well as trace metal concentrations from sequential leachates of OM-rich shales of the Cryogenian basal Datangpo Formation, Yangtze Platform (South China). Our study shows variable distribution of conservative and bioavailable trace metals as well as Cd isotope compositions between sequential leachates of carbonate, OM, sulphide, and silicate phases. We show that the Cd isotope compositions obtained from OM leachates can be used to calculate the ambient Cryogenian surface seawater of the restricted Nanhua Basin by applying mass balance calculations. By contrast, early diagenetic Mn carbonates and sulphides incorporated the residual Cd from dissolved organic matter that was in isotopic equilibrium with deep / pore waters of the Nanhua Basin. Our model suggests that the Cd isotopic composition of surface seawater at that time reached values of modern oxygenated surface oceans. However, the deep water Cd isotope composition was substantially heavier than that of modern fully oxygenated oceans and rather resembles deep waters with abundant sulphide precipitation typical for modern oxygen minimum zones. This argues for incomplete recycling of Cd and other bioavailable metals shortly after the Sturtian glaciation in the redox stratiﬁed Cryogenian Nanhua Basin. Our study highlights the importance of sequential leaching procedures when dealing with impure authigenic sediments such as OM-rich carbonates, mudstones, or shales to achieve reliable trace metal concentrations and Cd isotope compositions as proxies for (bio)geochemical metal cycling in past aquatic systems. the relative enrichment of the element dissolved in the ambient bottom / pore water represented by the carbonate leachate, (2) the relative enrichment in undissolved OM and sulphates, and (3) the relative enrichment of the element leached authigenic and / or diagenetic sulphides. Our results showed that relative depletion of Ba, Zn, and Mo together with a more than 10-fold relative enrichment of Mn into carbonate phases argues for a reducing depositional


The Cryogenian Datangpo Formation and the Emergence of Substantial Algae Growth
The Neoproterozoic is one of the geologically most fascinating eras in Earth's history as it witnessed profound paleoenvironmental changes such as the breakup of the Rodinia supercontinent [1], various alternations between icehouse and greenhouse conditions [2,3], deep ocean oxygenation [4,5], and the first appearance of macroscopic multicellular organisms [6][7][8][9]. Within the Neoproterozoic Era, the Cryogenian Period is of particular interest for (geo)biologists as sedimentary deposits during nonglacial intervals already showed the first geochemical evidence for a redox-stratified ocean with an oxidized surface layer [10,11]. This paleoenvironment is widely discussed as the habitat set up for the emergence of the first animals in the following Ediacaran Period [12][13][14]. Furthermore, biomarkers from Cryogenian strata indicate that demosponges [15] and algae [16] became important in Earth's oceans by that time. These two organisms presumably led to subsequent oxygenation of the deep oceans by filtration and recycling of sinking organic matter (OM) as well as by photosynthesis, and, in the case of algae, as efficient nutrient recyclers and main food source for the emerging animal classes in the late Neoproterozoic following the Marinoan glaciation [17]. It was possibly the combination of extreme environmental conditions during the Cryogenian (low pO 2 and periodic widespread ice coverage of the ocean, blocking the light/energy source for phototrophic organisms) that triggered evolutionary inventions which took place in the following Ediacaran and Cambrian periods. However, the research community still lacks a general understanding of if nutrient conditions in the Cryogenian ocean were generally unfavourable for life and which role primary producers played during interglacial periods as these processes may have led the foundation for the following 'explosion' of life.
The marine sediments of the Datangpo Formation (Fm.) were deposited during the nonglacial interval between the presumably globally occurring, low-latitude Sturtian and Marinoan glaciations in the intra-Yangtze Platform Nanhua Basin [18]. These two glaciations are also locally known as the Tiesi'ao and Nantuo glaciation, respectively. Magmatic zircons sampled in ash beds at the base and top of the interglacial Datangpo Fm. have been dated to 663 ± 4 Ma (isotope dilution thermal ionization mass spectrometry (ID TIMS) U-Pb age) and 654.5 ± 3.8 Ma (sensitive high-resolution ionprobe (SHRIMP) U-Pb age) [19]. Magmatic zircons of ash layers in the overlying glacial Nantuo Fm. have been dated between 636.3 ± 4.9 Ma (SHRIMP U-Pb age [19]) and 635.2 ± 0.6 (ID TIMS U-Pb age [20]). The Nanhua Basin is a restricted basin that developed on the southwest to northeast striking passive continental margin of the Yangtze Platform that was occasionally connected to the open ocean during the Cryogenian Period [21]. Extensive sedimentological and geochemical data have been obtained from shales of the Datangpo Fm. including iron speciation, trace metal, and rare earth and yttrium (REY) abundances as well as stable C, N, S, Fe, and Mo isotope compositions [10,11,[21][22][23][24][25][26]. The combined geochemical studies suggest that the interglacial Nanhua Basin environment is best described by a stratified water column with anoxic and limited euxinic deep waters. However, deep waters temporarily changed to suboxic conditions associated with widespread dissolution of Mn-oxy/hydroxides (positive Ce anomalies [24]) being deposited on the basin floor which subsequently were transformed into diagenetic Mn carbonates under reducing pore water conditions [22,23]. Reducing pore waters further resulted in a near-quantitative Fe-oxy/hydroxide transformation into sulphides leading to relatively heavier Fe isotope compositions in the remaining seawater [11]. Additionally, a large, dissolved organic carbon pool existed shortly after thawing of the Tiesi'ao glaciers along with high denitrification rates that prevailed throughout the low oxygen conditions of the Nanhua Basin's deep waters [25,26]. Molybdenum isotopes of Datangpo Fm. shales show uniformly high δ 98 Mo values in the lower part and combined low and variable δ 98 Mo values in the upper part, respectively, suggesting euxinic seawater conditions in the lower part of the formation and variable redox conditions in the stratigraphically upper units [21].

(Bio)Geochemical Cycling of Trace Elements in Precambrian Oceans
The abundance and availability of trace elements in the Precambrian oceans were greatly different from those of Phanerozoic and modern oceans [27]. Factors controlling the trace elemental budget in the oceans are their local source inputs (riverine, aeolian, benthic flow, and hydrothermal influx), ocean mixing time, and element valency, controlling the elements' chemical behaviour and their availability to form complexes with inorganic or organic ligands or their adsorption onto particle surfaces (cf. [28]).
While Ni, Zn, and Cd are commonly present as divalent cations in seawater and form complexes with chloride (Cl − ), carbonate (CO 3 2− ), or oxi/hydroxides (OH − ), other transition metals such as U, Fe, V, Mn, Mo, Co, Cr, and Se exist in a variety of oxidation states which influence their (bio)geochemical cycling in the oceans by variable reaction rates, i.e., their partitioning into authigenic phases such as sulphides, carbonates, or oxi/hydroxides or their availability for complexation with organic ligands [28]. The concept of mean ocean residence time (MORT) includes all processes involved in the modern ocean trace element cycling as the residence time of an element, i.e., the average time from source to sink, depending on how fast the element can be adsorbed on, or incorporated into, a phase that is available for sedimentary deposition [29]. However, in the Precambrian oceans the MORT of most of the elements is unknown and availability of metals and coincident (bio)geochemical cycling may have been very different from modern environments [30][31][32]. During the Proterozoic and Archean, the oceans had different source contributions [33,34] and the Precambrian oceans were most likely stratified in terms of their oxygen availability with an oxidized surface layer and anoxic to euxinic (containing free H 2 S) deep waters [12,35]. Pristine, i.e., nondiagenetic overprinted, nonbiogenic carbonates are believed to be reliable archives for the geochemical compositions of the ambient fluid they precipitate from ( [36]). In a fully oxygenated modern ocean, the ability of phototrophic organisms to take up trace elements using special transport enzymes generally causes a relative depletion of enzymatically required trace elements in the photic surface ocean. Together with a relative enrichment of these trace elements in the deeper water column (due to the oxidative dissolution of sinking dead biomass), this creates a nutrient-like element concentration profile. Under stratified ocean conditions the uptake mechanisms of phototrophic organisms in the surface waters should act in the same direction leaving the upper water column relatively depleted in vital trace elements. However, the deposition and preservation of carbonates and OM under reducing conditions can happen under variable environmental conditions influencing the ability of marine sediments to record pristine geochemical signals of the ambient seawater: • Manganous condition. Under manganous conditions the dissolution of Mn oxy-hydroxides together with a decreased oxidative OM dissolution may ultimately lead to positive Ce anomalies in shale normalized rare earth element (REE) patterns together with enriched Mn concentrations and negative stable C isotope compositions in carbonates [37][38][39]. Authigenic Mn carbonates, i.e., rhodochrosite, form within the sediment body [40]. • Anoxic condition. Under truly anoxic (pO 2 <10 µmol) conditions, OM will not be recycled anymore and deposited in OM-rich shales showing common enrichments of the redox sensitive elements U and Mo [41,42]. By contrast, carbonates will be generally enriched in Mn and Ce due to the breakdown of Mn-oxides [30]. Further, redox sensitive trace metals, such as U, Re, V, and Mo but also Cd, have been shown to be relatively enriched in authigenic carbonate under reducing and anoxic conditions possibly due to supressed oxidative biogeochemical metal recycling [43]. • Euxinic condition. Under euxinic (pO 2 <10 µmol, free H 2 S) conditions in the deep water/pore water space authigenic sulphides (e.g., framboidal pyrite grains) or early-diagenetic sulphides (e.g., idiomorphic pyrite grains) can form, binding crucial trace metals such as Sn, Zn, Cd, Ni, and Cu [44][45][46][47].
To qualify the relative enrichments of the above mentioned elements in OM-rich sediments, the analysed bulk rock element abundances are usually normalized to global crustal rock reference materials [48,49] or to the detrital input from the local hinterland [50]. By analogy, the relative enrichment of trace elements in nonbiogenic carbonates has been normalized to abundances in Phanerozoic carbonate reference materials [37,51]. These methods have shown to be reliable proxies to assess the relative enrichment of vital elements in marine lithologies and (bio)geochemical cycling of these elements in the ancient ocean. However, these proxies are prone to diagenetic modification, especially in the case of redox sensitive and fluid-mobile elements, which are strongly influenced by intermittent oxidizing conditions in the pore water space, remobilizing trace metals and redistributing them within the sediment body [52]. To achieve a better understanding of the element cycling within the water column and the redistribution within the sediment body a holistic view at the level of each individual phase within a mixed sediment lithology is useful.

Status Quo of Cd Isotopes in Marine Lithologies
Cadmium isotopes in marine environments are a promising new tool for studying (bio)geochemical cycling of micronutrients in the oceans [53][54][55] and in sedimentary archives, such as Fe-Mn crusts [56], carbonates [57,58], microbialites [30], mudstones, and black shales [57][58][59]. Cd concentrations in seawater show a nutrient-like behaviour, with surface depletion and deep water enrichments. However, the processes leading to this distinct behaviour of Cd in modern seawater are yet to be fully understood. Utilization of Cd by marine photosynthetic organisms is the most favoured explanation, leading to a depletion of the surface water Cd pool [55]. As a possible driver for the Cd uptake, Price and Morel (1990) [59] described a new variant of the enzyme carbonic anhydrase (Cd-CA), which is active in the resorption of carbon dioxide during photosynthetic processes. In this enzyme, the usual Zn atom used for the catalysis is replaced by Cd. However, Cd-CA is only known from a few diatom species [60] which first appeared in the Mesozoic.
Cadmium concentrations in the water column generally vary in the opposite direction to the fractionation of stable Cd isotopes [53,54]. Thus, Cd-depleted surface waters have isotopically 'heavy' Cd signatures, while deep waters exhibit higher Cd concentrations but 'lighter' isotopic compositions. This is considered to reflect surface removal of light isotopes by marine photosynthetic organisms and their sinking and remineralization in bottom waters [54,55]. The isotopic composition of seawater is set by the balance of Cd inputs (and their isotopic composition) and the fractionation on removal to sedimentary sinks. Utilization of Cd by phytoplankton may be a major process controlling modern seawater Cd budgets and creating a depth gradient in Cd concentrations and isotopic composition [54,55,59]. Under incomplete oxidative OM recycling, biological uptake may act as an effective Cd sink and authigenic minerals may store the ambient seawater Cd isotope composition. Consequently, stable Cd isotope compositions show broad variations in ancient sediments that may be linked to variations in paleo productivity and redox state [30,57,58,[61][62][63]. Many more processes govern the Cd isotopic compositions of sediments than in seawater. Therefore, when applying Cd isotopes in ancient marine sediments, additional fractionation processes with different fractionation factors have to be evaluated very carefully. Besides biological utilization, partitioning into calcite [64] or authigenic sulphides [45,46] and variable uptake under changing salinity [57,65,66] may play a role. Unlike the aforementioned processes, inorganic adsorption onto Mn oxides has been shown experimentally to have no significant isotope fractionation at long time scales (>24 h) [67]. This makes Mn oxides, and potentially diagenetic Mn carbonates formed after reductive oxide dissolution (assuming the second order kinetic fractionation factor into carbonates), potential archives for recording the Cd isotope composition of contemporaneous fluids.
We here used a sequential leaching approach targeting the systematics of bioavailable trace metals and Cd isotopes in various lithologies from the basal Datangpo Fm., Southwest China, to extract organic and inorganic fractions. We applied a slightly modified sequential leaching procedure [68] on five Cryogenian OM-rich marine sediment samples (one black shale and four overlying OM-rich shales with abundant Mn carbonates and sulphide grains) to study the (bio)geochemical metal cycling and Cd isotope fractionation in the Cryogenian Nanhua Basin seawater. To better understand the vertical metal cycling in the ambient seawater by phototropic life, we here also present for the first time Cd isotope compositions in different leaching fractions from Cryogenian mixed lithologies.

Geological Sampling Location and Mineralogy
Samples for this study, were collected from the Cryogenian basal Datangpo Fm. from drill core zk2303, which was drilled by the Geological Team 103 in the Daotuo area in NE Guizhou Province, South China ( Figure 1). The Datangpo Fm. conformably overlies diamictites of the Sturtian Tiesi'ao Fm. [69] and hosts a variety of OM-rich sediments such as black shales, mudstones, and Mn carbonates that were deposited in the semi-restricted Nanhua Basin on the SE slope of the Yangtze Platform with periodical exchange to open ocean [70]. The basal Datangpo Fm. comprises 22 m of OM-rich black shales (members 1 and 3) which are intercalated by a 6-m-thick OM-rich manganese carbonate bearing layer (member 2) ( Figure 1). The black shales are stratigraphically followed by sandy shales of the middle Datangpo Fm. (member 4) and grey siltstones of the upper member 5 (not shown in Figure 1) [25,26]. The Datangpo Fm. is topped by the diamictites of the Nantuo Fm. indicative for the start of the Marinoan ice age [2].
Fm. [69] and hosts a variety of OM-rich sediments such as black shales, mudstones, and Mn carbonates that were deposited in the semi-restricted Nanhua Basin on the SE slope of the Yangtze Platform with periodical exchange to open ocean [70]. The basal Datangpo Fm. comprises 22 m of OM-rich black shales (members 1 and 3) which are intercalated by a 6-m-thick OM-rich manganese carbonate bearing layer (member 2) (Figure1). The black shales are stratigraphically followed by sandy shales of the middle Datangpo Fm. (member 4) and grey siltstones of the upper member 5 (not shown in Figure 1) [25,26]. The Datangpo Fm. is topped by the diamictites of the Nantuo Fm. indicative for the start of the Marinoan ice age [2].
In this study, leachates of five representative samples of the lower Datangpo Fm. were analysed. The samples consisted of four OM-rich shales hosting abundant Mn-rich carbonate layers (member 2) that are overlying one black shale sample of member 1. Overall, 30 analytes derived from 5 individual leaching steps of five samples and one duplicate were measured for their trace metal systematics and Cd isotope compositions. The black shale at the base of the Datangpo Fm. bore detrital minerals, OM, and authigenic sulphides, while scanning electron backscatter diffraction microscopy (EBSD) together with energy-dispersive X-ray spectroscopy (EDX) performed on a JEOL JXA-8230 equipped with an Oxford X-MaxN 20 IE250 at the State Key Laboratory of Marine Geology, School of Ocean and Earth Sciences, Tongji University, Shanghai, and by Xiao et al. (2017) [24] on the OM-rich Mn carbonates obtained rhodochrosite and calico-rhodochrosite minerals, detrital quartz, clays, and feldspar as well as framboidal pyrites of ~10 µm diameter. Manganese carbonate minerals were elongated and showed two generations with globular shaped minerals (~5 µm diameter) with high Mn content (>25 wt.%) in the centre surrounded by diffuse Mn beds with Mn contents <25 wt.% (see Supplementary Materials Figure S1). These Mn carbonates are closely associated to framboidal sulphides with individual grain sized between ~5-10 µm (see S1) and, less commonly, idiomorphic bigger sulphide grains (~30 µm).  [26]. Mn carbonate samples marked in stratigraphic profile as red circles, black shale sample is marked as a red diamond.  [26]. Mn carbonate samples marked in stratigraphic profile as red circles, black shale sample is marked as a red diamond.
In this study, leachates of five representative samples of the lower Datangpo Fm. were analysed. The samples consisted of four OM-rich shales hosting abundant Mn-rich carbonate layers (member 2) that are overlying one black shale sample of member 1. Overall, 30 analytes derived from 5 individual leaching steps of five samples and one duplicate were measured for their trace metal systematics and Cd isotope compositions. The black shale at the base of the Datangpo Fm. bore detrital minerals, OM, and authigenic sulphides, while scanning electron backscatter diffraction microscopy (EBSD) together with energy-dispersive X-ray spectroscopy (EDX) performed on a JEOL JXA-8230 equipped with an Oxford X-MaxN 20 IE250 at the State Key Laboratory of Marine Geology, School of Ocean and Earth Sciences, Tongji University, Shanghai, and by Xiao et al. (2017) [24] on the OM-rich Mn carbonates obtained rhodochrosite and calico-rhodochrosite minerals, detrital quartz, clays, and feldspar as well as framboidal pyrites of~10 µm diameter. Manganese carbonate minerals were elongated and showed two generations with globular shaped minerals (~5 µm diameter) with high Mn content (>25 wt.%) in the centre surrounded by diffuse Mn beds with Mn contents <25 wt.% (see Supplementary Materials Figure S1). These Mn carbonates are closely associated to framboidal sulphides with individual grain sized between~5-10 µm (see S1) and, less commonly, idiomorphic bigger sulphide grains (~30 µm).

Sequential Leaching Protocol for Trace Metal Extraction from OM and Mn Carbonate-Rich Shales
OM-rich shales deposited in restricted basin environments host a variety of different authigenic and detrital phases. The authigenic phases comprise carbonates and OM, leachable in strong alkaline solutions [71], as well as hydroxides, sulphates, and authigenic sulphides. The detrital phases comprise silicate minerals and colloids (such as ashes, dust particles, etc.). For the geochemical analyses we chose five bulk powdered samples of approximately 1 g from a sample suite that has been described in detail by Wei et al. (2016) [26]. We then performed a sequential leaching procedure ( Figure 2) modified after Henrique-Pinto et al. (2016) [68] to extract trace metals associated with a variety of host phases within the powdered rocks. This method, in contrast to traditional Aqua Regia/hydrofluoric acid (HF) bulk rock digestion methods, allows a sequential extraction of moderately volatile elements at low temperature that are bound to organic or colloidal phases and usually provides a better recovery than the total digestion methods of organic rich rocks [68]. Further, sequential leaching methods gave us the possibility to study fractionation of trace metals and its isotopes between different authigenic phases rather than studying bulk rock geochemical signatures that may be corroborated by signals from local detrital sedimentation background into a restricted basin environment.

Sequential Leaching Protocol for Trace Metal Extraction from OM and Mn Carbonate-Rich Shales
OM-rich shales deposited in restricted basin environments host a variety of different authigenic and detrital phases. The authigenic phases comprise carbonates and OM, leachable in strong alkaline solutions [71], as well as hydroxides, sulphates, and authigenic sulphides. The detrital phases comprise silicate minerals and colloids (such as ashes, dust particles, etc.). For the geochemical analyses we chose five bulk powdered samples of approximately 1 g from a sample suite that has been described in detail by Wei et al. (2016) [26]. We then performed a sequential leaching procedure ( Figure 2) modified after Henrique-Pinto et al. (2016) [68] to extract trace metals associated with a variety of host phases within the powdered rocks. This method, in contrast to traditional Aqua Regia/hydrofluoric acid (HF) bulk rock digestion methods, allows a sequential extraction of moderately volatile elements at low temperature that are bound to organic or colloidal phases and usually provides a better recovery than the total digestion methods of organic rich rocks [68]. Further, sequential leaching methods gave us the possibility to study fractionation of trace metals and its isotopes between different authigenic phases rather than studying bulk rock geochemical signatures that may be corroborated by signals from local detrital sedimentation background into a restricted basin environment. Step 1: Twenty mL of 1 M double distilled HAc was added to the sample powders in 50 mL centrifuge vials, the solution was stirred for 30 min and reacted overnight. Then the samples were placed in an ultrasonic bath for 30 min, centrifuged at 3500 rpm for 15 min, and finally the leachate L1 was pipetted out and dried down. The residual solids of step 1 remained in the vials.
Step 2: Residual solids of step 1 were treated with 20 mL 2 M double distilled HCl, stirred for 30 min, and placed in an ultrasonic bath for 15 min. The leachate L2 then was centrifuged at 3500 rpm for 15 min, pipetted, and dried down. The residual solids of step 2 remained in the vials.
Step 3: Twenty mL 1 M Suprapur ® NaOH were added to the residual solids of step 2. The solution was then stirred for 30 min and reacted for 4 hrs in a 60 °C water bath. After placing the solution in an ultrasonic bath for 15 min, it was then centrifuged at 3500 rpm for 15 min. Finally, the leachate L3 was pipetted and dried down. The residual solids of step 3 remained in the c vials.
Step 4: The residual solids of step 3 were treated with 7 mL Aqua Regia (4 mL conc. double distilled HCl and 3 mL conc. double distilled HNO3). The solution reacted for 24 hrs before the reaction vessels were placed in a water bath at 100 °C with slightly opened caps. Then, the leachate L4 was diluted with 20 mL Milli-Q, centrifuged at 3500 rpm for 15 min, pipetted, and dried down. The residual solids of step 4 remained in the vials. Step 1: Twenty mL of 1 M double distilled HAc was added to the sample powders in 50 mL centrifuge vials, the solution was stirred for 30 min and reacted overnight. Then the samples were placed in an ultrasonic bath for 30 min, centrifuged at 3500 rpm for 15 min, and finally the leachate L1 was pipetted out and dried down. The residual solids of step 1 remained in the vials.
Step 2: Residual solids of step 1 were treated with 20 mL 2 M double distilled HCl, stirred for 30 min, and placed in an ultrasonic bath for 15 min. The leachate L2 then was centrifuged at 3500 rpm for 15 min, pipetted, and dried down. The residual solids of step 2 remained in the vials.
Step 3: Twenty mL 1 M Suprapur ® NaOH were added to the residual solids of step 2. The solution was then stirred for 30 min and reacted for 4 hrs in a 60 • C water bath. After placing the solution in an ultrasonic bath for 15 min, it was then centrifuged at 3500 rpm for 15 min. Finally, the leachate L3 was pipetted and dried down. The residual solids of step 3 remained in the c vials.
Step 4: The residual solids of step 3 were treated with 7 mL Aqua Regia (4 mL conc. double distilled HCl and 3 mL conc. double distilled HNO 3 ). The solution reacted for 24 hrs before the reaction vessels were placed in a water bath at 100 • C with slightly opened caps. Then, the leachate L4 was diluted with 20 mL Milli-Q, centrifuged at 3500 rpm for 15 min, pipetted, and dried down. The residual solids of step 4 remained in the vials.
Step 5: In the final step the residual solids of step 4 were transferred into perfluoroalkoxy alkane (PFA) beakers and reacted with 2 steps of 5 mL of conc. double distilled HF and 1 mL conc. double distilled HNO 3 , heated to 120 • C for 24 hrs with closed caps, and dried down afterwards. To break down the remaining fluorides, the leachate (full digestion of remaining silicates) L5 was treated with repeated dry-down steps of small volumes of conc. HNO 3 .
Finally, all dried down sequential leachates L1 to L5 were redissolved in double distilled 3% HNO 3 and an internal 500 ppb Rh standard was added. The solutions were analysed on a Thermo Finnigan ICP-MS Element 2, using internal standardization and external calibration methods. We obtained concentrations of potentially bioavailable and partly redox-sensitive trace metals Co, Zn, Ba, V, Cr, Mo, U, Ni, Cu, Mn, and Cd as well as conservative behaving trace metals such as Zr, Ti, Sr, Nb, and Th. Limits of detection (LOD) were calculated as the signal of the analytical blank + 3x standard deviation (SD) of the blank. Repeated measurements of bulk digestions of the GBW07107 (GSR-5) and the USGS SDO-1 shale reference materials proved the methods' analytical accuracy was between 5% and 30% relative standard deviation (RSD) (see bottom of Table 1) relative to the concentrations reported in Govindaraju (1994) [72]. However, for V and Th distinct deviations (>50%) from the reference values were observed. Cd concentrations obtained in bulk analyses of the GSR-5 standard were 92% accurate compared to compiled literature values, whereas its concentrations obtained in bulk analyses of the certified standard reference material (SRM) SDO-1 were about 70% too high. Therefore, we obtained more precise Cd concentrations by isotope dilution (ID) and from here on only Cd ID concentration data is discussed. A linear correlation between Cd concentrations obtained by ICP-MS standard calibration and ID TIMS measurements was obtained, with Cd concentrations obtained by ICP-MS being in average 37% higher (Supplementary Materials Figure S4).

Stable Cd Isotope Analysis
Cadmium stable isotope analyses were performed after adequate Cd purification from four leaching steps, L1, L2, L4, and L5 (total Cd concentrations in L3 were below 5 ng and, therefore, below our detection limit for accurate TIMS analyses). Therefore, we used aliquots containing approximately 50-100 ng of Cd mixed with an optimal amount of a 106 Cd-108 Cd double spike solution, [62] dried, and re-equilibrated the mixture in ultrapure 0.5M HBr. The 106 Cd-108 Cd double spike was originally calibrated against an in-house JMC Cd Plasma solution (Lot: 15922032) at the Max-Planck-Institut für Chemie in Mainz, assuming 110 Cd/ 112 Cd = 0.520089 [73] for unspiked cadmium. The solutions were then loaded onto BioRad ® Polyprep columns filled with 200 µL AG1-X8 anion-exchange resin (100-200 mesh) to retain the Cd. The resin was rinsed repeatedly with 1 N HCl (to remove traces of bromide and matrix) followed by elution of the Cd with 0.25 N HNO 3 . For further purification, the eluted Cd was acidified to~0.5 N HBr and passed a second time through the columns in order to achieve suitable yields of Cd. Then, 5-50 ng of Cd were loaded onto zone refined degassed Re filaments and measured at~1150 • C in static multi-collection mode on a Thermo Scientific Triton Thermal Ionisation Mass Spectrometer (TIMS) at the Max-Planck-Institut für Chemie in Mainz. The data reduction was achieved by using an in-house coded double spike algorithm assuming exponential natural and instrumental fractionation law. Our procedural blank for the Cd chemistry was 120 pg/g, obtained by isotope dilution and was insignificant compared to Cd concentrations being at least 50 to 500 times higher in the samples. Following [74], variations in the Cd isotope compositions are expressed as ε 112/110 Cd values (deviations of 112 Cd/ 110 Cd in parts per 10,000 from the zero reference material (RM) of the National Institute of Standardisation (NIST) SRM 3108: Long-term reproducibility of spiked NIST SRM 3108 Cd isotope measurements at MPI Mainz is 0.4 ε 112/110 Cd units.

Stable C and O Isotope Analyses in Carbonate
Stable O and C isotope composition analyses in the phosphoric acid leachable carbonate fraction followed the protocol described in [75]. Therefore, approximately 100 µg of carbonate powder were reacted with 30 µL of 100% phosphoric acid at 70 • C for ca. 1.5 h. The O and C isotope ratios were then analysed in the liberated CO 2 gas using a Thermo Finnigan Gasbench II coupled with a Thermo Scientific MAT-253 mass spectrometer at the State Key Laboratory for Mineral Deposits Research at Nanjing University. As a reference gas, CO 2 calibrated against the Vienna Pee Dee Belemnite (V-PDB) standard was used. The reported isotope ratios were all given in the delta notation, as per mill deviation from the V-PDB standard [76]. The analytical precision, based on the reproducibility of our in-house CaCO 3 carbonate standard reference material (SRM) GBW04405, was better than 0.1% (standard deviation, 1SD) for both O and C isotopes.

Trace Metal Concentrations of Sequential Leachates
Trace element concentrations of five different leaching steps from the studied lithologies of the basal Datangpo Fm. showed little variation between each sample but varied systematically between the performed leaching steps and between Mn carbonate and shale lithologies (Table 1, Figure 3A). We here report the absolute obtained concentrations of trace metals that were considered important in supporting phototrophic life: Cd, Co, Cu, Ni, Zn, Ba, Mo, Ba, and Mn [28] ( Figure 3A). These concentrations were obtained on five sequential leaching steps (L1 to L5) and the relative distribution of each element between the different leachates is given in Figure 3B. Finally, relative concentrations of the bioavailable metals in L1, L2, L3, and L4 were calculated relative to concentrations in L5 (total digestion of remaining detrital silicates) representing the detrital nondissolved weathering input into the basin ( Figure 3C). The equation used to normalize the relative trace element concentration in each leachate to L5 is given below: norm. tr. metal = conc. Li tr. metal conc. L5 tr. metal (2) with i = 1 to 4.
Geosciences 2020, 1, x FOR PEER REVIEW 3 of 31 The element distribution in Cryogenian seawater in the Nanhua Basin can be represented by leachates of OM-rich shale and Mn carbonates building up the basal Datangpo Fm. which most likely precipitated under manganous to slightly euxinic bottom water conditions and under the presence of a massive DOC pool [11,[20][21][22][23]. Under reducing conditions many vital trace elements (TE), although potentially abundant in Nanhua Basin seawater (given the overlap in elemental concentration of L5 with UCC), were adsorbed onto sinking OM or bound in sedimentary sulphides or thiomolybdates (Zn, Ba, Mo) and not recycled in the water column. However, relative enrichments in Co, Cu, Ni, Cd, and Mn in the carbonate leachates reveal that some critical trace metals may have had the chance to be recycled and to remain in the water column and have been available for marine microbial life during formation of the basal Datangpo Fm. In leaching step one (L1, leaching of carbonates), average Cd and Mo concentrations were lower by an order of magnitude compared to Co, Ba, Zn, Cu, and Ni, while Mn was enriched by three orders of magnitude with an average concentration of 2723 µg/g (±1373 SD) and making up to 84% of the total Mn budget, which was to be expected for Mn-rich carbonates being leached in L1. Besides the Ti and Mn concentrations, which were five times higher in the Mn carbonates than in the studied shale sample (black diamonds in Figure 3B,C), all other trace element concentrations did not differ between the two lithologies. Cd concentrations in this leachate accounted for 14% to 51% of the total Cd budget ( Figure 3B). Normalizing the obtained concentrations to the L5 revealed clear relative enrichments of Cd, Ni, and Mn in the leached phases, while Ba, Zn, and Mo were relatively depleted ( Figure 3C). While in leaching step two (L2, leaching of organic matter), Cd, Co, Cu, and Mo as well as Ni showed similar low average concentrations (0.068 µg/g (±0.068), 2.4 µg/g (±1.9), 5.3 µg/g (±2.7), 8.8 µg/g (±5.8), and 3.9 µg/g (±21.8) (1SD)) with Cd concentrations making up 5% to 28% of the total Cd budget ( Figure 3B). The average concentrations of Zn, Ba, and Mn were significantly higher (23.4 µg/g (±9.8), 69.3 µg/g (±63.3), 260 8.8 µg/g (±109) (1SD)). For most elements no significant variations between lithologies were observed. Ba and Mn showed the lowest concentrations in the shale. Mn and Mo were relatively enriched against the local detrital background in the leached phases while Zn and Ba were depleted and Ni, Co, Cd, and Cu showed a ubiquitous enrichment behaviour ( Figure 3C). Leaching Step three (L3, leaching of sulphates) did not yield any Ni concentrations above the detection limit and also average Cd concentrations obtained by ICP-MS were the lowest within these leached phases (27 ng /g (±24, 1SD)). Cd ID TIMS data was not obtained on this leachate due to low concentrations. Average Ba concentrations were the highest in this leachate phase due to the dissolution of Ba sulphates with 1306 µg/g (±5193, 1SD). Ba concentration in L2 made up 17% to 75% of total Ba budgets ( Figure 3B). Concentrations of Co and Cu were low while Mn, Mo, and Zn showed higher abundances, with the latter three also revealing detrital background normalized enrichments ( Figure 3C). No significant trace metal concentration variations between lithologies were observed.
Leaching step four (L4, leaching of sulphides) revealed relatively high average Cd concentrations of (0.19 µg/g (±0.04. 1SD), which accounted for 36% to 66% of the total Cd budget ( Figure 3B). However, absolute Cd concentrations were about two orders of magnitude lower than the concentrations of Mo, Co, Zn, Ni, Ba, and Cu in the leached phases and up to three orders of magnitude lower than Mn. Concentrations of Ba were extremely low in the shale sample (310 µg/g) but rose to over 6000 µg/g at the top of the Mn carbonates, while Ti concentrations were the highest in the L4 of the shale lithology. Despite Zn and Ba, all other trace metals were relatively enriched over the local detrital background ( Figure 3C).
In leaching step five (L5, total digestion of the remaining detrital silicates) increasing average concentrations were obtained from Mo, Co, Ni, Cu, Zn, and Mn to Ba, matching concentrations reported for the average upper continental crust (dashed coloured lines [77]). Average obtained Cd concentrations for two samples were (0.03 µg/g (±0.02, 1SD)) denoting for 4% to 16% of total Cd budgets ( Figure 3B), which was lower than values reported for the upper continental crust UCC but in the same range as the SRM GSR-5, arguing for the leachate of detrital minerals.

Stable C, O, and Cd Isotope Compositions
The δ 13 C carb values in the basal Datangpo Fm. increased from −12.5 to −6.7% (Table 2) from base to top in the studied shale and Mn carbonates. The δ 13 C org (−31.8 to −31.5% ) values had almost no variation in these samples [26]. The δ 18 O carb value of the basal shale sample was -11.0% while the overlying Mn carbonates were between −9.3 and −08.4% (Table 2). We further obtained Cd isotopic composition in four out of five sequential leaching steps, excluding L3 due to very low Cd concentrations which was not suitable for precise analyses on the TIMS. All Cd isotopic compositions are expressed here as ε 112/110 Cd (Equation (1)). The 18 performed analyses (5 × L1, 5 × L2, 5 × L3 + 1 duplicate, and 2 × L5) showed a total variation of 7.5 ε units (Table 3, Figure 3A). Variations in ε 112/110 Cd of L1 ranged from −4.3 to +0.66, with the highest total variation and the heaviest value being recorded in the shale lithology (diamonds in Figure 3A). Values obtained in L2 ranged from +0.13 to +2.9 (again with the heaviest value obtained in the shale). The L4 values ranged from +0.59 to +3.2 (+1.5 in the shale). Duplicate analyses of one L4 sample obtained ε 112/110 Cd values of +0.59 and +0.99, respectively, (overlapping within the long-term reproducibility of the spiked NIST SRM 3108 Cd isotope measurements of 0.4 ε 112/110 Cd units). Due to low Cd concentrations in L5, only two reliable stable Cd isotope compositions were obtained, the shale sample having a ε 112/110 Cd of 0.02 and 42.5 ng/g Cd and the Mn carbonate sample having a ε 112/110 Cd of −0.34 and 15.5 ng/g Cd. These values overlapped in error with average loess compositions reported in [56] (ε 112/110 Cd relative to NIST SRM 3108 = +0.04 ± 0.31) as an up-to-date best fit for the composition of the upper continental crust.

The Leaching of Mixed Marine Lithologies: Bioavailable Trace Metal Distributions Associated with (In) Organic Phases
This section highlights and investigates the distribution of bioavailable trace metals Ba, Mo, Cd, Zn, Co, Mn, Mo, Ni, and Cu as well as Zr and Ti that are almost immobile elements clearly associated with detrital aluminosilicates between the five operational leaching phases L1 to L5. We postulate that the leaching step 1 attacked abundant Mn carbonate minerals such as rhodochrosite and calico-rhodochrosite that have been reported for the mixed lithologies of the basal Datangpo Fm. at Daotuo [24,26] and resulted in high absolute Mn concentrations obtained in L1. Further, elevated average Ni, Cd, Co, and Cu concentrations in the leachate (making up 31%, 30%, 20%, and 15% of the total budget, respectively) may be explained by the substitution of these metals for Ca 2+ into carbonate lattices [63,76,77]. In L1, concentrations of Ti and Zr were exceptionally low and, in the case of Zr, below the detection limit, arguing against leaching of any detrital material. L2 was supposed to leach remaining carbonates as well as the labile metal compounds bound onto the fulvic acid compounds of OM that were extractable by strong acids. The metal concentrations and distributions obtained were usually in the range of L1 with the exceptions of Mo (up to 9 µg/g and 20%) which can be attributed to Mo affinity to organic complexes. L3 described leaching of sulphate components and showed low trace metal concentration levels (below 1 µg/g). Only Zn and Ba were enriched and made up 34% and 55% of the bulk budget, respectively. L4 broke down sulphides, which were abundant as framboidal authigenic pyrite grains in the shales and of the lower Datangpo Fm. (see EBSD pictures in Supplementary Materials Figures 1 and 2 and the existing literature on the lower Datangpo Fm. [11,22,24]). The metal concentrations in L4 were high; however, distributions varied distinctively between the bioavailable metals: While concentrations of up to 0.27 µg/g made up 50% of the whole rock Cd budget, Zn concentrations of 19 µg/g only accounted for 12% of the whole rock concentration bound in sulphides. Further, concentrations of Mo, Ni, and Cu were high in L4 and accounted for the highest or second highest proportion of the bulk metal concentration. L5 fully digested the remaining detrital silicates reported for the basal Datangpo Fm. [24]. High Ti concentrations of 1114 to 17,049 µg/g argued for the dissolution of the remaining silicate minerals. The obtained metal concentrations obtained in L5 were highly variable throughout the different sample set. High variations within the sample suite were observed in Zn, Cd, Ni, and Cu where concentrations of leaching step 5 in sample SMS-11 exceeded all other samples by up to one order of magnitude.
It is notable that average Cd, Mo, Zn, and Ba concentrations of the L5 matched the reported values for the upper continental crust (UCC) [78] while Cu, Zn, and Mn only overlapped in their total obtained range with UCC values. Cu and Zn revealed slightly lower concentrations in the silicate leach than anticipated for the average UCC composition. However, the differences were small and may reflect different hinterland compositions distributing variable detrital material into the Nanhua Basin or UCC values that were based on bulk rock analyses which comprised a suite of different minerals including nonsilica phases. We, therefore, argue that the L5 leachate represents the best representative for the local detrital background in the Nahua Basin at the beginning of the Datangpo Formation.

The Redox State and (Bio)Geochemical Metal Cycling in the Cryogenian Yangtze Ocean
All analysed trace metal concentrations in sequential leachates from the Cryogenian Datangpo OM-rich shale and Mn carbonates were normalized to the local detrital hinterland composition to diminish the impact of detrital material on the chemical sediment´s signal. Hence, redox-sensitive element enrichments gave us information about (1) the relative enrichment of the element dissolved in the ambient bottom/pore water represented by the carbonate leachate, (2) the relative enrichment in undissolved OM and sulphates, and (3) the relative enrichment of the element leached authigenic and/or diagenetic sulphides. Our results showed that relative depletion of Ba, Zn, and Mo together with a more than 10-fold relative enrichment of Mn into carbonate phases argues for a reducing depositional environment in which Mn-oxides dissolved, releasing excess Mn into the ambient seawater. Together with prominent enrichments in Mn in the carbonate leach but coincident depletions in the sulphide leachate of Mn this argues for anoxic and manganous conditions in the bottom/pore waters at the deposition of the studied samples. The leached OM was generally enriched in the vital trace metals Cu, Mo, Mn, and Cd, arguing for adsorption onto OM or nutrient-like uptake in phototrophic organisms. By contrast, sulphate leachate L3, i.e., sedimentary barites and (oxy)thiomolybdates, showed depletion in Cu, Co, Mn, and Cd but enrichment in Zn, Ba, and Mo due to their deposition under slightly euxinic conditions in the pore water space [41]. Mo had low affinity to be incorporated into CaCO 3 , while it can be adsorbed onto positively charged organic particles and incorporated into (oxy)thiomolybdates or sulphides [41]. Consequently, our carbonate leachate L1 showed low Mo concentrations while the sulphide leachate L3 showed higher Mo concentrations, further arguing for euxinic pore water/bottom water conditions. This finding is further supported by the high δ 98 Mo values recorded at the base of the Datangpo Fm. [21]. The enrichments of Cd, Cu, Ni, and Co in the sulphide leach argue for incorporation under euxinic conditions being present in the ambient fluids, supported by the presence of framboidal sulphide grains (EBSD pictures in Supplementary Materials Figure S2).
Iron speciation and S isotope analyses [10] suggest that the early Cryogenian Nanhua Basin can be best explained by a vertically stratified redox model. However, variable U and Mo enrichments in lower Datangpo Fm. shales in combination with shale-normalized REY patterns showing distinct positive Ce and negative Y anomalies argue for at least episodic oxic to suboxic conditions in the basin that allowed for the formation of Mn oxides [22,24] as well as Fe oxides [11] somewhere in the water column. While iron speciation and S isotope analyses [10] showed that the early Cryogenian Nanhua Basin can be best explained by a vertically stratified redox model, variable U and Mo enrichments in lower Datangpo Fm. shales as well as shale normalized REY patterns having positive Ce and negative Y anomalies argued for at least episodic oxic to suboxic episodes in the basin that allowed for the formation of Mn oxides [22,24] as well as Fe oxides [11] somewhere in the water column.
The element distribution in Cryogenian seawater in the Nanhua Basin can be represented by leachates of OM-rich shale and Mn carbonates building up the basal Datangpo Fm. which most likely precipitated under manganous to slightly euxinic bottom water conditions and under the presence of a massive DOC pool [11,[20][21][22][23]. Under reducing conditions many vital trace elements (TE), although potentially abundant in Nanhua Basin seawater (given the overlap in elemental concentration of L5 with UCC), were adsorbed onto sinking OM or bound in sedimentary sulphides or thiomolybdates (Zn, Ba, Mo) and not recycled in the water column. However, relative enrichments in Co, Cu, Ni, Cd, and Mn in the carbonate leachates reveal that some critical trace metals may have had the chance to be recycled and to remain in the water column and have been available for marine microbial life during formation of the basal Datangpo Fm.

The Reliability of Cd Isotopes in Leachates from Mixed Lithology Samples
The use of authigenic mineral phases as archives for ancient environmental conditions strongly relies on the assumption that the pristine elemental budget of the mineral of interest is not masked by detrital contamination and/or post-depositional alterations during fluid-rock interactions. Cadmium isotopic compositions vs. Cd concentrations ( Figure 4) and O isotopic compositions (Supplementary Materials Figure S5) in the leached phases showed neither correlations with highly immobile elements (e.g., Al, Zr, Ti) that were clearly associated with detrital material nor with fluid mobile elements (e.g., Sr, O). Cd/Ti ratios found in nondetrital phases were three orders of magnitude higher than values reported in loess as a representative for the UCC composition (Figure 4a,f). Cd isotopic composition of the HF leaching step further overlapped with the upper continental crust (Figure 4f), indicating that the results of this leaching step most closely represented the local detrital input from the continental hinterland into the Nanhua Basin. However, both Cd/Zn ratios and Cd isotope compositions of carbonate, OM, and sulphide leachates substantially differed from the siliciclastic leachates and average UCC values, suggesting that detrital material had negligible impact on the elemental and isotopic budget of the carbonates. Furthermore, the lack of positive correlations between Cd concentrations and the isotopic compositions of all leached mineral phases with highly fluid mobile Sr and O isotopes (Supplementary Materials Figure S5) strongly suggest that post-depositional alteration during diagenesis or metamorphic processes did not affect the Cd budget of the here-studied mineral phases. Thus, carbonate fractions reliably traced the Cd budget from the fluid of which they were precipitated and can be used as prime geochemical archives to reconstruct (bio)geochemical metal cycling processes within the Nanhua Basin. Figure S5) in the leached phases showed neither correlations with highly immobile elements (e.g., Al, Zr, Ti) that were clearly associated with detrital material nor with fluid mobile elements (e.g., Sr, O). Cd/Ti ratios found in nondetrital phases were three orders of magnitude higher than values reported in loess as a representative for the UCC composition (Figure 4a,f). Cd isotopic composition of the HF leaching step further overlapped with the upper continental crust (Figure 4f), indicating that the results of this leaching step most closely represented the local detrital input from the continental hinterland into the Nanhua Basin. However, both Cd/Zn ratios and Cd isotope compositions of carbonate, OM, and sulphide leachates substantially differed from the siliciclastic leachates and average UCC values, suggesting that detrital material had negligible impact on the elemental and isotopic budget of the carbonates. Furthermore, the lack of positive correlations between Cd concentrations and the isotopic compositions of all leached mineral phases with highly fluid mobile Sr and O isotopes (Supplementary Materials Figure S5) strongly suggest that postdepositional alteration during diagenesis or metamorphic processes did not affect the Cd budget of the here-studied mineral phases. Thus, carbonate fractions reliably traced the Cd budget from the fluid of which they were precipitated and can be used as prime geochemical archives to reconstruct (bio)geochemical metal cycling processes within the Nanhua Basin.  [56], surface ocean [79][80][81], deep ocean [54], and hydro genetic Mn crusts [82]. Long-term reproducibility of spiked Cd isotope reference material NIST SRM 3108 is 0. 4ε units (2SD).

(Bio)Geochemical Cycling of Cd Isotopes in the Nanhua Basin
Cross-plots of Cd isotope compositions vs. Cd concentrations over trace element concentration plots ( Figure 4) show a systematic behaviour of isotope fractionation and element incorporation between the studied leachates and comparable marine and crustal rock standard reference materials. We here discuss that three general processes, i.e., (1) authigenic Cd enrichment trend, (2) Cd adsorption trend, and (3) biologic Cd enrichment trend, control incorporation, and fractionation of Cd and its isotopes in paleo-marine mixed sediment records. In Figure 4a  Loess data [56], surface ocean [79][80][81], deep ocean [54], and hydro genetic Mn crusts [82]. Long-term reproducibility of spiked Cd isotope reference material NIST SRM 3108 is 0. 4ε units (2SD).

(Bio)Geochemical Cycling of Cd Isotopes in the Nanhua Basin
Cross-plots of Cd isotope compositions vs. Cd concentrations over trace element concentration plots ( Figure 4) show a systematic behaviour of isotope fractionation and element incorporation between the studied leachates and comparable marine and crustal rock standard reference materials. We here discuss that three general processes, i.e., (1) authigenic Cd enrichment trend, (2) Cd adsorption trend, and (3) biologic Cd enrichment trend, control incorporation, and fractionation of Cd and its isotopes in paleo-marine mixed sediment records. In Figure 4a, a linear authigenic Cd enrichment trend (I) can be seen between Cd/Zn ratios and log [Cd/Ti] ratios. The trend steadily increases from a loess composition analysed by   [56] that we used as the best representative for a upper continental crust (UCC) via marine ferro-manganese nodules [82] and the Phanerozoic marine carbonate reference material CAL-S [30] towards our values obtained in individual leachates L1 to L5 from the basal Datangpo formation. We argue that this authigenic Cd enrichment trend (1) was mainly governed by the incorporation of Cd into OM in the water column and accumulation in authigenic phases in marine sediments [52] relative to the UCC and Mn nodules (hence authigenic nondetrital phases showed lower Ti concentrations). Two further trends can be observed: The Cd adsorption trend (2) is defined by a straight increase of Cd/Zn ratios with no relative variation in Cd over Mn, Cu, Ni, and Mo rations, respectively (Figure 4b-e). Even though Zn tended to be adsorbed onto ferro-manganese crusts under oxic conditions, its adsorption and enrichment factor relative to ambient seawater was less than those of the other bioavailable trace metals [83]. The positive-correlated biologic Cd enrichment trend (3) between Cd and other bioavailable trace metals, which can be seen for authigenic formed marine minerals such as sulphides and carbonates as well as for OM, was most likely controlled via Cd uptake. Further, adsorption of Cd onto clay minerals such as illite has been described [84], that may denote for the L5 leachate element ratios plotting on the Cd enrichment trend as well. Cd isotope compositions help us to further understand the processes of Cd enrichment in the leachates from OM-rich mixed marine sediments. Figure 4f-j show similar trends between Cd isotope compositions and Cd over trace metal ratios. In Figure 4f a positive correlation between ε 112/110 Cd vs. Cd/Ti indicates that an isotopic fractionation took place during the uptake of Cd into OM that further enriched the Cd relative to Ti in the biomass. The surface water in which phototrophic organism utilized the Cd as a result was extremely depleted in Cd and had a heavy Cd isotope composition of up to 4 to 6 ε 112/110 Cd values (light blue shade in Figure 4f-j) cf. [54,79,85,86]. The oxidative decay of sinking dead biomass led to a replenishment of Cd in deep water and further Cd fractionation into marine authigenic minerals such as (Mn) carbonates and sulphides precipitated under reducing and euxinic conditions, respectively [45,46,85,87]. The secondary formed minerals tended to have the same Cd/Ti ratios, however, they showed distinct lighter Cd isotope compositions than the (deep) ocean waters they formed in. This argues for second order fractionation of light isotopes into authigenic minerals that has been described for carbonates [64] and sulphides [45], respectively. While Figure 4h,i show the same trend of ε 112/110 Cd vs. stagnant Cd/Cu and Cd/Ni ratios (arguing for nearly identical incorporation behaviour of the metals into authigenic minerals), Figure 4g shows a negative correlation of ε 112/110 Cd with Cd/Mn, which can be attributed to the formation of Mn carbonates under manganous deep/pore water conditions in the Nanhua Basin during the deposition of the basal Datangpo Fm. Finally, Figure 4j shows the opposite trend with increasing Cd/Mo vs. decreasing ε 112/110 Cd, which may be attributed to the low affinity of Mo to be incorporated into carbonate lattice.

Assessing the Cd Isotope Composition of Cryogenian Nanhua Basin Surface and Deep Seawater
The Cryogenian Yangtze Ocean was situated at the banks of the Yangtze Craton which had a passive continental margin with a SW facing slope towards the open ocean. On this continental shelf environment, several protracted basins developed in the terminal Proterozoic that were periodically restricted from the open ocean via NE to SW striking tectonic barriers. While Cryogenian oceans already witnessed the rise of phototrophic algae [17], the temporally restricted exchange of waters in the Nanhua Basin led to stratified redox conditions with oxic to suboxic surface waters overlying anoxic to euxinic deep waters [10,25]. Under these circumstances, a large organic carbon pool built up in the deep ocean which was only periodically and partially oxidized [25], limiting large amounts of OM-associated bio-essential trace metals for nutrient recycling. Examples for modern analogues to such an environment of high phototrophic primary productivity with redox stratification are the Black Sea, some Norwegian fjords, the Peruvian Margin, and the Cariaco Trench. Although scarce, few existing Cd isotope datasets in these environments hint to combined effects of the above described biologic Cd uptake, causing Cd depleted 'heavy' ε 112/110 Cd SSW surface seawaters (SSW) and Cd enriched 'light' ε 112/110 Cd DW deep waters (DW) due to Cd release via sinking OM decay. The Cd uptake in modern closed and semi-closed marine systems such as the Southern Ocean and the Black Sea is described to follow simple Rayleigh-distillation models with fractionation factors α 112/110 Cd SSW−OM between 1.0001 to 1.0006 [79,85,88]. Further, laboratory experiments with Cd stock solutions and freshwater phytoplankton Chlamydomonas reinhardtii and Chlorella spirulina by Lacan et al. (2006) [55] supported fractionation factors in the same order. Figure 5a shows a simplified model for modern oxygenated oceans with constant fractionation of a the Cd pool in the surface waters under semi-closed to closed system Rayleigh-type biologic fractionation (i.e., with or without Cd replenishment) [88]. Here an averaged fractionation factor of α SSW−OM = 1.00035 was used to illustrate fractionation of Cd into biomass via biological uptake [79]. An additional fractionation factor of α SSW−carb of 1.000227 [64] for the kinetic uptake of the remaining Cd in the surface waters into biogenic carbonates, i.e., foraminifera tests, was assumed. For example, in the Cd supply-limited modern Southern Ocean's surface water the Cd isotope compositions were in the order of +4 to +6 ε 112/110 Cd units, while deep water compositions ranged between +1 and +2 ε 112/110 Cd [79]. Together, Cd-depleted surface waters and Cd-enriched deep waters mimic the typical seawater macro-nutrient phosphorous' depth concentration gradient (purple line Figure 5a) [89]. These patterns can be best explained by oxidative Cd release from sinking degrading OM, further leading to the recycling of Cd and potentially other bio-essential trace metals [80]. However, negative deviations of the analysed Cd content in deep waters from the expected concentrations have been linked to oxygen depletion zones (ODZ) (<100 µmol/L O 2 ) and the presence of free H 2 S [46]. Under sulphidic conditions, additional Cd isotope fractionation into Cd sulphides has been described in a H 2 S plume near the Peruvian Margin [85] as well as in the euxinic deep waters of the Black Sea [87]. Thus, these individual and combined fractionation factors have the potential to shift the typical "nutrient-like" ε 112/110 Cd depth gradient to heavier values in modern oxygen-depleted deep water environments (Figure 5a, dashed purple line).
under semi-closed to closed system Rayleigh-type biologic fractionation (i.e., with or without Cd replenishment) [88]. Here an averaged fractionation factor of = 1.00035 was used to illustrate fractionation of Cd into biomass via biological uptake [79]. An additional fractionation factor of of 1.000227 [64] for the kinetic uptake of the remaining Cd in the surface waters into biogenic carbonates, i.e., foraminifera tests, was assumed. For example, in the Cd supply-limited modern Southern Ocean's surface water the Cd isotope compositions were in the order of +4 to +6 ε 112/110 Cd units, while deep water compositions ranged between +1 and +2 ε 112/110 Cd [79]. Together, Cd-depleted surface waters and Cd-enriched deep waters mimic the typical seawater macro-nutrient phosphorous' depth concentration gradient (purple line Figure 5a) [89]. These patterns can be best explained by oxidative Cd release from sinking degrading OM, further leading to the recycling of Cd and potentially other bio-essential trace metals [80]. However, negative deviations of the analysed Cd content in deep waters from the expected concentrations have been linked to oxygen depletion zones (ODZ) (<100 µmol/L O2) and the presence of free H2S [46]. Under sulphidic conditions, additional Cd isotope fractionation into Cd sulphides has been described in a H2S plume near the Peruvian Margin [85] as well as in the euxinic deep waters of the Black Sea [87]. Thus, these individual and combined fractionation factors have the potential to shift the typical "nutrient-like" ε 112/110 Cd depth gradient to heavier values in modern oxygen-depleted deep water environments (Figure 5a, dashed purple line). (a) Schematic model for the Cd isotope 'nutrient-like' depth profile in the modern oceans. Cd isotopes are fractionated from the continental derived input in the photic zone via biological uptake leading to progressively 'heavy' Cd-enriched surface waters. This effect may be accompanied by fractionation of the Cd pool into carbonate tests with a further positive fractionation factor. Sinking dead biomass within the water column can be oxidized and the Cd gets recycled, causing modern deep waters to have 'lighter' Cd isotope compositions. In oxygen minimum zones, the formation of authigenic Cd (a) Schematic model for the Cd isotope 'nutrient-like' depth profile in the modern oceans. Cd isotopes are fractionated from the continental derived input in the photic zone via biological uptake leading to progressively 'heavy' Cd-enriched surface waters. This effect may be accompanied by fractionation of the Cd pool into carbonate tests with a further positive fractionation factor. Sinking dead biomass within the water column can be oxidized and the Cd gets recycled, causing modern deep waters to have 'lighter' Cd isotope compositions. In oxygen minimum zones, the formation of authigenic Cd sulphides under the presence of free H 2 S can cause further isotope fractionation and a shift to 'heavier' deep water values. (b) Schematic model of a theoretical Cd isotope profile in the redox stratified Cryogenian Nanhua Basin. Only Cd utilization by algae may fractionate the surface waters' Cd pool in the photic zone while carbonate shells were not existent. In the anoxic deep and pore waters of the Nanhua Basin the formation of authigenic carbonates and sulphides (represented by the obtained values in the L1 and L3 leachates) shifted deep water Cd isotope compositions that were in isotopic equilibrium with the DOC (recorded in the obtained values from OM leachate L2) to 'heavier' values that are above the values recorded in modern oxygenated deep oceans. In order to model the surface seawater for the stratified Cryogenian ocean at the time of the basal Datangpo Fm. deposition we assumed that the OM leachate L2 recorded the Cd isotope composition of the phototrophic biomass that fractionated the Cd pool in the Nanhua Basin photic zone. Perhaps the simplest approach would be to use the average value of the obtained ε 112/110 Cd values from L2 and add a constant average isotopic offset ∆ 112/110 Cd SSW-OM = +3.5 derived from average α 112/110 Cd SSW−OM obtained in modern surface ocean water. This calculation would place the ε 112/110 Cd SSW of the Nanhua Basin at~+4.5 which is in the order of the modern oxygenated surface ocean [54,79,81,85]. Figure 5b shows a whisker plot of the Cd isotope compositions obtained in three different operational leaching phases L1, L2, and L4 representing carbonate, OM, and sulphide leachates. Following [90], authigenic or early diagenetic (Mn) carbonates and sulphides in mixed OM-rich lithologies have the potential to record deep water/pore water compositions that are in equilibrium with dissolved organic carbon (DOC). From this carbon pool, and under the reductive dissolution of Mn oxy-hydroxides, inorganic carbonates such as Mn and Fe carbonates can precipitate under manganous/ferruginous conditions and, in the presence of free H 2 S, authigenic sulphides form within the sediment under euxinic conditions. EBSD pictures from the studied lithologies presume that thin horizons of Mn carbonates have indeed been formed in early diagenetic stage from colloidal Mn oxy-hydroxides (Supplementary Materials Figure S3). Further, the presence of Mn oxy-hydroxides as precursors for the studied Mn carbonates is underlined by rare earth element studies by [24] showing that Cryogenian Mn ores have similar REY patters but slightly enriched total REY concentrations comparable to modern ferromanganese crusts and nodules. The abundant framboidal sulphides (Supplementary Materials Figures S1 and S2) are closely associated to the Mn carbonate horizons arguing for their formation under progressively reducing/sulfidic conditions in the sediment. If the degrading OM was in isotopic equilibrium with the pore waters, instantaneous formed carbonates and sulphides would then be represented by an offset of ∆Cd fluid-carb = 2.27 and ∆Cd fluid-CdS = 1.6, respectively. This gave us the chance to estimate the Cryogenian ε 112/110 Cd DW of the Nanhua Basin at the time of deposition which was in the order of ε 112/110 Cd dw =~2.5. This value was significantly heavier than those values reported for modern oxygenated open ocean [79,81] although only few Cd isotope datasets from redox stratified modern oceanic environments exist. We here, therefore, present a modelling approach for calculation of theoretic surface seawater as well as deep water values for the Cryogenian Nanhua Basin that can be easily adapted to modern redox stratified marine environments using simple mass balance equations following John et al. (2019) [58] using the following two assumptions as prerequisite: (1) There is one dominant source of Cd into the oceans which is continental weathering leading to a homogenous input into the oceans of ε 112/110 Cd UCC =~0 [56] via riverine runoff and aeolian dust [91], and (2) three main Cd sinks exist in the Nanhua Basin which are OM, as well as authigenic and early diagenetic (Mn) carbonates and cadmium sulphides (CdS). It has been previously shown that the overall fractionation of Cd into a calcite/dolomite lattice under saline open marine conditions is in the order of ∆ 112/110 Cd fluid-carbonate = -2.27 [64] and insensitive to the cooperation of other major cations like Mg 2+ . Guinoiseau et al. (2018) [45] showed that Cd isotope fractionation into CdS under saline fully marine conditions is ∆ 112/110 Cd fluid-CdS = -1.6. Finally, the fractionation of Cd from the surface seawaters into phototrophic biomass is unknown and has to be assessed: The heaviest observed Cd isotope value in the Datangpo Fm. was about +2.9 Cd ε 112/110 Cd obtained in L2, arguing for an overall fractionation of Cd isotopes in the same order into OM. This value is only slightly smaller than the observed Cd fractionation in modern surface oceans (presumably governed by nonquantitative uptake into phototrophic organisms) which is in the order of 3.5 to 6.5 Cd ε 112/110 Cd [79]. This information leads to the calculation of the Cd isotope composition of Cryogenian surface seawater compositions at time of deposition using Equation (3) (4) where ƒ CdS + ƒ carbonate = 1. Given this equation, ε 112/110 Cd SSW (at a given amount of Cd uptake by phototrophic biomass) with ε 112/110 Cd DW under variable sulphide and carbonate formation in the deep water can be calculated. Figure 6A shows the evolution of the theoretical Nanhua Basin surface water vs. deep water Cd isotope compositions in coloured variable-sized circles. Increasing primary productivity shifts both the surface and the deep water ε 112/110 Cd to heavier values (colour change from red to green in Figure 6A). While increased carbonate over CdS precipitation (big vs. small circles in Figure 6A) shifted both water masses, Cd compositions to even higher values due to the bigger fractionation factor were associated with Cd incorporation into the carbonate lattice. Green dashed and dash/dotted lines in Figure 6a represent possible evolution of the instantaneous formed biomass ε 112/110 Cd OM values with α 112/110 Cd SSW−OM = 1.0002 [79] and 1.00035 [55]. We proposed earlier that deep waters of the Nanhua Basin can be considered to be in isotopic equilibrium with the DOC as long as no secondary authigenic mineral fractionation occurred that would have shifted the deep waters' composition to heavier values (orange and blue arrow in Figure 6A). Therefore, we applied two possible ∆ 112/110 Cd DW-OM of +2 and +3.5 to calculate the compositions of surface waters without authigenic mineral formation. The results were plotted as green arrays of 1.5 ε units' length (difference between the two ∆ 112/110 Cd DW-OM offsets) in Figure 6A, which overlap with potential OM surface water compositions derived from phototrophic primary production. Figure 6B shows the Cd isotope composition vs. Cd concentration obtained in the L1 and L4 leachates of authigenic carbonate and sulphide, respectively (light shaded symbols), and the Cd composition of theoretic ambient fluid they precipitated off. For simplicity, theoretic Cd isotope compositions of the ambient deep/pore water from which authigenic or early diagenetic Cd sulphides and Mn carbonates instantaneously precipitated were calculated without concentration changes of the fluid during mineral precipitation: ε 112/110 Cd fluid = ε 112/110 Cd CdS + ∆Cd fluid-CdS (5) ε 112/110 Cd fluid = ε 112/110 Cd carb + ∆Cd fluid-carb (6) Although the plotted ambient fluid ε 112/110 Cd fluid vs. log [Cd] concentrations obtained in the respective authigenic minerals plot on a linear trend, i.e., are directly correlated, it is unlikely that this phenomenon represents a typical Rayleigh-type fractionation, because it is opposite to that observed in the modern oxygenated Southern Ocean [79] or ancient microbial carbonates [30]. In both cases, the heavy ε 112/110 Cd fluid correlated with high Cd concentrations arguing for closed-system Rayleigh-type fractionation via biomass Cd uptake leading to a Cd depletion in the ambient fluid. In the analysed authigenic Mn carbonate and sulphide leachates from the basal Datangpo Fm., however, heavier ε 112/110 Cd fluid values correlated with higher Cd concentrations, arguing for Cd incorporation into the mineral lattices with nonlimited Cd supply via OM decomposition. In essence, our sequential leaching approach showed that it is possible to decipher the paleo deep and surface water Cd partitioning in mixed marine lithologies by analysing the authigenic phases presumably formed in variable depths of the water column. We have shown that OM in a reducing environment likely preserved the Cd isotope composition of phototrophic organisms taking up light Cd as a nutrient in the surface waters and leaving the surface water isotopically heavy. Assuming that the fractionation factor was in the order of modern ocean phytoplankton primary production, the restricted Nanhua Basin surface seawater ε 112/110 CdSSW = +4.5 matches with values of the modern surface open ocean [79] and the assumed surface water composition in the early Cambrian at the same location [90]. By contrast, the leachates of early diagenetic Mn carbonates and framboidal sulphides from the basal Datangpo Fm. represent the Cd isotopic composition of Cryogenian deep water for the restricted Nanhua Basin (ε 112/110 CdDW = +2.5), which was strikingly heavier than that of modern deep waters of the open ocean [54]. Applying the idea of Hohl et al. [90], we suggest that the decoupling of the DOC's Cd isotope composition and that of presumable ambient seawater was likely the result of fluctuating and progressively increasing deep water oxygenation in the early Cambrian that ultimately led to ε 112/110 Cd values in the deep and surface waters that are in the range of modern oceans. That period of extreme ocean redox instability led to perturbations in the nutrient cycle at the Precambrian Cambrian boundary that were necessary for the evolution of animal life. In essence, our sequential leaching approach showed that it is possible to decipher the paleo deep and surface water Cd partitioning in mixed marine lithologies by analysing the authigenic phases presumably formed in variable depths of the water column. We have shown that OM in a reducing environment likely preserved the Cd isotope composition of phototrophic organisms taking up light Cd as a nutrient in the surface waters and leaving the surface water isotopically heavy. Assuming that the fractionation factor was in the order of modern ocean phytoplankton primary production, the restricted Nanhua Basin surface seawater ε 112/110 Cd SSW = +4.5 matches with values of the modern surface open ocean [79] and the assumed surface water composition in the early Cambrian at the same location [90]. By contrast, the leachates of early diagenetic Mn carbonates and framboidal sulphides from the basal Datangpo Fm. represent the Cd isotopic composition of Cryogenian deep water for the restricted Nanhua Basin (ε 112/110 Cd DW = +2.5), which was strikingly heavier than that of modern deep waters of the open ocean [54]. Applying the idea of Hohl et al. [90], we suggest that the decoupling of the DOC's Cd isotope composition and that of presumable ambient seawater was likely the result of fluctuating and progressively increasing deep water oxygenation in the early Cambrian that ultimately led to ε 112/110 Cd values in the deep and surface waters that are in the range of modern oceans. That period of extreme ocean redox instability led to perturbations in the nutrient cycle at the Precambrian Cambrian boundary that were necessary for the evolution of animal life.

Conclusions
The paleo-environmental reconstruction of the Nanhua Basin during the Cryogenian Period yields surface water Cd isotope compositions obtained in OM leachates of ε 112/110 Cd SSW = +4.5 that overlaps with modern surface ocean and the assumed Nanhua Basin surface water composition in the early Cambrian. Deep water, represented by authigenic carbonate and sulphide phases, were in equilibrium with bottom and pore waters and show ε 112/110 Cd DW = +2.5, which is heavier than that of modern, oxidized, deep ocean water masses. In the Cryogenian Nanhua Basin, the decoupling of OM and deep-water compositions was yet not achieved and the deep waters' Cd isotope composition was most likely in an equilibrium with that of deposited and partly decomposed OM. As a result, the deep water and the surface water gradient in ε 112/110 Cd space was much smaller than in the Phanerozoic oceans and the supressed nutrient recycling may have ultimately led to the postponed evolution of animal life at the end of the Neoproterozoic.
Future research in redox-stratified Precambrian ocean sediments needs to keep these processes in mind and researchers should apply sequential leaching methods in order to obtain Cd isotope compositions of individual marine authigenic phases and the fluids that were in equilibrium with them.