• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Soil-plant co-stimulation during forest vegetation restoration in a subtropical area of southern China

    2020-10-20 08:13:38ChanChenXiFangWenhuaXiangPifengLeiShuaiOuyangandYakovKuzyakov
    Forest Ecosystems 2020年3期

    Chan Chen,Xi Fang,2*,Wenhua Xiang,2,Pifeng Lei,2,Shuai Ouyang,2 and Yakov Kuzyakov,3,4

    Abstract

    Keywords: Vegetation restoration, Soil physicochemical properties, Soil organic carbon, Vegetation features, Driving factors

    Background

    Forest vegetation restoration has become a priority study area in efforts to solve global environmental problems,as highlighted by the Bonn Challenge, a global effort to restore 150 million hectares of degraded land and deforested forests by 2020 (Crouzeilles et al. 2016). Establishing the mechanisms of plant communities in the process of recovery has concentrated mainly on species composition,and their quantitative characteristics and spatial distribution. While these factors are relatively clear (Xiang et al. 2013; Chen et al. 2019), there is still a lack of indepth research on the feedback relationships between plant and soil,and the succession processes and regulation mechanisms of plant communities (Hu et al. 2017; Wang et al. 2018a). The feedback relationship between vegetation and soil has a great impact on the plant community,soil nutrient cycling,and soil and water conservation during vegetation restoration (Demenois et al. 2018).Insights into vegetation-soil feedback relationships are instrumental in predicting future scenarios under varying environmental conditions (van der Putten et al. 2013), as well as in designing measures for vegetation restoration at different succession stages(Huang et al.2018).

    The interactive effects of soil and vegetation suggest that both are always co-evolving and developing, which are recognized as an important mechanism for forest succession and development (van der Putten et al. 2013).The association between soil and aboveground vegetation may shift over the course of restoration (Huang et al. 2015). In the early stage of vegetation restoration,soil resources are the main limiting factors(van Der Maarel and Franklin 2013). Research has shown that the enrichment, spatial distribution, and redistribution of soil nutrients significantly affect the growth, reproduction, distribution,succession,and net primary productivity of plants(Alday et al. 2012). In particular, soil nutrients and water are the key factors in regulating vegetation development,as confirmed by the results of some fertilization experiments(Chang and Turner 2019) and different forest succession series (Huang et al. 2017). In turn, vegetation development can drive changes in the development and maintenance of soil(Huang et al.2018).Especially in the late stage of vegetation restoration, the accumulation of plant biomass leads to an increase in the return of soil organic carbon (SOC)and nutrients (Gu et al. 2019). Furthermore, soil nutrient storage reflects the balance of the main ecological processes, including nutrients stored in aboveground biomass,nutrients decomposed and returned to soil, and nutrient leaching, these mixed results may cause the complexity of the interaction between soil and vegetation (Huang et al.2018). Therefore, knowledge of how soil, vegetation and their interaction act on vegetation restoration is of particular importance for predicting future ecological restoration and development.

    Subtropical forest covers an extensive area and supports a high level of biodiversity and a global carbon store, particularly in China which has 71%of the current total forest area in the subtropics according to the MODIS landcover layer for 2012, with abundant rainfall and abundant forest resources(Corlett and Hughes 2015).However,long-term severe human disturbance has a serious effect on subtropical forest ecosystems, with complex topography and climate change resulting in fewer climax forests and a decrease in the functioning of an ecological security barrier(Huang et al.2018).The Chinese government initiated a series of state-funded forestry ecological projects,including programs to protect natural forests, the Grain to Green program,and the construction of shelterbelts in the middle and upper reaches of the Yangtze River. Consequently,forest vegetation has been rapidly restored,forming a series of secondary vegetation communities at different stages of restoration in this area (Ouyang et al.2016).During vegetation restoration,aboveground vegetation and soil physicochemical properties gradually change(Zhang et al. 2019). Changes in plant development and soil variables during vegetation restoration have been demonstrated in several studies (Ayma-Romay and Bown 2019; Wang et al. 2018a; Zhang et al. 2019), but the restorative effect of soil or vegetation has rarely been explored, and there is little information on how soil physicochemical properties and vegetation act together to affect vegetation restoration(Chang and Turner 2019).To our knowledge, no studies have addressed the question of the relative importance of the effects of soil, vegetation and their synergism on promoting vegetation restoration.It has therefore become a burning issue to elucidate the coordinated control effect of vegetation restoration, soil,and water on vegetation ecology and restoration ecology(Chang and Turner 2019).

    In this study,we followed the succession process of subtropical forest communities,and selected four distinct restoration periods (i.e. 4-5, 10-12, 45-46 years and old growth forest), which represent the four main stages of vegetation restoration in the subtropics of China. We selected permanent plots and determined soil physicochemical properties and vegetation features; i.e. species diversity, biomass, height (H), and diameter at breast height (DBH). Our objective was to investigate how soil physicochemical properties and vegetation features change and how soil and vegetation stimulate vegetation restoration individually and collectively.We formulated two hypotheses: (1) that vegetation restoration would have an obvious positive effect on soil physicochemical properties and vegetation features; and (2) that soil properties and vegetation features would collectively promote vegetation restoration, especially would have a significant impact on biomass. In addition, the main impact factors of biomass would be different in different restoration periods.

    Methods

    Study site

    As shown in Fig. 1, the study site was located in Changsha County (28°23′-28°24′ N, 113°17′-113°27′ E), situated in the middle of Hunan Province, China. The topography features a typical low hilly landscape, at an altitude of 55-260 m above sea level with an average slope of 18°-25°. The climate is characterized by southeast monsoon and a mid-subtropical humid climate with an annual average precipitation of 1416.4 mm (primarily between April and August) and an annual mean air temperature of 17.3°C. minimum and maximum air temperatures are 10.3°C in January and 39.8°C in July and August, respectively. The soils are mainly composed of red earth, which developed from slate and shale and are categorized as Alliti-Udic Ferrosols in the Chinese Soil Taxonomy, corresponding to Acrisol in the World Reference Base for Soil Resources (IUSS Working Group WRB 2006). Evergreen broadleaved forests are the climax and primary vegetation, but have been disturbed in varying degrees by human activities such as firewood collection. Natural forest protection programs in the past two decades have resulted in a variety of vegetation communities at different restoration stages in this area.

    Vegetation sampling

    In October 2015, four adjacent vegetation communities,with basically similar environmental conditions (site,slope, soil and climate) as showed in Table 1 were selected to represent a vegetation restoration gradient(using the method of space-for-time substitution). These communities were:

    (1) 4–5 yrs. restoration period. Controlled burns and site preparation were carried out in native evergreen broadleaved forest in the winter of 1965.A Pinus massoniana plantation was established in 1966 without any fertilization during this operation and then clear-cut in 1990. The woodlands were repeatedly cut until 2012.Since that time the vegetation has naturally recovered. The community is dominated by well-grown herbs, presently accompanied by some young shrubs,and belongs to the early stage of restoration according to the succession process of subtropical evergreen broadleaved forest (Xiang et al.2016).

    Table 1 Stand characteristics of the four forest types

    (2) 10–12 yrs. restoration period.Native evergreen broadleaved forest underwent a prescribed burn in 1965 and deforested to establish a Cunninghamia lanceolata plantation in 1966. This C. lanceolata plantation was clear-cut in 1989. The woodlands were logged every 3 to 5 years until 2004. The vegetation has naturally recovered to form a shrub community with well-grown shrubs and belongs to midrestoration stage according to the succession process of subtropical evergreen broadleaved forest(Xiang et al. 2016). However, the shrub community has no obvious arbor layers and herbaceous plant is relatively infrequent.

    (3) 45–46 yrs. restoration period.This period represents the secondary stage of mid-restoration. Native evergreen broadleaved forest was deforested in the early 1970s, and then naturally recovered to coniferous and broadleaved mixed forest. The communities are now about 45-50 years old, and have abundant seedlings and saplings, with larger plant density.However, the proportion of large diameter individuals is relatively low.

    (4) Old growth forest (representing the late stage of restoration).Native evergreen broadleaved forest has been well protected against human disturbances. According to the survey with local residents, this forest has been more than 90 years.

    In October 2015, we randomly established 4 fixed sample plots for long-term observation in each restoration period (Fig. 1). In the 4-5 and 10-12 years restoration periods, the plots were set at 20 m×20 m. In the 45-46 years restoration periods and old growth forest, the plots were established at 30 m×30 m. The 4 fixed plots in each restoration period were set in different mountains as far as possible,and the space distance between the two plots was more than 1000 m. To investigate the floristic components and tree spatial patterns of the forests, each plot (20 m×20 m) in the 4-5 and 10-12 years restoration periods was subdivided into four subplots (10 m×10 m), and each plot (30 m×30 m) in the 45-46 years restoration periods and old growth forest was subdivided into nine subplots (10 m×10 m).

    Species diversity measurement

    Species identities were recorded and measurements were taken of total H, the lowest live branch and crown width, and DBH for all individuals with DBH >1 cm in each plot. The data were used to calculate vegetation structural parameters of the different restoration periods; i.e. density of main tree species, average DBH, and average H. The Shannon-Wiener index (SD) was used to quantify the diversity of woody plants species in each plant community with the equation below (Madonsela et al. 2018).

    In Eq. 1, n represents the total number of species in the community, and Pirepresents the relative frequency of species i in the community. Table 1 summarizes the characteristics and site factors of each community.

    Biomass measurement

    Based on community surveys, biomass was measured by the harvest method and calculated by establishing relative growth equations of the organic biomass of the main tree species. For the 4-5 years restoration period, we collected all vegetation (shrubs,vines, herbs) in 2 m×2 m quadrats which were on plot peripheries and then classified the same plants according to the following criteria: shrubs were composed of fruit, leaf, branch, stem, and root;vines were composed of fruit, leaf, stem, and root;and herbs were composed of aboveground and underground parts. A 1 m×1 m quadrat was set up at the center of each 10 m×10 m subplot to determine litter biomass. All litter was collected from the ground in these quadrats and transported to the laboratory. Determined samples were freshly weighed and then oven-dried at 85°C to a constant weight to measure their dry mass for estimating biomass per plot area.

    For the 10-12 years restoration period, according to the average DBH and average H of the shrub (>1.5 m), and with the aim of ensuring that at least 3 average sample trees per dominant tree species were collected, 3 sample trees were selected and collected for each dominant tree in each plot periphery to determine fresh weight. Shrub samples were composed of fruit, leaf, branch, stem, and root. After oven-drying at 85°C to a constant weight, we determined moisture content and calculated each biomass component of each tree species, establishing their relative growth equations to calculate biomass per shrub plant(Table 2). The biomass determination of shrubs(below 1.5 m), vines, herbaceous layers, and litter layer used the same method as the 4-5 years restoration period. Finally, estimated biomass per plot area was based on data from community surveys.

    For the 45-46 years restoration period, 3 sample trees were selected for each dominant tree in each plot periphery according to average DBH and average H, with the same aim as that for the 10-12 years restoration period. Stratified samples (1.3 m, 3.6 m) were collected for the aboveground part and complete samples were excavated for the underground part (within 1.5 m of the tree stump) to measure fresh weight.Tree samples were composed of leaf, branch, stem,and root, in which root included fine root (<0.2 cm),rootlet (0.2-0.5 cm), thick root (0.5-2.0 cm), large root (>2.0 cm), and root apex. After determining fresh weight, samples were oven-dried at 85°C to a constant weight to calculate moisture content. We then estimated each biomass component of each tree species, established their relative growth equations and then calculated the biomass per tree plant (Table 2). The same methods as above were used to determine the biomass of shrubs, vines, herbaceous layers,and litter layer. Estimated biomass per plot area was based on data from community surveys. For the old growth forest, the relative growth equations for the main tree species in the tree layer were established using a similar method as the 45-46 years restoration period. However, because of the ban on logging in the old growth forest, the general growth equations of Cyclobalanopsis, deciduous broadleaf, evergreen broadleaf, and C. lanceolata, which were established by Ouyang et al. (2016) and Liu et al. (2010), were also used to estimate the biomass in the tree layer(Table 2).

    Soil sampling and analysis

    Each permanent plot was divided into 3 equal grids of cells along the diagonal for soil sampling. In each cell, soil profile characteristics were surveyed in 2015 to illustrate the consistency and comparability of soil background in different vegetation restoration periods, as shown in Table 3. Soil samples were taken by using cylindrical cores with a volume of 200 cm3collected at depths of 0-10, 10-20, 20-30, and 30-40 cm in December 2015, and in April, June, and October 2016.

    Soil samples from three cells at the same depth within a plot were mixed into a composite sample.Plant roots, debris, and gravels were cleared. Soil samples were air-dried and sieved through a 2-mm mesh for soil pH, available phosphorus (AP), and available potassium (AK); through a 1-mm mesh for soil alkaline hydrolysis nitrogen (AN); and through a 0.25-mm mesh for soil SOC, TN, total phosphorus(TP), total potassium (TK), total calcium (Ca), and total magnesium (Mg) determinations. The following properties were determined in the soil samples:

    (1) Bulk density (BD) was calculated using weights of the dried soil sample from the known cylindrical core volume. (2) pH value was analyzed in a soil-towater (deionized) ratio of 1:2.5 using a pH meter(FE20, Mettler Toledo, Switzerland). (3) SOC content was determined by the K2Cr2O7-H2SO4oxidation method. (4) TN content was determined using a semi-micro Kjeldahl method (Bremner 1996). (5) TP,TK, Ca, and Mg were extracted via aqua regia and 1:1 HCl. After extraction, TP was determined by spectrophotometry and TK, Ca, and Mg by atomic emission spectrometry with inductively coupled plasma(ICP-OES) using a Perkin Elmer Optima 7300DV optical emission spectrometer (Nicia et al. 2018). (6) For

    AN and AK, we used the alkaline diffusion method and the ammonium acetate extraction flame spectrophotometer method (ISSCAS 1978). (7) For AP, we used the Olsen method (Olsen et al. 1983).

    Table 2 Relative growth equations of different biomass components of the main tree species

    Table 3 Structural characteristics of soil profile

    Statistical analysis

    For data processing, we used the Microsoft Excel package(Office 2010).All statistical analyses were conducted using the R statistical software package (R Development Core Team 2016). In order to reflect the annual average situation of the soil properties, the arithmetic mean of four seasons in the same soil layer of each plot was calculated.At the same time, taking into account the great changes between soil layers of each variable, a weighted average of four soil layers was carried out. The parameter content of a soil layer as a percentage of the sum of four soil layers(fi) was calculated using Eq. 2, and the weighted average(X0)was calculated using Eq.3.

    In Eqs. 2 and 3, n represents the number of soil layers;Lirepresents the parameter content of a soil layer; and Xirepresents the parameter content of a soil layer.

    The response rate was used to determine the effects of restoration periods on soil properties and vegetation features, calculated by Eq. 4.

    In Eq.4,X1represents one of the soil properties or vegetation features in the 4-5, 10-12, and 45-46 years restoration periods, and X2represents one of that in the old growth forest. In this study, only X1in the 4-5 years restoration period is selected to reflect the extent of variables variation over the whole vegetation restoration.A positive value indicates an increase,a negative value indicates a decrease,and greater absolute values indicate greater change.Figure 2 was drawn by the geom_histogram function of ggplot2 package in the R statistical software. Before drawing, the values were normalized to a proportion of maximum value (= 1) and by min-max normalization to keep a common scale ranging from 0 to 1(Jain et al.2005).The min-max normalization was calculated by Eq.5.

    Principal component analysis(PCA)was used to determine the main factors in soil properties and vegetation features influencing vegetation restoration, and the correlations between soil properties and vegetation features.The PCA was implemented using the prcomp function and drawn by the ggplot2 package of the R software.The selection criteria for principal components included a cumulative contribution rate over 85% and eigenvalues greater than 1. Indicator whose absolute value of a loading matrix was greater than 0.7 was selected as the dominant factors (Armstrong 1967) for vegetation restoration. The cosine values of the angles between variables indicate relationship strength; angles ranging from 0° to 90° indicate variables have positive correlations,and 90° to 180° indicate negative correlations.

    Based on the results of PCA, we used linear function analysis to further examine the significant correlations of soil properties and vegetation features. Before fitting the linear function, data were normalized by min-max normalization for unifying dimensions, and also calculated using Eq. 5. It was assumed that the relation between soil properties and vegetation features can be expressed by Eq. 6, where k represents slope, and b represents a constant.

    Figure 3 was produced via the lm function and plot function in R. Variation partitioning analysis (VPA) was performed to quantify the relative contributions of soil factors, vegetation factors and their joint action to changes in biomass by the varpart function of vegan package. Before VPA, the suitably independent variables with the variance inflation factor (VIF)<3 were selected by using the car package (Yang et al. 2017), and then the factor analysis (FA) were used to reduced soil factors and vegetation factors into a common factor respectively using psych package in R. Figure 4 was drawn by the geom_bar function of ggplot2 package in R. Following the results of PCA and linear function analysis, we conducted a stepwise regression analysis (SRA) to analyze,screen, and eliminate variables that cause multicollinearity, and to determine the leading impact factors of biomass per restoration period. The SRA was performed by the step function in R.

    Results

    Changes in vegetation features and soil physicochemical properties during vegetation restoration

    Vegetation features and soil physicochemical properties varied in the regularity of change according to the different restoration periods (Fig. 2). Vegetation features(species diversity, biomass, DBH, and H) increased remarkably with vegetation restoration, and the response rates increased by 86.36%, 2906.52%, 128.11%, and 596.97% respectively. Specifically, the highest values of species diversity and biomass were observed in the old growth forest, and the highest values of DBH and H were observed in the 45-46 years restoration period.The change trends of biomass, DBH, and H were basically the same (Fig. 2a). The maximum values of soil BD,pH value, Mg content, and AK content occurred in the 10-12 years restoration period. BD, pH, and AK content showed a decreasing trend whereas Mg showed an increasing trend with vegetation restoration (Fig. 2b). The response rates of BD, pH, and AK were negative but changed slightly. However, the contents of SOC, TN,TP, TK, Ca, AN, and AP increased with vegetation restoration, and their maximum values were recorded in the old growth forest except for TK (Fig. 2c and d). The response rates of SOC, TN, TP, TK, Ca, AN, and AP ranged from 10.63% to 544.35%, with AN having the highest response rate of 544.65%, followed by TN(222.48%) and SOC (140.76%).

    Factors of soil properties and vegetation features influencing vegetation restoration and their relationships

    The results of PCA showed that soil properties and vegetation features explained 81.54% of the variations(PC1=48.70%; PC2=20.21%; PC3=12.63%), revealing three main correlated variable groups of vegetation restoration (Fig. 5). There was a strong positive correlation between PC1 and SOC, TN, AN, AP, biomass, DBH,and H, and a negative correlation between PC1 and soil pH. As shown in Fig. 5, the successful discrimination of the 45-46 years restoration periods and old growth forest from other periods were strongly influenced by PC1.In the selection criteria, PC2 was correlated positively with Mg, and PC3 with TK. Figure 5 also shows that the successful discrimination of the 10-12 years restoration period from the 4-5 years restoration period was highly influenced by PC2 and PC3. Therefore, the key factors influencing vegetation restoration can be summarized as soil water and fertilizer conservation capacity (pH), organic matter (SOC), macro nutrients (TN, TK), medium nutrients (Mg), available nutrients (AN, AP), and the plant community growth situation (biomass, DBH, H).

    The results of PCA also showed that biomass,DBH,and H had significant correlations with each other, while species diversity was weakly correlated with them. Biomass,DBH, and H had similar relationships with soil factors(Fig.5).Specifically,the order of high positive correlations with soil factors was SOC >TN >AN >TP >AP, whereas a high negative correlation was with soil pH.The order of factors with high positive correlations with species diversity was Ca >AP >AN >TN(Fig.5).

    As shown in Fig.3,the results of linear function analysis revealed that as SOC,TN,AN,and AP increased,biomass,DBH, and H significantly increased (p <0.05). However,biomass decreased remarkably with the increased in pH(p <0.001). With the increasing of Ca and AP, species diversity showed a great increase trend(p <0.01).

    Effects of soil properties and vegetation features (DBH and H) on biomass variation

    The VPA results showed that the combination of soil properties and vegetation features explained 90.51% of the variation in biomass in the whole restoration process, and explained 83.44%, 99.99%, 99.99%, and 98.15% of the variation in 4-5, 10-12, 45-46 years vegetation periods and old growth forest,respectively(Fig.4).The interaction of soil properties and vegetation features all had the highest explanation for the variation in biomass, ranging from 55.55%-72.32%. The soil properties alone explained 3.30%-31.44% of the variation, and the vegetation features alone explained 5.09%-24.32%,among which soil properties had higher individual explanation than vegetation features except the 45-46 years restoration period.

    The results of SRA (Table 4) indicated that the factors influencing biomass in the whole restoration process included DBH and SOC. The fitting equation was: ybio-mass=7151.27xDBH+7595.62xSOC(R2=0.914, p=0.000).However, there were different factors influencing the biomass per restoration periods. In the 4-5 years restoration period, SOC was the only dominant factor, and the fitting equation was: ybiomass=-966.94xSOC(R2=0.903,p=0.050).H,pH,and AP were the main influential factors in the 10-12 years restoration period.The fitting equation was: ybiomass=15,620.74xH-1.00xpH-3484.06xAP(R2=0.990, p=0.000). H and pH were the main factors in the 45-46 years restoration period (ybiomass=-10,432.46xH+14,071.07xpH;R2=0.990,p=0.000).In the old growth forest, SOC, TN and AP were the impacting factors (ybio-mass=45,060.13xSOC+18,771.33xTN+26,287.80xAP; R2=0.990,p=0.000).In all periods,AN was not screened into the regression equation.

    Discussion

    Soil physicochemical properties during vegetation restoration

    Our results showed that soil BD decreased, and the contents of SOC, TN, TP, TK, Ca, AN, and AP increased with vegetation restoration (Fig. 2), indicating that soil physicochemical properties improved significantly.These results are partially consistent with our hypothesis and with the results of Zhang et al. (2019).

    The rapid recovery of SOC at our study site has been proven to be affected by plant biomass and soil nutrients(Gu et al. 2019). The response rate of SOC (140.76%) in this research was higher than the results under semi-arid conditions (71%) recorded by Boix-Fayos et al. (2009),which may be due to the more humid conditions in subtropical regions. Consistent with the rapid accumulation of SOC, the rates of change in TN and AN were greater than the SOC change. This result differs from the results of studies in the same subtropical area of southwest China (Xu et al. 2018), which may be due to differences in the degree of degradation and type of vegetation system. Additionally, soil N is also input from other N sources, such as atmospheric N deposition, and symbiotic N fixation by legumes (Alday et al. 2012). This explains why the recovery rates of TN and AN were greater than SOC. Our results for the increase in TP and AP contents are consistent with the results of Zhang et al. (2019), who proposed that soil TP and AP contents gradually increase with the composition of tree species,annual litter yield, and SOC content along with the development of a forest’s second succession. This is also supported by the significant positive correlation of soil TP and AP contents with the species diversity, biomass,and contents of SOC, TN and AN observed in this study(Fig. 5), which suggests that the accumulation of SOC improves soil nutrients during vegetation restoration(Zhang et al. 2019).

    The variation ranges of BD and pH in the subtropical regions of China are 0.97-1.47 g·cm-3and 4.5-6.0, respectively (Hunan Provincial Department of Agriculture1989). Our results were in the variation ranges. All the pH samples in our study indicate that soil pH (4.54-4.96) was lower than the results of Takoutsing et al.(2016), being formed by a moderate ferrallitic effect under high temperature and high humidity conditions in subtropical regions(Li et al.2012b).Meanwhile,decreasing soil BD and pH has also been attributed to the accumulation of organic matter, which is conducive to the formation of soil aggregates and the improvement of soil microbial activity (Bienes et al. 2016). This in turn releases a large number of small molecular organic acids during the decomposition of organic matter (van Breemen et al. 1984), resulting in a decline in soil BD and pH values. The SOC in our study increased and showed negative relationships with BD and pH during vegetation restoration (Fig. 5). In addition, this study shows that biomass stimulated the decrease in soil pH during vegetation restoration (Fig. 3). The accumulation of biomass led to increased biomass in the roots, almost certainly reflecting the development of the vegetation community from annual plant species to perennial plants, which is more conducive to the release and accumulation of the various acid exudates (Pang et al. 2018). Although BD and pH showed a general declining trend, their values reached a peak in the 10-12 years restoration period(Fig. 2). These results may be caused by a combination of factors (i.e. soil texture, vegetation types and soil acid-base equilibrium). Firstly, as herbs developed into shrubs in our study site, the erosion effect of rainwater on soil silt and clay particles resulted in a high proportion of sand particles in the 10-12 years restoration period, reflecting the transformation of soil texture to sandy soil with high BD (Wang et al. 2018b). Secondly,changes in vegetation types could be a major driver behind the difference in cations absorption of the vegetation and consequent variation in the proportions of soil cations (Gu et al. 2019). From soil acid-base equilibrium, the increase in cations (especially Mg and AK)suggests that the soil H+was replaced by increased alkaline ions (Berthrong et al. 2009). Due to the similar soil parent materials at different restoration stages, soil Ca,Mg, and K contents, which are all derived from parent rock materials, change little in response rates during vegetation restoration(Takoutsing et al. 2016).

    Table 4 Stepwise regression of corresponding factors for biomass(n=16)

    Vegetation development during restoration periods

    In our study, species diversity increased with an 86.36%recovery rate as restoration progressed, and these results are consistent with the results of Wang et al. (2018a).The amount of biomass increased significantly with the greatest recovery rate (2906.52%) over the old growth forest, followed by H (596%) and DBH (128%) in the 45-46 years period. These results are partially consistent with our hypothesis and are very similar to the results of Hu et al. (2017).

    Improvements to the soil environment can provide community habitat quality which then promotes the enrichment of community diversity (Huang et al. 2015). Ca content had a significant positive effect on species diversity (Fig. 3), reflected in the following mechanisms.Firstly, Ca2+has the function of maintaining the homeostasis of intracellular ions, especially in acidic soil where higher Ca2+content can counterbalance the toxicity of aluminum ions for plants, further improving plant resistance to adversity and being conducive to the improvement of community diversity (Roem et al. 2002).Secondly, the increase in soil Ca content alongside vegetation restoration can be instrumental in the coexistence of species with different Ca requirements and the settlement of calcium-loving species (Hooper 1998). Additionally, soil P determines the species composition of a vegetation community (Huang et al. 2015); thus, soil AP content was considered as another major factor determining species diversity increase (Fig. 3).

    In our study, biomass, DBH, and H had basically the same changing trend, and were all significantly affected by SOC, TN, AN, and AP contents (Figs. 4 and 5). This is consistent with Brandies et al. (2009), who demonstrated that there are significant positive growth rates and similar effect factors between biomass, DBH, and H in a general case. Data analysis of our study site confirmed that the percentages of individual trees with DBH greater than 8 cm and H greater than 5 m were larger in the 45-46 years restoration period (54% and 77% respectively) than in the old growth forest (41% and 63%respectively) (Chen et al. 2019). The greatest values of DBH and H in the 45-46 years restoration period may be because Pinus massoniana, as the dominant species,is a fast-growing heliophilous plant that gets more light by increasing vertical growth (H)(Cheng et al. 2011).

    Soil SOC, TN, AN, and AP contents were leading factors in stimulating the increase in biomass, DBH, and H(Fig. 3). As the environmental basis for vegetation survival, improving the soil provides a better habitat and essential nutrients for vegetation growth (Huang et al.2018), ultimately promoting the positive succession of vegetation (Liang et al. 2010). The accumulation of SOC affects biomass, DBH, and H mainly by decomposing and releasing large amounts of nutrients to meet plant growth needs, and by improving soil texture and promoting microbial activity which provide a better growing environment for vegetation (Alday et al. 2012). Moreover, the increase in N content promotes growth of the leaf area and improves plant photosynthesis, providing sufficient energy for the growth of individual plants. P is the nutrient that most limits productivity and species richness (Huang et al. 2015), and also controls leaf litter decomposition (Zeng et al. 2016). In addition, soil P changes the structure of the root system, promotes the formation and growth of fine roots, lateral roots and secretions of root exudates,and thereby stimulates plants to make more efficient use of soil nutrients(Li et al.2017).

    Key factors affecting vegetation restoration

    Soil factors (pH, SOC, TN, TK, Mg, AN, and AP) and vegetation features (biomass,DBH, and H) were the main factors influencing vegetation restoration at our study site.This is consistent with the finding that the recovery of degraded ecosystems not only relies on soil rehabilitation,but also on the reconstruction, productivity, and function of vegetation(Liang et al.2010;Peng et al.2012).

    The soil properties and vegetation features can be viewed as three distinct groups. The first group can be summarized as soil pH, SOC, TN, AN, AP, biomass,DBH, and H across the vegetation restoration periods.The roles of soil pH, SOC, TN, AN, and AP have been analyzed above. Specifically, soil resource is the main limiting factor in the early period of vegetation restoration. However, in the later period of vegetation restoration, the change in community characteristics leads to light conditions becoming a limiting factor (van Der Maarel and Franklin 2013). With the accumulation of biomass, a complex community structure reduces the understory light transmittance, controlling the vegetation in the understory including the growth and mortality of tree seedlings and saplings (Montgomery and Chazdon 2001). Therefore, the shade tolerant species are successively established,increasing understory vegetation richness. On the other hand, heliophilous species are shaped by increasing H and diameter to gain more light by adapting to strong interspecific competition. At increasingly larger H and DBH, light transmittance could further influence a species’ light-capturing ability and distribution (Cheng et al. 2011). The limitation of light conditions for vegetation growth and performance in the late vegetation restoration period means that increases in biomass, DBH, and H are the key growth factors, determining restoration success.

    The second group of soil variables includes Mg. The increase in Mg during the restoration periods was accompanied by a series of improvements in the plants’physiological processes, such as photosynthetic efficiency, carbohydrate metabolism, and synergistic absorption with P (Unger 2010). The third group showed that the vegetation restoration development was determined by TK. Besides N and P, K is the limiting nutrient with a significant influence on vegetation growth and development (Pang et al. 2018), mainly reflected in the impact on plant photosynthesis and respiration by controlling the regulation of stomata opening (Unger 2010), even though here the effect of TK on vegetation development was not significant.

    Previous studies have suggested that species diversity was the dominant vegetation factor for vegetation restoration in a large scale (Crouzeilles et al. 2016), because higher species richness can enhance ecosystem stability and increase nutrient use efficiency (Hu et al. 2017).However, species diversity was not considered to be an influential factor for vegetation restoration in our study.The difference could be due to the non-significance of the relationships between species diversity and the main soil physicochemical properties or biomass, indicating that species diversity had no significant effect on the recovery of soil fertility and plant communities at our study site.In addition, species diversity showed a decreasing trend in the 45-46 years restoration period (Fig. 2), in which dominant species transformed from simple shrubs and herbs to pioneer species such as Pinus massoniana.In fact, needles of some Pinus species have been reported as a hindering factor which influences the regeneration of native plants and increases in species diversity (Navarro-Cano et al.2010).It is reasonable that species diversity has no significant effect on vegetation restoration in specific study area,but further research is needed.

    Soil and vegetation factors affecting biomass

    The variation of biomass was one of the important indexes reflecting vegetation restoration (Mansourian et al.2005). Therefore, the relative importance of soil properties and vegetation features in driving biomass development can reflect the degree of their individual and joint influence on vegetation restoration.

    Our study revealed that the change in biomass was strongly influenced by the interaction of soil properties and vegetation features, which explained 55.55%-72.32%of the biomass variation (Fig. 4). This dominant contribution by joint influence to biomass may be explained by the close interaction between vegetation and soil(Liang et al. 2010). As we discussed above, there was a clear co-evolutionary relationship between soil factors(pH, SOC, TN, AN, and TP) and vegetation features(DBH and H) across the restoration periods. This result suggests that the variations in key soil factors (pH, SOC,TN, AN, and TP) were likely to promote the growth of plant and the restoration of vegetation structure (Alday et al. 2012). In turn, vegetation features (DBH, and H)could influence improvements in the soil environment(Fig. 3). These results also offer the further evidence for the hypothesis that the mechanisms of plant and soil promote vegetation restoration synergistically.

    This study also found that soil properties explained 3.30%-31.44% of the variation in biomass,which was basically higher than explanation of vegetation (5.09%-24.32%). This result provides evidence that the importance of soil properties in driving the changes observed in biomass is more than that of vegetation features in the study region, which is most likely because the advantage of hydrothermal conditions in the subtropical region accelerates the material circulation, and promotes the enrichment of soil organic matter (Corlett and Hughes 2015); thus, providing a fertile environment for plant growth. The regulation mechanism of soil properties on biomass development had been discussed previously. With vegetation restoration, the increased in plant species has intensified the competition of aboveground parts for light resources and underground roots for soil resources (Cheng et al. 2011; Li et al. 2017),which further induces the variations of individual growth and morphological structure of trees (DBH and H). As DBH and H increased, more fine materials and litters can be intercepted and accumulated by plants, further enhancing the accumulation of biomass (Li et al. 2017).

    The biomass development at our study site was influenced by different soil and vegetation factors in different restoration periods. In the early restoration period (4-5 years), SOC was the major influential factor (Table 4).The possibility is that SOC is the main source of most nutrients, and that the accumulation of SOC promotes improvements in other soil factors, such as TN, AN, and AP, which have a notable effect on vegetation growth and development (Alday et al. 2012). In the 4-5 years restoration period, SOC content was low (Fig. 2), which is not conducive to the improvement of soil structure or the accumulation of nutrients (Bienes et al. 2016).Therefore, the low SOC not only limits the growth of plant roots, but also intensifies the contradiction between the demand of plant growth for water and nutrients and the supply of soil water and nutrients, resulting in hindrances to plant growth.

    H, pH, and AP were the main factors driving biomass development in the 10-12 years restoration period. This could be attributed to the competition of shrubs for light,which would drive the increasing of H to adapt to interspecific competition (Cheng et al. 2011). Additionally, the accumulation of biomass impels plants to need more Nand P-rich substances (such as enzymes, transport proteins, and amino acids) to participate in metabolic activities (Qin et al. 2016). Therefore, shrubs need to absorb more N and P for growth than do herbs.In particular,P is an important limited factor in red soil area of south China(Gao et al. 2014).However,in the 10-12 years restoration period, the increasing of pH affected the availability of P(Duan et al. 2008), suggesting that the role of AP may intensity the inequity of competition among plants, rather than promote the accumulation of biomass.

    Biomass in the 45-46 years restoration period was conditioned by the synergistic effect of H and pH.The significant effects of H and pH may be caused by a combination of two factors. Firstly, the dominant tree species (Pinus massoniana) of 45-46 years restoration period obtains more light by increasing H and canopy density (Cheng et al.2011),resulting in lower density of woody plants(Table 1);thus,H had a negative effect on biomass.Secondly,low soil pH is beneficial to improve soil permeability, aggregates and porosity(such as BD),and the accumulation of soil nutrients (such as SOC, N and P) (Ramírez et al. 2015), and enhances the availability of P, K, Ca, and Mg (Duan et al.2008).Meanwhile,soil pH decreased with vegetation restoration,and the bioaccumulation and material circulation increased with advantageous hydrothermal conditions(Corlett and Hughes 2015),which were beneficial to the increment of soil nutrient content;thus restoration stimulates the increase of biomass.

    In the old growth forest (sub-climax community), the structure of plant community has reached a state of stable (Peng et al. 2012), which means that the development of vegetation features (DBH and H) has entered a slow growth stage and has less of an impact on biomass.Instead, as a nutrient bank and soil health indicator(Bienes et al. 2016), SOC continues to influence biomass growth. In addition, evergreen trees with a long leaf life need to accumulate more organic substances (such as lignin) to construct defensive structures, and require higher N and P content to maintain normal growth and metabolism (Zeng et al. 2016). Therefore, the supply capacity of soil N and P largely determines the effectiveness of vegetation restoration(Li et al. 2012a).

    Conclusions

    The present work has shown that vegetation restoration can improve significantly soil texture and fertility (especially N, P,and SOC) and vegetation features(species diversity, biomass, DBH, and H). The study showed a clear coupling relationship between some soil factors (pH,SOC, TN, AN, and TP) and vegetation development and structural components (biomass, DBH, and H). At the same time, soil properties and vegetation features had a strongly cooperative influence on the variation of biomass, which suggested that the successful restoration of a degraded forest was driven mainly by their synergistic effect. The individual effect of soil factors on biomass development was greater than that of vegetation factors.Notably, the controlling factors of biomass had differed in the different restoration periods.

    Abbreviations

    AK: Available potassium; AN: Alkaline hydrolysis nitrogen; AP: Available phosphorus; BD:Bulk density; Ca:Total calcium; DBH: Diameter at breast height; K: Potassium; Mg: Total magnesium; MODIS: Moderate Resolution Imaging Spectroradiometer; N: Nitrogen; P: Phosphorus; PCA:Principal component analysis; R2: Adjustment decision coefficient; SOC: Soil organic carbon; SRA: Stepwise regression analysis; H:Height; TK: Total potassium;TN: Total nitrogen;TP: Total phosphorus

    Acknowledgments

    We thank the administrative staff at the Dashanchong Forest Farm,Changsha County, Hunan Province, for their support.

    Authors’ contributions

    XF and CC designed the idea and study, and coordinated the manuscript preparation. CC, XF, WX, PL, SO, and YK processed the data and analyzed the results. CC, XF, WX, and YK contributed to the manuscript writing and editing. All authors read and approved the final manuscript.

    Funding

    This work was supported by the National Forestry Public Welfare Industry Research Project (grant no. 201504411) and the National Natural Science Foundation of China (grant nos. 31570447 and 31300524).

    Availability of data and materials

    The datasets generated and/or analyzed during the current study are not publicly available due[the data is a part of the author’s graduation thesis]but are available from the corresponding author on reasonable request.

    Ethics approval and consent to participate

    Not applicable.

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1Faculty of Life Science and Technology, Central South University of Forestry and Technology, Changsha 410004,China.2Huitong National Field Station for Scientific Observation and Research of Chinese Fir Plantation Ecosystem in Hunan Province, Huitong 438107, China.3Department of Soil Science of Temperate Ecosystems, Georg-August University of Gottingen, 37077 G?ttingen, Germany.4Department of Agricultural Soil Science, Georg-August University of Gottingen, 37077 G?ttingen, Germany.

    Received: 18 November 2019 Accepted: 17 April 2020

    菩萨蛮人人尽说江南好唐韦庄| 亚洲成色77777| 精品人妻1区二区| 国产在线观看jvid| tube8黄色片| 国产真人三级小视频在线观看| 国产成人系列免费观看| 狠狠婷婷综合久久久久久88av| 一区在线观看完整版| www.熟女人妻精品国产| 久久久久网色| 少妇被粗大的猛进出69影院| 又大又爽又粗| 亚洲欧洲精品一区二区精品久久久| 午夜91福利影院| 亚洲国产成人一精品久久久| 日本黄色日本黄色录像| 香蕉国产在线看| 亚洲欧洲精品一区二区精品久久久| a 毛片基地| 男人添女人高潮全过程视频| 国产福利在线免费观看视频| 可以免费在线观看a视频的电影网站| 大片免费播放器 马上看| 国产成人欧美| 久久毛片免费看一区二区三区| 成年人免费黄色播放视频| 亚洲黑人精品在线| 桃花免费在线播放| 在现免费观看毛片| 美女中出高潮动态图| 一级毛片黄色毛片免费观看视频| 大码成人一级视频| 夜夜骑夜夜射夜夜干| 国产视频一区二区在线看| 三上悠亚av全集在线观看| 老汉色∧v一级毛片| 一级黄色大片毛片| 免费看av在线观看网站| 久久精品国产a三级三级三级| 国产日韩欧美在线精品| 午夜免费男女啪啪视频观看| 大香蕉久久网| 最近最新中文字幕大全免费视频 | 婷婷色综合www| 久久久久国产精品人妻一区二区| 纯流量卡能插随身wifi吗| av在线老鸭窝| 天堂中文最新版在线下载| 精品熟女少妇八av免费久了| 午夜日韩欧美国产| 人人妻,人人澡人人爽秒播 | 男女高潮啪啪啪动态图| 国产视频一区二区在线看| 亚洲精品一二三| 麻豆av在线久日| √禁漫天堂资源中文www| 亚洲午夜精品一区,二区,三区| 国产精品久久久av美女十八| 天天操日日干夜夜撸| 777米奇影视久久| 99国产精品一区二区三区| 无遮挡黄片免费观看| 人体艺术视频欧美日本| 日韩中文字幕视频在线看片| 男人爽女人下面视频在线观看| 天天影视国产精品| 永久免费av网站大全| 美女扒开内裤让男人捅视频| 久热爱精品视频在线9| 一区二区日韩欧美中文字幕| 电影成人av| 亚洲欧洲日产国产| 韩国高清视频一区二区三区| 亚洲成人国产一区在线观看 | 嫩草影视91久久| 人人妻人人澡人人爽人人夜夜| 国产一区二区在线观看av| 99久久综合免费| 老熟女久久久| 国产精品av久久久久免费| 久久这里只有精品19| 啦啦啦在线免费观看视频4| 中文字幕制服av| 国产欧美亚洲国产| 成人国产av品久久久| 久久久国产一区二区| 在现免费观看毛片| 国产一卡二卡三卡精品| 最黄视频免费看| 丁香六月欧美| 亚洲精品一区蜜桃| cao死你这个sao货| 亚洲情色 制服丝袜| 久久久久久免费高清国产稀缺| 尾随美女入室| 99香蕉大伊视频| 久久国产亚洲av麻豆专区| 久久人人97超碰香蕉20202| 成人免费观看视频高清| 亚洲成国产人片在线观看| 99国产精品一区二区三区| 午夜福利乱码中文字幕| 国产精品偷伦视频观看了| 看免费av毛片| 狠狠精品人妻久久久久久综合| 成人免费观看视频高清| 中文字幕色久视频| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品.久久久| 亚洲欧美成人综合另类久久久| 国产成人精品久久二区二区免费| 日韩视频在线欧美| 久久久亚洲精品成人影院| 一本—道久久a久久精品蜜桃钙片| 亚洲欧美日韩高清在线视频 | 国产色视频综合| 成人国产av品久久久| 黄片小视频在线播放| 久久中文字幕一级| 男人操女人黄网站| 精品久久蜜臀av无| 国产免费视频播放在线视频| 侵犯人妻中文字幕一二三四区| 国产日韩欧美亚洲二区| 91精品三级在线观看| 亚洲av欧美aⅴ国产| e午夜精品久久久久久久| 啦啦啦在线免费观看视频4| 人人妻,人人澡人人爽秒播 | 女性被躁到高潮视频| 午夜免费观看性视频| 又大又爽又粗| 欧美乱码精品一区二区三区| 18禁国产床啪视频网站| 午夜老司机福利片| 欧美日韩av久久| 午夜福利,免费看| 婷婷色综合大香蕉| av电影中文网址| 国产淫语在线视频| 免费人妻精品一区二区三区视频| 午夜福利视频精品| 欧美在线黄色| 欧美精品一区二区免费开放| 亚洲欧美清纯卡通| 午夜福利,免费看| 久久这里只有精品19| 在线观看国产h片| 一级a爱视频在线免费观看| 久久毛片免费看一区二区三区| 一区在线观看完整版| 欧美日韩国产mv在线观看视频| 久久国产精品男人的天堂亚洲| 天天躁夜夜躁狠狠久久av| 青青草视频在线视频观看| 亚洲av成人不卡在线观看播放网 | 欧美日韩亚洲国产一区二区在线观看 | 新久久久久国产一级毛片| 亚洲七黄色美女视频| 亚洲精品国产色婷婷电影| 在线亚洲精品国产二区图片欧美| 天堂中文最新版在线下载| 欧美另类一区| 老熟女久久久| 中国美女看黄片| 亚洲 欧美一区二区三区| 成年av动漫网址| 叶爱在线成人免费视频播放| 国产高清videossex| 亚洲国产欧美日韩在线播放| 中文字幕人妻丝袜制服| 丰满迷人的少妇在线观看| 三上悠亚av全集在线观看| 欧美xxⅹ黑人| 91精品伊人久久大香线蕉| 国产一区有黄有色的免费视频| 亚洲欧美一区二区三区黑人| 校园人妻丝袜中文字幕| 亚洲国产看品久久| 国产成人av教育| www日本在线高清视频| 考比视频在线观看| 首页视频小说图片口味搜索 | 亚洲国产av影院在线观看| 老司机在亚洲福利影院| 天天添夜夜摸| 国产精品一区二区免费欧美 | 丝袜美腿诱惑在线| 91老司机精品| 青春草视频在线免费观看| 久久精品亚洲av国产电影网| av电影中文网址| 宅男免费午夜| av又黄又爽大尺度在线免费看| 国产日韩欧美视频二区| 肉色欧美久久久久久久蜜桃| av国产精品久久久久影院| tube8黄色片| 91国产中文字幕| 亚洲国产欧美网| 精品人妻在线不人妻| 肉色欧美久久久久久久蜜桃| 黄频高清免费视频| 晚上一个人看的免费电影| h视频一区二区三区| 一级,二级,三级黄色视频| 亚洲av电影在线观看一区二区三区| 侵犯人妻中文字幕一二三四区| 国产不卡av网站在线观看| 最新的欧美精品一区二区| 日韩伦理黄色片| 国产激情久久老熟女| 午夜日韩欧美国产| 伊人久久大香线蕉亚洲五| 亚洲国产日韩一区二区| 欧美人与性动交α欧美精品济南到| 我的亚洲天堂| 每晚都被弄得嗷嗷叫到高潮| 男女床上黄色一级片免费看| 国产精品一国产av| 国产av国产精品国产| 亚洲三区欧美一区| 大香蕉久久网| 国产视频首页在线观看| 国产成人系列免费观看| 亚洲国产av新网站| 另类精品久久| 婷婷色麻豆天堂久久| 精品国产国语对白av| 国产精品国产av在线观看| 操美女的视频在线观看| 飞空精品影院首页| 一区二区三区精品91| 欧美性长视频在线观看| 90打野战视频偷拍视频| 天堂8中文在线网| 日本91视频免费播放| 超色免费av| 国产精品久久久久成人av| 一本一本久久a久久精品综合妖精| 精品久久蜜臀av无| 女人高潮潮喷娇喘18禁视频| 青草久久国产| 国产黄色视频一区二区在线观看| 你懂的网址亚洲精品在线观看| 一级毛片女人18水好多 | 亚洲av男天堂| 中文字幕制服av| 亚洲欧美精品自产自拍| 亚洲精品国产一区二区精华液| 啦啦啦视频在线资源免费观看| 午夜老司机福利片| 十分钟在线观看高清视频www| 免费久久久久久久精品成人欧美视频| 国产99久久九九免费精品| 亚洲专区国产一区二区| 久久av网站| 免费人妻精品一区二区三区视频| 黄色a级毛片大全视频| www.自偷自拍.com| 亚洲欧美精品自产自拍| 国产精品三级大全| 欧美成人午夜精品| www.999成人在线观看| 婷婷色综合大香蕉| 亚洲av片天天在线观看| videos熟女内射| 男女之事视频高清在线观看 | 午夜日韩欧美国产| av有码第一页| 黄色a级毛片大全视频| 国产淫语在线视频| 亚洲精品乱久久久久久| 精品亚洲成a人片在线观看| 免费观看a级毛片全部| a 毛片基地| 久久人妻福利社区极品人妻图片 | 中文字幕另类日韩欧美亚洲嫩草| 亚洲色图综合在线观看| av视频免费观看在线观看| 最新在线观看一区二区三区 | 免费高清在线观看视频在线观看| a 毛片基地| 你懂的网址亚洲精品在线观看| 啦啦啦视频在线资源免费观看| 91麻豆精品激情在线观看国产 | 这个男人来自地球电影免费观看| 精品国产国语对白av| 精品人妻熟女毛片av久久网站| 亚洲国产中文字幕在线视频| 国产精品久久久av美女十八| 少妇的丰满在线观看| 亚洲精品国产av蜜桃| 丝袜人妻中文字幕| 国产视频一区二区在线看| 大型av网站在线播放| 国产在视频线精品| 国产有黄有色有爽视频| 国产成人精品久久二区二区免费| 亚洲精品一二三| 亚洲国产中文字幕在线视频| 成年动漫av网址| 一级毛片 在线播放| 少妇粗大呻吟视频| 精品熟女少妇八av免费久了| 亚洲国产成人一精品久久久| 中文字幕人妻熟女乱码| 国产一级毛片在线| 一区在线观看完整版| 一边摸一边做爽爽视频免费| 天天躁夜夜躁狠狠久久av| 欧美性长视频在线观看| 美女脱内裤让男人舔精品视频| 91精品国产国语对白视频| 亚洲男人天堂网一区| www.精华液| 国产真人三级小视频在线观看| 女人精品久久久久毛片| 久久久国产欧美日韩av| 操出白浆在线播放| 中文字幕亚洲精品专区| 日本午夜av视频| 在线观看www视频免费| 亚洲av国产av综合av卡| 国产日韩欧美视频二区| 日韩人妻精品一区2区三区| 国产老妇伦熟女老妇高清| 黑丝袜美女国产一区| 精品高清国产在线一区| www.av在线官网国产| 高潮久久久久久久久久久不卡| 午夜免费男女啪啪视频观看| 成人影院久久| 91麻豆精品激情在线观看国产 | 亚洲三区欧美一区| 国产精品成人在线| 99热全是精品| 久久精品亚洲av国产电影网| 国产无遮挡羞羞视频在线观看| 免费高清在线观看视频在线观看| 天堂俺去俺来也www色官网| 视频在线观看一区二区三区| 久久久久国产精品人妻一区二区| 18在线观看网站| a 毛片基地| 天天操日日干夜夜撸| 女性被躁到高潮视频| 欧美xxⅹ黑人| 9热在线视频观看99| 51午夜福利影视在线观看| 中文字幕精品免费在线观看视频| 大话2 男鬼变身卡| 国产精品久久久久久精品电影小说| 久久人妻福利社区极品人妻图片 | 这个男人来自地球电影免费观看| 午夜两性在线视频| 亚洲色图 男人天堂 中文字幕| 日韩一区二区三区影片| 国产精品国产三级专区第一集| 国产精品99久久99久久久不卡| 亚洲一码二码三码区别大吗| 亚洲欧美一区二区三区国产| 久久久亚洲精品成人影院| 老司机在亚洲福利影院| av网站免费在线观看视频| 亚洲欧美中文字幕日韩二区| 美女大奶头黄色视频| 水蜜桃什么品种好| 亚洲精品自拍成人| 日韩一区二区三区影片| 只有这里有精品99| 欧美 日韩 精品 国产| 国产熟女午夜一区二区三区| 久久久久久久大尺度免费视频| 婷婷色综合大香蕉| 国产深夜福利视频在线观看| 日韩一卡2卡3卡4卡2021年| av天堂在线播放| 精品一区二区三区四区五区乱码 | 如日韩欧美国产精品一区二区三区| 性色av一级| 后天国语完整版免费观看| 国产成人欧美| 国产麻豆69| 亚洲图色成人| 最新在线观看一区二区三区 | 人人妻人人澡人人看| 中文字幕亚洲精品专区| 日本vs欧美在线观看视频| 国产亚洲精品第一综合不卡| 成人三级做爰电影| 久久精品aⅴ一区二区三区四区| 免费av中文字幕在线| 精品人妻在线不人妻| 亚洲欧洲精品一区二区精品久久久| a级片在线免费高清观看视频| 国产人伦9x9x在线观看| 国产欧美日韩精品亚洲av| 高潮久久久久久久久久久不卡| 欧美日本中文国产一区发布| 在线观看免费高清a一片| 夫妻性生交免费视频一级片| 一本综合久久免费| 亚洲av日韩在线播放| 国产老妇伦熟女老妇高清| 老鸭窝网址在线观看| 日韩 亚洲 欧美在线| 日韩大码丰满熟妇| 亚洲av综合色区一区| 精品国产乱码久久久久久小说| 亚洲图色成人| 精品欧美一区二区三区在线| 午夜免费观看性视频| 欧美精品av麻豆av| 亚洲国产日韩一区二区| 天天躁日日躁夜夜躁夜夜| 欧美xxⅹ黑人| 精品少妇一区二区三区视频日本电影| 免费观看a级毛片全部| 自线自在国产av| 亚洲伊人久久精品综合| 一个人免费看片子| 久久久精品94久久精品| 人人妻人人爽人人添夜夜欢视频| 少妇猛男粗大的猛烈进出视频| 亚洲成色77777| 无限看片的www在线观看| 最近最新中文字幕大全免费视频 | 日本一区二区免费在线视频| 97在线人人人人妻| 国产精品免费大片| 亚洲国产精品999| 校园人妻丝袜中文字幕| 色网站视频免费| 无限看片的www在线观看| 人成视频在线观看免费观看| 一边摸一边做爽爽视频免费| 欧美日韩一级在线毛片| 亚洲少妇的诱惑av| √禁漫天堂资源中文www| 免费在线观看完整版高清| 少妇 在线观看| 欧美人与善性xxx| 欧美日韩亚洲国产一区二区在线观看 | 国产精品一区二区免费欧美 | 少妇被粗大的猛进出69影院| 一级毛片 在线播放| 国产一级毛片在线| 亚洲美女黄色视频免费看| 丰满少妇做爰视频| 宅男免费午夜| 自拍欧美九色日韩亚洲蝌蚪91| 国产一级毛片在线| 又大又爽又粗| 建设人人有责人人尽责人人享有的| 欧美成狂野欧美在线观看| 婷婷色av中文字幕| 两人在一起打扑克的视频| 久久国产精品男人的天堂亚洲| 亚洲欧洲精品一区二区精品久久久| 日本黄色日本黄色录像| 夜夜骑夜夜射夜夜干| 日韩人妻精品一区2区三区| 国产精品久久久久久精品古装| 久久精品久久久久久噜噜老黄| 亚洲成人免费av在线播放| 2018国产大陆天天弄谢| 尾随美女入室| 久久精品国产综合久久久| cao死你这个sao货| 美女视频免费永久观看网站| 色婷婷av一区二区三区视频| 久久久精品免费免费高清| 一级毛片女人18水好多 | 少妇 在线观看| 后天国语完整版免费观看| 极品少妇高潮喷水抽搐| www.熟女人妻精品国产| 五月开心婷婷网| 婷婷色综合www| 亚洲精品日韩在线中文字幕| 在线亚洲精品国产二区图片欧美| 久久久久国产一级毛片高清牌| kizo精华| 国产免费一区二区三区四区乱码| 母亲3免费完整高清在线观看| 国产精品久久久久久精品古装| 操美女的视频在线观看| 国产精品国产av在线观看| 日日摸夜夜添夜夜爱| 免费观看av网站的网址| 国产黄色视频一区二区在线观看| 国产欧美日韩综合在线一区二区| 欧美精品高潮呻吟av久久| 国产精品一区二区在线不卡| 日韩熟女老妇一区二区性免费视频| 老司机亚洲免费影院| 青草久久国产| 午夜福利视频在线观看免费| 国产高清videossex| 性高湖久久久久久久久免费观看| www日本在线高清视频| 久久性视频一级片| 女人高潮潮喷娇喘18禁视频| 国产高清国产精品国产三级| 性高湖久久久久久久久免费观看| 9热在线视频观看99| 国产熟女午夜一区二区三区| 天天影视国产精品| 高潮久久久久久久久久久不卡| 一级毛片黄色毛片免费观看视频| 亚洲三区欧美一区| 国产97色在线日韩免费| 熟女av电影| 大片免费播放器 马上看| 亚洲国产欧美一区二区综合| 中文字幕亚洲精品专区| 2018国产大陆天天弄谢| 岛国毛片在线播放| 久久久久网色| 日韩伦理黄色片| 国产淫语在线视频| 国产成人免费无遮挡视频| 欧美在线一区亚洲| 狠狠精品人妻久久久久久综合| 欧美日韩亚洲国产一区二区在线观看 | 久久九九热精品免费| 丁香六月天网| 中国国产av一级| 亚洲av电影在线进入| 一级黄色大片毛片| 免费不卡黄色视频| 亚洲,欧美精品.| 一本综合久久免费| 欧美精品一区二区免费开放| 国产成人系列免费观看| 看免费av毛片| 精品福利永久在线观看| 久久人人97超碰香蕉20202| 夫妻午夜视频| 免费高清在线观看日韩| 亚洲熟女毛片儿| 成人国语在线视频| 亚洲精品美女久久久久99蜜臀 | 久久99精品国语久久久| 丝袜喷水一区| 久久久久久久精品精品| 亚洲欧美一区二区三区黑人| 少妇粗大呻吟视频| 国产精品欧美亚洲77777| 久久精品人人爽人人爽视色| 亚洲欧美日韩另类电影网站| 最新在线观看一区二区三区 | 久久久精品94久久精品| 欧美日韩综合久久久久久| 日韩中文字幕欧美一区二区 | 午夜福利,免费看| 亚洲自偷自拍图片 自拍| 成人黄色视频免费在线看| 青春草视频在线免费观看| 免费一级毛片在线播放高清视频 | 大陆偷拍与自拍| 日韩中文字幕欧美一区二区 | 亚洲男人天堂网一区| 超碰成人久久| 精品国产乱码久久久久久小说| 两个人免费观看高清视频| av国产久精品久网站免费入址| 在线观看免费日韩欧美大片| 中文字幕亚洲精品专区| 操美女的视频在线观看| 人妻一区二区av| 精品少妇久久久久久888优播| 久久国产精品影院| 亚洲,欧美精品.| 亚洲av电影在线观看一区二区三区| 亚洲成av片中文字幕在线观看| 国产av国产精品国产| 五月开心婷婷网| 91九色精品人成在线观看| 少妇裸体淫交视频免费看高清 | 国产精品一区二区在线不卡| 中文字幕人妻丝袜制服| 午夜视频精品福利| 久久久久精品国产欧美久久久 | 又大又爽又粗| 国产一区有黄有色的免费视频| 亚洲,一卡二卡三卡| 久久久久久久久久久久大奶| 精品人妻一区二区三区麻豆| 欧美精品av麻豆av| 天堂中文最新版在线下载| 国产欧美日韩一区二区三区在线| 伦理电影免费视频| 欧美xxⅹ黑人| www.自偷自拍.com| 久久久久久久久免费视频了| 亚洲黑人精品在线| 亚洲少妇的诱惑av| 脱女人内裤的视频| 又大又爽又粗| 精品一区二区三卡| 黄色毛片三级朝国网站| 日韩av在线免费看完整版不卡| tube8黄色片| 美女福利国产在线| 国产色视频综合| 日韩欧美一区视频在线观看| 最近手机中文字幕大全| 欧美日韩福利视频一区二区| 国产午夜精品一二区理论片| av不卡在线播放| 亚洲国产av新网站| 捣出白浆h1v1|