The importance of non-carbonate mineral weathering as a soil formation mechanism within a karst weathering profile in the SPECTRA Critical Zone Observatory, Guizhou Province, China

Soil degradation, including rocky desertification, of the karst regions in China is severe. Karst landscapes are especially sensitive to soil degradation as carbonate rocks are nutrient-poor and easily eroded. Understanding the balance between soil formation and soil erosion is critical for long-term soil sustainability, yet little is known about the initial soil forming processes on karst terrain. Herein we examine the initial weathering processes of several types of carbonate bedrock containing varying amounts of non-carbonate minerals in the SPECTRA Critical Zone Observatory, Guizhou Province, Southwest China. We compared the weathering mechanisms of the bedrock to the mass transfer of mineral nutrients in a soil profile developed on these rocks and found that soil formation and nutrient contents are strongly dependent upon the weathering of interbedded layers of more silicate-rich bedrock (marls). Atmospheric inputs from dust were also detected.


Introduction
Carbonate rocks are widely distributed across China, covering a total area of *3,430,000 km 2 with the exposed area of the carbonate making up *13% of China's territory (Chinese Academy of Sciences 1979; Wang et al. 2014). The soils that form atop carbonate bedrock tend to be thin and highly susceptible to degradation, especially through subsurface erosion in the case of karst terrain. In the absence of significant aeolian input, formation of soils on carbonate bedrock is dependent on non-carbonate mineral ''impurities' ' [i.e. silicate and metal (hydr)oxide minerals] because (Ca,Mg)CO 3 dissolves completely during chemical weathering, leaving behind no residual mineral matrix to support soils. In a study of carbonate and non-carbonate rocks and associated soils in Guizhou Province, Southwest China, Chen and Bi (2011) found that non-carbonate rocks and soils have three times the trace element content of the carbonate rocks and soils and that higher value crops are grown in the non-carbonate regions.
Karst landscapes occupy 61.9% of the 30,240 km 2 of Guizhou Province (Chen and Bi 2011) and 25.2% of the province is affected by rocky desertification (Wang et al. 2014), creating serious economic and ecological problems (Huang and Cai 2007). Soil loss is the primary cause of karst rocky desertification (e.g. Wang 2003;Wang et al. 2014). The preferential loss of the most biologically active surface soil may significantly impact soil fertility, even in fertilized agricultural systems. The long-term sustainability of soils on a landscape is dependent on a balance between soil formation via weathering or deposition and soil loss 11th International Symposium on Geochemistry of the Earth's Surface.
& Heather L. Buss h.buss@bristol.ac.uk 1 via erosion. Despite this, the initial soil forming mechanisms in karst landscapes are poorly understood. Here we examine the chemical weathering mechanisms of bedrock in a karst watershed in Guizhou Province and the loss and gain of nutrient elements in the soil to determine the controls on soil formation and nutrient production in this sensitive environment.

Field site and sampling
The SPECTRA Critical Zone observatory (SCZO) is located within Puding County, Guizhou, Southwest China. The area has a humid sub-tropical climate, with a mean annual air temperature of 15.1°C (Wang et al. 2014).
Rainfall is monsoonal with a mean annual precipitation of 1378 mm. Within Puding County there is approximately 863.7 km 2 of exposed carbonate rocks, accounting for 79.2% of the county's total area. The rocks are mostly dolomite, limestone and contain some argillaceous minerals of Triassic and Permian age (Wang et al. 2014). Sampling was conducted in June 2016 in the Chenqi catchment in an abandoned pear orchard (105°46 0 17 00 E, 26°15 0 48 00 N WGS84 datum). A soil pit was dug to 70 cm, where bedrock was reached, and bulk soil samples were collected from each horizon. Due to a lack of accessibility to the different lithologies at the bottom of the soil pit, rock samples were taken from a nearby outcrop from each of three visibly different beds, including limestone and two marls (carbonate-rich mudstone, one with a shale-like, platy cleavage) both fresh and weathered (Fig. 1). The outcrop lies just below the pit, stratigraphically.

Laboratory analyses
Polished 30 lm thick petrographic thin sections were made of each bedrock lithology, with the samples only exposed to an anhydrous cooling fluid during preparation. A Hitachi S-3500N scanning electron microscope (SEM) and a Thermo Scientific 10 mm 2 Silicon Drift Detector were used to couple backscattered electron (BSE) imaging with 1024 9 1024 pixel resolution elemental maps, measuring 2.8 9 2.1 mm. Elemental maps were produced at 20 kV, rastering the electron beam whilst the SDD records X-ray counts at different excitation energies corresponding to different elements. The elemental X-ray maps were produced with 40 frames per map (frame time = 100 s). To determine mineral modal abundances, phase analysis was conducted on the X-ray maps using the Noran System Seven V3.2 software, by grouping together pixels of the scan area that display similar X-ray spectra (i.e., are chemically similar). Following lithium metaborate fusion, bulk rock samples were analyzed by inductively coupled plasma-optical emission spectrometry (ICP-OES) for solidstate elemental concentrations to support the results of the X-ray maps. Both the X-ray and ICP-OES results were normalized to exclude CO 3 . The solid-state elemental concentrations of bulk soils were also measured by ICP-OES. Mass transfer within the soil was determined by normalizing the measured solidstate concentrations of mobile elements to an assumed immobile element. The net mass loss or gain of the mobile component, relative to the parent material can then be determined by calculating the mass transfer coefficient, tau (Anderson et al. 2002): where (s) = mass transfer coefficient, tau; i = immobile component; j = mobile component; C = mass fraction, p = parent material and w = weathered material. When s i;j ¼ 0; no mobilization of j has occurred relative to the parent composition. When s i;j ¼ À1; component j has been completely lost and when s i;j [ 0; component j has been gained relative to the parent. Therefore, tau can be used to track the progress of weathering through a critical zone profile. Here we use the deepest sample in the soil pit as the parent material and assume that Ti is immobile (e.g. White et al. 1998).
3 Results and discussion

Bedrock composition
The limestone bedrock is composed primarily of calcite with some dolomite and minor quartz ( Table 1). The unweathered Marl A bed contains *75 vol% calcite, with significant proportions of K-rich and Fe-rich glauconites (Table 1). Marl A also contains minor sphene and apatite. The Marl B bed contains the lowest proportion of calcite (*67%) of all the beds at the outcrop and also contains K-rich and Fe-rich glauconite. Both marls show relative increases of quartz with weathering, suggesting it is conserved, and less of both glauconites (  Mineral abundances were determined using phase analysis on X-ray maps. Errors represent one standard error of the mean (n = 5) samples, with Marl A and B having similar concentrations ( Table 2). The limestone contains small amounts of Mg and Si due to the presence of dolomite and quartz (Table 1)

Soil composition
The soil profile has a mean bulk density of 1.16 gÁcm -3 , with most elements displaying partial depletion relative to the deepest sample (Fig. 2). Ca is totally depleted above *25 cm depth (Fig. 2a), supporting the observation from the bedrock analysis of a preferential loss of Ca from the system. To a lesser extent, Mg displays a similar trend to Ca, with *75% loss at the soil surface (Fig. 2a). Al, K and Na show initial losses of *50%, and Fe of 20%, in the deepest 20 cm, with little change in the overlying soil ( Fig. 2b-c). Mn is the only element in the profile that shows an increase relative to the parent material ([100%) towards the land surface (Fig. 2d). P shows a loss (*30%) in the deepest 20 cm followed by a slight increase, so that the shallow soil is only depleted relative to the deepest sample by *15% (Fig. 2d). Increases of P towards the surface could be due to (1) biolifting by vegetation, (2) application of phosphate fertilizer, or (3) atmospheric dust deposition. A study of dust in Jinzhong town, Guizhou Province, located 140 km away from Chenqi, attributed elevated Mn and P with local phosphorus mining activities (Yang et al. 2015), which may also be the source of Mn and P enrichments in our profile. Similarly, Tang et al. (2016) found that the dust in the Maolan Forest (located 241 km from Chenqi) was of local provenance mixed with anthropogenic sources.

Implications for soil formation and nutrient sources
The congruent dissolution of carbonate minerals leaves the weathered rocks comparatively enriched in silicate minerals, which are themselves relatively enriched in Al and Si due to preferential loss of Ca and Mg during incongruent dissolution. Ca and Mg are also preferentially lost during continued weathering in the soil profile. The residual silicate minerals form a structural matrix with which to support soil formation. We can use Al, a major component of silicate minerals that is not easily leached from soils, to estimate the potential of the different bedrock lithologies to produce a 70 cm thick soil. Using the rock and soil bulk densities (gÁcm -3 ) plus the solid-state Al concentrations (gÁg -1 ;  Elemental abundances, expressed as oxides, were determined using X-ray maps and compared to ICP-OES data to confirm the results, both are normalized to exclude CO 3 a Errors of the X-ray data represent one standard error of the mean (n = 5) and reflect heterogeneity of the rock 130 cm 3 , respectively, again with no erosive loss. Exposed rock in Puding County is *79% carbonate (Wang et al. 2014); assuming that the marls analysed here are representative of the remaining *21%, it is clear that although significantly less abundant than limestone in the catchment, the marl beds are disproportionately important in terms of soil formation. Similarly, natural sources of the nutrient elements in the profile, with the exception of Ca, must derive from weathering of the marl bedrock or dust input. The depletion profiles presented here (Fig. 2) are consistent with in situ soil formation from bedrock, although transport and deposition of soil and weathering products in other parts of the catchment are likely.

Conclusions
Both Ca and Mg are preferentially lost during weathering of all rock types and in the soil profile, leaving both the weathered rocks and soils relatively enriched in non-carbonate (silicate and oxide) minerals. Similarly, the low nutrient content of carbonate rock, which is due to the scarcity of non-carbonate minerals, means that the spatial variation in ecosystem productivity and environmental resilience may be strongly controlled by the distribution of relatively silicate-rich rocks. Although the limestone bedrock contains some of these minerals, they are likely insufficient to support significant soil formation and nutrient production. Therefore, the marl layers in the catchment represent an integral and critical part of the local karst ecosystem.