Characteristics of soil profile CO2 concentrations in karst areas and their significance for global carbon cycles and climate change

CO2 concentrations of 21 soil profiles were measured in Zhaotong City, Yunnan Province. The varying characteristics of soil profile CO2 concentrations are distinguishable between carbonate and noncarbonate areas. In noncarbonate areas, soil profile CO2 concentrations increase and show significant positive correlations with soil depth. In carbonate areas, however, deep-soil CO2 concentrations decrease and have no significant correlations with soil depth. Soil organic carbon is negatively correlated with soil CO2 concentrations in noncarbonate areas. In carbonate areas, such relationships are not clear. This means that the special geological process in carbonate areas – carbonate corrosion – absorbs part of the deep-soil-profile CO2. Isotope and soil pH data also support such a process. A mathematical model simulating soil profile CO2 concentration was proposed. In noncarbonate areas, the measured and the simulated values are almost equal, while the measured CO2 concentrations of deep soils are less than the simulated in carbonate areas. Such results also indicate the occurrence of carbonate corrosion and the consuming of deep-soil CO2 in carbonate areas. The decreased CO2 concentration was roughly evaluated based on stratigraphic unit and farming activities. Soil pH and the purity of CaCO3 in carbonate bedrock deeply affect the corrosion. The corrosion in carbonate areas decreases deep-soil CO2 greatly (accounting for 5.2 %– 66.3 % with average of 36 %) and naturally affects the soil CO2 released into the atmosphere. Knowledge of this process is important for karst carbon cycles and global climate changes and it may be a part of the “missing carbon sink”.


Introduction
In recent years, there has been increasing worldwide concern about carbon exchange among the atmosphere, the ocean, and terrestrial ecosystems.Specifically, there have been ongoing questions regarding the problem of carbon flux, or carbon source versus carbon sink.The missing carbon sink has puzzled scientists since Callendar (1938) presented the imbalance of absorbed and released CO 2 .The missing carbon sink reaches as much as 1.7 Pg, accounting for as much as 24 % of total carbon (Sundquist, 1993).There are differing viewpoints regarding the spatial distribution and absorption strength of the missing carbon sink in terrestrial ecosystems (Fan et al., 1998;Potter and Klooster, 1999).The carbon cycle in karst areas has attracted great interest due to the absorbed and released CO 2 via carbonate corrosion and its share in regulating atmospheric CO 2 (Li and Yuan, 1995;Martin et al., 2013).Therefore, some scholars have looked for the "missing sink" within the absorbed and released carbon in karst systems, and the estimated values reach a dominating part (almost one-third) of the missing sink (Jiang and Yuan, 1999).
Soil carbon, with storage of 1300-2000 Pg C, and as much as 2-3 times that of vegetation storage, plays an important role in maintaining carbon balance (Fearnside, 2018) so that a slight change imposes a great effect on the atmospheric CO 2 concentration.Several factors affecting soil CO 2 concentration, such as environmental factors (soil temperature, moisture, water content, etc.) and human activities, have been widely discussed (Bajracharya et al., 2000;Dai et al., 2004;Owens et al., 2018;Fearnside, 2018).In karst areas, however, the important geological process, carbonate corrosion, has been largely ignored in discussions of soil CO 2 levels and there is no documents detailing the soil CO 2 concentration and its relationship with global climate change in karst areas.Several problems puzzle us.Is there any difference between soil profile CO 2 concentrations in carbonate areas and those in noncarbonate areas?If so, is the difference caused by carbonate corrosion?By how much is it affected?Moreover, studies have revealed that there is a CO 2 imbalance between carbon released into atmosphere and that produced by organic matter in carbonate areas (Jiang and Yuan et al., 1999;Pan et al., 2000) but there is no reasonable explanation.Lack of research work on these questions restricts our understanding of soil CO 2 transfer, limits further study of the mechanisms, and impedes learning of its significance for the carbon cycle.
In order to understand the varying characteristics of soil CO 2 concentrations in karst areas and their potential effect on the global carbon cycle, soil profile CO 2 was measured and samples of soils and rocks were gathered in the typical karst area of Zhaotong City, Yunan Province, China.The objectives of this paper are to (1) analyze comparatively the varying characteristics of the soil profile CO 2 concentration in carbonate and noncarbonate areas; (2) discuss the relationship between soil CO 2 concentration and other parameters, and clarify the effect of carbonate corrosion on soil CO 2 ; and (3) develop a mathematical model of soil CO 2 transfer and quantitatively evaluate the effect scale of carbonate corrosion on soil CO 2 concentration, and discuss its significance for the global carbon cycle and climate change.

Study area
The study area including Zhenxiong County and Weixin County in Zhaotong City, northern Yunnan Province, China, was selected.The area contains high mountains and steep gorges.Many of the mountain peaks tower above 2000 m and there are many different natural watersheds.The area is subtropical and humid.It has a plateau climate with an average annual temperature of 11.7 • C and an average precipitation of 1200 mm.Monthly precipitation is above 100 mm and vertical climate belts with four seasons are clearly demarcated.The soil types include mainly yellow, black and brown earth with a wide thickness range (from a few centimeters up to 70-80 cm).The flora is dominated by grass, shrubs, and partly by secondary forest.
The bedrock is composed predominantly of Mesozoic limestone and dolomite with flysch and associated sedimentary rocks.The widely exposed strata include mainly Ordovician, Permian, Triassic, Jurassic, and Quaternary units.Devonian strata are not present and Precambrian, Cambrian, and Silurian strata occur in limited outcrops or as inclusions among other strata.Ordovician, Permian, and Triassic rocks are mainly marine carbonate deposits and Jurassic and Quaternary units are mainly composed of terrestrial clastic deposits.

Sampling and analyzing methods
In order to comprehensively reveal characteristics of soil CO 2 concentration in the karst area, soil profiles of different stratigraphic units and vegetation types were selected.And profiles in carbonate or noncarbonate areas were both involved.In total, CO 2 concentrations of 21 soil profiles and organic carbon of 12 soil profiles were analyzed.The profile sites are shown in Fig. 1, and among these profiles in carbonate areas include the Lower Ordovician Meitan Formation (O 1 ), the Middle and Upper Ordovician Baota Formation (O 2−3 ), the Lower Permian Xixia and Maokou formations (P 1 m(q)), the Upper Permian Changxing Formation (P 2 c), and the Middle Triassic Guanling Formation (T 2 g).Sites in noncarbonate areas include Middle Permian basalt (P 2 β), shale in the Upper Permian Longtan Formation (P 2 l), mudstone in the Lower Triassic Feixianguan Formation (T 1 f), and siltstone intercalated with shale in the Upper Triassic Xujiahe Formation (T 3 x).
CO 2 concentration within the soil pores was measured every 10 cm from the surface down to the rock-soil interface using a GASTEC 801 instrument and 2LL or 2L CO 2 detector tube (GASTEC Co., Japan).The profile soil samples were of one-to-one correspondence with the gas samples and also taken every 10 cm.
The starting samples were air-dried naturally and then pulverized (particle diameter < 150 µm).Soil organic carbon was determined using the potassium dichromate volumetric method.Soil pH was measured in distilled water at a solid/solution ratio of 1/5, with the instrument model PHS-2.Water contents of soils were synchronously measured by a cutting ring.CaO and MgO contents of rocks were determined by inductively coupled plasma atomic emission spectrometry (ICP-AES) with a charge injection detector (CID), model TJA IRIS/AP.The standard materials (GBW07401, GBW07408) were used for quality control with relative deviation less than 5 %.

Varying CO 2 concentration characteristics of soil profiles
Figure 2 shows soil profile CO 2 concentrations varying with soil depths in seven noncarbonate areas.The data show a distinct tendency of increasing CO 2 concentration with soil depth with R 2 = 0.8-0.92(Table 1).The reasons may be the higher soil bulk density, more condensed soil pores, and difficulty of CO 2 diffusion in the deeper soil.In fact, soil profile CO 2 has been widely reported to be correlated with soil depth by previous research (Rustad et al., 2000;Dai et al., 2004;Malak et al., 2018) and even the following linear equation has been developed (James and George, 1991).Mean CO 2 = 0.035 + 0.0015 (Depth) (R 2 = 0.99, P <0.05).Our observations in noncarbonate areas are concordant with these reports and support soil profile CO 2 increases with soil depth in noncarbonate areas.
Fourteen soil profile CO 2 concentrations with soil depth in carbonate areas were gained (Fig. 3).The results show a complex and inverse relationship between soil CO 2 and soil depth in carbonate areas.Most soil profile CO 2 increases with soil depth in the upper sections, such as the Mangbu O 2−3 grassy profile (Fig. 3a), Mangbu O 2−3 shrub profile (Fig. 3b), Mangbu O 2−3 farmland profile (Fig. 3d), Banqiao O 1 m grassy profile (Fig. 3e), and Banqiao P 2 c grassy profile (Fig. 3j).CO 2 concentrations decrease with soil depth when they increase from surface to a certain depth in Mangbu O 2−3 farmland profile (Fig. 3i), Banqiao P 1 m(q) grassy profile (Fig. 3k), Gaotian T 2 g grassy profile (Fig. 3m), and Mangbu O 1 m grassy profile (Fig 3n).Those of Banqiao O 1 m farmland profile (Fig. 3f) and Banqiao P 1 m(q) shrub profile (Fig. 3l) even decrease all along with soil depth; two farmland profiles of Mangbu O 2−3 (Fig. 3c) and Tangfang P 2 c (Fig. 3h) fluctuate and have no regularity due to the effect of human farming activities.Generally, except Mangbu O 2−3 farmland profile (Fig. 3c) and Tangfang P 2 c farmland profile (Fig. 3h), which are disturbed by farming, CO 2 concentrations of other profiles in carbonate areas all decrease with soil depth at the rock-soil interface (Fig. 3b, e, j).Moreover, there is no correlation of soil CO 2 concentration with soil depth because sequestration of deep-soil CO 2 concentration occurs in carbonate areas.Why does the sequestration only take place in carbonate areas but not in noncarbonate ones?Naturally, the particular carbonate process -carbonate corrosion -is considered, i.e., part of deep-soil CO 2 is consumed and CO 2 sequestration occurs and there is no linear relationship between CO 2 concentration and soil depth in carbonate areas.In fact, Buyannovsky and Wagner (1983), Solomon and Cerling (1987), and Xu and He (1996) all reported that soil CO 2 concentrations reach a peak at a certain depth and then decrease with soil depth in carbonate areas.CO 2 concentrations in Banqiao O 1 m farmland profile (Fig. 3f) and Banqiao P 1 m(q) shrub profile (Fig. 3l) continue to decrease with depth through the integral profile and they also had the highest concentration at the 10 cm layer.Instances of CO 2 concentrations in surface layers higher than those in bottom layers are scarcely documented in carbonate areas.

Relationship between soil profile CO 2 concentrations and soil organic carbon
Soil organic carbon (SOC) was analyzed in a part of the profiles, corresponding with CO 2 concentration.Results are given in Fig. 4, among which panels 4a-h indicate profiles in carbonate areas and panels 4i-l indicate those in noncarbonate (shale) areas.
Correlation analysis of soil profile CO 2 concentration and SOC in shale areas is listed in Table 2. CO 2 shows a decreasing tendency with increasing SOC with high regression coefficients (R 2 = 0.67-0.85).An exception of 0.29 occurs in Wufeng P 2 l secondary forest, which possibly is caused by stronger root respiration and a higher ratio of CO 2 generated by the roots.Therefore, SOC is directly affected by the release of soil CO 2 and the key problem for soil carbon storage is to slow down the renewing of soil organic matter (Chen et al., 2002).The reason for nonsignificance (P >0.05) may be that soil CO 2 concentration is related not only to SOC but also to soil respiration and microbe activities.However, there is no such tendency in carbonate areas as that in shale areas (Table 3) and even those of Banqiao O 1 m farmland profile and Banqiao P 1 m(q) shrub profile show an increasing tendency.Previous studies in carbonate areas such as Shilin, Lunan City, and Guizhou Plateau also showed no correlation between CO 2 concentration and SOC (Liang et al., 2003).
What is the reason for the poor relationship between soil CO 2 and SOC in carbonate areas?A possible answer may be carbonate corrosion.By means of corrosion, deep-soil CO 2 is partly consumed and its level decreases.Consequently, the relationship becomes poor.In addition, varying characteristics of SOC cannot explain well the decrease in deep-soil CO 2 levels in carbonate areas.

Varying characteristics of profile soil pH
Soil pH curves varying with soil depth are drawn in Fig. 5ah, indicating carbonate profiles, and Fig. 6i-l, indicating noncarbonate (shale) profiles.In noncarbonate areas, there is a complex relationship between pH and depth; pH increases obviously at the rock-soil interface, whereas pH nonsignificantly varied with soil CO 2 and SOC.Conversely, in carbonate areas, pH generally increases with soil depth in the surface layer except in the Banqiao O 1 m farmland profile.Moreover, from Figs. 3 and 5 it is evident that soil CO 2 concentrations decrease where soil pH decreases too, and even CO 2 levels in the Banqiao O 1 m farmland profile decrease from the surface to the bottom with soil pH through the entire profile.These observations imply that the decrease in deep- farmland, 7 is Tangfang P 2 c grass, 8 is Tangfang P 2 c farmland, 9 is Wufeng P 2 l shrub, 10 is Wufeng P 2 l second growth, 11 is Wufeng P 2 l grass, 12 is Mohei P 2 l farmland, 13 is Mangbu O 2−3 farmland, 14 is Banqiao P 2 c grass, 15 is Banqiao P 1 m(q) grass, 16 is Banqiao P 1 m(q) shrub, 17 is Tangfang P 2 β grass, 18 is Shuitian T 2 g shrub, 19 is Jiucheng T 3 x shrub, 20 is Mangbu O 1 m grass, and 21 is Zhenxiong T 1 f grass).
By means of this reaction, deep-soil CO 2 is consumed by the corrosion of the underlying carbonate rock and pH decreases synchronously.This reaction cannot take place in soil over areas with noncarbonate bedrock so here the deep-soil CO 2 concentration does not decrease, but increases.

Carbonate corrosion and the global carbon cycle
Many studies have observed that soil CO 2 concentration in carbonate areas decreases with depth when it reaches a maximum at a certain soil depth in carbonate areas (Buyannovsky and Wagner, 1983;Li and Yuan, 1995;Xu and He, 1996;Liang et al., 2003).There has, however, been no reasonable explanation for the observations.= 0.67 to 0.85, although significance is not clear because soil CO 2 is determined not only by organic matter but also by other factors such as root respiration and microbe activities.By contrast, such a correlation in carbonate areas is poor, which was concluded also by Li and Yuan (1995) and Liang et al. (2003) from experiments in carbonate areas.Soil CO 2 of carbonate areas, in every depth at different sites, is negatively correlated with SOC and the relationship becomes worse with increasing soil depth.This observation means that SOC content cannot explain well the decreased CO 2 concentration of deep soil in carbonate areas but rather may be related to carbonate corrosion.Soil pH in carbonate areas always decreases with soil CO 2 and this may imply that H + generated by carbonate corrosion mixes into the deep soil increasing soil acidity.
Previous work has determined the imbalance between soil CO 2 produced and released in carbonate areas.Pan et al. (2000) observed and simulated field data in Yaji, Guangxi Province, concluding that CO 2 produced by decomposition of organic matter is more than that released into the air.This confirms that the rock and the soil have an obviously "absorbing effect" for CO 2 .The data account for an absorbing coefficient of 22-130 g m −2 a −1 .
Isotopes can effectively trace the carbon source of soil CO 2 .Figure 6 reflects the δ 13 C value of soil CO 2 and SOC overlying different bedrock types according to data from Li et al. (2001).It shows that in deep soil, CO 2 has a higher δ 13 C value than the SOC in limestone and dolomite areas, whereas the isotope ratios are more equivalent in clay-stone areas.Such an observation may support the conclusion that deep-soil CO 2 in clay-stone areas is mainly or completely from soil organic matter and that in limestone and dolomite www.earth-syst-dynam.net/10/525/2019/Earth Syst.Dynam., 10, 525-538, 2019 areas there must be an additional carbon source whose δ 13 C should be more than −14 ‰ .CaCO 3 in carbonate has δ 13 C values of −3 ‰ to ∼ +1 ‰.It must, therefore, be recognized that carbon in CaCO 3 of carbonate bedrock mixes into soil CO 2 since the corrosion reaction is reversible.
It has been examined that the karst carbon cycle is an important trace for the global carbon cycle and that further study is important to the hunt for the missing carbon sink (Jiang and Yuan, 1999).From what is presented above, with a focus on the process of carbonate corrosion and comparison of different parameters in carbonate and noncarbonate areas, it is logical to conclude that carbonate corrosion causes the decreased CO 2 concentration at the rock-soil interface in carbonate areas.As a result, the decreased CO 2 level caused by corrosion will, of course, impose effects on atmospheric CO 2 and the karst carbon cycle.This is significant for the potential fixation of carbon, the study of global carbon cycle balance, and the hunt for the missing carbon sink.

Mathematical model of soil profile CO 2 transfer
In this model, only the molecular diffusion of CO 2 is considered, neglecting other processes such as viscous flow and Knudsen diffusion in karst soil because of the weak air pressure gradient.Moreover, density gradient was regarded as the dominant dynamic of CO 2 diffusion and temperature gradient was neglected because of its low contribution (0.2 %-0.4 %) to CO 2 flow.Therefore, the transport of soil CO 2 can be described by the following one-dimensional diffusion equation according to Fick's second law and laws of conservation of mass (Zeng and Zheng, 2002), assuming horizontal homogeneity:  Here, θ a is the air content, θ w is the water content, C a is the gaseous CO 2 concentration, J da is the gaseous CO 2 flow due to diffusion, J dw is the dissolution CO 2 flow due to diffusion, J ca is the gaseous CO 2 flow due to convection, J cw is the dissolution CO 2 flow due to convection, S is the carbon source, Q is the water absorbed by roots, t is time, and z is the space coordinate.
Such an equation can be gained according to Fick's first law: where D a is the gaseous CO 2 diffusion coefficient in soil substrate, D w is the dissolution CO 2 diffusion coefficient in soil substrate, q a is the soil air transference amount, and q w is the soil water transference amount.Equation ( 3) can be deduced from Eqs. ( 1) and (2) if it is assumed that soil water is stable and gaseous and dissolution CO 2 flows are not considered: Previous studies were referenced when the parameters were determined and all the parameters should be gained in win-ter of the same working period: q w = τ exp(− z δ ), presented by Yoyam and Dani (1993), and Q = τ exp(− z δ ) δ by Warren and Michael (1984), and in winter q w = 0 and Q = 0; D a = D 0 a ( θ a θ w )(θ a )( T T 0 ) 1.823 , by Collin and Rasmuson (1988).Here, D 0 a is the CO 2 diffusion coefficient in air at the reference temperature T 0 .
For the carbon source, the rate of CO 2 produced by root respiration and microbes can be expressed as follows: S(z) = S 0 exp(−z/z s ), where S(z) is the soil profile CO 2 at depth of z, S 0 is the CO 2 concentration in the surface soil, z is the soil depth, and z s is the depth gradient.It also considered the CO 2 produced by organic matter expressed as follows: 3.3∂z 2 .Then Eq. ( 4) is achieved: where ∂C a ∂t and ∂C w ∂t are stable when they are from the same time and soil profile.
Based on the studies above, the soil profile CO 2 concentration varying with soil depth can be expressed by the following equation: According to Tailor's formula and it can be roughly expressed like the following equation when x < 1: When Eq. ( 7) is applied to Eq. ( 5), Eq. ( 8) can be gained to express profile CO 2 concentration (C a ) varying with soil depth (z): Here, a, b, and c are uncertain parameters, which vary with θ a , θ w , S 0 , T , and D a of different profiles.
That means it can be expressed as a linear or parabolic relationship of soil profile CO 2 concentration and soil depth.Actually, many observations and simulations also confirmed the same results (James and George, 1991;Zeng and Zheng, 2002;Malak et al., 2018).Therefore, it seems reasonable to express a linear or parabolic relationship of soil profile CO 2 concentration and soil depth.

3.6
The rough evaluation of CO 2 decreased by corrosion SPSS software was used to simulate the curve of measured soil CO 2 concentration and soil depth in noncarbonate areas (Fig. 7 and Table 4), resulting in parabolas with multiple regression coefficients R 2 = 0.8 to 1. Multiple regression coefficient of P 2 c secondary forest profile shows the lowest level at 0.79, which may be due to the different root respiration and the absorbed water at different depths.The simulation evidence that the model is reliable and can be used to roughly reveal the laws of soil profile CO 2 concentration.
In carbonate areas, however, there is no linear or parabolic relationship between soil profile CO 2 concentration and soil depth and the measured values are inconsistent with the simulated ones.Linear or parabolic relationships can be found in the surface soil.Since it is carbonate corrosion that decreases the CO 2 concentration in the deep soil of carbonate areas, the CO 2 concentration in the surface layer can be used to predict the CO 2 concentration of deep soil based on the developed model.The predicting equation and results are listed in Fig. 8.It shows that there is a strong difference between the measured and the predicted values and that all the pre-  dicted values are greater than the measured ones in deep soil.It can also be deduced that deep-soil CO 2 is consumed by carbonate corrosion.
The method of subtraction of predicted and measured values can be used to evaluate the decreased CO 2 concentration in carbonate areas caused by carbonate corrosion, and the results are listed in Table 5.If synthesis factors, such as vegetation types and soil types, were considered, the rough evaluation of the decreased CO 2 concentration of every stratigraphic unit can be gained by taking the average (Fig. 9).

3.7
The main controlling factors of decreased CO 2 concentration Figure 9 shows great dissimilarity in the decreased CO 2 concentration with different stratigraphic units in the following order: T 2 g > P 1 m(q) > O 1 m > O 2−3 > P 2 c. Figure 10 shows the calculated results of the decreased CO 2 concentration in farmland and natural soil (grass and shrub, respectively) of the same stratigraphic unit.CO 2 concentration in T 2 g and P 1 m(q) farmland is lacking but the comparative analysis of O 1 m , O 2−3 , and P 2 c can demonstrate that the decrease in CO 2 in natural soil profiles is obviously less than that in farmland profiles.It is clear that corrosion was strengthened by farming activities and more CO 2 was consumed in the deep soil, which may be due to higher CO 2 levels and acidity caused by farming.Therefore, the decreased CO 2 concentrations of T 2 g and P 1 m(q) should be more than the calculated values, when farming activities are considered.The decreased CO 2 concentration in different farmland profiles is remarkably distinguishable at different sites, even on their same stratigraphic units (Table 5).It seems that the degree of human activity and the quantities of imported or exported energy determine the corrosion to some degree.
Several parameters, such as CaO and MgO contents of carbonate, water content, and pH of the overlying soil, were determined to address some natural factors affecting decreased CO 2 concentration.The parameters are shown in Fig. 11.Deep-soil pH is negatively correlated with decreased CO 2 concentration and the stronger the soil acidity, the more the decreased CO 2 concentration.Water content of deep soil does not impose an effect on corrosion.CaO content of carbonate is positively correlated with the decreased CO 2 con-  centration and the purer the CaCO 3 in carbonate rock, the stronger is the corrosion.MgO content of carbonate is not correlated with corrosion, which indicates that it is CaCO 3 corrosion and not that of MgCO 3 that is consuming soil CO 2 .Simulation by SPSS software results in an equation (y = −3E −08x 2 +0.0002x +6.976) of decreased CO 2 concentration and soil pH with a multiple regression coefficient R 2 = 0.9779 and a second equation (y = 0.0012x + 17.857) of decreased CO 2 level and CaO content of carbonate with a multiple regression coefficient R 2 = 0.4191 (Fig. 12).A field experiment of carbonate corrosion in the southern part of Guizhou (Nie et al., 1984), a laboratory simulation using citric acid to corrode limestone (Cao et al., 2001), and an experimental study on the stability of CaCO 3 and MgCO 3 under acid rain conditions (Teir et al., 2006) led to the conclusion that corrosion is related closely with soil acidity and carbonate purity.The calculated results can support the same conclusion and accord well with their studies.

Discussion and conclusions
It is not surprising that soil CO 2 concentrations decrease in the deep layers over carbonate bedrock areas, especially at the bottom of soil profiles, as has been observed by many experiments (Buyannovsky and Wagner, 1983;Li and Yuan, 1995;Xu and He, 1996;Liang et al., 2003), and is now supported by this paper.The explanation by some studies (Li and Yuan, 1995) that decreased CO 2 is caused by decreased microbe or root respiration in deep soil is challenged by our data.At first, one important reason leading to the earlier conclusion lies, perhaps, in the lack of comparative analyses of soil CO 2 levels in carbonate and noncarbonate areas.The underlying foundation of soluble carbonate in carbonate areas was not taken into consideration, and, most importantly, there was no proof or data to support this idea.Secondly, there is no decrease in CO 2 in soil profiles of noncarbonate areas (mudstone, basalt, shale, or siltstone areas); also, it seems reasonable to expect CO 2 decrease by lower microbe and root respiration rates in deep soil layers of both carbonate or noncarbonate areas.Thirdly, decrease in soil CO 2 takes place in 20-30 cm soil layers, and even from the soil surface in some profiles, so it may be unreasonable to attribute CO 2 decrease to microbe respiration in such shallow occurrences.
Additionally, soil profile CO 2 only decreases in carbonate areas and SOC content is positively correlated with soil CO 2 concentration in noncarbonate areas (R 2 = 0.67-0.85),although there is no significant correlation in some profiles because soil CO 2 is not only related with organic carbon but also with other factors such as root respiration.Soil CO 2 and organic carbon in different depths of carbonate areas are positively correlated with low correlation coefficients but not in soil profiles of these carbonate areas.This means that organic carbon cannot be responsible for the decreased CO 2 concentrations.Furthermore, CO 2 consumed by carbonate corrosion leads to an uncorrelated relationship between soil CO 2 and organic carbon levels in carbonate areas.Soil profile pH in carbonate areas always suddenly and sharply decreases at the depth of CO 2 decrease and this can be explained well by carbonate corrosion.Analysis of δ 13 C isotope, which mixes into a CO 2 in deep soil layers of carbonate bedrock areas (dolomite or limestone) also demonstrates that there is another carbon source, whose δ 13 C level is more than −14 ‰.In soil of clay-stone areas, however, soil CO 2 and soil organic carbon have the same δ 13 C value.This provides strong evidence that carbonate corrosion occurs, and thus deep-soil CO 2 is consumed in carbonate areas.Simply stated, our work strongly indicates that carbonate corrosion leads to the decrease in soil profile CO 2 concentration in areas with carbonate bedrock.
Further, a mathematical model of soil CO 2 transfer was developed, showing that soil CO 2 concentration can be roughly expressed as a linear or parabolic increase with soil depth.The linear or parabolic increase can be demonstrated, strongly supported by both field data and the models.Soil  firmation.In the deep soil of carbonate areas, however, especially at the rock-soil interface, the simulated values are always higher than the field measurements.All of these points may also indicate that carbonate corrosion occurs in the deep soil and that a part of soil CO 2 is consumed by carbonate corrosion.In addition, the decreased CO 2 concentration caused by carbonate corrosion can be evaluated by the subtraction of measured and simulated CO 2 .The decreased CO 2 concentration is related closely to deep-soil pH and CaO content of carbonate rock (correlation coefficients, respectively, R 2 = 0.97 and 0.41), together with farming activities, but not with deep-soil water content and MgO content of carbonate.These results and conclusions can be supported by experiments and are widely accepted by karst scholars, who add validity to our results and conclusions.
The carbon cycle in karst areas has attracted great attention because of the imbalance of the global carbon cycle and in recent years there has been a search to resolve the missing carbon sink related to the absorbing and releasing of carbon in CaCO 3 systems (Jiang and Yuan, 1999).Experiments and calculations indicate that 1.774×10 7 t of carbon are absorbed by karstification in China and that 2.2×10 8 ∼ 6.08×10 8 t of carbon are drawn back from the atmosphere worldwide every year (Jiang and Yuan, 1999).It is obviously significant in regards to the increasing atmospheric temperature.Soil, as an important carbon storage area, is of great importance to atmospheric CO 2 concentration and slight variations may impose great effects on the global carbon cycle.Several factors affecting soil CO 2 concentration have been discussed, such as environmental ones (soil temperature, moisture, water content, etc.), microbe activities, and human activities, but no published details about the effect of carbonate corrosion on soil CO 2 concentration can be found.Our study argues that deep-soil CO 2 concentrations in carbonate areas are obviously decreased, especially at the rock-soil interface, and that this is mainly caused by carbonate corrosion.If this conclusion is correct, then naturally the atmospheric CO 2 levels in carbonate areas should be affected by the corrosion and this should be very significant in the hunt for the missing carbon sink.

Figure 3 .
Figure 3. Varying characteristics of soil profile CO 2 concentration in carbonate areas (profile no. in brackets).

Figure 4 .
Figure 4. Varying characteristics of soil profile organic carbon in karst areas (profile no. in brackets).

Figure 5 .
Figure 5. Soil pH of different profiles in karst area (a-h indicate those in carbonate areas and h-l indicate those in shale areas; profile no. in brackets).

Figure 6 .
Figure 6.Varying δ 13 C of soil CO 2 and soil organic carbon with soil depth overlying different bedrock types (data are after Li et al., 2001).

Figure 7 .
Figure 7.The measured and the simulated CO 2 concentrations of soil profiles in noncarbonate areas.

Figure 8 .
Figure 8.The measured and predicted soil profile CO 2 concentrations in carbonate areas.

Figure 9 .
Figure 9.The evaluation of the decreased CO 2 concentration caused by carbonate corrosion based on stratigraphic units.

Figure 10 .
Figure 10.The decreased CO 2 concentration in farmland and natural soil of the same stratigraphic unit.

Figure 11 .
Figure 11.Relationship of the decreased CO 2 concentration and deep-soil pH, water content, CaO, and MgO contents of carbonate.

Figure 12 .
Figure 12.Correlation analysis of soil pH, CaO of carbonate, and decreased CO 2 concentration.

Table 1 .
Regression analysis of soil CO 2 concentration and profile depth in noncarbonate areas.
Note: numbers in parenthesis give the profile no., * means significant correlation at 0.05 level.

Table 2 .
Correlation analysis of soil CO 2 and soil organic carbon in shale areas of karst.No decrease in soil CO 2 in noncarbonate areas is found, and, furthermore, the depths with decreasing CO 2 concentrations were distinguishable in different profiles, even at only 20-30 cm depths.The decreased CO 2 concentration could be attributed to decreased microbe numbers or root respiration at such depths.By comparative analysis of soil CO 2 concentration in areas of carbonate and noncarbonate bedrock, it should be suggested that the explanation is due to the special geological process of carbonate corrosion.Soil CO 2 and SOC in noncarbonate areas have a good negative correlation with correlation coefficients R 2 P Wufeng P 2 l shrub (9) y = −618.67x+4199.6 0.67 0.387 Wufeng P 2 l second growth (10) y = −766.39x+7548.9 0.29 0.239 Wufeng P 2 l grass (11) y = −13093x + 69 890 0.74 0.351 Mohei P 2 l farmland (12) y = −8646.2x+49 490 0.85 0.077in the deep soil, but there are no scientifically observed data to support this idea and it remains only a hypothesis.

Table 3 .
Correlation analysis of soil CO 2 and soil organic carbon in carbonate areas of karst.

Table 5 .
The evaluated results of the decreased CO 2 concentration in carbonate areas caused by carbonate corrosion.