Soil organic matter (SOM), as a kind of natural polymers, affects the migration and transport of ions and particles in soil system due to its surface characteristics and interaction and then causes significant changes in soil quality such as soil fertility loss and pollutant transfer. It is of great importance to study the temporal and spatial evolution of soil organic matter and its driving mechanism for soil health management. This study aims to fully reveal the evolution characteristics and driving mechanism of soil organic carbon (SOC) in farmland of the Beijing plain based on a six-year site monitoring. According to the research results, there is a significant difference in the overall soil organic content during the 6-year period. The temporal stability of SOC is moderate, and it is inversely proportional to SOC content in terms of spatial distribution. SOC content increases as organic fertilizer input rises, and an extra unit (15 kg·ha−1) of organic fertilizer input leads to an increase of 0.057 g·kg−1 in SOC content. The soil with higher clay content exhibits higher SOC content. The organic carbon content in different soil texture types increases with time, and there is a significant difference between the increases in medium loam and light loam. The grain field plantation system exhibited declining SOC content, while the protected vegetable fields, open vegetable fields, and orchards all showed an increase in SOC content. According to our results, the SOC content of farmland in the plain areas of Beijing is largely dependent on the input of organic carbon if other conditions remain unchanged or exhibit insignificant changes.

1. Introduction

With the development of modern macromolecule chemical industry and the intensification of agricultural production, many synthetic organic compounds have been discharged into soil directly or indirectly through the atmosphere or water body in modern human production and life. Microplastics are widely used in remediation of contaminated soil due to their small size and large specific surface area [1, 2]. However, due to the specific environmental behavior and potential environmental risks of engineering microplastics, many researchers believe that the use of engineering nanomaterials should be prohibited in in situ repair [3]. Therefore, more and more researchers focus on the nature and application of natural soil microplastics. Soil organic matter (SOM), as an important source for natural microplastic in soil, participates in various ecological processes including element cycling and pollutant migration [4], because it has had large specific surface area and is easy to adsorb dissolved organic matter in soil-water. SOM can improve soil adsorption, buffering, and fertilizer retention and makes soil more loose, noncohesive, and easy to cultivate. In addition, it also can store carbon from the atmosphere and other pollutants and alleviates the biological hazards of chemical pesticides, heavy metals, and other pollutants in soil [57]. For example, if humic acid particles, a kind of important ingredient of soil organic matter, are adsorbed on the surface of microplastics, the steric hindrance of microplastics will increase, thus affecting the migration of microplastics [8, 9].

Therefore, it is of great significance and urgent to study the spatiotemporal changes and driving mechanisms of soil organic matter and to reasonably and effectively control the content of soil organic matter for maintaining soil fertility, reducing soil pollution, and promoting the sustainable utilization of soil resources.

In the natural environment, adsorption will affect the migration, transformation, and biodegradation of pollutants [10]. Surfactants are used to treat refractory organic pollutants in soil-water environment to promote special desorption of soil and increase apparent water phase concentration. Several studies have been carried out on the application of surfactants to strengthen remediation of specially polluted soil, including soil cleaning application and bioremediation [11, 12]. SOC content has obvious characteristics of temporal evolution, and many studies have been conducted on characterizing the evolution of farmland SOC content [1316]. Carbon simulation models have been adopted in some analyses of evolution features of regional SOC content. At present, the carbon models that are frequently used worldwide include DNDC [17], RothC [18], CENTURY [19], and APSIM [20]. All these models have been applied to simulate changes in farmland SOC content at both national and regional levels [1316, 21]. The method of mathematical modeling has more advantages in simulating SOC dynamics of farmland in larger spatial and temporal scales. But the simulation results are not representative of the actual situation. Therefore, the study on SOC content change based on data from long-term monitoring at fixed locations is needed.

Research on the driving mechanism of SOC change based on sample points has been carried out. However, most of the factors used to analyze the driving mechanism are static, e.g., elevation and soil type which mainly affects the spatial distribution pattern of SOC content. The dynamic development of farmland SOC is also driven by the natural environment and artificial management. Human activities impose significant effects on SOC. According to Davidson and Ackerman [22], the organic carbon in the 30 cm surface layer of soil could be reduced by as much as 40% on average after natural vegetation is reclaimed into farmland. According to some studies [23, 24], improvements in agricultural management, such as minimum and no tillage, and the increase of chemical and organic fertilizers, can alleviate or reverse the loss of SOC. As a result, it is necessary to explore the impact of human activities on the change of soil organic carbon content.

In short, this study, taking the farmland of Beijing plain as the study area, aims to (1) reveal the evolution characteristics of SOC during 2010 to 2015 and (2) explore the influence of human activity like organic manure and planting system and soil properties such as soil texture on the change of SOC based on the long-term location monitoring data.

2. Data and Methods

2.1. Overview of the Research Zones

According to the digital elevation of Beijing and the spatial distribution of farmland, the 100 m contour was set as the research boundary. The research area is located in the eastern plain of Beijing, with a total area of about 7779.12 km2 (Figure 1). The research zones were divided into three functional areas, including cover the water source conservation area, the agricultural security area, and the urban life area. The area of light loam is the largest, followed by sandy loam and middle loam (Figure 2). And the international standard soil texture classification standard was adopted in this study.

2.2. Data Processing and Analysis

Long-term monitoring points have been established according to the distribution of farmland. There were 149, 137, 129, 123, 104, and 91 sampling points from 2010 to 2015 (Figure 1). To ensure the representativeness of the samples, the cross-section soil samples were all collected from the central points and the two neighboring angles at each sampling point. Furthermore, the soil samples collected at the same depth were mixed prior to testing (Figure 1). Potassium dichromate titration was adopted to assess SOC (Johnson, 2010). Aside from the data collected at the monitoring points, the digital elevation (resolution: 25 m×25 m), the 1:50000 soil texture map, and the input of organic fertilizer were obtained via the Beijing municipal digital soil system and field investigation.

SPSS software (SPSS, Inc., Chicago, IL, USA, 2012) was applied to conduct descriptive statistical analysis and variance analysis. The temporal stability, the spatial prediction, and spatial correlation analysis were conducted in the ArcGIS 10.2 (ESRI, 2010).

2.3. Study Methods
2.3.1. Empirical Bayesian Kriging

The traditional linear kriging methods have strong smoothing effects and require a relatively large number of sampling points. Meanwhile, the data also needs to satisfy the spatial stationarity hypothesis. The empirical Bayesian kriging (EBK) method can measure unstable variables and is more accurate than other kriging methods for small data sets [25, 26]. EBK estimates basic semivariation functions to describe the introduced errors. The input data is used to simulate several semivariation functions to describe the uncertainty of semivariation function estimation, and the predicted standard error should be smaller [2528].

Additionally, we also validated whether the different number of samples in the above six years can meet the research requirements. Referring to the Technical Specification for Soil Environmental Monitoring, the variability coefficient and relative deviation are used to determine whether the sample points meet the research requirements. The formula is as follows.

where n is the number of samples, t is the value of t at a certain degree of freedom at a given confidence level (generally 95% for soil monitoring), CV is the coefficient of variation (%) obtained from the collected data, m is the acceptable relative deviation (%), and soil monitoring is generally limited to 10-20%.

In this study, the coefficient of variation (CV) of SOC in the above six years was 39.31%, 33.87%, 47.02%, 40.62%, and 48.57%, respectively. According to the analysis results, the average coefficient of variation (%) of all monitoring indicators was 42.03% and the acceptable relative deviation m was taken as the minimum value of 10%. Through calculation, we need to set up sampling points 40. Based on the above results, the minimum sampling points of 91 in 2015 can meet the research requirement. As a result, the number of samples in different years can meet the analysis requirements.

2.3.2. Overall Temporal Stability Coefficient of SOC

The SOC temporal stability in the 6-year study period was quantified by the variation coefficients of measurements obtained from repeated sampling. This study adopts the variation coefficient calculation method first proposed by Blackmore [29]. The equation for the temporal stability coefficient (TSCVti) is as follows:where SOCit refers to the measured value of soil property obtained at the “ith” sampling point at the “tth” sampling time, and “n” is the frequency of sampling (n=6).

The temporal stability coefficients were divided into five levels in three types, namely, “Stable” (<10; S), “Medium Stability” (10-25; MS), and “Unstable” (>25; US). “Medium Stability” was further divided into three levels, namely, MS1 (10-15), MS2 (15-20), and MS3 (20-25) [29].

2.3.3. Quantitative Analysis of Relationships between SOC Content and Influencing Factors

Simple spatial correlation analysis (R) can be used to explore the relationship between two specific variables.where Rxy is the spatial correlation coefficient of x and y, xi is the SOC content in the “ith” year, yi is the input of organic fertilizer in the “ith” year, is the mean SOC content of six years, is the mean organic fertilizer input of six years, and n is the total number of years (n=6).

3. Results and Discussion

3.1. Overall Evolution Features of SOC Content

The organic matter in the soil interacts with the engineering polymer to a certain extent, and the degradation of the artificial polymer is induced where the organic matter exists. Organic matter is one of the important reasons for the degradation of artificial macromolecules. The microplastics discharged into the environment will interact with the organic matter in the soil. The addition of humus in the soil will increase the steric hindrance between the microplastics themselves, making the particle size of the microplastics significantly smaller than that in the presence of humic acid. Therefore, when humic acid exists, the penetration ability of the microplastics is enhanced, and its retention in the medium is reduced. Understanding soc content in Beijing Plain is helpful for remediation of soil macromolecular pollution in the later period.

We calculated the statistical eigenvalues of SOC content in different years. The Shapiro-Wilk test (the W method) is used for normal distribution test.

The mean value of SOC content increased with time. According to variance analysis, there was a significant difference (F=3.042, df=727, p⩽0.001) in SOC content during the six years. As the mean value of SOC content increased, the standard deviation and range also increased; i.e., the regional difference became increasingly significant with the human interference. The different interannual SOC contents were consistent with the linear relationship, and the fitting effect was significant (df=4, R2 = 0.8733 > = 0.88114). According to the normal distribution test of the W method, the SOC contents from different years were consistent with the normal distribution (p>0.05). The cumulative variation ratio of farmland SOC content in the plain areas of Beijing was 12.89%, with an average annual variation ratio of 2.58% (Table 1). In the literature conducted by Liang et al. [30] found that the average farmland SOC content of North China was 17.70 g·kg−1, which is higher than the maximum observed in this study. The above results indicate that there is still huge potential for increases in the farmland SOC content.

The EBK method was then applied to obtain the SOC temporal stability variation coefficient figure (Figure 3), which shows the time-based variation trend of SOC content. The subset of the empirical Bayesian method was 100; the number of overlap factors was 3; the number of simulations was 100; there were 4 search sections; the search radius was 29.55 km2; and the output grid was 100 m×100 m (the size of lattice was 1 ha2).

According to Figure 3, the average of SOC temporal stability coefficients in the research zones was 20.25. The SOC temporal stability mainly stayed in MS3, between stability and instability. The areas of MS1, MS2, and MS3 were 42.29 km2, 1027.14 km2, and 718.11 km2, respectively. The unstable areas were mainly located in Haitian District, Chaoyang District, Sheringham District, and Changing District in the urban life area. The temporal stability of SOC content was totally different from the SOC content in terms of spatial distribution, and regions with higher SOC content showed a lower temporal stability.

As far as functional area is concerned, SOC content was the least stable in the urban life area and the peripheral circle with a radius of 35 km. For over 80% of the urban life area, temporal stability stayed in US. In contrast, the agricultural security area and the water source conservation exhibited opposite patterns, and these areas mainly stayed in MS.

3.2. Analysis of Influencing Factors for the Evolution of SOC Content

This study selects factors that reflect farmland management such as the input of organic fertilizer and the plantation system as well as the soil texture of the natural environment to analyze the driving mechanism behind the evolution of farmland SOC content.

3.3. Input of Organic Fertilizer

From 2010 to 2015, the average annual fertilizer input was 38.24 kg·ha-1, 38.16 kg·ha-1, 38.81 kg·ha-1, 40.43 kg·ha-1, 42.25 kg·ha-1, 44.20 kg·ha-1, respectively. The correlation analysis between the input of organic fertilizer (scalar quantity, kg·ha−1) and the mean SOC content of corresponding years is carried out (Figure 4(a)). R at the annual scale was employed to calculate the spatial correlation (Figure 4(b)).

Figure 4(a) shows the quantitative fitting relationship between the mean SOC contents and the input of organic fertilizer (scalar quantity) during the 6 years. The average SOC content and the input of fertilizer were consistent with a linear equation, and the fit is significant (df=4, R2 = 0.8433 < = 0.88114). An increase in the input of organic fertilizer led to an increase in SOC content, and an extra unit (1 kg·ha−1) of organic fertilizer led to an increase of 0.057 g·kg−1 in SOC content.

The average spatial correlation coefficient between SOC content and the input of organic fertilizer over several years was 0.49 (Figure 4(b)), and there was an overall positive spatial correlation between these two variables. The regions with a positive correlation had a total area of 6120.5 km2, accounting for 78.69% of the total study area. Specifically, 19.72% and 31.30% of the total area exhibited significance at the p<0.05 level and the p<0.01 level, respectively. The input of organic fertilizer can effectively increase SOC content. The above conclusions are consistent with the achievements of relevant studies (Dendoncker et al., 2011). According to the studies of Shen et al. [31], Chen et al. [32], Wu et al. [33], Dendoncker (2011), and Liang et al. [34], the average carbon fixation rate of the farmland soil in Belgium nearly doubled with the use of organic fertilizer. Shen et al. [31] also argued that the long-term use of organic fertilizer or a mixture of organic fertilizer and NPK mineral fertilizer helped increase total SOC content. Chen et al. [32] found that long-term use of organic fertilizer could promote the continual growth of soil organic content and that the extent of growth rose with the increase in the input of organic fertilizer. An experiment by Wu et al. [33] indicated that organic fertilizer, especially the mixture of chemical fertilizer and organic fertilizer, could enhance SOC content.

3.4. Soil Texture

To quantify significant differences in SOC content among different years, 6 groups (6 years) of SOC content were used for the variance analysis. According to Levine’s homogeneity of variance test, the mean values of SOC content were all consistent with the homogeneity of variance. Therefore, the least significant difference method was adopted.

Throughout the 6 years, the content of organic carbon in different soil texture is medium loam, light loam, sandy loam, and sandy soil from high to low (Table 2). There was a significant difference among different types of soil texture in the same year. There was an extremely significant difference in SOC content between sandy soil and medium loam as well as other types of soil texture (p⩽0.01). This indicates that soil with greater adhesion has higher organic carbon (organic) content. This result is consistent with the conclusions of previous studies [3537].

The SOC content of medium and light loams generally increased with time, and there was a significant interannual difference in the SOC content (F=0.08, p<0.05). The average extent of increase in the SOC content of medium loam was 8.26%, and the SOC exceeded 15 g·kg−1 by 2015. The SOC content of medium loam increased by 2.39% on average, and the SOC approached 14 g·kg−1 by 2015. The SOC content of was relatively low at approximately 10 g·kg−1. There was no interannual difference in the SOC for sandy loam and sandy soil, and the SOC content of sandy loam declined overall. The standard deviation of SOC for all soil textures exhibited an increasing trend. This indicates that the regional SOC content depends on human and natural effects.

3.5. Plantation System

Pesticides play an important role in agricultural production, but due to the toxicity of pesticides, such as ecological damage, water, soil, and air pollution. Pesticide residues refer to various substances such as pesticide protoplasts, toxic metabolites, and degradants that remain in agricultural products after the use of pesticides. Pesticides in farmland belong to macromolecule pollutants, and organic matter, especially humus, can degrade pesticides.

The plantation systems deployed in our study area includes vegetable fields (protected and open ones), grain fields, and orchards, and these substantially differ in terms of agricultural management.

Figure 5 shows the average soil organic carbon content in different crop cropping systems and its change over time. In general, the mean SOC content of greenhouse vegetable field, open vegetable field, grain field, and orchard field were 13.03, 12.26, 11.02, and 10.06 g·kg−1, respectively during 2010 to 2015. And we can see that soil organic carbon content in different crop systems in the same year had significant differences except 2010 and 2013. On the other hand, the SOC content of the same crop planting system varies with time. Firstly, the mean values of SOC content for the greenhouse vegetable fields from 2010 to 2015 were (12.39±5.98) a, (12.21±4.03) a, (11.51±6.25) a, (12.88±5.51) a, (14.59±7.31) a, and (15.66±8.50) a g·kg-1, respectively. With the year of 2013 as the inflection point, the average SOC content of protected vegetable fields first declined and then increased throughout the 6 years. The trend was consistent with the linear equation, and the correlation coefficient (R2) was 0.7025. Secondly, the mean value of SOC content of open vegetable fields from 2010 to 2015 were (11.23±3.85) a, (11.91±3.80) a, (12.73±5.35) a, (11.71±4.76) a, (13.30±3.62) a, and (13.93±4.06) a g·kg−1, respectively. The trend was consistent with the linear equation as a whole, and the R2 was 0.7430. Thirdly, the mean values of SOC content of orchards from 2010 to 2015 were (9.46±3.35) a, (11.07±4.43) a, (10.82±6.53) a, (11.81±6.37) a, (11.24±3.46) a, and (11.69±5.33) a g·kg−1, respectively. The SOC content of orchards showed an overall increasing trend. This trend was consistent with the linear fitting equation, with a correlation coefficient of 0.6734. Compared with other plantation systems, the SOC variation of orchards was relatively stable, and the SOC content remains approximately 11 g·kg-1. Lastly, the mean values of SOC content in the grain fields from 2010 to 2015 were (10.54±3.80) a, (9.92±2.84) a, (9.32±2.31) a, (10.54±3.38) a, (9.84±2.87) a, and (10.14±3.08) a g·kg-1, respectively. The variation of average SOC content of grain fields during the 6 years was relatively irregular and showed and overall trend of decline. The trend fitted the linear equation, showing a negative correlation.

The mean value of several years and interannual variation of the SOC content of different plantation systems showed the following pattern: protected vegetable fields > open vegetable fields > orchards > grain fields. The variation of SOC content can be roughly divided into two types. The first type showed a tendency of increase, including protected vegetable fields, open vegetable fields, and orchards, and the second type showed an overall tendency of decline, i.e., grain fields. Protected vegetable fields feature high input, high output, and multiple-cropping systems, and this is significantly different from the ecological environment of uncultivated soil. Grain fields are characterized with a low level of inputs, and most of the fertilizers were inorganic. Orchards feature high but singular input, and the management mode remained largely unchanged. In recent years, Beijing has published such documents as the Beijing Municipal Plan for the Development of Modern Agriculture during the 12th “Five-Year Plan” Period (2011-2015) and Opinions on Promoting the Development of Protected Agriculture to increase the input for protected agriculture and orchards. By 2020, the area of grain fields, vegetable fields, and orchards in Beijing will reach 53,333ha, 46,667ha, and 66,667ha, respectively, showing a “Three-Pillar” pattern. As the area of grain fields declines and the area of vegetable fields and orchards increases, corresponding inputs should be enhanced.

4. Conclusions

In this study, long-term monitoring at fixed locations in conjunction with the temporal stability coefficient, the regression slope, the EBK, spatial correlation analysis, and variance analysis were used to reveal the evolution character of the SOC content and relevant driving mechanisms in the plain areas of Beijing. The main conclusions are listed as follows:

(1) There was significant variance in the mean SOC content during the 6-year period and in general, the SOC content increases with time. (2)As a whole, the temporal stability of SOC was moderately stable. However, it was completely different from the spatial distribution of SOC. Specifically, the temporal stability of SOC is low in regions with high SOC content. (3) The application of organic fertilizer can significantly increase the soil organic carbon content. This study showed that the organic carbon content increased by 0.057 g·kg−1 for every 1 unit of organic fertilizer. (4) The grain field plantation system exhibited declining SOC content, but protected vegetable fields, open vegetable fields, and orchards all showed an increase in SOC content. (5) Regional SOC content depends on human and natural effects with human factor playing a dominant role, although soil with a higher clay content exhibits higher SOC content. (6) As a kind of polymer pollutant, engineering polymer has influence on soil properties. Clarifying the influence and mechanism of organic matter on the adsorption and migration behavior of engineering polymers in water-soil multimedia environment, and the distribution of SOC in Beijing plain area are helpful for remediation of farmland pollution in the later period.

Data Availability

The authors declare that the data supporting the findings of this study are available within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.


This work was supported by the National Key Research and Development Program of China (2016YFD0300801).