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

    Effect of thinning intensity on understory herbaceous diversity and biomass in mixed coniferous and broad-leaved forests of Changbai Mountain

    2021-02-28 09:11:58GerongWangYueSunMoZhouNaiqianGuanYuwenWangRunhuaJiangZhiyuLiuMengjiaWuandFucaiXia
    Forest Ecosystems 2021年4期

    Gerong Wang,Yue Sun,Mo Zhou,Naiqian Guan,Yuwen Wang,Runhua Jiang,Zhiyu Liu,Mengjia Wu and Fucai Xia*

    Abstract

    Keywords: Selective thinning, Mixed coniferous and broad-leaved forests, Herbaceous plants, Diversity, Biomass

    Background

    Understory herbaceous plants are an important part of the species diversity of forest ecosystems (Han et al. 2016). In temperate forests, more than 80% of vascular plants are herbaceous plants (Whigham 2004), and herbaceous plants are of great significance for maintaining the species diversity of forest ecosystems. Compared with trees and shrubs, understory herbs are more sensitive to disturbances and are the most adaptable to the environment. If the disturbance severely destroys the original forest ecosystem,herbaceous plants will usually recover fast,thereby quickly improving the soil site conditions,reducing the intensity of solar radiation on the ground,and providing important protection for the regeneration of trees and shrubs and the growth of seedlings (Huang et al. 2016).Therefore, the species diversity and biomass of the herb layer can reflect the vegetation restoration status after forest disturbance.

    Forest thinning is an important disturbance in the forest managements. Thinning can have an important impact on the stability of the forest ecosystem, and even destroy the balance of the forest ecosystem (Lei et al. 2005). As an important disturbance in the forest managements,thinning has long attracted people’s attention (Zhou et al. 2012). Reasonable thinning can maintain the stability and structure of the forest ecosystem while obtaining timber,and achieve sustainable forest management (Zhou et al. 2015). Thinning can affect forest growth,understory vegetation regeneration,forest spatial structure,forest species diversity(Rong et al.2014;You et al.2015;Gong et al.2015;Huang 2016).Some studies have reported that thinning can not only promote the growth of tree height and diameter at breast height,increase forest stock and improve wood quality (An et al. 2012), but also improve the understory environment,reduce understory competition among species,and promote the increase in the growth and diversity of understory species, especially the understory herb layer(Li et al.2012).

    Due to long-term human disturbance, the original broad-leaved Korean pine forests in Northeast China have been severely damaged, and now there are a large number of secondary coniferous and broad-leaved mixed forests (Zhang et al. 2009). Reasonable thinning can accelerate the succession rate of this forest type, realize its sustainable management, and play an important role in maintaining the stability of timber yield and species diversity (Liu et al. 2014). Through thinning of different intensities in the near-mature secondary coniferous and broad-leaved mixed forest in the State-owned Forest Protection Center of the Jilin Province Forestry Experimental Zone, this study investigated the species diversity and the seasonal changes of biomass of the understory herb layer after thinning, and evaluated the influence of thinning on herbaceous plants. The study aims to provide a theoretical basis for reasonable tending and thinning for secondary coniferous and broad-leaved mixed forests in Northeast China.

    Materials and methods

    Study area

    The study site was located at Dapo Forest Farm (43°57′ N,127°43′ E), the state-owned forest protection center of Jilin Province Forestry Experimental Zone. In July 2011, 4 permanent monitoring plots with an area of 1 ha (100 m ×100 m)were set up in the near-mature forest.Each plot was divided into 25 plots of 20 m × 20 m and the forest stand survey was conducted. In December 2012, 4 plots were thinned with actual intensity of 0,17.2%,32.7%,and 50.8%respectively, which were defined as unthinning (CK), light thinning (LT), moderate thinning (MT) and heavy thinning(HT)(Liu et al.2019).

    Herbaceous plant diversity

    At the beginning of May 2015,4 herbal plots of 1 m ×1 m were randomly set up in each 20 m×20 m plot,and a total of 100 herbal plots were set up. From early May to early October 2015, the species name, abundance, average height,and coverage of the herb plots in the four plots were investigated,and the total coverage was also recorded.The relative abundance(RA),relative coverage(RC),relative height(RH),and importance value(IV)of herbaceous plants based on the obtained data were calculated. Furthermore, the Margalef index (R), Shannon-Wiener index (H), Simpson index (D)and Pielou index(Jsw)were also calculated as follows:

    where IV is the importance value, RA is the relative abundance,RC is the relative coverage, RH is the relative height, S is the total number of species in the sample plot, N is the total individuals of species in the sample plot, and Niis the individual of the i-th species Number,Piis the relative value of i, and the calculation formula of Piis:

    Herbaceous plant biomass

    In early May, mid-July, and early September 2015, the understory herbaceous biomass was sampled and surveyed three times. Specifically, two 1 m × 1 m subplots were randomly set up in each 20 m × 20 m plot. We investigated all species names, abundance, height, coverage, total coverage, and finally harvested the whole plant for all the herbs in the subplot. The samples were taken back to the laboratory to clean the soil at the roots with water. The above- and belowground parts were cut with scissors, and weighed their fresh weight (MA for aboveground part and MB for belowground part). Then the sample was first quenched at 105°C for 30 min, adjusted to 85 °C and dried to a constant weight. After drying,weigh the dry weight of the aboveground and belowground parts and record them as WA and WB. The herb biomass of the entire subplot was calculated based on the coverage, the dry weight of the herbaceous plants,and the area of the entire subplot (Sun et al. 2017).

    Data analysis

    Correlation analysis and one-way analysis of variance(one-way ANOVA) were performed to test thinning effects on herbaceous diversity and biomass. A total of 100 1 m × 1 m herbal plots in each thinning treatment was defined as replicates. Multiple comparisons were performed using the least significant difference method(LSD) once the thinning effects were significant. All data analysis was performed using SPSS18.0 software. The significance level was set at α <0.05.

    Results

    Effect of thinning on plant diversity

    Species composition

    During the whole growing season, a total of 87 species of herbaceous plants were investigated in 4 plots,belonging to 37 families and 74 genera. From the distribution of each family, the number of species belonging to Ranunculaceae, Liliaceae, Umbelliferae and Compositae was the largest, with 11, 7, 7 and 5 species respectively. There were 19 families containing only one species and 8 families including 2 species. Only 2 families (Rosaceae and Papaveraceae) included 3 species. Four families(Cyperaceae, Rubiaceae, Caryophyllaceae and Urticaceae)included 4 species, accounting for 10.8% of all families.From the distribution of the genera, there are 64 genera with 1 species, accounting for 86.5% of the total number of genera. There are 10 genera with more than 2 species, accounting for 13.5% of the total number of genera. Among them, there are 7 genera (Lilium,Viola, Geranium, Hippochaete, Cimicifuga, Arisaema and Anemone) containing 2 species and 3 genera(Rubia, Carex, and Urtica) containing 3 species. In summary, a single family and a single genus had an absolute advantage in the survey plots, and the phenomenon of few families and few genera was obvious, which was a typical species distribution of temperate forest herbs.

    Compared with the control, thinning increased the number of species and genera of herbaceous plants.Light thinning, moderate thinning and heavy thinning increased 2 species and 3 genera, 8 species and 9 genera,1 species and 1 genus, respectively. For the family, moderate thinning increased 1 family, while light thinning and heavy thinning both decreased by 1 family compared with the control (Table 1).

    The number of species, genera and families of herbaceous plants in each thinning treatment showed seasonal variations. In May, the number of species, genera and families in the thinning were higher than those of the control, and the number of species, genera and families in the heavy thinning plot was the largest. In June,the number of species, genera and families in the thinning were all low compared with the control, and light thinning had the smallest number of species, genera and families. In July, August and September, moderate thinning had the largest number of species, genera and families while heavy thinning was opposite. In October, the number of herbaceous species in each plot decreased sharply,with little difference among thinning treatments.The number of species in the control and moderate thinning was slightly higher than that of light thinning and heavy thinning.

    Seasonal variation of importance values

    The importance values of herbaceous plants in different thinning in each month are shown in Table 2. In May,the top 5 importance values of herbaceous plants in the control were Anemone raddeana, Hylomecon japonica,Aegopodium alpestre, Meehania urticifolia and Corydalis ambigua. The sum of the importance values of these plants was 51.38%, dominating in the plot. Among them, there were 15 species with importance values greater than 1%, and the sum of their importance values accounted for 93.54% of all herb species. Carex pilosa,

    Anemone raddeana, Filipendula palmata, Meehania urticifolia, and Cardamine leucantha were the top 5 importance values of herbaceous plants in the light thinning, with the sum of the importance values of 51.21%. There are 14 species of herbs with importance value greater than 1%, and the sum of their importance values accounted for 88.35% of all herb species. Carexpilosa, Anemone raddeana, Filipendula palmata, Meehania urticifolia and Poa annua were the top 5 importance values of herbaceous plants in the moderate thinning. The sum of the importance values of these plants was 49.36%, and there were 18 species with the importance value greater than 1% and the sum of their importance values accounted for 89.56% of all herbaceous species. Filipendula palmata, Aegopodium alpestre, Anemone raddeana, Anemone amurensis, and Meehania urticifolia were the top 5 herbaceous plants in the heavy thinning plots. The sum of the important values of these plants was 46.09%. There were 17 species of herbaceous plants with the importance value greater than 1%, and the sum of their importance values accounted for 89.93% of all herb species. The proportion of important value of early spring herbs (including

    Table 1 Changes in the number of families,genera and species of herbaceous plants in each thinning treatment in different months

    Anemone raddeana, Anemone amurensis, Hylomecon japonica, Gagea lutea, Corydalis ambigua, Adonis amurensis, etc.) in the control was relatively high(42.64%), while that of light thinning, moderate thinning and heavy thinning were 24.90%,18.64%and 22.09%respectively, which was significantly lower than that of the control, indicating that the phenological period of herbs after thinning was advanced. Some early spring herbs in early May had completed the whole life cycle and died.

    In June, the top 5 important values of herbaceous plants in the control were Carex pilosa, Poa annua,Meehania urticifolia, Aegopodium alpestre and Filipendula palmata. The sum of the importance values of these plants was 50.88%. There were 18 species with the importance value greater than 1, and the sum of their importance values accounted for 92.57% of all herbaceous species. The top 5 important values of herbaceous plants in the light thinning were Carex pilosa,Athyrium multidentatum, Aegopodium alpestre, Filipendula palmata and Brachybotrys paridiformis. The sum of the importance values was 59.92%, and there were 16 species with the importance value greater than 1. The top 5 importance values of herbaceous plants in the moderate thinning and heavy thinning plots were Carex pilosa, Meehania urticifolia, Carex remotiuscula, Aegopodium alpestre and Filipendula palmata, but their importance values were different. The order of importance values in the moderate thinning was Carex pilosa >Meehania urticifolia >Carex remotiuscula >Aegopodium alpestre >Filipendula palmata. The importance values of the heavy thinning were Carex pilosa >Filipendula palmata >Aegopodium alpestre >Carex remotiuscula >Meehania urticifolia. In the moderate thinning and heavy thinning, the sum of the importance values of these plants was 58.10 % and 56.22%, respectively. Herbs with importance values greater than 1% were 20 and 17,and accounted for 91.30%and 91.49%of all herb species,respectively. The number of species in the control was large, because the poor light conditions in the control led to late plant phenology, and there were more early spring herbs in the control, which led to the existence of early summer herbs and early spring herbs at the same time,resulting in a large number of species; while in the thinning plot,due to the better light conditions caused by thinning,most of the early spring herbs had died,resulting in the number of species has declined.

    In July, the top 5 importance values of herbaceous plants in the control were Carex pilosa, Laportea bulbifera, Filipendula palmata, Cardamine leucantha, and Meehania urticifolia. The sum of the importance values of these plants was 56.33%, there were 16 species with the importance value greater than 1%, and the sum of their importance values accounted for 90.20%of all herb species. The top 5 important values of herbaceous plants in the light thinning were Carex pilosa, Filipendula palmata, Aegopodium alpestre, Matteuccia struthiopteris and Anemone amurensis. The sum of the importance values of these plants was 60.67%, and there were 19 species with the importance value greater than 1%. The sum of their importance values accounted for 90.52% of all herb species. The top 5 important values of herbaceous plants in the moderate thinning and heavy thinning plots were Carex pilosa, Meehania urticifolia,

    Carex remotiuscula, Aegopodium alpestre and Filipendula palmata. The order of importance values in the moderate thinning was Carex pilosa >Filipendula palmata >Carex remotiuscula >Meehania urticifolia >Aegopodium alpestre, and the importance values of the plots with heavy thinning were Carex pilosa >Filipendula palmata >Aegopodium alpestre >Meehania urticifolia >Carex remotiuscula. In the moderate and heavy thinning, the sum of the importance values of these plants was 54.94 % and 60.25%, respectively. The herb plants with importance values greater than 1% were 21 and 22 species and accounted for 92.47 % and 93.50% of all herb species, respectively. 38, 39, 46 and 36 species of herbs were found in control, light thinning, moderate thinning and heavy thinning respectively. Compared with the control, moderate thinning significantly increased the number of herbs. This is due to the good light under the forest and the invasion of light loving herbs species after moderate thinning, and there are more shade environment under the forest and more shade loving herbs are retained, which makes the light loving and shade loving species co-exist, resulting in a relatively large number of species. After heavy thinning,the shade environment under the forest is less and the light condition is excellent, which is especially suitable for the survival of photophilous herbs, but the loss of photophilous species.

    Table 2 The top 5 herbs with relative importance value in different months

    In August, the top 5 importance values of herbaceous plants in the control were Milium effusum, Meehania urticifolia, Carex pilosa, Aegopodium alpestre and Filipendula palmata. The sum of the importance values was 60.17%, and there were 15 species with the importance value greater than 1, accounting for 90.59% of all herb species. The top 5 importance values of herbaceous plants in the light thinning were Carex pilosa, Aegopodium alpestre, Filipendula palmata, Brachybotrys paridiformis and Meehania urticifolia, with the sum of the important values of 58.26%. There were 19 species with the importance value greater than 1%, accounting for 94.23% of all herb species. The top 5 important values of herbaceous plants in the moderate thinning were Carex pilosa, Meehania urticifolia, Aegopodium alpestre, Filipendula palmata and Brachybotrys paridiformis. The sum of the important values was 59.42%,and there were 16 species with the importance value greater than 1%, accounting for 87.55% of all herb species. The top 5 important values of herbaceous plants in the heavy thinning were Carex pilosa, Aegopodium alpestre, Filipendula palmata, Carex callitrichos and Meehania urticifolia. The sum of the importance values was 54.16%, and there were 21 species with the importance value greater than 1%, accounting for 95.51% of all herb species. Among the top 5 species of herbs in different thinning plots, there are 4 common species, namely Carex pilosa, Meehania urticifolia, Aegopodium alpestre and Filipendula palmata. From the important values of these 4 species, the important values of Carex pilosa increased significantly after thinning, while that of Meehania urticifolia decreased. Aegopodium alpestre and Filipendula palmata were relatively consistent in different plots, and there was no significant change.

    In September, the top 5 importance values of herbaceous plants in the control were Carex pilosa, Aegopodium alpestre, Milium effusum, Meehania urticifolia,and Cardamine leucantha, the sum of the important values was 67.97%. There were 11 species of herbs with importance values more than 1%, accounting for 94.00% of all herb species. The top 5 important values of herbaceous plants in the light thinning were Carex pilosa, Aegopodium alpestre, Filipendula palmata,Athyrium multidentatum and Carex remotiuscula. The sum of the important values was 69.67%. There were 15species with importance values greater than 1%, accounting for 96.65% of all herb species. The top 5 importance values of herbaceous plants in the moderate thinning were Carex pilosa, Meehania urticifolia, Aegopodium alpestre, Filipendula palmata and Athyrium multidentatum. The sum of the important values was 65.76%. There were 16 species with importance value more than 1%, accounting for 92.82% of all herb species. The top 5 importance values of herbaceous plants in the heavy thinning were Carex pilosa, Filipendula palmata, Aegopodium alpestre, Urtica angustifolia and Meehania urticifolia. The sum of the important values was 78.24%. There were 13 species with the importance value greater than 1%, accounting for 99.14% of all herb species. 23, 21, 25 and 14 species of herbs were found in the control, light thinning, moderate thinning and heavy thinning, respectively. This is due to the fact that more photophilous species, such as Impatiens noli-tangere,have invaded the plots after heavy thinning. The earlier invasive species will occupy a wide area in the sample plot, and prevent the invasion of other species. Consequently, the number of species will decrease after heavy thinning. With the coming of autumn, the number of photophilous species in the heavy thinning plots in September was lower because of their higher individuals and more easily withered and died.

    Table 2 The top 5 herbs with relative importance value in different months (Continued)

    In October, the top 5 importance values of herbaceous plants in the control were Carex pilosa, Meehania urticifolia, Aegopodium alpestre, Milium effusum and Cardamine leucantha, the sum of the importance value was 92.65%. There are 8 species with importance value greater than 1%, accounting for 98.60% of all herb species. The top 5 importance values of herbaceous plants in the light thinning were Carex pilosa, Aegopodium alpestre, Meehania urticifolia, Carex callitrichos and Filipendula palmata. The sum of the important values was 94.52%, and there were 7 species with the importance value greater than 1%, accounting for 99.32% of all herb species. The top 5 importance values of herbaceous plants in the moderate thinning were Carex pilosa,Aegopodium alpestre, Carex remotiuscula, Meehania urticifolia and Carex callitrichos. The sum of the important values of was 94.91%. There were 6 species with importance value greater than 1%, accounting for 97.38% of all herb species. The top 5 importance values of herbaceous plants in the heavy thinning were Carex pilosa, Aegopodium alpestre, Carex callitrichos, Meehania urticifolia, and Carex remotiuscula.The sum of the important values was 88.10%. There were 8 species with important value greater than 1%,accounting for 99.34% of all herb species. There were 10, 8, 10 and 9 species of herbs in the control, light thinning, moderate thinning and heavy thinning, respectively. There was no significant difference in the number of species in different plots, and the species in different plots were almost the same.

    Seasonal changes in diversity index

    Calculation of diversity index showed that there were significant variations (P <0.05) in different thinning. Except for the control, which showed an increasing trend from May to June, the Margalef index showed a similar trend in all plots in other months (Fig. 1a). It gradually decreased from May to July, increased from July to August, and then fell rapidly from August to October. The Margalef index of different treatments showed various patterns in different months with moderate thinning >heavy thinning >light thinning >control in May, and control >moderate thinning >light thinning >heavy thinning from June to August.From August to October, the Margalef index of the four plots all showed a significant decreasing trend with control >light thinning >moderate thinning >heavy thinning.

    In the four treatments, the Shannon-Wiener index and Simpson index had similar patterns (Fig. 1b and c). The Shannon-Wiener index and Simpson index gradually decreased from May to July, increased gradually from July to August, and then decreased sharply. In June, July and September, the Shannon-Wiener index and Simpson index in the control were higher than the other three thinning treatments. In May, August and October, the Shannon-Wiener index and Simpson index of the four plots were relatively the same.

    Throughout the growing season, the Pielou index was roughly in the shape of an inverted “N” (Fig. 1d), namely it had a high value at the beginning of growth, a low value in June or July, and the second peak in August,and then it decreased sharply. The Pielou index of the control was generally higher than that of the other thinning treatments. The control plot and the heavy thinning showed the first decrease in June, while the light and moderate thinning showed the first decrease in July.

    Effect of thinning on biomass

    Seasonal changes in biomass

    In May, thinning can promote the increase of total biomass, above-ground biomass and belowground biomass of herbaceous plants, but there were no significant differences among thinning treatments. Total biomass and above-ground biomass were the largest in light thinning,and the belowground biomass was the largest in moderate thinning. In terms of root-to-shoot ratio, the moderate thinning was significantly higher than the other three treatments, and the light and heavy thinning were not significantly different from the control. Compared with the control, the number of herbaceous plants decreased after light thinning, but the difference was not obvious.Moderate and heavy thinning significantly increased the number of herbaceous plants. In terms of average height and total coverage, moderate thinning was significantly higher than that of light thinning.

    Fig.1 Monthly change of herbaceous plant diversity index

    The biomass of herbaceous plants increased in July compared with May. Compared with the control, thinning promoted the increase of total biomass, above- and belowground biomass. Among them, the total biomass,above- and belowground biomass were significantly higher in heavy thinning than those of the control. The difference between moderate thinning and the control was not significant. In terms of root-shoot ratio, the difference between the four treatments was not significant,indicating that thinning intensity did not change the biomass allocation between the above- and belowground parts. Thinning increased the number of plants, and the number increased with increasing thinning intensity.The average height of herbaceous plants after thinning decreased, but the difference was not significant. In addition, thinning increased the total coverage of herbaceous plants. The coverage was significantly higher in the heavy thinning than that in the control, while the light and moderate thinning were not significantly different from the control.

    Compared with the control, thinning significantly increased the total herb biomass and belowground biomass in September. Heavy thinning rather than light thinning and moderate thinning significantly increased above-ground biomass. After thinning, the root-shoot ratio increased significantly, but the difference in the root-shoot ratio was not significant among light, moderate and heavy thinning, indicating that thinning changed the distribution of above- and belowground biomass of herbaceous plants. The nutrients were mostly allocated to the belowground part, but the thinning intensity did not change the distribution of plant biomass between above- and belowground parts. Thinning increased the number of herbaceous plants, and the number of herbaceous plants in the moderate and the heavy thinning was significantly higher than that in the control. The average height of herbs decreased after thinning, which was evident in light thinning. Thinning significantly increased the total coverage of herbs, but there were no significant differences among light, moderate and heavy thinning.

    The correlation analysis of the biomass in different months is shown in Fig. 2. In May, the total biomass significantly correlated with the number, average height,and the total coverage, but it did not significantly correlate with the root-shoot ratio. Above-ground biomass was positively correlated with the number, average height and total coverage, but not correlated with the root-shoot ratio. Belowground biomass was positively correlated with root-shoot ratio, the number and average height, but did not change with the total coverage.

    Fig.2 Correlation of various biomasses of herbaceous plants in different months

    In July, the number and total coverage of herbaceous plants were positively correlated with total biomass,above- and belowground biomass. The average height of herbaceous plants did not correlate with total biomass,above- and belowground biomass. The root-shoot ratio significantly correlated with total biomass and belowground biomass.

    In September, the root-shoot ratio positively correlated with total biomass and belowground biomass. The number of herb plants had a significantly and positively correlation with total biomass and belowground biomass. The average height of herb plants was significantly and negatively correlated with the total biomass and above-ground biomass.The total coverage of herb plants was positively related to the total biomass, above- and belowground biomass.

    It was found that the herb biomass was different over time (Table 3). In the control, the total biomass, aboveand belowground biomass all showed their maximum in July. After thinning, the total biomass, above- and belowground biomass all gradually increased over time and reached their maximum in September. In the thinning treatments, the total biomass, above- and belowground biomass in September were significantly higher than those in May.

    Correlations between diversity and biomass

    The correlation analysis between herb biomass and diversity index in different months is shown in Fig. 3.In May, total herb biomass, above- and below ground biomass were significantly and positively correlated with Shannon-Wiener index, Simpson index,and Pielou index. The correlation between herb biomass and diversity index in July and September was not significant. Compared with May, the relationship between species diversity and biomass were changed,indicating that seasonal changes will change the relationships.

    Discussion

    Importance values variation

    Forest thinning, as an important forest management practice, is an important way to change the composition and structure of the forest, adjust forest density, and improve the understory environmental conditions. It will directly affect the composition and diversity of understory vegetation, and influence tree growth and community stability (Zhang and Tang 2008). The importance values of each species varied significantly amongdifferent thinning treatments. The reason may be that the growth of overstorey changes with season, which in turn causes changes in the understory environment. The importance values of species in different treatments showed similar patterns with seasonal changes. The importance values of early spring plant were greater in May, and the top 5 importance values of various treatments were Anemone amurensis, Hylomecon japonica,Anemone raddeana, and Corydalis ambigua. Furthermore, these plants constitute the main body of herbaceous plants under the broad-leaved Korean pine forest in early spring. At this time, the overstorey has not yet reached canopy closure, the understory has sufficient sunlight, and the understory temperature continues to rise. In addition to early spring plants, other herbaceous plants have also sprouted in large numbers, so there are more species in May (Chen and Matter 2017; Kleinn et al. 2020). In June, the trees grow fast, but the overstorey is not completely closed. The light intensity of understory gradually decreases, and the plants wither and die in early spring. Except for the control, the number of species significantly decreased in other thinning treatments. At this time, Carex pilosa, Aegopodium alpestre, Meehania urticifolia, and Filipendula palmata replaced early spring plants as the dominant species of understory herbs, and ferns also showed greater advantages due to their vigorous growth (Gao et al. 2019). The overstory grows fast in July. The canopy density is relatively high, and the understory light intensity is low. At this time, the understory is mainly covered by shadetolerant species such as Carex pilosa, Carex remotiuscula, Carex callitrichos, Carex siderosticta, Meehania urticifolia, Laportea bulbifera, Urtica angustifolia. As precipitation begins to increase, understory wet plants grow quickly, leading to an increase in the number of species (Hulík and Douda 2017). The overstory reaches its growth peak in August, and the canopy closure continues to increase, resulting in a continuous decrease in understory light intensity. Early spring plants and some early summer plants withered and died, and the number of understory herbaceous species declined significantly,but the dominant species did not change (Gong et al.2016).The temperature drops significantly in September,and the overstory stops growing, some tall trees begin to defoliate, and the canopy closure is reduced. At this time, most of the herbaceous plants also stop the growth of the above-ground parts and enter the nutrient reserve period, which causes a significant decrease in the number of species. Understory herbaceous plants still consist of dominant species such as Carex pilosa, Aegopodium alpestre and Meehania urticifolia (Wu et al. 2018). In October, the Changbai Mountain area has entered the end of autumn. Although most of the trees in the upper layer of the canopy have fallen leaves, there is sufficient sunlight understory. But due to the low temperature,most of the herbaceous plants have withered and died,causing the number of species to drop to the lowest value in the growing season (Dai et al. 2013). Based on the data of the entire growing season, the top herbaceous plants in the various treatments are mainly composed of Carex pilosa, Aegopodium alpestre, Meehania urticifolia and Filipendula palmata, indicating that these species have strong adaptability and a wide ecological range.

    Table 3 Biomass of herbaceous plants in different months

    Herbaceous plant species variation

    Fig.3 Correlation between herb biomass and diversity index in different months

    The monthly changes in the number of herbaceous plant species in different thinning treatments indicate that the number of species in the thinning in May is higher than that in the control. This may be due to the increased understory light intensity after thinning, which leads to the advance of plant phenology (Ye et al. 2015). The number of species in the thinning treatments decreased in June. This is due to the advance of the phenological period after thinning that causes some early spring short-lived plants to wither and die, resulting in a decrease in the number of species in the thinned treatments (Li et al. 2006). From July to September, the number of species in the four treatments gradually decreased. This is because July is the rainy season in this area. The soil under the forest is moist and the temperature is high, which is most suitable for plant growth. Therefore, the number of species in July reached the maximum. Although the enhancement of understory light intensity might have caused the death of some shade-tolerant species in moderate thinning, heliophile will quickly invade the site, leading to species reaching the maximum in the moderate thinning during this period. Although light thinning may not significantly increase light intensity, the thinning will form gaps, which will lead to the death of some rare shade-tolerant herbaceous plants, and the heliophile cannot invade the plot,resulting in the decrease of species numbers. In the heavy thinning treatment, the understory sunlight is significantly enhanced and there are more gaps, causing the death of most of the shade-tolerant herbaceous plants.Although the heliophile can invade the site, they cannot compensate for the death number of shade-tolerant plant, resulting in a significant decrease in species numbers (Ding et al. 2016). The understory temperature and light intensity of the four treatments are very low in October, and the understory environmental conditions are roughly the same. Most of the herbaceous plants withered and died, and only a few herbs with strong cold tolerance could survive(Small and McCarthy 2002).

    Herbaceous plant individuals number variation

    Throughout the growing season, the number of individual plants in the heavy thinning is the largest, indicating that thinning can promote the increase in the number of individual plants. This is due to the increase in the intensity of thinning, which leads to the increase of understory light intensity, subsequently promoting the increase of understory herb plants (Westerband and Horvitz 2017; Didion 2020). The number of individuals in the moderate thinning was higher than that in the heavy thinning in July. This is due to the presence of large herb plants in the heavy thinning, occupying a larger space under the forest, resulting in the death of some small one (Zhang et al. 2019). The number of individual plants in the control showed a decreasing trend.This was because most of the control in May were early spring plants. Although the early spring plants are small,there are more individuals on the same stem, so the number of individuals was the largest. Early spring plants gradually die from June to August, early summer plants and late summer plants dominant the herbaceous plants. These plants are relatively large, occupying more space understory, resulting in a decrease in the number of individual plants. The overstory defoliate in September, the increase in understory light intensity increases the ground temperature, and the secondary growth of some plants leads to an increase in the number of individuals. The temperature gradually decreases in October,and most understory herbs are not cold-resistant, and will gradually wither and die, leading to a decrease in the number of individuals (Grant et al. 2019). The increase in the number of individual plants in May and June after thinning is due to the increase in understory light and the increase in understory temperature due to thinning, which advances the phenological period of plants. There are two types of plants including early spring plants and early summer plants in June (especially in heavy thinning). At this time, the plant individuals are small, and the understory space can meet the needs of plant growth. Therefore, the number of individual plants tends to increase.After June, the size of individual plants gradually increases after June, resulting in insufficient understory space and the competition causing the death of some individuals, which led to a decrease in the number of individuals(Li et al. 2017).

    Herbaceous plant biomass variation

    The total biomass, aboveground biomass and belowground biomass of the thinning treatments were greater than those in the control, indicating that thinning can promote the biomass growth of understory herbaceous plants. The greater the thinning intensity, the faster the plant biomass will grow. This is because thinning can significantly improve the understory light environment,thereby promoting biomass growth of understory herbs(Tinya and ódor 2016). The plant heights of the light thinning and the control showed a peak curve, and the plant heights of the heavy and moderate thinning showed a double-peak curve. This is because July is the rainy season and the temperature is high, and the environmental conditions are moderate. Both overstory and understory entered the vigorous growth period, and the average height of the herbs reached the maximum. In August, as the stand canopy density increase, the understory light of the light thinning and the control was relatively poor, some tall herbs gradually withered and died due to insufficient growth. The light environment of the moderate and the heavy thinning improved significantly after thinning. After the highest canopy closure time in August, understory herbaceous plants have a significant high growth. For example, the heights of Carex pilosa,ferns, and Impatiens noli-tangere, Veratrum nigrum and other plants in the moderate and heavy thinning were significantly higher than those in the light thinning and the control(Sun et al. 2017).

    The diversity of herbaceous plants is an important indicator of the succession and development of forests.This study only studied four-intensity thinning in nearmature forests. In order to better reflect the effects of thinning on the understory environment, herbaceous plant diversity and biomass, therefore, it is necessary to further carry out relevant experiments in forest with different development stages and stand characteristics. It is also a necessity to analyze the impact of thinning on forest species diversity and biomass, guide forest management and accelerate the development of secondary forests to primary forests. In addition, this study also demonstrates that the species diversity and biomass of shrubs will also be affected by thinning. In the future, research can be conducted on the impact of thinning on the species diversity and biomass of shrubs.

    Conclusions

    The changes of biodiversity in the herb layer under the coniferous and broad-leaved mixed forest in Changbai Mountain were mainly affected by seasons. However,thinning changed the light environment in the forest,which could improve the plant diversity and biomass of herb layer in a short time. Moderate thinning had the most significant effect on the number of species, and the number of individual plants, average height and total coverage of plants increased significantly in all months of the growing season. With the increase of thinning intensity, the biomass of herb layer increased accordingly.Total herb biomass, above- and belowground biomass showed positive correlations with Shannon-Winner index, Simpson index and Pielou index in May. The findings provide theoretical basis and reference for reasonable thinning and tending in coniferous and broadleaved mixed forests.

    Acknowledgements

    We would like to thank the two anonymous reviewers and the associate editor for their constructive comments that greatly helped improving the manuscript.

    Authors’contributions

    Gerong Wang designed the study and developed the main hypotheses, Mo Zhou and Mengjia Wu collected and prepared the data, Naiqian Guan,Yuwen Wang, Runhua Jiang and Zhiyu Liu analyzed the data, Fucai Xia and Yue Sun drafted the manuscript. All authors commented preliminary versions of the manuscript and contributed to improve the final version. The author(s) read and approved the final manuscript.

    Funding

    This work was financially supported by a grant from the National Key Research and Development Program of China (2017YFC0504102) and the Science and Technology Program of Jilin Provincial Education Department(JJKH20180349KJ).

    Availability of data and materials

    The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

    Declarations

    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

    1Key Laboratory of State Forestry Administration on Conservation and Efficient Utilization of Precious and Rare Forest Resource in Changbai Mountain, Forestry College,Beihua University,132013 Jilin, China.2West Campus of Hebei Agricultural University, 2596 Lekai South Street, Lianchi District, 071000 Baoding, China.

    Received: 26 December 2020 Accepted: 14 May 2021

    国产在线精品亚洲第一网站| 91国产中文字幕| 午夜视频精品福利| 亚洲真实伦在线观看| 在线观看一区二区三区| 精品国产美女av久久久久小说| 老汉色av国产亚洲站长工具| 亚洲av成人av| 欧美日韩乱码在线| 久久午夜亚洲精品久久| 国产亚洲欧美98| 淫秽高清视频在线观看| 欧美不卡视频在线免费观看 | 丰满人妻熟妇乱又伦精品不卡| 一本精品99久久精品77| 黄色毛片三级朝国网站| 757午夜福利合集在线观看| 国产成人欧美| 日韩免费av在线播放| 久久亚洲精品不卡| 国产av不卡久久| 欧美黑人巨大hd| 人人妻人人澡人人看| 久久草成人影院| 精品久久久久久久久久久久久 | 国产又爽黄色视频| 在线天堂中文资源库| 大型av网站在线播放| 欧美在线一区亚洲| 精品国产乱码久久久久久男人| 欧美日韩亚洲综合一区二区三区_| 久久精品国产亚洲av高清一级| 99在线人妻在线中文字幕| 欧美又色又爽又黄视频| 亚洲自拍偷在线| 女同久久另类99精品国产91| 国产成人欧美| 无遮挡黄片免费观看| 他把我摸到了高潮在线观看| 色综合亚洲欧美另类图片| 999久久久国产精品视频| 一区二区三区高清视频在线| 国产成人欧美| 亚洲精品在线美女| 亚洲国产日韩欧美精品在线观看 | 午夜福利成人在线免费观看| 国产在线精品亚洲第一网站| av电影中文网址| 精品国产超薄肉色丝袜足j| 免费看a级黄色片| 久久99热这里只有精品18| 久久精品人妻少妇| 在线观看66精品国产| 变态另类成人亚洲欧美熟女| 韩国精品一区二区三区| 亚洲精品在线观看二区| 一进一出好大好爽视频| 日韩欧美免费精品| 十八禁人妻一区二区| 免费在线观看亚洲国产| 成人18禁在线播放| 一本综合久久免费| 精品国产亚洲在线| 欧美在线黄色| 亚洲av日韩精品久久久久久密| 午夜免费激情av| 久久久精品国产亚洲av高清涩受| 日韩欧美国产在线观看| 中文字幕精品免费在线观看视频| 国产午夜福利久久久久久| 国产精品影院久久| 亚洲久久久国产精品| 精品国产美女av久久久久小说| 91大片在线观看| 亚洲av成人一区二区三| 国产精品亚洲美女久久久| 久久青草综合色| 搞女人的毛片| 亚洲真实伦在线观看| 国产精品日韩av在线免费观看| 后天国语完整版免费观看| 人妻久久中文字幕网| 高潮久久久久久久久久久不卡| 中文亚洲av片在线观看爽| 黄色女人牲交| 黄色视频,在线免费观看| 观看免费一级毛片| 中文资源天堂在线| 亚洲一区二区三区色噜噜| 亚洲自拍偷在线| 一区二区日韩欧美中文字幕| 男女下面进入的视频免费午夜 | 变态另类成人亚洲欧美熟女| 欧美成人免费av一区二区三区| 亚洲天堂国产精品一区在线| 国产精品久久电影中文字幕| 色婷婷久久久亚洲欧美| 亚洲美女黄片视频| 国产精品av久久久久免费| 欧美色视频一区免费| svipshipincom国产片| 婷婷亚洲欧美| 精品电影一区二区在线| 他把我摸到了高潮在线观看| 精品久久久久久久毛片微露脸| 久久久久国内视频| 欧美黄色淫秽网站| 黑人欧美特级aaaaaa片| 人人妻人人澡欧美一区二区| 精品欧美国产一区二区三| 国产精品美女特级片免费视频播放器 | 精品国内亚洲2022精品成人| www日本在线高清视频| 婷婷精品国产亚洲av| 天天躁夜夜躁狠狠躁躁| 日韩免费av在线播放| 色精品久久人妻99蜜桃| 美女高潮喷水抽搐中文字幕| 少妇裸体淫交视频免费看高清 | 午夜两性在线视频| 此物有八面人人有两片| 久久久久久久午夜电影| 国产成人av教育| 亚洲 欧美一区二区三区| 亚洲三区欧美一区| 亚洲久久久国产精品| 黄片大片在线免费观看| 亚洲精品国产一区二区精华液| 日日摸夜夜添夜夜添小说| 丝袜在线中文字幕| 国产私拍福利视频在线观看| 女性生殖器流出的白浆| 亚洲精品久久国产高清桃花| 久久久国产成人精品二区| 午夜免费鲁丝| 日韩精品中文字幕看吧| 自线自在国产av| 欧美黄色片欧美黄色片| 久久久精品欧美日韩精品| 操出白浆在线播放| 50天的宝宝边吃奶边哭怎么回事| av电影中文网址| 亚洲电影在线观看av| 一区福利在线观看| 日韩欧美免费精品| 国产麻豆成人av免费视频| 又大又爽又粗| ponron亚洲| 巨乳人妻的诱惑在线观看| 在线天堂中文资源库| 母亲3免费完整高清在线观看| 日韩有码中文字幕| 免费一级毛片在线播放高清视频| 黄网站色视频无遮挡免费观看| av欧美777| 男女视频在线观看网站免费 | 夜夜爽天天搞| 很黄的视频免费| 国产成人啪精品午夜网站| 亚洲在线自拍视频| 久久久久免费精品人妻一区二区 | 久久香蕉国产精品| 啪啪无遮挡十八禁网站| 久久草成人影院| 国产伦在线观看视频一区| 999久久久精品免费观看国产| 精品国产一区二区三区四区第35| 黄色女人牲交| 欧美日韩黄片免| 色综合亚洲欧美另类图片| 十八禁网站免费在线| 午夜免费观看网址| 亚洲真实伦在线观看| 国产黄a三级三级三级人| 黄色视频,在线免费观看| 91国产中文字幕| 欧美丝袜亚洲另类 | 午夜福利高清视频| 变态另类丝袜制服| 久久香蕉激情| 一级毛片女人18水好多| 久久精品91蜜桃| 免费电影在线观看免费观看| 中亚洲国语对白在线视频| 中文在线观看免费www的网站 | 亚洲专区中文字幕在线| 超碰成人久久| 黄色 视频免费看| 国产私拍福利视频在线观看| 满18在线观看网站| 中文字幕人妻熟女乱码| 真人做人爱边吃奶动态| 欧美亚洲日本最大视频资源| 国产高清激情床上av| 久热爱精品视频在线9| 9191精品国产免费久久| 一进一出抽搐gif免费好疼| 亚洲激情在线av| 国产精品二区激情视频| 一边摸一边抽搐一进一小说| 国产99久久九九免费精品| 午夜福利在线在线| 欧美zozozo另类| 国产高清videossex| 久久香蕉精品热| 宅男免费午夜| 一级片免费观看大全| 久久人妻av系列| 少妇粗大呻吟视频| 波多野结衣高清作品| 日韩大码丰满熟妇| 老司机福利观看| 老熟妇乱子伦视频在线观看| 国产又色又爽无遮挡免费看| 男女视频在线观看网站免费 | 国产在线观看jvid| 色老头精品视频在线观看| 欧美av亚洲av综合av国产av| 欧美在线一区亚洲| 国产av一区在线观看免费| 午夜激情av网站| 在线免费观看的www视频| 真人一进一出gif抽搐免费| 国产激情欧美一区二区| 亚洲精品中文字幕在线视频| 18禁国产床啪视频网站| 中亚洲国语对白在线视频| 色尼玛亚洲综合影院| 国产精品乱码一区二三区的特点| 国内揄拍国产精品人妻在线 | a级毛片a级免费在线| av视频在线观看入口| 美女高潮到喷水免费观看| 国产色视频综合| 我的亚洲天堂| 激情在线观看视频在线高清| 国产精品香港三级国产av潘金莲| 欧美在线一区亚洲| 熟女少妇亚洲综合色aaa.| 美国免费a级毛片| av中文乱码字幕在线| 亚洲狠狠婷婷综合久久图片| 欧美日韩亚洲国产一区二区在线观看| 欧美激情久久久久久爽电影| 国产亚洲精品第一综合不卡| 国产伦人伦偷精品视频| 亚洲av电影不卡..在线观看| 免费在线观看视频国产中文字幕亚洲| 亚洲国产毛片av蜜桃av| 亚洲精品色激情综合| 丝袜人妻中文字幕| 狂野欧美激情性xxxx| 成人精品一区二区免费| 99re在线观看精品视频| 黄片大片在线免费观看| 麻豆一二三区av精品| 久久国产乱子伦精品免费另类| 免费看美女性在线毛片视频| 亚洲一区高清亚洲精品| 757午夜福利合集在线观看| 欧美av亚洲av综合av国产av| 亚洲九九香蕉| 免费在线观看影片大全网站| 久久婷婷成人综合色麻豆| 男女下面进入的视频免费午夜 | 亚洲九九香蕉| 99精品久久久久人妻精品| xxx96com| 国产亚洲精品av在线| 午夜精品久久久久久毛片777| 亚洲av电影不卡..在线观看| 国产免费男女视频| 男女下面进入的视频免费午夜 | 国内精品久久久久久久电影| 无限看片的www在线观看| 婷婷精品国产亚洲av在线| 亚洲性夜色夜夜综合| 中文字幕另类日韩欧美亚洲嫩草| 18禁黄网站禁片免费观看直播| 制服丝袜大香蕉在线| 色综合欧美亚洲国产小说| 岛国在线观看网站| 国产av在哪里看| 国产亚洲精品久久久久5区| 99精品久久久久人妻精品| 麻豆一二三区av精品| 黄色 视频免费看| 亚洲激情在线av| 精品免费久久久久久久清纯| 在线av久久热| 韩国av一区二区三区四区| 国产爱豆传媒在线观看 | 亚洲精品色激情综合| 51午夜福利影视在线观看| 美女高潮喷水抽搐中文字幕| 欧美激情高清一区二区三区| 成在线人永久免费视频| 久久久精品国产亚洲av高清涩受| 日日夜夜操网爽| 国产野战对白在线观看| 两性夫妻黄色片| 欧美日韩瑟瑟在线播放| 天堂√8在线中文| 久久久久国产一级毛片高清牌| 夜夜看夜夜爽夜夜摸| 久久精品亚洲精品国产色婷小说| 熟女电影av网| 久久精品国产清高在天天线| 夜夜爽天天搞| 丰满人妻熟妇乱又伦精品不卡| 妹子高潮喷水视频| 亚洲狠狠婷婷综合久久图片| 91大片在线观看| 精品久久久久久,| 国产欧美日韩一区二区精品| 色婷婷久久久亚洲欧美| 午夜福利欧美成人| 人人妻人人澡人人看| 亚洲av五月六月丁香网| 日本精品一区二区三区蜜桃| 亚洲一码二码三码区别大吗| 国产三级黄色录像| 一本一本综合久久| 欧美最黄视频在线播放免费| 亚洲国产毛片av蜜桃av| 亚洲国产欧美一区二区综合| 窝窝影院91人妻| 久久欧美精品欧美久久欧美| 免费女性裸体啪啪无遮挡网站| 色播亚洲综合网| 精品乱码久久久久久99久播| 男人的好看免费观看在线视频 | 久久伊人香网站| 一区二区三区激情视频| 久久中文字幕一级| 国产精品永久免费网站| 国产高清视频在线播放一区| 亚洲欧美一区二区三区黑人| 校园春色视频在线观看| 国产爱豆传媒在线观看 | 一边摸一边抽搐一进一小说| 人人澡人人妻人| 中文亚洲av片在线观看爽| 亚洲激情在线av| 一区二区三区精品91| 中文在线观看免费www的网站 | 日本免费一区二区三区高清不卡| 国产高清有码在线观看视频 | 天堂动漫精品| 18禁国产床啪视频网站| 欧美色欧美亚洲另类二区| 欧美不卡视频在线免费观看 | 此物有八面人人有两片| 91国产中文字幕| 99国产精品一区二区三区| 18禁国产床啪视频网站| 精品无人区乱码1区二区| 日本熟妇午夜| 中文字幕人妻丝袜一区二区| 久久久久久亚洲精品国产蜜桃av| 午夜两性在线视频| 黄色视频,在线免费观看| 两个人看的免费小视频| 国产又爽黄色视频| 很黄的视频免费| 色哟哟哟哟哟哟| 男女之事视频高清在线观看| 最新美女视频免费是黄的| 美女免费视频网站| av福利片在线| 国产精品1区2区在线观看.| 国产三级黄色录像| 欧美中文综合在线视频| 亚洲成人久久性| 欧美日韩乱码在线| 久久伊人香网站| 午夜视频精品福利| 亚洲七黄色美女视频| 精品国内亚洲2022精品成人| 91成人精品电影| 久久国产精品男人的天堂亚洲| 村上凉子中文字幕在线| 麻豆av在线久日| 日日摸夜夜添夜夜添小说| 精华霜和精华液先用哪个| 老鸭窝网址在线观看| 欧美乱色亚洲激情| 日本精品一区二区三区蜜桃| 最新在线观看一区二区三区| 国产亚洲av高清不卡| 黄色视频,在线免费观看| 一级作爱视频免费观看| 天堂√8在线中文| 18禁黄网站禁片午夜丰满| 搞女人的毛片| 免费电影在线观看免费观看| 国产亚洲欧美98| 国产高清视频在线播放一区| 女生性感内裤真人,穿戴方法视频| 91麻豆av在线| 国产精品综合久久久久久久免费| 国产成人精品无人区| 欧美一级a爱片免费观看看 | 日韩精品中文字幕看吧| 亚洲激情在线av| bbb黄色大片| 精品国产亚洲在线| 免费观看人在逋| 成人国语在线视频| 欧美成人免费av一区二区三区| 午夜两性在线视频| 免费人成视频x8x8入口观看| 黑人欧美特级aaaaaa片| 国产激情欧美一区二区| 久热这里只有精品99| 男人操女人黄网站| 久久久久国产精品人妻aⅴ院| 宅男免费午夜| 伦理电影免费视频| avwww免费| 国产日本99.免费观看| 精品国产亚洲在线| 51午夜福利影视在线观看| 国产精品野战在线观看| 国产一区在线观看成人免费| 黄色片一级片一级黄色片| 侵犯人妻中文字幕一二三四区| 久久香蕉国产精品| 夜夜躁狠狠躁天天躁| 亚洲第一青青草原| 成人18禁在线播放| 男人舔女人下体高潮全视频| 国产97色在线日韩免费| 亚洲全国av大片| 高清毛片免费观看视频网站| 欧美激情久久久久久爽电影| 亚洲欧美一区二区三区黑人| 香蕉久久夜色| 女性生殖器流出的白浆| 99精品在免费线老司机午夜| 欧美黄色片欧美黄色片| 久久午夜综合久久蜜桃| 村上凉子中文字幕在线| 国产国语露脸激情在线看| 久久精品91蜜桃| 性色av乱码一区二区三区2| 日韩视频一区二区在线观看| 人妻久久中文字幕网| 欧美成人午夜精品| 91成年电影在线观看| 两个人看的免费小视频| 久久99热这里只有精品18| 国产精品亚洲av一区麻豆| 免费在线观看日本一区| 窝窝影院91人妻| 黄色女人牲交| www国产在线视频色| 97碰自拍视频| 这个男人来自地球电影免费观看| 久久精品国产综合久久久| 亚洲,欧美精品.| 老司机福利观看| av中文乱码字幕在线| 国产亚洲av嫩草精品影院| 一区二区三区激情视频| 亚洲最大成人中文| 亚洲国产欧洲综合997久久, | 午夜福利18| 成年女人毛片免费观看观看9| 国产精品久久久av美女十八| 狂野欧美激情性xxxx| 波多野结衣av一区二区av| 成人国语在线视频| 97人妻精品一区二区三区麻豆 | 天天一区二区日本电影三级| 老汉色av国产亚洲站长工具| 黑人欧美特级aaaaaa片| 国产精品一区二区三区四区久久 | 韩国精品一区二区三区| 久久午夜亚洲精品久久| 国产亚洲精品久久久久久毛片| 国产精品 国内视频| 超碰成人久久| or卡值多少钱| 国产成人av激情在线播放| 亚洲国产毛片av蜜桃av| 麻豆国产av国片精品| 国产精品1区2区在线观看.| 欧美性猛交╳xxx乱大交人| 日日夜夜操网爽| 波多野结衣高清作品| 50天的宝宝边吃奶边哭怎么回事| 色综合亚洲欧美另类图片| 国产欧美日韩一区二区精品| 桃红色精品国产亚洲av| 欧美精品啪啪一区二区三区| 免费女性裸体啪啪无遮挡网站| 热re99久久国产66热| 亚洲国产欧洲综合997久久, | 人人妻人人澡欧美一区二区| 久久久久久亚洲精品国产蜜桃av| 亚洲 国产 在线| 女警被强在线播放| 999久久久国产精品视频| 女警被强在线播放| 韩国精品一区二区三区| 非洲黑人性xxxx精品又粗又长| 天天躁夜夜躁狠狠躁躁| 在线观看免费午夜福利视频| 此物有八面人人有两片| 无人区码免费观看不卡| 成人av一区二区三区在线看| 国产麻豆成人av免费视频| 人人妻人人澡欧美一区二区| 中出人妻视频一区二区| 国产高清激情床上av| 精品午夜福利视频在线观看一区| 国产在线观看jvid| 亚洲一区高清亚洲精品| 中文字幕人成人乱码亚洲影| 国产激情久久老熟女| 国产91精品成人一区二区三区| 精品第一国产精品| 亚洲第一青青草原| 最近最新中文字幕大全电影3 | 国产精品久久久久久人妻精品电影| 成人永久免费在线观看视频| 欧美日韩中文字幕国产精品一区二区三区| 一级毛片女人18水好多| 亚洲熟妇熟女久久| 国内精品久久久久精免费| 国产三级黄色录像| 在线国产一区二区在线| 制服丝袜大香蕉在线| 搞女人的毛片| 国产成人精品久久二区二区91| 亚洲午夜精品一区,二区,三区| 最近最新免费中文字幕在线| 成人特级黄色片久久久久久久| 欧美色欧美亚洲另类二区| 国产私拍福利视频在线观看| 国产亚洲精品一区二区www| 日韩大码丰满熟妇| 午夜影院日韩av| 国产精品一区二区免费欧美| 欧美日韩乱码在线| 久久婷婷成人综合色麻豆| 欧美中文日本在线观看视频| 亚洲熟妇中文字幕五十中出| 不卡av一区二区三区| 两性夫妻黄色片| 精华霜和精华液先用哪个| 国产99久久九九免费精品| 欧美在线一区亚洲| 韩国精品一区二区三区| 日韩欧美国产在线观看| 免费观看人在逋| 老汉色∧v一级毛片| 午夜福利成人在线免费观看| 久久久久久久午夜电影| 女同久久另类99精品国产91| 最新在线观看一区二区三区| 黄片小视频在线播放| 欧美 亚洲 国产 日韩一| 国产熟女午夜一区二区三区| 日本 av在线| 波多野结衣巨乳人妻| 这个男人来自地球电影免费观看| 色在线成人网| 99精品在免费线老司机午夜| av超薄肉色丝袜交足视频| 婷婷精品国产亚洲av在线| xxx96com| 大型av网站在线播放| 韩国精品一区二区三区| 久久国产精品人妻蜜桃| 久久精品国产亚洲av高清一级| 国内少妇人妻偷人精品xxx网站 | 狂野欧美激情性xxxx| 国产黄a三级三级三级人| 亚洲欧美精品综合一区二区三区| 亚洲成a人片在线一区二区| 中文字幕人妻丝袜一区二区| 少妇熟女aⅴ在线视频| 女警被强在线播放| 亚洲成人国产一区在线观看| 不卡av一区二区三区| 国产精品美女特级片免费视频播放器 | 国产免费男女视频| 在线观看日韩欧美| 亚洲激情在线av| 日本熟妇午夜| 三级毛片av免费| 日韩欧美一区视频在线观看| 久久热在线av| 成人亚洲精品一区在线观看| 欧美性长视频在线观看| 中亚洲国语对白在线视频| 国产一区二区三区视频了| 成年免费大片在线观看| www日本黄色视频网| 中文在线观看免费www的网站 | 国产精品美女特级片免费视频播放器 | 久热爱精品视频在线9| 最新在线观看一区二区三区| 99在线人妻在线中文字幕| 久久国产精品影院| 亚洲第一欧美日韩一区二区三区| 色综合站精品国产| 精品一区二区三区视频在线观看免费| 亚洲avbb在线观看| 叶爱在线成人免费视频播放| 老熟妇仑乱视频hdxx| 丰满人妻熟妇乱又伦精品不卡| 国产伦人伦偷精品视频|