Effects of Shelter Forests on Soil Organic Carbon of Irrigated Soils in the Taklimakan Desert

: An experiment was carried out to test the effects of artiﬁcial vegetation on soil organic carbon in sandy soil. The soils were collected from the Taklimakan desert highway shelter forests with different stand age (the stand ages are 5, 7, 10, 13, 16, respectively, and three shrubs named Calligonum mongolicunl , Tamarix chinensis and Haloxylon ammodendron were planted along the highway) in Xinjiang, northwest of China. The soil organic carbon stock in soil vertical layers were calculated. We measured four passive organic carbons (i.e., micro-aggregates organic carbon, humic organic carbon, acid-resistant organic carbon and antioxidant organic carbon). Furthermore, we analyzed the correlations and ratios among the different passive organic carbons. Finally, the chemical composition of humus was detected and the relative contents of C=O and CH groups were determined. The main results showed that, (1) the soil organic carbon and organic carbon stock were decreased with the increase of depth, mainly in 0–50 cm. (2) With the increase of stand age, only in Tamarix chinensis forest, the total soil organic carbon stock increased a little. (3) Total soil organic carbon had more closely correlation with contents of micro-aggregate organic carbon and humic organic carbon. (4) C=O/C-O-C increased a little after 10 years; CH/C-O-C had no obvious change with stand age; CH 2 /CH 3 did not change obviously after 13 years. The Tamarix chinensis forest is the most helpful for carbon sequestration in sandy soil and stabilization in surface layer than Calligonum mongolicunl and Haloxylon ammodendron .


Introduction
Buildup of an artificial shelter forest can bring extensive changes to the local soil properties [1]. SOM is closely associated with soil development [2,3]. SOM accumulation in surface soil indicates early soil development [4]. In a desert environment, after an artificial shelter forest was planted, the litter began to enter into surface soil and can provide nutrition for soil microorganisms, promoting the decomposition of litter [5,6]. SOC is a major component of SOM [7]. Therefore, research on changes of SOC, along with the increase of plantation time of artificial vegetation, can preliminarily clarify the OC dynamics of sandy soil under the influence of artificial vegetation.
SOC is contained in plant, animal and microbial residues in all stages of decomposition; SOC can be divided into active, intermediate and passive pools based on the turnover rate of different compounds [8][9][10][11]. Turnover rate of carbon is most often determined by two methods: (a) decomposition studies, (b) natural labeling of SOM using stable 13 C during the development of sandy soil. Thus, we selected several shelter forests along the Taklimakan desert highway with different stand age. By measuring soil passive carbon fractions, we would reveal the functional characteristics of soil humus and calculate soil carbon stock. At last, we would provide information for estimating the carbon balance for arid areas.

Area Description
The study area is located in the Tarim Oil Field of the Taklimakan Desert in Xinjiang, northwest of China. The annual average air temperature has been 12.4 • C for the last 20 years. The temperature on the surface of sand can reach 75.3 • C in summer. The average annual rainfall is 24.6 mm and the average annual potential evaporation is 3638.6 mm. Severe wind and sand hazards occur more than 130 days for one year. The soil chemical and physical properties before vegetation planted were shown in Table 1.

Soil Sampling
Soil samples were collected from the Taklimakan shelter forests with different stand ages in May 2011. The species of tree planted in the shelter forest are three shrubs named C.m., T.c. and H.a. which are resistance to drought, salt and alkali. The plant spacing is 1 m × 1 m. The width of the shelter forest is 72-78 m. The type of the soil is sandy soil (Arenosols) [38]. The shelter forest was drip-irrigated with saline water. The shelter forest was irrigated once every 15 days from March to May and September to October, once every 10 days from June to August and no irrigation from November to February. Overall, 30 mm per time was irrigated.
The five plots were comprised of shrubs with different stand ages (5, 7, 10, 13, and 16 with the plantation time of 2006, 2004, 2001, 1998, and 1995, respectively). The geographic coordinate and altitude are showed in Table 2. At each plot, as for each plant species, we selected three plants randomly which looked like eachother, and can mirrorthe characteristics of the samples in all the plants. The sampling location was 50 cm far from the roots of plant. Furthermore, 1 kg of sub-samples were collected from the soil layers of 0-10, 10-20, 20-35, 35-50, 50-70, 70-100 cm, respectively. The auger was used to soil sample Finally, one sample was mixed from six soil layers. Three sampling locations were collected in one plot. The number of total soil samples was 120. At the same time, soil BD was measured with a drying method in the layers of 0-10, 10-20, 20-35, 35-50, 50-70, 70-100 cm, respectively. The soil samples were air-dried for laboratory analysis. The method for soil aggregate separation was the dry sieving method. The sieving was done manually to achieve aggregate separation by moving the sieve up and down, 50 times in 2 min. First, each soil sample of air-dried was sieved to pass the 2 mm. Then, each soil was divided into two parts. One was for determining the soil TOC, (no fractionates with particle size bigger than 2 mm), HCl hydrolysis, and NaClO oxidation; the other was for next aggregation separation. Next, the soil was sieved to pass 0.25 mm. Soil particle size less than 0.25 mm was micro-aggregate (d < 0.25 mm). OC MIA was determined by the dichromate oxidation method [39].

Soil OC HS , OC NaClO and OC HCl
OC HS was determined by a NaOH-Na 4 P 2 O 7 extraction-separation technique according to Zhang et al. [29]. Moreover, 10 g soil (d < 0.25 mm) was extracted by NaOH (0.1 M)-Na 4 P 2 O 7 (0.1 M) alkali mixture solution (10 g soil/200 mL extraction) for 20 h, and filtered. The filtrate was then divided into two parts. One part was freeze-dried for 72 h and then mixed with KBr (1:120). The characteristics of the chemical structure of humus was determined by Fourier Transform infrared spectroscopy (Bruker Tensor 27, Leipzig Germany). The other part was dried at a constant temperature of 60 • C for three weeks and then determined by the dichromate oxidation method [39].
The OC NaClO was determined according to Siregar et al. [25], with minor modifications. In short, 20 g soil samples (d < 2 mm) were mixed with 50 mL 6% NaOCl solution (pH = 8). After extraction for 24 h at 25 • C, the soil residue was washed by 50 mL 1M NaCl followed by deionized water, then centrifuged (10,000× g rpm for 5 min). The procedure was repeated twice and the residue was then dried at 60 • C for 72 h for the determination of antioxidant organic carbon. The OC HCl was extracted according to Pual et al. [40], with minor modification, by changing the extraction time from 18 h to 3 h. The procedure was similar to that for extraction of OC NaClO , except the extractant was 6 M HCl.
The hierarchical structure of SOC fractions is shown in Figure 1.
We used the following formula (Equation (1)) for the calculation of OC stock in the i layer: OCi, Stock = OCi × BDi × Thicknessi × 10 (1) OCi was the concentration of organic carbon (g kg −1 ) passed through 0.25 mm sieves in the bulk soil of the i layer, BD was the BD of the i layer (g cm −3 ), and thickness was the soil layer thickness (cm) of the i layer. OCStock (g m −2 ) is the sum of all the SOC in each layer (Equation (2)).

Statistical Analysis
Effects of stand age on soil organic stock and passive fractions were verified statistically using one-way ANOVA analysis, by the statistical softness of SPSS 16.0. Means were compared using new multiple range test (Duncan), and the relationships between different organic carbon fractions were analyzed with Pearson correlation coefficient. Difference was declared significantly at α < 0.05. All the figures were conducted by the softness of Origin 8.0.

SOC and OCstock at Different Stand Ages and Soil Depth
The changes of SOC with soil depth and stand age can be seen from Table 3
We used the following formula (Equation (1)) for the calculation of OC stock in the i layer: OC i was the concentration of organic carbon (g kg −1 ) passed through 0.25 mm sieves in the bulk soil of the i layer, BD was the BD of the i layer (g cm −3 ), and thickness was the soil layer thickness (cm) of the i layer. OC Stock (g m −2 ) is the sum of all the SOC in each layer (Equation (2)).

Statistical Analysis
Effects of stand age on soil organic stock and passive fractions were verified statistically using one-way ANOVA analysis, by the statistical softness of SPSS 16.0. Means were compared using new multiple range test (Duncan), and the relationships between different organic carbon fractions were analyzed with Pearson correlation coefficient. Difference was declared significantly at α < 0.05. All the figures were conducted by the softness of Origin 8.0.

SOC and OC stock at Different Stand Ages and Soil Depth
The changes of SOC with soil depth and stand age can be seen from Table 3. Mainly in 0-100 cm, the contents of H.a. (in 5 years), T.c. (in 10 years) and C.m. (in 10 and 16 years) decreased with soil depths, and did not change obviously in those layers. The content of SOC in surface layer was the highest. As for three kinds of plants, C.m have more influence on the content of the soil carbon in a vertical space than those of H.a. and T.c. Table 3. The changes of SOC (g kg −1 ) with soil depth in the shelter forest under different stand ages. Data were mean ± standard error (n = 3). Different lowercase letters indicate significant differences among soil depth under the same shelter forest and the same year at p < 0.01 level, and different capital letters indicate significant differences among stand age under the same shelter forest and the same soil depth at p < 0.05 level. The value of OC stock in 0-100 cm was calculated in Table 4. The results showed that the content of OC stock was higher under T.c. than those under H.a. (except 10 years), and C.m. (except for 10 years and 13 years). This indicates that in early times during the development of the sandy soil, T.c. has more effects on accumulation of OC stock in BD than those of H.a. and C.m. Table 4. Changes of soil OC stock (g m −2 ) in different stand ages in the soil layer of 0-100 cm. Data were mean ± standard error. Different lowercase letters indicate significant differences among plant species at p < 0.01 level, and different capital letters indicate significant differences among stand age at p < 0.05 level.

Dynamic Changes of Different Carbon Fractions under Different Stand Age
Dynamic changes of different carbon fractions under different stand age can be seen from Figure 2. Only the total SOC stock in T.c. forest increased with the stand age. This indicates that the stand ages in T.c. forest have great influence on the sequestration of soil carbon.
With the increase of stand age, all the contents of TOC, OC MIA , OC HS , OC HCl and OC NaClO increased ( Figure 2). The contents of TOC, OC MIA , OC HS increased year by year. The contents of OC HCl and OC NaClO increased mainly in the first few years (except OC NaClO in T.c. forest and OC HCl in C.m. forest). The contents of OC HCl and OC NaClO in H.a. forest increased faster in 5-7 years. the content of OC NaClO in C.m. forest increased faster in 5-10 years. the content of OC HCl in T.c. forest increased faster in 5-10 years. The order of five carbon content for three kinds of forest was: TOC > OC MIA (except 5 years in T.c. forest), OC NaClO > OC HS > OC HCl (Table 5). This indicates that the T.c. is more helpful for passive carbon accumulation in surface soil than those of H.a. and C.m.  (Table 5). This indicates that the T.c. is more helpful for passive carbon accumulation in surface soil than those of H.a. and C.m.

The Relationships among Carbon Fractions
The correlation coefficients between OCMIA, OCHS and TOC were higher than those of OCHCl, OCNaClO and TOC. This indicates that the humic and aggregate organic carbon changed in line with TOC during the early development stage of sandy soil, and had a closer relationship with TOC.
Correlation coefficients between the various passive carbon fractions showed that OCHS and OCMIA were the highest (0.954, 0.981, 0.954), followed by OCNaClO and OCMIA (0.909, 0.834, 0.656) ( Table 6). This indicates that the relationship between OCHS and OCMIA was closer than those between other passive carbon fractions.

The Relationships among Carbon Fractions
The correlation coefficients between OC MIA , OC HS and TOC were higher than those of OC HCl , OC NaClO and TOC. This indicates that the humic and aggregate organic carbon changed in line with TOC during the early development stage of sandy soil, and had a closer relationship with TOC.
Correlation coefficients between the various passive carbon fractions showed that OC HS and OC MIA were the highest (0.954, 0.981, 0.954), followed by OC NaClO and OC MIA (0.909, 0.834, 0.656) ( Table 6). This indicates that the relationship between OC HS and OC MIA was closer than those between other passive carbon fractions.
The values of OC MIA /TOC, OC HS /TOC, OC HCl /TOC and OC NaClO /TOC reflected the relative percentage of passive carbon accounting for soil TOC (Table 6). OC HCl /TOC and OC NaClO /TOC in the different forests almost decreased with the increase of stand age. In H.a. forest, OC HS /TOC was enhanced with the increase of stand age. In T.c. and C.m. forests, the value of OC HS /TOC was less than the values of OC HCl /TOC and OC NaClO /TOC. This indicates that, with the increase of stand age, the role of humification may become a more important factor in carbon sequestration in a chemical way. The percentage size of four fractions in TOC was OC MIA /TOC > OC NaClO /TOC > OC HS /TOC > OC HCl /TOC. The relative proportion of OC MIA and OC NaClO in TOC is higher than those of OC HCl and OC HS (Table 7). This indicates that aggregation and resistance to oxidation may play more important role than humification and resistance to hydrolysis in soil carbon sequestration during early times in the soil forming process; aggregation played the most important role in the sequestration of soil carbon. The role of resistance to hydrolysis was the weakest. The relative content of passive organic carbon fractions of T.c. forest was the highest. This indicates that the role of T.c. in carbon stabilization was higher than those of C.m. and H.a.

The Chemical Structure of the Humic Molecular
The baseline correction of original infrared spectra on shapes of humus was conducted first. The infrared spectra were obtained by taking the absorbance as peak height, and wave number as peak position (Figure 3) [29,33,35,41].
est. This indicates that the role of T.c. in carbon stabilization was higher than those of C.m. and H.a.

The Chemical Structure of the Humic Molecular
The baseline correction of original infrared spectra on shapes of humus was conducted first. The infrared spectra were obtained by taking the absorbance as peak height, and wave number as peak position (Figure 3) [29,33,35,41]. It showed that the C=O group appeared at 1670 cm −1 , hydroxyl hydrogen bonding associated a peak at 3400 cm −1 , and the vibrations of O-H surface deformation and O-H stretching appeared at 920 and 2491 cm −1 , respectively. The peaks of aromatic C=C appeared at 1637 cm −1 and the bands appeared around 1456 cm −1 due to aliphatic C−H. Two bands appeared respectively at 2897 and 2839 cm −1 due to the vibration of aliphatic C−H. No absorption peaks of CH3 and CH2 groups appeared in 5 years.
The FTIR spectra were half quantitatively analyzed at two absorption bands. In this study, we take the original peak height of C-O-C group (at 1134 cm −1 ) as a, and the original peak height of C=O group (at 1670 cm −1 ) as b, and the value of b/a represents the amount of C=O groups relative to the C-O-C groups, providing information about the relative content of the C=O group in SOM fractions [32,33]. The calculation results showed that the value of b/a increased within 10-16 years ( Table 8). The values of c/a and d/a represented CH3/C-O-C and CH2/C-O-C and provide information about the relative content of CH group (CH3 and CH2, respectively) in the organic matter molecule [32,33]. Furthermore, the value of d/c (CH2/CH3) indicates the ratio of branched chain to carbon chain [35]. With the increase of stand age, the values of c/a and d/a did not change noticeably. The value of d/c had no obvious change with stand age. The ratios of C=O, CH2 and CH3 to C-O-C were: C=O > CH2 > CH3. It showed that the C=O group appeared at 1670 cm −1 , hydroxyl hydrogen bonding associated a peak at 3400 cm −1 , and the vibrations of O-H surface deformation and O-H stretching appeared at 920 and 2491 cm −1 , respectively. The peaks of aromatic C=C appeared at 1637 cm −1 and the bands appeared around 1456 cm −1 due to aliphatic C−H. Two bands appeared respectively at 2897 and 2839 cm −1 due to the vibration of aliphatic C−H. No absorption peaks of CH 3 and CH 2 groups appeared in 5 years.
The FTIR spectra were half quantitatively analyzed at two absorption bands. In this study, we take the original peak height of C-O-C group (at 1134 cm −1 ) as a, and the original peak height of C=O group (at 1670 cm −1 ) as b, and the value of b/a represents the amount of C=O groups relative to the C-O-C groups, providing information about the relative content of the C=O group in SOM fractions [32,33]. The calculation results showed that the value of b/a increased within 10-16 years ( Table 8). The values of c/a and d/a represented CH 3 /C-O-C and CH 2 /C-O-C and provide information about the relative content of CH group (CH 3 and CH 2, respectively) in the organic matter molecule [32,33]. Furthermore, the value of d/c (CH 2 /CH 3 ) indicates the ratio of branched chain to carbon chain [35]. With the increase of stand age, the values of c/a and d/a did not change noticeably. The value of d/c had no obvious change with stand age. The ratios of C=O, CH 2 and CH 3 to C-O-C were: C=O > CH 2 > CH 3 .

The Changes of SOC and Soil Carbon Stock
The accumulation rate of SOC is associated with the amount of decomposable exogenous organic matter. Litter and dead plant roots are important exogenous organic matter. In our study, the increases of organic carbon and organic carbon stock with depth mainly appeared in the 0-50 cm soil layer. This is because litter is above this layer, and plant root is below this layer. The amount of litter may be greater than the amount of dead root. Litter decomposition rate is affected by their own characteristics, including contents of N, C/N and contents of lignin and cellulose in litter [42]. The litter characteristics of T.c. may be more conducive to the processes of litter decomposition and nutrient release, so the soil carbon content in T.c. forest is higher.

The Effects of Artificial Vegetation on Chemical Passive Organic Carbon
The content of TOC was 0.73-3.55 g kg −1 in the surface layer of sandy soil, which was the smallest compared with other soils developed from different parents, such as podzol 35.33 g kg −1 [16], Mollisols 5.5-8.5 g kg −1 [43], loess 5.5-8 g kg −1 [44], and fen 200.4 g kg −1 [45]. This indicates that the pedogenic process of arenosols may be the slowest.
In our study, the contents of TOC were 2.53, 2.80 and 3.61 g kg −1 and those of the OC MIA were 1.79, 2.47, 3.40 g kg −1 in 13 years, for different artificial vegetation ( Figure 2). This result is smaller than the findings of Tan et al. [37] who reported that the contents of TOC and OC MIA were 27.6 and 42.79 g kg −1 , respectively in 12 years. This could be because the dominant species and litter quantity in two sites are different. Thus, the quantity and chemical composition of plant tissues entering the soil should be different, which is related to the content of SOM [46]. The main species in this study is Calligonum mongolicunl, Tamarix chinensis and Haloxylon ammodendron, and the one in Tan's is Pinus sylvestris var. mongolica. Different litter quality and quantity may lead to different contents of TOC and OC MIA in two sites.
In our study, the content of TOC was higher than that of OC MIA, which is consistent with Yu et al. [47]. However, it is not consistent with Tan et al. [37]. This may be due to one of the stabilization mechanisms of soil OC: physicochemical leads to occlusion of organic matter into micro-and macro-aggregates by binding SOM and mineral particles [48]. This process may mainly occur in macro-aggregates but not micro-aggregates during the early period (about 10 years) of sandy soil development, and a similar phenomenon occured in previous studies [49,50]. After that, the process of occlusion of organic matter into micro-aggregates may play a more important role in soil OC stabilization, so the difference between TOC and OC MIA may be shortened and even TOC may smaller than OC MIA . The first stage may continue for a longer time in sandy loam soil and crop ecological environment.
The TOC and OC MIA were increasing with stand age (Figure 2). This is consistent with Jin [36]. However, it is not consistent with Tan [37]. This may because the amount of TOC and OC MIA may stop accumulating in sandy soil after a specific development time scale.
Humus is more resistant to decomposition than the non-humic material [51], and is generally fixed in the process of humification [15]. In our study, the content of OC HS increased with stand age, being consistent with the findings of Zhang [29,52]. The reason for this phenomenon could be attributed to the amount of litter that can be decomposed and changed into humus in the soil surface increasing with the increase of stand age.
The contents of OC HCl and OC NaClO increased slowly after 10 years. The OC HCl and OC NaClO simulate the oxidation and hydrolysis decomposition process of organic matter in nature [22,23]. They are two complex chemical reaction processes. Chemical reaction rate is affected by local climate, such as the air temperature rainfall and the characteristics of the plant species. The air temperature influences soil temperature. The rainfall influences soil humidity. The plant species, such as the content of litter from lignin and cellulose, affect the possibility of the reaction substrates. Therefore, the specific reason for these two kinds of passive fractions contents stop increasing the need for intensive research.

The Relation among Carbon Fractions
The correlations between OC MIA and TOC in our study are 0.949, 0.950, 0.900 (Table 7); it is close to the results of Tan [37] and Yu [47], who reported that the correlation between OC MIA and TOC is 0.979 and 0.964, respectively. This indicates that the correlation between OC MIA and TOC is relatively stable, although the soil development time and vegetation species are not the same.
The average value of OC MIA /TOC (84.8%) is similar to Tan. et al. (80.13%) [37], but lower than Yu et al. (97.29%) [47]. The development time of Tan's study is 39 years and Yu's is more than 100 years. Therefore, a longer development time may influence the value of OC MIA /TOC.
The humification plays a more important role in carbon stabilization by chemical way. The reason for this may be that NaClO attack and HCl hydrolysis are chemical activities carried out by enzymes and biodegradation. However, SOM's major part is bound to the mineral soil component, and mineral particles act as a surface for HS self-assembly [53], so humus exists in the soil with an inorganic-organic compound. Therefore, with the increase of stand age, humification may become more active.
The content of OC HCl , OC NaClO , OC HCl /TOC and OC NaClO /TOC decreased. This indicated that resistance to hydrolysis and oxidation were restrained. The reason for this phenomenon may be related to two SOM stabilization mechanisms: (1) the innate stability (recalcitrance) of the organic matter; and (2) the absence of the required suite of microbial decomposers [18].
The absence of the required suite of microbial decomposers can make the turnover time of OM longer [18]. In initial years, after the shelter forest had been planted, the amount of litter in the surface soil was very low, so the microbial available carbon and nitrogen source was limited, leading to the low microbial population [54]. Therefore, the exogenous organic matter was difficult to mineralize, but the possibility changed into passive organic carbon, such as OC HCl and OC NaClO increasing, so OC HCl /TOC and OC NaClO /TOC was very high. With the increase of stand age, exogenous organic matter accumulated in the surface soil. The situation appeared opposite, causing the decomposition rate of soil carbon pool to intensify. The possibility of active carbon pool transforming to passive carbon would decrease, while the TOC continued increasing, so the OC HCl /TOC and OC NaClO /TOC decreased.

The Influence of Artificial Vegetation on Chemical Composition of Humus
The main functional groups appeared in the humus molecular of sandy soil during humification process in Tarim desert is C=O, CH and C=C ( Figure 3). This is in line with Zhang [29]. The number of aliphatic groups decreased except for C-O-C in 5 years and C=O, CH in 13 years. This is in line with Zhang [29]. The stand age of Zhang is 0, 20, 32, 43, 51, and the major plant species were Caragana korshinskii, Artemisia or dosica, Hedysarium. This indicates that the stand age and species of artificial vegetation would have little effect on the aliphatic changes in the humic chemical structure during the development of sandy soil.
The year when CH group appeared is different to that of Zhang [29]. Our study took place in 7 years, but CH appeared until 51 years in Zhang. The reason for this inconsistency may also be explained as composition of the plant species is different in the two studies and could have impact on the development of humic sustance in soils. Thus, the different properties of litters in the two studies could lead to a difference in CH groups of humic sustance. This indicated that the evolution rate of CH groups may be different under different artificial vegetation. Similar results can be seen in different ecological systems [33,55].
Zhang [29] showed in Tengger Desert that, the changes of the chemical structure of an aromatic compound of a humic substance developed constantly over the years. The reason may be because the time scale in two studies is different. The first age of the former study is 5 years, but the latter is 20 years. The former increased in 13 years, and may continue to increase, so the number changes of aromatic groups in humus molecular may increase with stand age.
The content of C=O (C=O/C-O-C) in soil organic matter increased from 10 years to 16 years (Table 8). This means that the ability of soil to form a mineral-organic compound to stabilize carbon is increased. This may be related to the soil clay particles contents which are affected by stand age. SOM interacts with soil minerals mainly via functional C=O groups [8,31]. With the stand age increased, the quality of soil increased, and the soil clay particles contents increased, which made soil humification intensity strengthen and form more clay particles in soil [56]. However, the process is very slow at the beginning of the development of sandy soil, so the C=O/C-O-C increased a little.
The contents of C=O, CH 2 and CH 3 after 7 years followed the order: C=O > CH 2 > CH 3 . SOM interacted with soil minerals mainly via functional C=O groups [8,34]. This means that the mineral-organic compound may play more important roles in humification process than alkyl carbon in early times of soil development. SOM was stabilized by interaction with mineral components [18,55], so the C=O > CH. Additionally, humus is a kind of high polymer formed by the partial decomposition of polymers from plant residues, followed by condensation and polymerization reactions [17]. CH 2 can be a structure monomer of polymethylene [28], but CH 3 cannot. Therefore, the humic molecular tends to form more CH 2 than CH 3 to strengthen its own stabilization.

Conclusions
The content of SOC reduced with the increase of soil depth, mainly in 0-100 cm. The T.c. forest played the most important role in soil carbon accumulation in unit mass and C.m. play the most in unit area. The space variation caused by C.m. is the biggest.
In terms of the contents, the passive organic carbon increased with stand age, and the OC MIA and OC HS increased years by years, while the OC NaClO and OC HCl increased remarkably in the first few years, then had no obvious change. Aggregation and resistance to oxidation played a more important role than humification and resistance to hydrolysis in carbon sequestration during the early soil development stage. With the increase of stand age, the humification may play a more important role.
The T.c. forest is the most helpful for carbon and sequestration carbon stabilization in sandy soil. The artificial vegetation influenced the molecular stability of the humic substance. Mineral-organic compounds may play more important roles in the humification process than alkyl carbon in the early times of soil development.
Some conclusions need further study to prove their accuracy. First, the vegetation species may influence: (1) the changes of TOC and OC MIA ; (2) the speed of CH group appeared in humic molecular. Second, the macro-aggregates may appear before microaggregates in sandy soil.
The stabilization mechanism of carbon in sandy soil needs more investigation to enclose the chemical structure of humus, and the chemical bonds of CH 2 and CH 3 . The time scale that the studies chose may influence the results in studying aromatic groups.