Yan Wu
Qiong Wang
Huimei Wang3
Shijie Han- 1Department of Biological Engineering, Da Qing Normal University, Daqing, China
- 2Department of Life Science, Henan University, Kaifeng, China
- 3Key Laboratory of Forest Plant Ecology, Northeast Forestry University, Harbin, China
The influence of shelterbelt afforestation on soils in different-depth profiles and possible interaction with climatic conditions is important for evaluating ecological effects of large-scale afforestation programs. In the Songnen Plain, northeastern China, 720 soil samples were collected from five different soil layers (0–20, 20–40, 40–60, 60–80, and 80–100 cm) in shelterbelt poplar forests and neighboring farmlands. Soil physiochemical properties [pH, electrical conductivity (EC), soil porosity, soil moisture and bulk density], soil carbon and nutrients [soil organic carbon (SOC), N, alkaline-hydrolyzed N, P, available P, K and available K], forest characteristics [tree height, diameter at breast height (DBH), and density], climatic conditions [mean annual temperature (MAT), mean annual precipitation (MAP), and aridity index (ARID)], and soil texture (percentage of silt, clay, and sand) were measured. We found that the effects of shelterbelt afforestation on bulk density, porosity, available K, and total P were observed up to 100 cm deep; while the changes in available K and P were several-fold higher in the 0–20 cm soil layer than that in deeper layers (p < 0.05). For other parameters (soil pH and EC), shelterbelt-influences were mainly observed in surface soils, e.g., EC was 14.7% lower in shelterbelt plantations than that in farmlands in the 0–20 cm layer, about 2.5–3.5-fold higher than 60–100 cm soil inclusion. For soil moisture, shelterbelt afforestation decreased soil water by 7.3–8.7% in deep soils (p < 0.05), while no significant change was in 0–20 cm soil. For SOC and N, no significant differences between shelterbelt and farmlands were found in all five-depth soil profiles. Large inter-site variations were found for all shelterbelt-induced soil changes (p < 0.05) except for total K in the 0–20 cm layer. MAT and silt content provided the greatest explanation powers for inter-site variations in shelterbelt-induced soil properties changes. However, in deeper soils, water (ARID and MAP) explained more of the variation than that in surface soils. Therefore, shelterbelt afforestation in northeastern China could affect aspects of soil properties down to 100 cm deep, with inter-site variations mainly controlled by climate and soil texture, and greater contribution from water characteristics in deeper soils.
Introduction
Globally, ecological shelterbelt engineering projects, such as the Great Plains Shelterbelt Project (Roosevelt Engineering) in the USA, the Great Plan for the Transformation of Nature in the former Soviet Union, forestry and water conservation projects in Japan, the Green Dam Engineering Project in the five countries of North Africa, and the Three-North Shelterbelt Program in China, have increased the scientific study of shelterbelt forests (Zhang et al., 2016). There are numerous forest plantations worldwide, many of which were planted in degraded or abandoned farmlands and are used as agricultural protection forests or bioenergy forests in China (Wang G. Y. et al., 2018; Zhang et al., 2018) and worldwide (Deniz and Paletto, 2018; Jha, 2018). The area of shelterbelt forests used for protecting soil and water increased to 330 million ha globally by 2010, accounting for 8% of all forest areas. The largest proportion of shelterbelt forests is in Asia (26%), 33% of which are in East Asia, and China's shelterbelt forests account for most of that area (60 million ha of the total 83 million ha) (Obschatko et al., 2010). There are approximately 6.67 million ha of poplar plantations that are widely distributed in China. The large shelterbelt forest area in China makes it a good example for studying the ecological functions of shelterbelt forests, and underground soil changes are an important issue to fully understand the functions of forests (Zhu, 2013; Wang et al., 2015, 2017b; Wu and Wang, 2016; Wang Q. et al., 2017; Zhong et al., 2017; Nan et al., 2018).
Black soils in northeastern China are mainly located in the Songnen Plain and Sanjiang Plain, which contain one of the three global black soil belts, and over 45% of the total grain output in northeastern China is produced in this region (Wang et al., 2009b). Although the black soils in northeastern China contain abundant soil organic matter and have high fertility compared with other soils (Cas, 1980; Hljtr, 1992), excessive historical reclamation has led to sharp decreases in soil fertility since the establishment of the People's Republic of China in 1949 (Wang et al., 1996) and nearly half of the nitrogen and soil organic matter has been lost from the black soils in northeastern China (Ding and Liu, 1980; Wang, 2002; Wang et al., 2011b). Several studies have shown that afforestation in cultivated farmland soils induced changes in most soil properties and soil fertility, contributing to soil improvement in different cases (Li and Cui, 2000; Wang Q. et al., 2014; Wang et al., 2017a). Shelterbelts of different ages and tree species could effectively reduce nitrate nitrogen by 22–60% (Jaskulska and Jaskulska, 2017), and also regulate soil physiochemical properties, fertility, and carbon sequestration (Wang et al., 2017a). In addition, soil physical properties could be altered from afforestation practices, including increases in soil bulk density and decreases in total porosity, water retention, and ventilation capacity (Wang, 2002; Wang et al., 2011b, 2017a). However, other studies also found that fast-growing plantations, such as larch, poplar, or eucalyptus, require more soil nutrients, and water (Chen, 1998; Mendham et al., 2003; Merino et al., 2004; Zhang et al., 2004; Li Y. et al., 2018), which is possibly induced by deep soil changes in various soil properties (Wang H. M. et al., 2014; Wang W. J. et al., 2014). Most of these studies have been undertaken in surface soils <40 cm deep with the assumption of neglectable changes in deep soils relative to surface soils. However, other studies have found that deep soils can sensitively react to land use changes (Fontaine et al., 2010; Rukshana et al., 2011), especially for tree species that have relatively longer roots compared to crops (Wang Q. et al., 2014; Wang S. et al., 2017).
Songnen Plain was named after the Songhua and Nenjiang Rivers running through this region. This plain has been recognized as the northern-most region of the Three-North Shelterbelt Program (Wu and Wang, 2016). Songnen Plain is about 18.28 million ha and locates in the transitional region between the semi-moist and semi-arid region, featured as saline-alkalinization and heavy farmland soil degradation (Li, 2000; Wang et al., 2011a) as well as natural forest degradation (Dai et al., 2018). Our previous study has shown that poplar shelterbelt afforestation in northeastern China slightly changed SOC sequestration and N nutrients in the surface (20 cm) soils, with sharp decreases in bulk density (Wu et al., 2018), with no consideration in deep soils (>20 cm). Moreover, glomalin-related soil carbon sequestration was higher in deep soils than that in surface soils, with more response to climatic changes in the farmlands of this region (Wang et al., 2017b). Most poplar roots concentrated in the 0–60 cm soils, and the influence of vegetation growth and microbial activities on soils may extend over the depth of the roots (Jobbágy and Jackson, 2000). Annual precipitation in Songnen Plain ranges from 300 to 500 mm, with a 2–3-fold higher annual evaporation (1,000–1,500 mm) (Li, 2000). This natural background, heavy pressure from farming and grazing, and fast saline-alkalinization in soil are important challenges for social development and livelihood in this region (Li, 2000; Wang et al., 2011a). The evaluation of shelterbelt afforestation on underground soils in this region must fully consider the variations in the widespread plain, and fully understanding of the underlying mechanisms needs more consideration on forest characteristics, climatic conditions including the aridity index (ARID), soil texture both at surface and deep soils (Wang W. J. et al., 2014; Wu et al., 2018).
In the present study, we alleged that deep soils at 100 cm depth should be included in the evaluation of various soil changes in poplar afforestation, and large inter-site variation in the shelterbelt forest-induced soil changes were related to local climatic differences, soil texture, and forest growth. We posed several research questions as follows: (1) Should deep soil layers be included in the evaluations of soil improvements from degraded farmlands to poplar forests and did these improvements differ in different soil parameters? (2) How great a difference among locations occurred in the shelterbelt-induced soil changes, and which factors of climatic condition, soil texture, and tree growth parameter were responsible for these variations? By evaluating the shelterbelt-induced soil changes in various properties in different soil layers, our data assisted the evaluation of underground soil changes in large scale shelterbelt programs, such as the Three-North Shelterbelt Program, particularly the quantification of the importance of deep soils for afforestation practices in degraded farmlands.
Materials and Methods
Study Sites and Sample Collection
The Three-North Shelterbelt Program established tree plantations around farmlands in northern China, northwestern China, and northeastern China in 1978 (Zhu, 2013). The general design was to plant 4–10 rows of poplars around 500 × 500 m of farmland, and large areas of shelterbelts around farmlands are found everywhere throughout the Songnen Plain of northeastern China (Figure 1). Nowadays, the most-used poplar variety in northeast regions (young forests) is Yinzhong poplar (Populus alba × Populus berolinensis), while historically, the most-planted poplars were Populus simonii, Populus × xiaohei, and Populus deltoides × P. canadensis, etc. (Wu et al., 2018).
Figure 1. Six study sites in the Songnen Plain, Northeastern China, and a typical poplar shelterbelt-farmland paired site. Parts of this figures was adapted from our previous publication (Wu et al., 2018).
Six study sites (Dumeng, Fuyu, Lanling, Mingshui, Zhaodong, and Zhaozhou) distributed in the Songnen Plain in the middle of northeastern China were selected as study sites (Figure 1). The soil types in the study region are typical black soils, including Chernozem (Fuyu, Lanling), Phaeozem (Mingshui), Cambosols (Dumeng), and some degraded soil, such as Solonetz (Zhaozhou, Zhaodong). This region has a continental monsoon climate, with MAT of 2.9–4.4°C, MAP of 350–500 mm, and ARID of 0.4–0.7.
Soil samples were collected from 72 paired shelterbelt plantations and farmland plots in the six study sites. Five soil profiles were collected from each paired plot. After the exclusion of the A0 layer, we sampled 100 cm of soil from 0 to 20, 20 to 40, 40 to 60, 60 to 80, and 80 to 100 cm depths. Additionally, we obtained a composite sample by mixing five samples from the same soil layers from each of the five soil profiles. In total, 720 soil samples (6 regions × 2 (farmland and shelterbelt) × 12 sites/region × 5 depths/site = 720 samples) were collected.
Determination of Soil Parameters
Sample preparation details and some of the soil parameters analysis (e.g., bulk density, soil moisture, SOC, total N, available N, total K, available K, total P, available P, and soil texture) have been described previously by Wu et al. (2018). Soil porosity was calculated by the following formula: soil porosity = (1−bulk density/specific gravity) × 100%. The pH of the soil solution (one-part soil to five-parts water) was measured with an acidity meter (Sartorius PT-21, Shanghai, China). Soil electrical conductivity (EC) was determined with an EC meter (DDS-307, Shanghai Precision Scientific Instruments Co., Ltd., Shanghai, China) (Bao, 2000). Soil carbon or nutrient storage were computed as:
Farmland soil carbon or nutrient storage = αf × ρf × 0.2 × (1–Vgravel)
Poplar soil carbon or nutrient storage = αp × ρp × 0.2 × ρf/ρp × (1–Vgravel)
where, αf and αp are the concentrations of farmland and poplar SOC (g kg−1); ρf and ρp are farmland and poplar soil bulk densities (Mg m−3), respectively; 0.2 is the soil thickness (0.2 m); and Vgravel is the proportion of gravel. Details regarding the bulk density correction can be found in Wuest (2009) and Wu et al. (2018).
Forest Characteristics, Soil Texture, and Climatic Data Collection
Poplar forest characteristics of tree density, tree height, and DBH were measured at each plot site. Regarding the distance between forest and farmland, 42% (30 plots) of plots were <3.4 m from neighboring farmland, while 47% (34 plots) were 3.4 to 6.7 m away. To reduce the influence of roots on neighboring farmlands, ditches of about 2 m in width and 2 m in depth were excavated between shelterbelt and farmland by local farmers. Forest characteristics data of 72 plots can be found in Wu et al. (2018).
The soil texture was the relative amount of sand, silt, and clay in the bulk soil, measured using a rapid and simple method described by Kettler et al. (2001) and Wu et al. (2018).
MAT and MAP at the six sites were obtained from the meteorological scientific data sharing service network of China (http://cdc.cma.gov.cn/) and the ARID was computed as the MAP over the mean annual reference evapotranspiration (Huo et al., 2013).
Data Analysis
Calculation of storage at different soils depth, such as 0–20 cm, 0–40 cm, 0–60 cm, 0–80 cm, and 1 m soil profiles, is a general rule for many previous studies for ease of comparison among studies (Wang et al., 2011b; Wei et al., 2014; Wang H. M. et al., 2017a; Deng et al., 2018). In this paper, in order to compare soil carbon or nutrient storage with other studies, we have re-grouped our data from 0 to 20, 20 to 40, 40 to 60, 60 to 80, and 80 to 100 cm into 0 to 20, 0 to 40, 0 to 60, 0 to 80, and 0 to 100 cm by combining the corresponding soil layer's measured data. For example, 0–20 and 20–40 cm were combined into one set of 0–40 cm by average of two data; Similarly, combining 0–40 and 40–60 cm into 0–60 cm, combining 0–60 and 60–80 cm into 0–80 cm, and combining 0–60 and 80–100 cm into 0–100 cm during dada analysis in order to evaluate the effects of afforestation on soil properties in five soil-depth profiles.
Multivariate analysis of variance (MANOVA) was used to determine the influence of land use type (shelterbelt forest and neighboring farmland), sampling location (Dumeng, Fuyu, Lanling, Mingshui, Zhaodong, and Zhaozhou), and their interaction on various soil parameters. The 19 parameters (soil bulk density, soil porosity, soil moisture, pH, EC, SOC concentration, total N concentration, alkaline hydrolyzed N concentration, total P concentration, available P concentration, total K concentration, available K concentration, SOC storage, total N storage, alkaline hydrolyzed storage, total P storage, available P storage, total K storage, and available K storage) were used as dependent variables.
A paired t-test was used to determine the difference in soil properties between shelterbelt plantations and farmlands at different soil depths, and the Duncan's test was used for multiple comparisons among different soil-depth profiles for all shelterbelt-induced soil changes. In the present study, the relative change [(forest—farmland)/farmland] of each soil parameter was treated as a dependent variable for the following analysis.
Redundancy analysis (RDA) was conducted to ordinate the complex associations between shelterbelt-induced variations in various soil properties and climatic conditions, soil texture, and forest characteristics (Canoco 5.0 software program). Ordination was performed in all five soil-depth profile, as we wanted to find the differences in different soil depths. Conditional term effects (excluding collinear effects among different dependent parameters) were derived from the RDA, and the possible factors contributing to the dependent variables (e.g., shelterbelt-forest-induced changes in soil properties as a whole). In conditional term effects, the significant factor with the highest explanation percentage showed the strongest contribution to the variation of soil properties. Details explaining the RDA ordination can be found in previous studies (Wang et al., 2017b, 2018b).
Stepwise regression analysis was used to explore the factors responsible for the poplar-induced changes in various soil properties in the five soil layers. Statistical significance was evaluated at p = 0.05, unless otherwise stated. Three groups of parameters, forest characteristics (tree height, DBH, and tree density), soil texture (sand, silt, and clay), and climatic conditions (MAT, MAP, and ARID), were tested as independent parameters. The more entering times, and more frequent parameter into the stepwise models indicate the stronger influence from these parameters for explaining the variations of soil properties from shelterbelt afforestation. By using this criteria, stepwise regression models were analyzed for simplifying the presentation of the data and facilitate data interpretation.
Results
Land Use Type and Sampling Location Affect All Soil Parameters: Manova Results
Table 1 showed the influences of soil use type, sampling location, and their interaction on the soil parameters at five depths. Significant land use effects (farmland and poplar forest) were observed in bulk density, porosity, and total P storage for all five depths. Significant differences between the two land uses on pH, EC, available K concentration, and available K storage were found in the 0–20 cm depth, whereas others, such as soil moisture, total K, total P, available P concentration, and available P storage were statistically different among the two land uses in the deeper soil layers (>20 cm).
Table 1. Shelterbelt plantation establishment, sampling regions, influences various soil parameters, and possible interacts at different depths.
Compared with land use differences, there were even larger significant location-related differences among all parameters in all five soil depths. Moreover, significant interactions existed among the influence of land use and sampling location on some soil parameters in different layers (Table 1). For example, the influence of shelterbelt plantations on porosity, total K concentration, and available P concentration significantly interacted with location in the five depths, indicating that these shelterbelt-induced changes significantly differed among the six locations in all soil layers.
Changes in Soil Properties Between Shelterbelt Plantations and Farmlands at Five Soil-Depth Profiles: Overall Patterns
The effects of shelterbelt construction on soil properties in the five soil-depth profiles and the differences in shelterbelt-induced soil changes in various parameters among the five profiles were shown in Table 2.
Table 2. A comparison in soil properties between shelterbelt plantation and farmland at five depths and the differences in shelterbelt-induced soil changes among five profiles.
The effects of shelterbelt construction on soil properties varied in five soil depth layers. Some indicators, such as bulk density, porosity, available K concentration, and total P storage, had significant difference at five soil depth. However, some indicators, such as SOC concentration (storage), total N concentration (storage), and available N concentration (storage) had no significant change on five profiles following shelterbelt establishment. Moreover, the effects of afforestation on the surface soils were more obvious than at depth for pH and EC. On the contrary, the effect of soil moisture was seen mainly in the deeper soil profiles (Table 2).
The significances among the five soil-depth profiles were distinct for different shelterbelt-induced soil properties changes. First, opposite trends (p < 0.05) were observed in shelterbelt-induced soil moisture and EC changes between surface and deeper soil layers. A 6.2% increase in soil moisture (poplar compared with farmland) was observed in the 0–20 cm layer, whereas there was a 7.3–8.7% decrease in the deep soil profiles. Contrary to soil moisture, EC was 14.7% lower in shelterbelt plantations than that in farmlands in the surface layer, but was 4.1, 6.2, and 4.2% higher in the 0–60, 0–80, 0–100 cm layers, respectively (Table 2). Second, the changes in available K and P were several-fold higher in the 0–20 cm soil profiles than that in the deeper profiles (p < 0.05). For example, a 117.4% increase in available K concentration was observed in the surface layer, whereas only a 24.6–39.2% increase was observed in the deeper layers. A 28.3% increase in available P concentration in shelterbelt plantations was observed in the surface layer, whereas a −0.2–6.1% change was found in the deeper layers. Third, no significant changes (p > 0.05) were found among the five soil profiles for the other properties.
Shelterbelt-Induced Soil Parameter Changes: Large Inter-site Variations Differed With Soil Depth
All soil parameters except total K concentration and total K storage showed marked location-related differences in the five depths among the different sites (Table 1). The vertical pattern and magnitude of differences are shown in Figure 2 and Table A1.
Figure 2. Changes of various soil properties in poplar forests compared with neighbor farmland, and differences at 1 m profiles. Dash line in the figure showed the zero line, indicating that no changes relative to neighbor farmland. Statistics of the mean values had shown in Table 2 and Table A1. (A) Soil bulk density change %; (B) Soil water change %; (C) Soil EC change %; (D) Soil porosity change %; (E) Soil pH change %; (F) Soil organic carbon concentration change %; (G) Soil total N concentration change %; (H) Soil alkaline hydrolyzed N concentration change %; (I) Soil total K concentration change %; (J) Soil alkaline K concentration change %; (K) Soil total P concentration change %; (L) Soil alkaline P concentration change %; (M) Soil organic carbon storage change %; (N) Soil total N storage change %; (O) Soil alkaline hydrolyzed N storage change %; (P) Soil total K storage change %; (Q) Soil alkaline K storage change %; (R) Soil total P storage change %; (S) Soil alkaline P storage change %.
The location differences were soil-depth dependent, i.e., in most cases, the surface soil layer showed much larger location differences in the poplar forest-induced changes in various soil properties (Figure 2). For example, inter-site differences in soil bulk density and soil water were −10 to 5% and −25 to 60% respectively, whereas those in deep soils (0–100 cm) were, respectively −9 to −1% and −20 to 0%. However, for other parameters, similar inter-site differences were found among surface and deep soils, with even larger variation in deep soils, for example, soil available P and SOC for both concentration and storage (Figure 2).
Depth-induced significant differences were observed in soil moisture, pH, EC, available K concentration and storage, and available P concentration and storage (Table 2), with large inter-site variations found among the different sites (Figure 2 and Table A1). There was a higher amount of available K concentration in the surface layer than that in the other four soil layers (Table 2). This trend mainly occurred in Dumeng, Fuyu, Mingshui, and Zhaodong (p < 0.05), with a 221.8% increase in the 0–20 cm layer and an average 0–20 cm layer and an average 46.8% increase in the other four profiles at Dumeng, whereas there was a 79.9% increase in the surface layer and an average 15.9% increase in the deeper profiles at Mingshui (Table A1). Moreover, although no significant differences were observed in three of the soil parameters (porosity, total K concentration, and total K storage) among the five soil layers (Table 2), inter-site differences were found among the different sites. For example, there were significant differences (p < 0.05) among the five soil profiles in total K concentration and storage at Dumeng, Fuyu, and Mingshui, with a 24.1% decrease in the 0–20 cm layer and an average 1.8% increase in the deeper profiles at Dumeng, and a 42.8% decrease in the surface and an average 17.5% decrease in the deeper profiles at Mingshui (Table A1).
RDA Ordination: Climatic, Soil Texture and Forest Controls on the Inter-site Variations and Differences Between Surface and Deep Soils
As shown in Table 3, in general, climatic conditions provided the largest explaining power for the inter-site variations of shelterbelt-induced soil changes. Moreover, at different soil layers, MAT was the most influential parameter, providing the highest explanation percentage. For example, MAT explained 4.9, 5.1, 4.4, 3.0, and 3.5% of the forest-induced soil variations for 20, 40, 60, 80, and 100 cm depth soils, respectively (Table 3). In deeper soils, ARID and MAP explained much more of the variation than that in the surface soil layer. For example, ARID in 20 cm, 40 cm soils did not show significant explanation for the variations, while in 60, 80, and 100 cm soils, ARID showed significant explaining powers ranged from 3.0–4.6% (p < 0.05); and MAP explained 6.9% of the variations for 100 cm soil layers (p < 0.01) (Table 3).
Table 3. Comparison on the explaining power from climatic condition, soil texture, and forest characteristics for the forest-induced soil changes at different locations from the RDA ordination-related conditional term effects excluding their collinear effects.
Following the climatic conditions, soil texture gave the next largest explanation power for the location-related variations and silt showed significant explanation percentage in four out of five soil layers (p < 0.05) (Table 3). In general, the deep soils, the less explaining power from soil texture of silt percentage. In 20 and 40 cm soils, silt's explaining percentage was 4.5–7.8% (p < 0.001), and this percentage was 3.3–3.5% in 60 and 80 cm soils (p < 0.05), and no significant explaining percentage was found in 100 cm soil (p = 0.124) (Table 3).
In addition, tree growth traits (DBH and height) also significantly explained the shelterbelt-induced soil variations at different locations, and their explaining powers ranged from 2.5 to 3.6% at different soil layers (p < 0.05) (Table 3).
Stepwise Regression Statistics: Factors Related to the Inter-site Variations and Differences Between Surface and Deep Soils
As a further step to decouple the association, stepwise regression analysis was used to determine the most possible parameters for the large inter-site variations (Table A2).
For each soil properties, we found different associations with climatic conditions. For example, positive correlations were observed between shelterbelt-induced soil bulk density change and MAT in all layers expect for the 0–40 cm depth. Although soil texture and tree growth significantly accounted for shelterbelt-induced porosity change in the 0–20 cm layer, MAT was the leading factor that determined the porosity change in the other four soil depths (r2 = 0.10–0.18, p < 0.05). The pH change was significantly related to MAT in the 0–20 and 0–40 cm depth layers, and the shelterbelt-induced EC decrease was accompanied with higher MAT in all five soil profiles. Significant positive correlations were found between SOC storage change, total K concentration (storage) changes, and MAT (p < 0.05). In the deep soil layers (0–80 and 0–100 cm), ARID was the significant affecting factor (r2 = 0.32, p < 0.001) for bulk density changes compared with the surface soil layer. Similar significant negative correlations were found between available P (concentration and storage) and ARID in the deeper soil profiles. Higher ARID accounted for the poplar-induced total K decrease in the five soil profiles (r2 = 0.15–0.46, p < 0.001), whereas the available K (concentration and storage) changes in the 0–100 cm depth could be explained by ARID and MAP (r2 = 0.39 and r2 = 0.40, respectively) (Table A2).
For soil silt percentage, significant positive correlations were found among silt percentage and EC, available P concentration, whereas marked negative correlations were found in the shelterbelt-induced differences in seven soil parameters in the different soil depths (including porosity, total N concentration, alkaline hydrolyzed N concentration, available K concentration, total N storage, alkaline hydrolyzed N storage, and available K storage) (Table A2). For tree growth parameters, DBH, Tree height, and Tree density have been found in different stepwise regression models in different soil layers; However, their appearances were not as often as that of climatic parameters and soil textures (Table A2).
By counting the entering times for each tested parameter observed in all stepwise regression models, we want to confirm the findings in RDA ordination, and the basic criteria is that the more entering times mean the more influences on soils from this parameter (Figure 3). Comparison among climate, soil texture and tree growth, we found the most entering times from climatic factors (10–12 entering times), followed by soil texture (3–7 entering times), and tree growth factors (0–7 entering times); This is the similar to those observed in RDA ordination (Table 3). In the case of different climatic parameters, we found that MAT showed the most influences (5–8 entering times), followed by ARID (2–5 entering times) and MAP (1–2 entering times) (Figure 3). At the vertical soil profiles, MAT's influences decreased from surface to deep soils, as shown by eight entering times in 20 cm soils and five entering times in 100 cm soils. However, ARID's influences showed a contrary pattern, i.e., lower influences were at surface soils (two entering times in 20 cm), while much stronger influences were in deep soils (five entering times) in 80 cm and 100 cm soils (Figure 3).
Figure 3. Differences of entering times of climatic factors (upper), soil texture (middle), and forest traits (lower) observed in the all stepwise regression models at different depths. Inset tables are the most observed parameters (MAT, ARID, MAP, Silt, and DBH) and their entering times all stepwise models at different soil depths. All the stepwise regression models were shown in Table A2. The more inclusion of the parameters into the stepwise models indicates their stronger contribution at that soil layers for explaining the forest-farmland differences in the studied soil properties.
Soil texture showed the similar entering times at the surface (5–7 entering times) and deep (3–7 entering times) soils, and 2–4 times entering into stepwise models were observed from silt percentage at five depths soil (Figure 3). In the case of DBH, two entering times were found in 100 cm soil, while 0–1 time entering to the stepwise models was found in other soils (Figure 3).
Discussion
Higher Water Consumption in Deep Soils, Saline-Alkalinization in Surface Soils, and Better Physical Structure in the Entire 100 cm Soil Profile Following Shelterbelt Plantations
Our study highlighted significant water consumption in deeper soil layers in poplar forests in northeast China. There were significant decreases (7.3–8.7%) in 0–40 cm and deeper soils (p < 0.05), while a slight increase was found in 0–20 cm layer (Table 2). Artificial afforestation can decrease soil moisture because of leaf interception and root uptake (Jin et al., 2011). Divergent hydrological response to large-scale afforestation has been found on a national scale throughout China (Li Y. et al., 2018). Previous studies have shown that soil moisture differed significantly between traditional farmland and introduced woody vegetation (Wang et al., 2009a, 2011b; Liu et al., 2010; Li Y. et al., 2018). In north and southeast China, the increased precipitation and increased forest area were not statistically significant, and had only a weak influence on soil moisture content (Li Y. et al., 2018). In southwest China, however, the afforestation practices have been shown to significantly reduce soil moisture in combination with decreased precipitation (Li Y. et al., 2018). In northeastern China (the same region as in the present study), soil moisture has been shown to be significantly decreased by −8.1 mm decade−1 (Li Y. et al., 2018). High water consumption following afforestation has been reported as an important feature of fast-growth tree species afforestation (Yang et al., 2012; Jia et al., 2017; Liang et al., 2018). Songnen Plain (the central part of northeastern China) is characterized as experiencing land degradation with a shortage of precipitation (Li, 2000), and the high water consumption from poplar shelterbelt plantations possibly intensifies the degree of drought in the deeper soil depths, with an average precipitation of 400–500 mm and a large area of saline-alkalinization land (Zhang et al., 2013). Currently, measures used by local people to prevent this water consumption include digging a root-cutting ditch to hinder root invasion into farmland. Possible other measures to counteract the over-consumption of water by plantations forests have been proposed (Ferraz et al., 2013); for example, the selection of suitable tree species with low water utilization, such as Picea spp. (e.g., P. jezoensis), plays an important role in the reduction and regulation of water use (Wang et al., 2017a), and an increased proportion of native forests and mosaic management could also stabilize water flow across plantation landscapes.
Saline-alkalinization is an important feature of local land degradation in the Songnen Plain (Li, 2000) and different methods have been invented for soil improvement and afforestation practices in this region (Wang et al., 2011a). Our study found that, compared to farmland, the establishment of shelterbelts increased soil pH in the 0–20 cm and 0–40 cm soil layers (p < 0.001) (Table 2), whereas there were no significant differences in the deeper soil layers. Therefore, poplar afforestation resulted in surface soil saline-alkalinization in the Songnen Plain of China. Previous meta-analyses have found site-scale soil acidification globally (Berthrong et al., 2009) or afforestation-induced soil neutralizing pH that favored acidifying alkaline soils (Hong et al., 2018). Our findings were different from these meta-analyses, which may be related to the following. First, the addition of plant residues can increase, decrease, or have little effect on soil pH (Tang and Yu, 1999; Marschner and Noble, 2000; Xu et al., 2006; Rukshana et al., 2011), and are mainly dependent on the amount of returning organic materials. In poplar shelterbelt plantations, the relatively smaller area (i.e., several rows) could result in limited litter decomposition and rhizospheric processes following afforestation. Second, the increases in evapotranspiration (Yao et al., 2016) caused by poplar afforestation reduces the leaching loss of base cations (Slessarev et al., 2016), and thus increases soil pH. Third, the upward vertical movement of water by deep-rooted trees (compared with crops) could generally induce the upward movement of soluble salts from the deep soils to the surface soils, resulting in soil saline-alkalinization (Li, 2000; Lu et al., 2017; Wang et al., 2017a).
Planting shelterbelts with fast-growing species such as poplar causes soil bulk density to be significantly reduced in the 0–20 cm soil layer (Wu et al., 2018), and similar significant improvement in soil physics (e.g., bulk density decreasing and porosity increasing) were observed over the entire 100 cm soil profile including 0–40, 0–60, 0–80, and 0–100 cm soil layers (Tables 1, 2, p < 0.05). Soil physical structure is very important for soil function (Han et al., 2018) and previous studies have found surface soil improvements, for example, Marta and Halina (2008) observed that total porosity, on average, in the entire 20 cm horizon of the studied afforested soils was 1.08 and 1.12 times higher (for soils of young and older stands, respectively) than that in the arable soils. Our results are in agreeance with previous findings, emphasizing the importance of shelterbelt afforestation to deep soil layers. The long-term farmland cultivation in northeastern China has seriously degraded black soils, with one important aspect being the degradation of the soil physics (Li, 2000; Wang, 2002). Our results clearly show that shelterbelt afforestation could strongly improve soil physics and suggests a possible way for local soil improvement, such as returning degraded farmland to forests, with such policy being implemented in China over the past years (Wang et al., 2011b).
Non-accrual of SOC Both in Surface and Deep Soil Layers Following Shelterbelt Forest Establishment
Results from various studies on the effects of afforestation on SOC are inconsistent. Some studies have found that afforestation increased SOC accumulation (Lemma et al., 2006; Wang et al., 2011b; Wei et al., 2012; Cukor et al., 2017), whereas other studies have shown that afforestation decreased SOC (Farley et al., 2004; Mao et al., 2010) or there was more initial loss than SOC gain (Paul et al., 2002; Wang et al., 2006; Ritter, 2007). In this paper, we did not find any significant changes between poplar shelterbelt plantations and neighboring farmlands.
According to our survey (data not shown here), 28 mg cm−3 poplar root system was in 1 m depth, and 95% of the root was distributed in 0–60 cm soil layer, especially in 20–40 cm (57%). Previous studies have also highlighted the possible differences in different soil layers, owing to the root differences between trees (long roots) and crops (short roots) (Wang H. M. et al., 2014). For example, Hooker and Compton (2003) found that SOC linearly accumulated in the subsoil (20–70 cm), but did not differ in the top 20 cm after afforestation. Wang H. M. et al. (2014) reported that, in larch forest plantations, the rate of change in SOM in the surface soil was 262.1 g kg−1year−1; however, a different trend in deeper soils resulted in no evident changes in the overall 80 cm soil profile. Our previous paper found no shelterbelt-induced SOC accumulation in the 0–20 cm soil layer (Wu et al., 2018). In the present study, we confirmed a similar finding (i.e., no significant SOC changes) in the 0–40, 0–60, 0–80, and 0–100 cm soil layers (Tables 1, 2, Table A1).
Possible reasons for the above-mentioned patterns include the following. First, the shelterbelt poplar planting area was generally 4–6 rows of poplars around large farmlands of ~25 ha in size (Figure 1). In this type of shelterbelt forest, the canopy litter is usually deposited on both farmland and forest simultaneously, which reduces the influence of shelterbelt poplar to forest soils with reference to neighboring farmlands. Second, high productive crops and tillage practices diminish the differences between farmland and forests. It is generally assumed that forests can improve soil carbon sequestration; however, different crops and tillage practices change this sequestration. In the present study, the high productive C4 crop (maize) was the main crop in this region and soybean was the second largest crop with high N-fixation ability (this N fixation favors SOC sequestration) (Lian et al., 2017; You et al., 2017). Proper chemical fertilizer utilization together with straw returning, which has been strongly implemented by the local government, improves the stabilization and accumulation of SOC (Li and Han, 2016). No chemical fertilization or organic manure were applied in the management of the shelterbelt forest, which was different from the neighboring farmlands. Third, microbial priming-induced SOC loss possibly also contributed to the patterns (Li L. J. et al., 2018). The root-exudate inputs in the deep soil could stimulate the decomposition of SOC by priming soil microbial activity (Marie-Anne et al., 2014), and the SOC mineralization might stimulate loss of the deeper SOC pool (Fontaine et al., 2007), which possibly resulted in the non-accrual of SOC storage in the deeper profiles in the poplar shelterbelts.
Non-evident Changes in all Nutrients Except Available K Recovery in 5 Soil Layers and Total P Depletion in Deep Soils Following Shelterbelt Afforestation
Shelterbelt-induced available K recovery was found in all five soil depths (p < 0.05) (Table 2), showing that K accumulation was not only in the surface soils but was also in the deeper soils. Returning farmland to forest could rehabilitate the soil K fertility in different areas worldwide using different tree species (Likens et al., 1994; Romanowicz et al., 1996). In China, Jiao et al. (2012) observed that available K was significantly higher in afforested sites than that in degraded croplands in the Loess Plateau. In the case of farmland fertilization practices, more K fertilizer together with N (the favorite fertilizer of local people) should be applied to ensure soil nutrient supply for crop productivity.
In the present study, shelterbelt-induced total P depletion was observed mainly in the deeper soil layers (p < 0.05) (Table 2). In large areas of larch plantations in northeastern China, Wang W. J. et al. (2014) observed that more P was stored in deeper soil layers, and >70% of P (total and extractable) was found in deeper soil layers (20–80 cm) during larch reforestation. The development of larch plantations could result in a general uplifting of SOM, N, and P based on vertical distribution data and this redistribution was accompanied by the depletion of N and P. The soil nutrient depletion could be related to the biological uplifting function and possible absorption related to tree growth (Jobbágy and Jackson, 2004). In addition, the differences in microbial decomposability between deep and surface soils might strengthen the depletion in deep soil layers (Fontaine et al., 2007; Xiang et al., 2008).
All tested nutrients including N, available N, available P, and total K did not significantly change throughout the entire 100 cm soil profile in shelterbelt poplar forests with reference to neighboring farmlands in a large-scale field in this paper (Table 2). A general observation from field surveys is that the growth of crops near the shelterbelt poplar is smaller and general assumptions are that shelterbelt afforestation can decrease soil nutrients owing to nutrient competition between poplar and crops. By the entire soil profile measurements in the present study, we updated this assumption and found that the most likely nutrient depletion was that of P depletion. However, for almost all other nutrients, such depletion was not found either in surface soils or deeper soils. This should be taken into consideration in future shelterbelt forest evaluations.
Large Inter-site Variations Closely Associated With Climatic Conditions and Deep Soils Showed Greater Dependence on Arid and Map: Implications
To determine the differences between plantation forest and neighboring farmlands, many previous studies have looked at a relatively smaller region to minimize the inferences of inter-site variation on the forest effects (Mao et al., 2010; Wang et al., 2011b; Wei et al., 2012; Cukor et al., 2017). In the present study, large-scale sampling was undertaken in Songnen Plain (at least 33,000 km2) to determine the general soil change patterns in a 100 cm soil profile. Large inter-site variations in the shelterbelt-induced soil changes were found (Figure 2 and Table A1). Currently, global climatic changes strongly affect local development including natural processes in farmlands, pastures, and forests (Li, 2000; Li Y. et al., 2018) and new developments has been reported in tree inventory methods (Wang et al., 2018a) and complex association analysis (Lv et al., 2018; Wang et al., 2018b; Yang et al., 2019). Decoupling the contribution of different components on inter-site variations of the shelterbelt-induced soil changes will favor the mechanical understanding of underlying processes, and is a possible strategy for large scale evaluation of shelterbelt poplar ecological functions (Wang et al., 2018b; Yang et al., 2019).
To identify the possible contributions from climatic conditions, soil texture, and tree growth that affected soil property changes following afforestation, statistical methods including RDA and stepwise regression analysis were applied, which have been proven as beneficial for determining the causal relationship between patterns (Eisenhauer et al., 2015; Wang et al., 2018b; Yang et al., 2019). Accordingly, we found that the inter-site variations in the shelterbelt-induced soil changes in surface soils were different from those in deeper soils.
In general, MAT was the most important parameter with the highest explanation powers, and changes in seven parameters (bulk density, porosity, soil moisture, pH, EC, SOC stock, and total K in different soil layers) were significantly accounted for by the variations in MAT (Figure 3, Table 3, Table A2). This is related to the fact that temperature is one of the main factors limiting soil nutrients following the planting of plantations (Wiesmeier et al., 2013). Similar to our conclusion, climate warming has been strongly associated with a decrease in the accumulation of glomalin proportion to total soil carbon in soils (Wang et al., 2017b). Moreover, RDA and stepwise regression analysis indicated that ARID and MAP gave greater explanation percentages in deeper soils than those in surface soils (Figure 3, Table 3, Table A2).
Potential evapotranspiration (PE) and aridity index (ARID) based on the observational data from 1961 to 2004 from 94 meteorological stations showed a general increasing trend in MAT, MAP, PE, and ARID (Zhao et al., 2007). Moreover, increases were more significant in MAT and PE than that in MAP and ARID (Zhao et al., 2007), showing that northeastern China is the most serious region experiencing global changes, especially for temperature warming. The ARID range of the six locations in Songnen plain is 0.49–0.67 and the semi-arid climate possibly affects the shelter-induced soil properties changes. Our findings indicate that global changes will shift the shelterbelt-induced soil changes with reference to neighboring farmlands. Moreover, different changes have been found between surface and deep soils. Compared with surface soils, the drying trends will give more influences in deep soils owing to the large explanation powers compared with the surface soils (Figure 3, Table 3, Table A2). The Three-North Shelterbelt Program has been evaluated as the most important natural environmental rehabilitation program in China (Bryan et al., 2018), and our findings highlight that soil changes should be carefully considered during any evaluation and divergent responses to climatic changes should be included in risk assessments.
Conclusion
By analyzing 720 soil samples from 72 paired sites of poplar shelterbelts and farmlands in Songnen Plain in northeastern China, we concluded the following: (1) Shelterbelt poplar plantations significantly improved soil physical properties by decreasing bulk density and increasing porosity down into the 100 cm depth; however, higher water consumption was mainly found in the deep soils and soil saline-alkalinization was mainly in the surface soils. (2) There were no evident changes in all nutrients except for available K recovery following shelterbelt afforestation in all five soil depths and shelterbelt-induced total P depletion occurred mainly in deep soils. (3) Large inter-site variations were found for all shelterbelt-induced soil changes (p < 0.05) except for total K in the 0–20 cm soil layer, and MAT and soil texture were the largest explanation powers for soil property changes in the different soil layers. However, in deeper soils, soil drought (ARID and MAP) gave more explanation percentages than that in surface soils. Our findings highlight that shelterbelt poplar plantations could divergently change different soil properties in different soil depths, and inter-site variation was strongly associated with climatic changes. Our findings favor shelterbelt poplar forest evaluation and the underlying reasons for the large inter-site variation could help find suitable parameters to reduce the uncertainty of future evaluations.
Author Contributions
WW, SH, YW, QW, and HW conceived and designed the experiments. WW and YW contributed reagents, materials, and analysis tools. QW, YW performed the experiments. WW, SH, QW, and YW analyzed the data. WW and YW wrote the paper. All authors approved the submitted and final versions.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
This study was supported financially by China's National Foundation of Natural Sciences (41877324, 41730641), Fundamental Research Funds for the Central Universities (2572017DG04, 2572017EA03), 13–5 Key Research and Development Project from China Ministry of Science and Technology (2016YFA0600802), Heilongjiang Province for Distinguished Young Scholars (JC201401), and Doctoral Scientific Fund Project of Da Qing Normal University (17ZR03). Thanks are due to Mr. Zhong Zhaoliang and other students (Mr. Du Hongju and Ms. Zhang Jianyu) for their help in laboratory analysis and field sampling.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00220/full#supplementary-material
References
Bao, S. D. (2000). The Method of the Soil and Agriculture Chemical Analysis. Beijing: China Agriculture Press.
Berthrong, S. T., Jobbagy, E. G., and Jackson, R. B. (2009). A global meta-analysis of soil exchangeable cations, pH, carbon, and nitrogen with afforestation. Ecol. Appli. 19, 2228–2241. doi: 10.1890/08-1730.1
Bryan, B. A., Gao, L., Ye, Y., Sun, X., and Connor, J. D. (2018). China's response to a national land-system sustainability emergency. Nature 559, 193–204. doi: 10.1038/s41586-018-0280-2
Cukor, J., Vacek, Z., Linda, R., and Bílek, L. (2017). Carbon sequestration in soil following afforestation of former agricultural land in the Czech Republic. Central Eur. Forest. J. 63, 74–104. doi: 10.1515/forj-2017-0011
Dai, L. M., Li, S. L., Lewis, B. J., Wu, J., Yu, D. P., et al. (2018). The influence of land use change on the spatial–temporal variability of habitat quality between 1990 and 2010 in Northeast China. J. Forest. Res. 1–10. doi: 10.1007/s11676-018-0771-x
Deng, L., Wang, K. B., Zhu, G. Y., Liu, Y. L., Chen, L., et al. (2018). Changes of soil carbon in five land use stages following 10 years of vegetation succession on the Loess Plateau, China. Catena 171, 185–192. doi: 10.1016/j.catena.2018.07.014
Deniz, T., and Paletto, A. (2018). Effects of bioenergy production on environmental sustainability: a preliminary study based on expert opinions in Italy and Turkey. J. Forest. Res. 29, 1611–1626. doi: 10.1007/s11676-018-0596-7
Ding, R. X., and Liu, S. T. (1980). A study on the fertility of black soil after reclamation. Acta Pedol. Sin. 1, 20–32.
Eisenhauer, N., Bowker, M. A., Grace, J. B., and Powell, J. R. (2015). From patterns to causal understanding: structural equation modeling (SEM) in soil ecology. Pedobiol. J. Soil Ecol. 58, 65–72. doi: 10.1016/j.pedobi.2015.03.002
Farley, K. A., Kelly, E. F., and Hofstede, R. G. M. (2004). Soil organic carbon and water retention following conversion of grasslands to pine plantations in the Ecuadorian Andes. Ecosystems 7, 729–739. doi: 10.1007/s10021-004-0047-5
Ferraz, S. F. B., Lima, W. D. P., and Rodrigues, C. B. (2013). Managing forest plantation landscapes for water conservation. Forest Ecol. Manage. 301, 58–66. doi: 10.1016/j.foreco.2012.10.015
Fontaine, S., Bardoux, G., Abbadie, L., and Mariotti, A. (2010). Carbon input to soil may decrease soil carbon content. Ecol. Lett. 7, 314–320. doi: 10.1111/j.1461-0248.2004.00579.x
Fontaine, S., Barot, S., Barré, P., Bdioui, N., Mary, B., and Rumpel, C. (2007). Stability of organic carbon in deep soil layers controlled by fresh carbon supply. Nature 450, 277–280. doi: 10.1038/nature06275
Han, Q. L., Zhou, X. B., Liu, L., Zhao, Y. D., and Zhao, Y. (2018). Three-dimensional visualization of soil pore structure using computed tomography. J. Forest. Res. doi: 10.1007/s11676-018-0834-z
Hljtr (1992). Soil of Heilongjiang Province, P.R. China. Harbin: Science and Technology press of Heilongjiang Province.
Hong, S., Piao, S., Chen, A., Liu, Y., Liu, L., Peng, S., et al. (2018). Afforestation neutralizes soil pH. Nat. Commun. 9:520. doi: 10.1038/s41467-018-02970-1
Hooker, T. D., and Compton, J. E. (2003). Forest ecosystem carbon and nitrogen accumulation during the first century afterag ricultural abandonment. Ecol. Appli. 13, 299–313. doi: 10.1890/1051-0761(2003)013[0299:FECANA]2.0.CO;2
Huo, Z. L., Dai, X. Q., Feng, S. Y., Kang, S. Z., and Huang, G. H. (2013). Effect of climate change on reference evapotranspiration and aridity index in arid region of China. J. Hydrol. 492, 24–34. doi: 10.1016/j.jhydrol.2013.04.011
Jaskulska, R., and Jaskulska, J. (2017). Efficiency of old and young shelterbelts in reducing the contents of nutrients in Luvisols. Agric. Ecosyst. Environ. 240, 269–275. doi: 10.1016/j.agee.2017.02.033
Jha, K. K. (2018). Biomass production and carbon balance in two hybrid poplar (Populus euramericana) plantations raised with and without agriculture in southern France. J. Forest. Res. 29, 1689–1701. doi: 10.1007/s11676-018-0590-0
Jia, X. X., Shao, M. A., Zhu, Y. J., and Luo, Y. (2017). Soil moisture decline due to afforestation across the Loess Plateau, China. J. Hydrol. 546, 113–122. doi: 10.1016/j.jhydrol.2017.01.011
Jiao, J. Y., Zhang, Z. G., Bai, W. J., Jia, Y. F., and Wang, N. (2012). Assessing the ecological success of restoration by afforestation on the Chinese Loess Plateau. Restor. Ecol. 20, 240–249. doi: 10.1111/j.1526-100X.2010.00756.x
Jin, T. T., Fu, B. J., Liu, G. H., and Wang, Z. (2011). Hydrologic feasibility of artificial forestation in the semi-arid Loess Plateau of China. Hydrol. Earth Syst. Sci. Discuss. 8, 2519–2530. doi: 10.5194/hess-15-2519-2011
Jobbágy, E. G., and Jackson, R. B. (2000). The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Applic. 10, 423–436. doi: 10.1890/1051-0761(2000)010[0423:TVDOSO]2.0.CO;2
Jobbágy, E. G., and Jackson, R. B. (2004). The uplift of soil nutrients by plants: biogeochemical consequences across scales. Ecol. Evol. 85, 2380–2389. doi: 10.1890/03-0245
Kettler, T. A., Doran, J. W., and Gilbert, T. L. (2001). Simplified method for soil particle-size determination to accompany soil-quality analyses. Soil Sci. Soc. Am. J. 65, 849–852. doi: 10.2136/sssaj2001.653849x
Lemma, B., Dan, B. K., Nilsson, I., and Olsson, M. (2006). Soil carbon sequestration under different exotic tree species in the southwestern highlands of Ethiopia. Geoderma 136, 886–898. doi: 10.1016/j.geoderma.2006.06.008
Li, J. Q., and Cui, G. (2000). On nature forest protection and degraded ecosystem restoration in Northwest China. J. Beijing Forest. Univ. 22, 1–7. doi: 10.13332/j.1000-1522.2000.04.001
Li, L. J., and Han, X. Z. (2016). Changes of soil properties and carbon fractions after long-term application of organic amendments in Mollisols. Catena 143, 140–144. doi: 10.1016/j.catena.2016.04.007
Li, L. J., Xia, Z. B., Ye, R., Doane, T. A., and Horwath, W. R. (2018). Soil microbial biomass size and soil carbon influence the priming effect from carbon inputs depending on nitrogen availability. Soil Biol. Biochem. 119, 41–49. doi: 10.1016/j.soilbio.2018.01.003
Li, X. J. (2000). The alkili-saline land and agricultural sustainable development of the Western songnen plain in China. Sci. Geogr. Sin. 20, 51–55.
Li, Y., Piao, S., Li, L. Z. X., Chen, A., Wang, X., Ciais, P., et al. (2018). Divergent hydrological response to large-scale afforestation and vegetation greening in China. Sci. Adv. 4:eaar4182. doi: 10.1126/sciadv.aar4182
Lian, T. X., Jian, J., Wang, G. H., Tang, C. X., Yu, Z. H., Li, Y., et al. (2017). The fate of soybean residue-carbon links to changes of bacterial community composition in Mollisols differing in soil organic carbon. Soil Biol. Biochem. 109, 50–58. doi: 10.1016/j.soilbio.2017.01.026
Liang, H. B., Xue, Y. Y., Li, Z. S., Shuai, W., Xing, W., Gao, G., et al. (2018). Soil moisture decline following the plantation of Robinia pseudoacacia forests: evidence from the Loess Plateau. Forest Ecol. Manage. 412, 62–69. doi: 10.1016/j.foreco.2018.01.041
Likens, G. E., Driscoll, C. T., Buso, D. C., Siccama, T. G., Johnson, C. E., Lovett, G. M., et al. (1994). The biogeochemistry of potassium at Hubbard Brook. Biogeochemistry 25, 61–125. doi: 10.1007/BF00000881
Liu, W. Z., Zhang, X. C., Dang, T. H., Ouyang, Z., Li, Z., Wang, J., et al. (2010). Soil water dynamics and deep soil recharge in a record wet year in the southern Loess Plateau of China. Agric. Water Manage. 97, 1133–1138. doi: 10.1016/j.agwat.2010.01.001
Lu, J. L., Shen, G., Wang, Q., Pei, Z. X., Ren, M. L., et al. (2017). Larch, ash, Scots pine, and farmland-induced differences on 17 soil parameters and their comprehensive analyses. Acta Ecol. Sin. 37, 1–10. doi: 10.1016/j.chnaes.2016.09.003
Lv, H. L., Wang, W., He, X., Wei, C., Xiao, L., and Zhang, B (2018). Association of urban forest landscape characteristics with biomass and soil carbon stocks in Harbin City, Northeastern China. PeerJ 6:e5825. doi: 10.7717/peerj.5825
Mao, R., Zeng, D. H., Hu, Y. L., Li, L. J., and Yang, D. (2010). Soil organic carbon and nitrogen stocks in an age-sequence of poplar stands planted on marginal agricultural land in Northeast China. Plant Soil 332, 227–287. doi: 10.1007/s11104-010-0292-7
Marie-Anne, D. G., Jastrow, J. D., Gilette, S., Johns, A., and Wullschleger, S. D. (2014). Differential priming of soil carbon driven by soil depth and root impacts on carbon lability. Soil Biol. Biochem. 69, 147–156. doi: 10.1016/j.soilbio.2013.10.047
Marschner, B., and Noble, A. D. (2000). Chemical and biological processes leading to the neutralisation of acidity in soil incubated with litter materials. Soil Biol. Biochem. 32, 805–813. doi: 10.1016/S0038-0717(99)00209-6
Marta, O., and Halina, S. (2008). The effect of afforestation with Scots pine (Pinus silvestris L.) of sandy post-arable soils on their selected properties. I. Physical and sorptive properties. Plant Soil 305, 157–169. doi: 10.1007/s11104-008-9537-0
Mendham, D. S., O'connell, A. M., Grove, T. S., and Rance, S. J. (2003). Residue management effects on soil carbon and nutrient contents and growth of second rotation eucalypts. Forest Ecol. Manage. 181, 357–372. doi: 10.1016/S0378-1127(03)00007-0
Merino, A. N., Solla-Gullón, F., and Edeso, J. M. (2004). Soil changes and tree growth in intensively managed Pinus radiata in northern Spain. Forest Ecol. Manage. 196, 393–404. doi: 10.1016/j.foreco.2004.04.002
Nan, H., Liu, X. X., Jumpponen, A., Setälä, H., Kotze, D. J., Biktasheva, L., et al. (2018). Over twenty years farmland reforestation decreases fungal diversity of soils, but stimulates the return of ectomycorrhizal fungal communities. Plant Soil 427, 231–244. doi: 10.1007/s11104-018-3647-0
Obschatko, E. S. D., Foti, M. P., Román, M. E., Cuenca Capa, P. R., Chinchilla, O., et al. (2010). Global Forest Resources Assessment 2010: Main Report. Fao Forestry Paper.
Paul, K. I., Polglase, P. J., Nyakuengama, J. G., and Khanna, P. K. (2002). Change in soil carbon following afforestation. Forest Ecol. Manage. 168, 241–257. doi: 10.1016/S0378-1127(01)00740-X
Ritter, E. (2007). Carbon, nitrogen and phosphorus in volcanic soils following afforestation with native birch (Betula pubescens) and introduced larch (Larix sibirica) in Iceland. Plant Soil 295, 239–251. doi: 10.1007/s11104-007-9279-4
Romanowicz, R. B., Driscoll, C. T., Fahey, T. J., Johnson, C. E., Likens, G. E., Siccama, T. G., et al. (1996). Changes in the biogeochemistry of potassium following a whole-tree harvest. Soil Sci. Soc. Am. J. 60, 1664–1674. doi: 10.2136/sssaj1996.03615995006000060009x
Rukshana, F., Butterly, C. R., Baldock, J. A., and Tang, C. (2011). Model organic compounds differ in their effects on pH changes of two soils differing in initial pH. Biol. Fertil. Soils 47, 51–62. doi: 10.1007/s00374-010-0498-0
Slessarev, E. W., Lin, Y., Bingham, N. L., Johnson, J. E., Dai, Y., Schimel, J. P., et al. (2016). Water balance creates a threshold in soil pH at the global scale. Nature 540, 567–569. doi: 10.1038/nature20139
Tang, C., and Yu, Q. (1999). Impact of chemical composition of legume residues and initial soil pH on pH change of a soil after residue incorporation. Plant Soil 215, 29–38. doi: 10.1023/A:1004704018912
Wang, C. M., Ouyang, H., Shao, B., Tian, Y. Q., Zhao, J. G., and Xu, H.-Y. (2006). Soil carbon changes following afforestation with Olga Bay Larch (Larix olgensis Henry) in Northeastern China. J. Integr. Plant Biol. 48, 503–512. doi: 10.1111/j.1744-7909.2006.00264.x
Wang, G. Y., Ma, O. Z., Wang, L. G., Shrestha, A., Chen, B. Z., Mi, F., et al. (2018a). Local perceptions of the conversion of cropland to forestland program in Jiangxi, Shaanxi, and Sichuan, China. J. Forest. Res. doi: 10.1007/s11676-018-0870-8
Wang, H. M., Wang, W., Chen, H., Zhang, Z., Mao, Z., and Zu, Y. G. (2014). Temporal changes of soil physic-chemical properties at different soil depths during larch afforestation by multivariate analysis of covariance. Ecol. Evolu. 4, 1039–1048. doi: 10.1002/ece3.947
Wang, H. M., Wang, W. J., and Chang, S. (2017a). Sampling Method and Tree-Age Affect Soil Organic C and N Contents in Larch Plantations. Forests 8, 28. doi: 10.3390/f8010028
Wang, J. G., Liu, H. X., and Meng, K. (1996). Songnen Plain Agricultural Ecosystem Research. Harbin: Harbin Engineering University Press.
Wang, J. K. (2002). An approach to the changes of black soil quality (I)—changes of the indices of black soil with the year(s) of reclamation. J. Shenyang Agric. Univ. 33, 43–47.
Wang, Q., Wang, W., He, X., Zhang, W., Song, K., Han, S., et al. (2015). Role and variation of the amount and composition of glomalin in soil properties in farmland and adjacent plantations with reference to a primary forest in north-eastern China. PLoS ONE 10:e0139623. doi: 10.1371/journal.pone.0139623
Wang, Q., Wang, W., He, X., Zheng, Q., Wang, H., Wu, Y., et al. (2017). Changes in soil properties, X-ray-mineral diffractions and infrared-functional groups in bulk soil and fractions following afforestation of farmland, Northeast China. Sci. Rep. 7:12829. doi: 10.1038/s41598-017-12809-2
Wang, Q., Wu, Y., Wang, W., Zhong, Z., Pei, Z., Ren, J., et al. (2014). Spatial variations in concentration, compositions of glomalin related soil protein in poplar plantations in northeastern China, and possible relations with soil physicochemical properties. Sci. World J. 2014:160403. doi: 10.1155/2014/160403
Wang, S., Wang, Z., and Gu, J. (2017). Variation patterns of fine root biomass, production and turnover in Chinese forests. J. Forest. Res. 28, 1185–1194. doi: 10.1007/s11676-017-0386-7
Wang, W., Xiao, L., Zhang, J., Yang, Y., Tian, P., Wang, H., et al. (2018). Potential of Internet street-view images for measuring tree sizes in roadside forests. Urban Forest. Urban Green. 35, 211–220. doi: 10.1016/j.ufug.2018.09.008
Wang, W., Zhang, B., Xiao, L., Zhou, W., Wang, H., and He, X. (2018b). Decoupling forest characteristics and background conditions to explain urban-rural variations of multiple microclimate regulation from urban trees. PeerJ 6:e5450. doi: 10.7717/peerj.5450
Wang, W. J., He, H. S., Zu, Y. G., Guan, Y., Liu, Z. G., Zhang, Z.-H., et al. (2011a). Addition of HPMA affects seed germination, plant growth and properties of heavy saline-alkali soil in northeastern China: comparison with other agents and determination of the mechanism. Plant Soil 339, 177–191. doi: 10.1007/s11104-010-0565-1
Wang, W. J., Lu, J. L., Du, H. J., Wei, C. H., Wang, H. M., Fu, Y., et al. (2017a). Ranking thirteen tree species based on their impact on soil physiochemical properties, soil fertility, and carbon sequestration in Northeastern China. Forest Ecol. Manage. 404, 214–229. doi: 10.1016/j.foreco.2017.08.047
Wang, W. J., Qiu, L., Zu, Y. G., Su, D. X., An, J., Hong-Yan, W., et al. (2011b). Changes in soil organic carbon, nitrogen, pH and bulk density with the development of larch (Larix gmelinii) plantations in China. Glob. Change Biol. 17, 2657–2676. doi: 10.1111/j.1365-2486.2011.02447.x
Wang, W. J., Wang, H. M., and Zu, Y. G. (2014). Temporal changes in SOM, N, P, K, and their stoichiometric ratios during reforestation in China and interactions with soil depths: importance of deep-layer soil and management implications. Forest Ecol. Manage. 325, 8–17. doi: 10.1016/j.foreco.2014.03.023
Wang, W. J., Zhong, Z. L., Wang, Q., Wang, H., Fu, Y., He, X., et al. (2017b). Glomalin contributed more to carbon, nutrients in deeper soils, and differently associated with climates and soil properties in vertical profiles. Sci. Rep. 7:13003. doi: 10.1038/s41598-017-12731-7
Wang, Z. Q., Liu, B. Y., Gang, L., and Zhang, Y. X. (2009a). Soil water depletion depth by planted vegetation on the Loess Plateau. Sci. China 52, 835–842. doi: 10.1007/s11430-009-0087-y
Wang, Z. Q., Liu, B. Y., and Wang, X. (2009b). Erosion effect on the productivity of black soil in Northeast China. Sci. China 52, 1005–1021. doi: 10.1007/s11430-009-0093-0
Wei, X., Qiu, L., Shao, M. G., Zhang, X., and Gale, W. (2012). The accumulation of organic carbon in mineral soils by afforestation of abandoned farmland. PLoS ONE 7:e32054. doi: 10.1371/journal.pone.0032054
Wei, X., Shao, M., Gale, W., and Li, L. (2014). Global pattern of soil carbon losses due to the conversion of forests to agricultural land. Sci. Rep. 4:4062. doi: 10.1038/srep04062
Wiesmeier, M., Prietzel, J., Barthold, F., Spörlein, P., Geuß, U., Hangen, E., et al. (2013). Storage and drivers of organic carbon in forest soils of southeast Germany (Bavaria) – Implications for carbon sequestration. Forest Ecol. Manage. 295, 162–172. doi: 10.1016/j.foreco.2013.01.025
Wu, Y., Wang, W., Wang, Q., Zhong, Z., Pei, Z., Wang, H., et al. (2018). Impact of poplar shelterbelt plantations on surface soil properties in northeast China. Can. J. Forest Res. 48, 559–567. doi: 10.1139/cjfr-2017-0294
Wu, Y., and Wang, W. J. (2016). “Poplar forests in NE China and possible influences on soil properties: ecological importance and sustainable development,” in Poplars and Willows, Cultivation, Applications and Environmental Benefits, ed M. V. Desmond (Hauppauge, NY: Novapublishers), 1–29.
Wuest, S. B. (2009). Correction of bulk density and sampling method biases using soil mass per unit area. Soil Sci. Soc. Am. J. 73, 3647–3654. doi: 10.2136/sssaj2008.0063
Xiang, S., Doylec, A., Holdenad, P. A., and Schimelc, J. P. (2008). Drying and rewetting effects on C and N mineralization and microbial activity in surface and subsurface California grassland soils. Soil Biol. Biochem. 40, 2281–2289. doi: 10.1016/j.soilbio.2008.05.004
Xu, J. M., Tang, C., and Chen, Z. L. (2006). The role of plant residues in pH change of acid soils differing in initial pH. Soil Biol. Biochem. 38, 709–719. doi: 10.1016/j.soilbio.2005.06.022
Yang, L., Wei, W., Chen, L. D., and Mo, B. R. (2012). Response of deep soil moisture to land use and afforestation in the semi-arid Loess Plateau, China. J. Hydrol. 475, 111–122. doi: 10.1016/j.jhydrol.2012.09.041
Yang, Y. B., Lv, H. L., Fu, Y. J., He, X. Y., and Wang, W. J. (2019). Associations between road density, urban forest landscapes, and structural-taxonomic attributes in northeastern china: decoupling and implications. Forests 10:58. doi: 10.3390/f10010058
Yao, Y., Wang, X., Zeng, Z., Liu, Y., and Peng, S (2016). The effect of afforestation on soil moisture content in Northeastern China. PLoS ONE 11:e0160776. doi: 10.1371/journal.pone.0160776
You, M. Y., Li, N., Zou, W. X., Han, X., and Burger, M. (2017). Increase in soil organic carbon in a Mollisol following simulated initial development from parent material: dynamics of SOC stocks at initial stage of soil restoration. Eur. J. Soil Sci. 68, 39–47. doi: 10.1111/ejss.12400
Zhang, G. Q., Zhang, P., and Cao, Y. (2018). Ecosystem carbon and nitrogen storage following farmland afforestation with black locust (Robinia pseudoacacia) on the Loess Plateau, China. J. Forest. Res. 29, 761–771. doi: 10.1007/s11676-017-0479-3
Zhang, H. Q., Wang, L. X., Sun, G. Y., and Yang, Y. (2013). Evaluation of saling-alkali land resource and development potential in low songnen plains. Chin. J. Agric. Resour. Region. Plan. 34, 6–11.
Zhang, X. Q., Kirschbaum, M. U. F., Hou, Z., and Guo, Z. (2004). Carbon stock changes in successive rotations of Chinese fir (Cunninghamia lanceolata (lamb) hook) plantations. Forest Ecol. Manage. 202, 131–147. doi: 10.1016/j.foreco.2004.07.032
Zhang, Y., Peng, C. H., Li, W. Z., Tian, L. X., Zhu, Q. A., Chen, H., et al. (2016). Multiple afforestation programs accelerate the greenness in the ‘Three North’ region of China from 1982 to 2013. Ecol. Indicat. 61, 404–412. doi: 10.1016/j.ecolind.2015.09.041
Zhao, Q., Zeng, D. H., Lee, D. K., He, X. Y., Fan, Z. P., and Jin, Y. H. (2007). Effects of Pinus sylvestris var. mongolica afforestation on soil phosphorus status of the Keerqin Sandy Lands in China. J. Arid Environ. 69, 569–582. doi: 10.1016/j.jaridenv.2006.11.004
Zhong, Z. L., Wang, W. J., Wang, Q., Wu, Y., Wang, H. M., and Pei, Z. (2017). Glomalin amount and compositional variation, and their associations with soil properties in farmland, northeastern China. J. Plant Nutri. Soil Sci. 180, 1–13. doi: 10.1002/jpln.201600579
Keywords: poplar shelterbelt, farmlands, soil properties change, deep-layer soil, analysis of causes
Citation: Wu Y, Wang Q, Wang H, Wang W and Han S (2019) Shelterbelt Poplar Forests Induced Soil Changes in Deep Soil Profiles and Climates Contributed Their Inter-site Variations in Dryland Regions, Northeastern China. Front. Plant Sci. 10:220. doi: 10.3389/fpls.2019.00220
Received: 28 September 2018; Accepted: 08 February 2019;
Published: 05 March 2019.
Edited by:
Zhiyou Yuan, College of Forestry, Northwest A&F University, ChinaReviewed by:
Lu-Jun Li, Northeast Institute of Geography and Agroecology (CAS), ChinaGuofan Shao, Purdue University, United States
Copyright © 2019 Wu, Wang, Wang, Wang and Han. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wenjie Wang, wwj225@nefu.edu.cn
Shijie Han, hansj@iae.ac.cn