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

    Tree mortality and regeneration of Euphrates poplar riparian forests along the Tarim River,Northwest China

    2021-02-28 09:11:18AyjamalKeramHalikTayierjiangAishanMaierdangKeyimuKadeliyaJiapaerandGuoleiLi
    Forest Ecosystems 2021年4期

    Ayjamal Keram,ümüt Halik ,Tayierjiang Aishan,Maierdang Keyimu,Kadeliya Jiapaer and Guolei Li

    Abstract

    Keywords:Tree mortality,Regeneration strategy,Seedling and sapling recruitment,Gap makers,Riparian forest,Tarim River

    Highlights

    · Effect of hydrological alterations on tree mortality of Populus euphratica riparian forest was quantified.

    · The contribution of different habitats (canopy gap area, undercanopy area and riverbank area) on plant regenerations was evaluated.

    · Recommendations for regeneration and succession of riparian forests at the Tarim River were put forward.

    Introduction

    Riparian forests provide fundamental ecosystem services in arid regions, including maintenance of ecological stability and prevention of natural disasters such as sandstorms, heatwaves, and desertification (Song et al. 2000;Chen et al. 2013; Betz et al. 2015; Mamat et al. 2018;Halik et al. 2019). The existence of trees, bushes, and grass vegetation on the banks of the Tarim River, which constitute the natural barriers of river oases, is absolutely dependent on the river (Zeng et al. 2002, 2006;Halik et al. 2006; Aishan et al. 2018). In recent decades,floodplain vegetation has been threatened by increased water scarcity, and large areas of riparian forests including seedlings and saplings, have withered (Gries et al.2003; Foetzki 2004; Aishan 2016; Zeng et al. 2020). With China’s rapid economic development and the implementation of the “Ecological Civilization” and “One Belt One Road” initiatives, the ecological restoration of the floodplain ecosystem along the Tarim River has been designated as a high priority by the Chinese government(Halik et al. 2006, 2019). Conservation of the remaining riparian forests and recovery of the degraded parts of the ecosystem are crucial for further sustainable development of the region (Rumbaur et al. 2015; Deng 2016;China Green Foundation 2018), particularly, transport infrastructure, including railways, highways, and oil and gas pipelines in the Tarim Basin are in need to be protected (Deng 2009, 2016; Aishan 2016; Halik et al. 2019).Thus, it is necessary to preserve the functions and services of these forests through adaptive management.

    The riparian forest (also known as Tugai forests) at the Tarim River constitute a natural green belt at the northern edge of the Taklimakan Desert. The Euphrates poplar(Populus euphratica Oliv.)is a dominant tree species in the floodplain ecosystems. As a “Green Corridor”,P. euphratica forests have become increasingly important for preventing the unification of two neighboring sandy deserts, the Taklimakan and the Kuruktagh (Xu et al. 2006; Betz et al. 2015). These forests have important ecological functions in addition to socio-economic and touristic value, such as protection of biodiversity,regulation of the climate and hydrologic conditions in oases, fertilization of soils, and maintenance of regional ecosystem balance (Wang et al. 1995; Huang 2002;Huang and Pang 2011; Mamat et al. 2019). Nevertheless,the Tarim floodplain ecosystem, especially the riparian vegetation, has been seriously degraded, as the water conveyance in the main channel has rapidly decreased and many tributaries of the watershed have been disconnected from the main river (Zhao et al. 2011; Chen et al.2013).

    Many studies have demonstrated that water shortages in the Tarim are mainly caused by climate change coupled with intensive anthropogenic activities, i.e., tree felling, overgrazing, excessive land reclamation and disorderly expansion of cotton monoculture. Therefore,restoring degraded floodplain ecosystems has become a major focus of applied forest research (Deng 2009; Rumbaur et al. 2015; Thomas et al. 2017; Thomas and Lang 2020). Among all processes of natural forest dynamics,the formation of forest canopy gaps was considered a vital regeneration strategy for maintaining floodplain forest structures (Han et al. 2011; Keram et al. 2019).Gap disturbance is the main driving force of forest dynamics as it creates environmental heterogeneity (Zhu et al. 2003, 2007; Albanesi et al. 2005; Bottero et al.2011). Forest gaps thus create important habitat for the regeneration of plant species, which may otherwise be suppressed by the undercanopy (Han et al. 2011; Sharma et al. 2018). They also play a vital role in forest regeneration and succession, especially in the establishment and development of plant species that differ in ecological recruitment (Runkle 1998; Mountford 2006; Rentch et al.2010; Han et al. 2013; Zhu et al. 2014; Jankovska et al.2015).Nagel et al.(2010)reported a high presence of advanced regeneration in a gap regeneration study of a(mixed) beech virgin forests. Besides, it also has been reported that canopy gaps close when the height of regenerations reached 20 m, consistent with the definition used by Nagel and Svoboda (2008). Zhu et al. (2021)found that with the decreased of gap size, pioneer species became the sub-canopy layer, and with the aging of gap, the light conditions changed over time, which was conductive to the recruitment of shade-tolerant species;Keram et al. (2019) also revealed that hydrological conditions (groundwater, runoff and water consumption)are the main driving force of the gap-scale disturbance of desert riparian forests along the Tarim River. In floodplain forests, canopy gaps may not be filled with regenerations within a short period because of high tree mortality. This leads to the continual expansion of canopy gaps (Keram et al. 2019). The majority of earlier studies have focused on the mortality of gap makers and the diversity of regeneration species in other various climatic zones, such as tropical, subtropical, north temperate, and cold temperate regions (Runkle and Yetter 1987; Yamamoto 2000; Dorot? and Thomast 2008;Petritan et al. 2013; Popa et al. 2017; Zhu et al. 2018;Kitao et al. 2018). However, systematic research on tree mortality and plant regeneration of desert riparian forests under various hydrological scenarios are relatively lacking. Therefore, it is necessary to scientifically understand the response of P. euphratica mortality rate to hydrological dynamics at regional scales, and seedling and sapling regeneration in the Tarim River Basin.

    Forest canopy gaps may provide opportunities for species regeneration and have therefore been widely exploited in forest recovery programs (Coates and Philip 1997; Schliemann and Bockheim 2011; Kern et al. 2013,2017; Nagel et al. 2016; Lu et al. 2018). Thus, it is necessary to conduct studies to better understand the mortality of gap makers under long-term hydrological processes and seedling and sapling establishment in different habitats, such as canopy gap areas (CGAs), undercanopy areas (UCAs), and riverbank areas (RBAs) of floodplain, as well as to clarify the contribution of canopy gaps to plant regeneration in riparian forests. Such information may provide insights into the efficiency of gap usage and seedling planting in forest restoration activities along the Tarim River. In the present study, we used comprehensive field investigation data to describe tree mortality as well as to compare the density of seedlings and saplings between the CGAs, UCAs, and RBAs to evaluate whether the CGAs are important for the regeneration of P. euphratica riparian forests along the Tarim River.In this work we made three hypotheses:

    (H1)

    Based on the observation that canopy gap disturbances have frequently emerged in the desert riparian forests and its composition and structure were drastically different from other forest types (Han et al. 2011; Keram et al.2019). We hypothesized that characteristics of tree mortality (number, DBH) strongly responded to hydrological dynamics (such as stream flow) at local and regional scales in the Tarim River Basin, resulting in the emergence of special forest gap structures.

    (H2)

    Conclusions on the size development of canopy gaps,once formed they would continuously expand (Keram et al. 2019). The important process of gap dynamics,i.e. plant regeneration, in floodplain forests has not been quantified yet. To test the general hypothesis,that canopy gaps promote plant regeneration in riparian forests, we formulated the hypothesis that canopy gaps drive tree regeneration to some degree, due to canopy gap provides more adequate light conditions than the UCAs.

    (H3)

    In floodplain forests, due to the degradation of plant habitat, canopy gaps may not be filled up with other regenerations in a relatively short period. By comparison of the species diversity between gap makers and gap fillers in canopy gaps, we hypothesized that canopy gaps may not be filled up with other regenerations within a short time, (i) which may be due to the degraded habitat of plant growth, (ii) prolonged and unexpected flooding disturbance resulted in the degradation of P. euphratica trees, and (iii) its mortality did not cause significant changes in tree distribution but caused changes at tree level.

    We intend to verification of these hypotheses by compiling data from six permanent monitoring plots located at undisturbed natural forest sites along the middle reaches of the Tarim River. We believe that addressing these hypotheses will contribute to improve the understanding of the current state of the Tarim riparian ecosystem, where P. euphratica is now in urgent need of protection. Additionally, it will provide an ecological framework for the design and implementation of closeto-nature restoration techniques to extend P. euphratica forests.

    Methods

    Study site

    Field work was conducted in Yingbazha village (41.22°N,84.31°E, ASL 1000 m) in Tarim Huyanglin Nature Reserve (THNR) along the middle reaches of the Tarim River (Fig. 1), which was established in 1984 and upgraded to a National Nature Reserve in 2006. The region is an extremely arid warm temperate zone, with an annual mean precipitation of 75 mm and an annual mean temperature of 11.05°C (Keyimu et al. 2018). The annual potential evapotranspiration ranges from 2500 to 3500 mm (Rouzi et al. 2018). According to the USDA(United States Department of Agriculture) soil classification system, the soil of the Tarim River is a member of the Aridisol order, and the soil is silty loam (Hu et al.2009; Yu et al. 2009). The local groundwater system is recharged from the surface water through bank infiltration(Huang et al. 2010). The forest structure and composition in the study area is relatively simple, and the local flora is mainly composed of P. euphratica, Tamarix ramosissima and Phragmites australis (Table 1). However, a few rare plant species are also present, including Populus pruinosa,Tamarix hispida,Tamarix leptostachys,Glycyrrhiza glabra,Inula salsoloides, Karelinia caspica, Lycium ruthenicum,Alhagi sparsifolia, Apocynum venetum, Halimodendron halodendron, Poacynum hendersonii, Cirsium segetum,Cynanchum auriculatum, Aeluropus pungens, Taraxacum mongolicum, Salsola collina, Calamagrostis pseudophragmites, Myricaria platyphylla, Halogeton glomeratus, and Elaeagnus oxycarpa (Gries et al. 2003; Chen et al. 2004;Thevs 2006; Thevs et al. 2008; Huang et al. 2010; Lang et al.2016;Halik et al.2019).

    Fig.1 Sketch map of the study area

    Field investigation

    A complete field survey was carried out from May 15 to June 20, 2017. A total of 60 canopy gaps were investigated in six replicate grid plots (50 m×50 m) (Table 1).Canopy gaps were created by tree mortality (i.e., death)or by >50% loss of tree branches. The ratio of gap diameter to tree height on the border of the gap (RD/Hvalue)ranged from 0.4 to 2.2 (Keram et al. 2019), which is consistent with measurements reported by Zhu et al. (2015).For each canopy gap, all surrounding trees and gap makers were identified, and the species, diameter at breast height (DBH), and tree height (TH) or log length were recorded for each tree. Morphological characteristics were used to determine the degree of tree decomposition and the life expectancy of gap makers was estimated under the guidance of local experienced forest workers (Liu and Hytteborn 1991; Droessler and Von 2005; Thomas and Jurij 2006; Petritan et al. 2013; Wen 2016; Yang et al. 2017; Keram et al. 2019), and the relative ages of the gap makers were thus determined. In addition, we determined whether there were significant differences in the DBH distribution among living and dead trees. We regarded trees bordering canopy gaps as living trees and gap makers as dead trees. To investigate regeneration species, we selected five replicate plots in each habitat (CGAs, UCAs, and RBAs) and recorded the species and numbers of seedlings and saplings. We simultaneously compared the species diversity and densityin these plots. We defined a CGA as an area where the canopy gap >10 m2. For the UCA, we chose plots around the canopy gaps in the study area; RBA plots were located near the riverbank, which were covered with excessive water. For each habitat, 15 repeated sampling plots (10 m×10 m) were set, giving a total of 45 sampling plots. In each plot, we counted the number of species, for seedlings (≤10 cm tall) and saplings (10–100 cm tall) (Walker 2000).

    Table 1 Characteristics of six natural stands sampled for analysis of tree mortality

    Data processing

    Classifying the DBH structure of gap makers

    Over the entire area of canopy gaps, gap makers were classified into five categories according to their DBH and TH (Xu et al. 2016): Class I comprised juvenile trees(DBH <5 cm, TH <4 m), Class II comprised young trees(5 cm <DBH <10 cm, TH ≥4 m), Class III comprised young/middle-aged trees (10 cm <DBH <15 cm), Class IV comprised middle-aged trees (15 cm <DBH <20 cm),and Class V comprised old trees (DBH >20 cm).

    Dynamic index of gap maker structure

    The dynamic structure of gap makers was quantitatively analyzed. The dynamic index of population structure(VPi) can be represented by the change trend of gap makers in various structures. If 0 <VPi≤1, the quantity of gap makers increased with the VPivalue; If -1 <VPi≤0, the quantity of gap makers decreased with the VPivalue. The dynamic index of gap maker individuals (Vn)without external disturbance is described by the following equation:

    where Snis the number of individuals in n DBH class,Sn+1is the number of individuals in n+1 DBH class of gap makers, and (…)maxis the maximum value in the column. Below is the formula used to calculate the dynamic index of gap maker structure (VPi):

    where Snand Vnare the same as in Eq. (1), K is the quantity of gap makers in different age classes, and the range of VPiis consistent with Vn. We used the formula proposed by Chen (1998) to calculate the dynamic index value under completely random (Vpi’) and non-random(Vpi′′) disturbances.

    Statistical analysis

    We have calculated species diversity according to the Margalef richness index and Shannon Wiener index (Ma 1994; Ma and Liu 1994). A fitting optimization t-test method was used to statistically analyze the change trend of gap makers that died in different years and to evaluate the distribution patterns of the DBH values of gap makers. We used a chi-square goodness-of-fit test to reveal the DBH distribution of dead trees versus living trees. Further, we comparatively analyzed plant regeneration density between the different habitats (CGAs,UCAs, and RBAs) using the Margalef and Shannon-Wiener indices. Values of P <0.05 were considered to be statistically significant. The data were organized using Excel 2015 (Microsoft, Redmond, USA) and analyzed using OriginPro 2016(OriginLab Corp2016).

    Results

    Tree mortality

    Among the 245 gap makers identified, the mean DBH was 20.80±9.03 cm (range 10.10–49.50 cm) and there were significant variations in different years. The chisquare goodness-of-fit test showed that the DBH structure of gap makers presented a common curve distribution (Fig. 2). The DBH distribution of gap makers approximated a Gaussian fit (R2=0.762, 0.749, 0.829,0.896, 0.918, and 0.364, respectively; P <0.05). However,large young gap makers (DBH <15 cm) were more common in the study area, especially in the 2007–2016(Gaussian fit, R2=0.364, P <0.05). During the investigation, we found that some young trees were injured or crushed by old trees.

    In this study, we found that stream water was significantly negatively correlated with tree mortality rates over the last three decades. Based on the above data analysis,there was an ascending trend in the dynamic index of gap maker’s population structure during the periods of 1977–1986, 1987–1996, 1997–2006, and 2007–2016,while in 1957–1966 and 1967–1976 there was a descending trend (Table 2). Furthermore, there was a considerable negative interaction between the dynamic index value of gap maker population structure and annual runoff; the equation of the fitted model was y=-0.0215x+0.6146 (R2=0.912, P <0.01) without external disturbances (Vpi), while the equation of the fitted model was y=-0.0005x+0.0146 (R2=0.8405,P <0.01) under completely random disturbances (Vpi′)and y=0.0032x+0.112 (R2=0.356, P <0.01) under non-random disturbances (Vpi′′). Vpiand Vpi′ decreased markedly with increase of stream runoff in the study area(Fig. 3). The Vpi′ value became positive and tended to increase when the river dropped to 25.86×108m3.

    Comparison of the DBH distribution of dead trees and living trees showed significant variation with a high proportion of young dead trees (Fig. 4). As mentioned above, the trees surrounding CGAs and gap makers represented living trees and dead trees,respectively. Grid chart analysis showed that the median DBH of dead trees, especially those died before 1996, were larger than that of living trees. However, by comparing DBH, we found that the trees that died in the 1997–2016 were smaller than living trees.

    Fig.2 Diameter class distribution of P.euphratica gap makers in the study area

    Forest regeneration

    Riparian vegetation was distributed on both sides of the riverbank and was relatively simple and sparse in structure. In our field investigation, we counted 23 seedling and sapling species belonging to 12 families and 20 genera. In general, there were obvious differences in seedling and sapling regeneration among the habitats(CGAs,UCAs, and RBAs)(Table 3).Among the 14 plant species found in the CGAs, 7 were also found in the UCAs and 19 were also found in the RBAs; P. euphratica, T. ramosissima, T. hispida, T. leptostachys, G. uralensis and E. oxycarpa commonly appeared in all three habitats. P. pruinosa, P. hendersonii, C. segetum, C. auriculatum, A. pungens, T. mongolicum, M. platyphylla and L. ruthenicum were rare species, which only appeared in the RBAs. P. australis, H. halodendron and A. sparsifolia grew in the CGAs, but not in the UCAs or RBAs.Additionally, the RBAs contained a high number of P.euphratica species (30) compared to other habitats, and the number of P. euphratica and other species in the CGAs was higher than in the UCAs.

    The seedling and sapling density of P. euphratica was significantly higher in the RBAs than in the CGAs and UCAs. However, compared with the UCAs, the CGAs promote tree regeneration to some extent by providing favorable conditions for the survival and growth of seedling and saplings (Table 4). We also performed quantitative analysis of seedling and sapling regeneration, based on species diversity and richness (Fig. 5). According toanalysis of the Margalef and Shannon-Wiener indices,the species diversity of seedlings and saplings was remarkably higher in the RBAs than in the other two areas,and higher in the CGAs than in the UCAs(P <0.05).Further,though the density and diversity of seedlings and saplings in the CGAs was not as high as in the RBAs,the survival rate of young trees was higher in the CGAs than in the RBAs(Table 4).

    Table 2 Analysis of population structure dynamics of gap makers under different volumes of runoff water

    Fig.3 Relationship between population structure dynamics of gap makers and annual runoff

    Comparison of gap makers and gap fillers

    The statistical analysis showed that P. euphratica is the most common species of gap maker. However, it also existed as gap filler, along with T. ramosissima and P.australis (Fig. 6). Other gap fillers included P. pruinosa,A. venetum, H. halodendron, P. hendersonii, C. segetum,C. auriculatum, A. pungens, T. mongolicum, S. collina,C. pseudophragmites, M. platyphylla, H. glomeratu and E. oxycarpa, as well as some herbs such as G. glabra, I.salsoloides, K. caspica, L. ruthenicum, and A. sparsifolia,although they were very rare.

    Discussion

    Fig.4 Boxplots of the DBH of dead trees and living trees on the grid charts in canopy gaps that were created at different times.The boxencloses the middle 80%of observations. Median and mean values are indicated by a vertical line(-) and a star(★),respectively.The upper whisker indicates the largest DBH value,while the lower whisker shows the smallest.The plus sign(+)indicates the maximum DBH

    Table 3 Seedling and sapling recruitment in different habitats in the study area

    One or more gap makers, i.e., dead trees, create canopy gaps in forests (Zang et al.1999). Keram et al. (2019)found that the mortality of gap makers did not substantially contribute to increased gap size in floodplain forests along the Tarim River. This could be explained by the small crown and DBH of gap makers. Tree mortality is caused by different agents, including drought or flooding, strong wind, diseases, and insect pests, and affects forest dynamics (Krasny and Whitmore 1992; Battles and Fahey 2000). Host-specific disturbance was considered to be the most common agent of tree mortality in old-growth mixed beech forests in the Western Carpathians (Orman and Dorota 2017), while wind disturbance led to a disproportionate number of gap makers in old-growth red spruce–Fraser fir forests in North Carolina (White et al. 1985) and red spruce–balsam stands in New Hampshire (Foster and Reiners 1986;Worrall et al. 2005). However, the causes of tree mortality in riparian forests are different from those in other forest types, such as tropical and subtropical forests. The results of the present study showed that tree mortality in riparian forests along the Tarim River was mainly influenced by hydrological factors (surface water), although strong winds or insect pests may also play a role (Han et al. 2011).

    Water is the most important limiting factor for the maintenance and growth of P. euphratica in arid ecosystems, and P. euphratica depends on ground water that available to its roots (Yu et al. 2019). It is known that water is the factor with most influence on the survival,distribution, and development of vegetation in hyperarid regions. In arid ecosystems,the vitality of P. euphratica forests strongly depend on the groundwater level(Halik et al. 2006, 2009; Wang et al. 2008; Ginau et al.2013; Keyimu et al. 2017). Halik et al. (2019) reported that the vitality of P. euphratica diminished when the groundwater level was not sufficient for the development of tree stands. In addition, Xu et al. (2016) studied the correlation between groundwater and surface water using stable isotope technology. Their results indicated that the groundwater level is significantly affected by surface water. Chen et al. (2004) and Keilholz et al.(2015) also demonstrated that groundwater is often supplied by stream flow and infiltration from watercourses.The groundwater depth in the Tarim Basin has dropped continuously as a result of the desiccation of streams and decreases in stream flow since an embankment was built in 2000. Thus, the physical characteristics of plant habitats have been changed. Consequently, the growth of P. euphratica cannot be sustained as before in this region, and forests have highly degraded or died. According to the present study, hydrological alterations may be the main cause of tree death in desert riparian forests along the Tarim River. DBH is one of the most common forest inventory variables and can reflect the population structure and growth status of tree stands (Condit et al.2000; Di et al. 2014). Due to the unique climate and habitat conditions of desert floodplain areas, tree shape variables are obviously different from those in temperate or/and subtropical forests (Ling et al. 2015). In the present research, it was found that most young gap makers (DBH <15 cm) died in the most recent two decades, i.e., 1997–2006 and 2007–2016. In the early 1970s, water resources were affected on a regional scale by climate change and anthropogenic activities (increased water use for irrigation), and then lag effects led to a progressive decrease in surface water from 1984 and caused a reduction in groundwater from 2004 (Keram et al. 2019). This resulted in habitat deterioration and adecline in the survival of young P. euphratica individuals. Based on quantitative correlation analysis, the present study revealed that gap maker mortality and DBH structure responded negatively to variations in surface water. This may be due to long-lasting exogenous disturbances, such as a scarcity of surface water or a decrease in groundwater level, having a notable influence on vegetation through drought stress. Thus, the water requirement of riparian vegetation was not satisfied.Therefore, gap makers of low DBH appeared frequently in P. euphratica forests, which experienced an increase in more resilient species. This may be because the vitality of P. euphratica in this arid area was directly affected by the hydrological conditions, and only a small proportion of young trees survived to reach canopy height.However, there were a large number of larger, mature trees among the living trees that have developed root systems and are therefore not significantly affected by the reduced groundwater level.

    Table 4 Comparison of P.euphratica seedlings and saplings in different habitats in the middle reaches of the Tarim River

    Fig.5 Species diversity of seedlings and saplings in different habitats

    Fig.6 Comparison of gap maker and gap filler species in the middle reaches of the Tarim River

    By creating environmental heterogeneity in terms of light availability, gap disturbances play a key role in forest regeneration as well as in the establishment and development of tree species with different ecological recruitment patterns (Runkle 1989; Peterken 1996;Mountford 2006; Zhang et al. 2006; Li et al. 2013). Gapbased restoration provides a flexible system that emulates natural regeneration under different environments. Van and Dignan (2007) reported that seedling and sapling diversity and richness were promoted by light intensity. Zhu et al. (2014) stated in a review of published studies that compared to UCAs, the seedlings and saplings of shade-tolerant species have significantly higher density in CGAs. Overall, water, light availability,and seedbed conditions potentially affect seedling emergence, germination, and establishment in floodplain riparian forests. It was found that appearance of seedling is closely related to flooding patterns in this study. Besides, a higher density of seedlings and saplings is expected in gaps as they provide full light conditions. The results of the present observations showed that the seedling and sapling species diversity was higher in the CGAs than in the UCAs. Hence, forest canopy gaps playing a pivotal role in the long-term germination and regeneration of plant species in some degree. Oliver and Larson(1996) stated that canopy gaps simultaneously result in mortality in some individuals and establishment and growth in others, and the ongoing process of death and replacement has a profound effect on forest structure and composition. Pham et al. (2004) showed that Abies balsamea is the most frequent successor in Abies forest stands and that Picea mariana is the most likely replacement species in Picea forest stands, regardless of the species of gap makers. In the present study, canopy gaps were filled with shrubs and herbs. When we examined the probability of replacement in floodplain forests,there seemed to be a reciprocal replacement of P.euphratica by T. ramosissima. and P. australis (Fig. 6). However, the probability of P. euphratica replacement was lower when the mortality of young trees was particularly high.

    Conclusion and suggestion

    As the largest riparian forest ecosystem in Northwest China, Euphrates poplar forests are facing great challenges of climate change and human interference (Zhang et al. 2016; Halik et al. 2019). Monitoring of tree mortality and plant regeneration in different habitats is considered the most effective way of ensuring the sustainable regeneration and ecological restoration of floodplain forests along the Tarim River. Based on evaluation of the correlation between tree mortality and water resources in the study area, hydrological conditions (runoff volume) are the main reason for the death of P. euphratica forests. At the beginning of the 1970s, hydrological conditions were affected on a regional scale by anthropogenic activities. Consequently, riparian forests were threatened under water scarcity, which resulted in high mortality among Populus species. Our observations of seedling and sapling density indicated that the frequency of plant species was different among the three studied habitats; it was highest in the RBAs, intermediate in the CGAs, and lowest in the UCAs. It may have been affected by light availability and hydrological conditions during seedling and sapling recruitment.

    In conclusion, searching the tree mortality and plant regeneration are considered to be the most effective ways to clarify the development trend of desert riparian forest. Furthermore, the protection and maintenance of the floodplain vegetation along the Tarim River are crucial to sustainable regional development in the Tarim Basin. Our study area is located in the Yingbazha section, where there are fragment pits of stagnant water by flooding. Moreover, though the density and diversity of seedlings and saplings in the CGAs was not as high as in the RBAs, the survival rate of young trees was higher in the CGAs than in the RBAs. Accordingly, we suggest that making full use of limited resources of stagnant water pits to establish a “Natural Nursery of the Poplar”and appropriately transplanted into the gap, the main purpose of which is to achieve an improved regeneration efficiency of the Populus euphratica riparian forests.Additionally, there is a need for further applied experiment to establish test areas for flood water allocation and utilization strategies, which, if successful, should then be extended to the whole basin. Thus, we argue that, to some degree, canopy gaps can provide a favorable environment for seedlings to grow into saplings,thereby reaffirming the importance of gap creation for forest regeneration.

    Abbreviations

    DRB: Distance from river bank; GWL: Groundwater level; DBH: Diameter at the breast height; TH: Tree height; CGAs: Canopy gap areas;UCAs: Undercanopy areas; RBAs: Riverbank areas; THNR: Tarim Huyanglin nature reserve

    Acknowledgments

    Cordial thanks to MSc-Students (Elyar, Wang Wenjuan, Ramile, Mihray and Gao Qing) from the Xinjiang University for their intensive support during the field investigation. Special thanks to Xinjiang LIDAR Applied Engineering Technology Research Center for giving us free use of the Riegl VZ-1000 Terrestrial Laser Scanner. We thank the anonymous reviewers for their constructive comments that greatly helped us improve the quality of this manuscript.

    Authors’contributions

    Ayjamal Keram conceived, designed and performed the experiments,contributed the data collection, analyzed the data, prepared figures and tables, and authored or approved the final draft of the manuscript. ümüt Halik conceived and designed the experiments, contributed the data collection, reviewed the manuscript, and approved the final draft. Tayierjiang Aisahan, Maierdang Keyimu and Kadeliya Jiapaer participated in field work,contributed the data collection/validation and language proofreading. Guolei Li edited and validated the manuscript with critical comments and reviewed the manuscript draft.All authors checked and approved the final manuscript.

    Funding

    The work was funded by the National Natural Science Foundation of China(31860134, U1703102, 31700386).

    Availability of data and materials

    The data set generated for the study area is 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

    1College of Resources and Environmental Science, Key Laboratory for Oasis Ecology of Ministry of Education,Xinjiang University, ürümqi 830046, China.2College of Forestry, Key Laboratory for Silviculture and Conservation of Ministry of Education, Beijing Forestry University, Beijing 100083, China.3State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085,China.

    Received: 26 February 2021 Accepted: 9 June 2021

    伊人久久国产一区二区| 亚洲欧美日韩卡通动漫| 多毛熟女@视频| 久久人人爽av亚洲精品天堂| 午夜福利视频精品| 成年av动漫网址| 亚洲在久久综合| 国产成人精品无人区| 桃花免费在线播放| 午夜老司机福利剧场| 国产成人欧美| 亚洲精品美女久久av网站| 男女下面插进去视频免费观看 | 水蜜桃什么品种好| 午夜福利乱码中文字幕| 母亲3免费完整高清在线观看 | 精品一区在线观看国产| 亚洲欧美日韩另类电影网站| 色5月婷婷丁香| 边亲边吃奶的免费视频| 多毛熟女@视频| 国产精品久久久久久精品电影小说| 国产亚洲av片在线观看秒播厂| 人人妻人人澡人人爽人人夜夜| 日本色播在线视频| 亚洲,一卡二卡三卡| 热re99久久精品国产66热6| 国产av一区二区精品久久| 高清黄色对白视频在线免费看| 亚洲成国产人片在线观看| 亚洲婷婷狠狠爱综合网| 高清在线视频一区二区三区| 大陆偷拍与自拍| 一区二区av电影网| 咕卡用的链子| 少妇精品久久久久久久| 纵有疾风起免费观看全集完整版| 国产国拍精品亚洲av在线观看| 99久国产av精品国产电影| 在线观看人妻少妇| 日本黄大片高清| 最近最新中文字幕免费大全7| 9191精品国产免费久久| 亚洲欧美色中文字幕在线| 九九在线视频观看精品| 高清av免费在线| 色婷婷久久久亚洲欧美| 日本91视频免费播放| 久久久久久久亚洲中文字幕| 一区二区三区四区激情视频| 青青草视频在线视频观看| 欧美精品av麻豆av| xxxhd国产人妻xxx| 一级毛片我不卡| 国产日韩欧美视频二区| 内地一区二区视频在线| 乱人伦中国视频| 亚洲欧美成人精品一区二区| 亚洲成色77777| 亚洲精品视频女| 国产国语露脸激情在线看| 2018国产大陆天天弄谢| 一区二区三区四区激情视频| 侵犯人妻中文字幕一二三四区| 男女无遮挡免费网站观看| 日本-黄色视频高清免费观看| 爱豆传媒免费全集在线观看| 又黄又粗又硬又大视频| 老女人水多毛片| 免费观看在线日韩| 黑人猛操日本美女一级片| 2021少妇久久久久久久久久久| 久久久久久久亚洲中文字幕| 亚洲一码二码三码区别大吗| 国产精品一区二区在线不卡| 亚洲国产成人一精品久久久| 日本免费在线观看一区| 啦啦啦在线观看免费高清www| 国产老妇伦熟女老妇高清| 久久婷婷青草| 午夜av观看不卡| 久久精品久久精品一区二区三区| 国产亚洲精品第一综合不卡 | 国产精品久久久久久av不卡| 大陆偷拍与自拍| 看非洲黑人一级黄片| 啦啦啦视频在线资源免费观看| 狠狠精品人妻久久久久久综合| 十八禁网站网址无遮挡| 成人国语在线视频| 一边摸一边做爽爽视频免费| 亚洲精品美女久久av网站| 国产极品粉嫩免费观看在线| 国内精品宾馆在线| 熟女av电影| 欧美日韩视频精品一区| 亚洲综合精品二区| 成年女人在线观看亚洲视频| 久久99蜜桃精品久久| 国产日韩一区二区三区精品不卡| 国产深夜福利视频在线观看| 亚洲精品aⅴ在线观看| 九色成人免费人妻av| 久久久久国产精品人妻一区二区| www.av在线官网国产| 日韩不卡一区二区三区视频在线| 国语对白做爰xxxⅹ性视频网站| 在线观看人妻少妇| 日韩人妻精品一区2区三区| 久久影院123| 亚洲国产欧美日韩在线播放| 在线看a的网站| 有码 亚洲区| 亚洲欧美日韩另类电影网站| 国产成人精品在线电影| 精品久久国产蜜桃| 熟女人妻精品中文字幕| 久久精品国产a三级三级三级| tube8黄色片| 精品久久久久久电影网| 国产亚洲欧美精品永久| 少妇猛男粗大的猛烈进出视频| 亚洲人与动物交配视频| 午夜影院在线不卡| 欧美最新免费一区二区三区| 大香蕉久久成人网| 亚洲av福利一区| 18禁在线无遮挡免费观看视频| 一边摸一边做爽爽视频免费| 精品国产一区二区三区四区第35| 国产免费现黄频在线看| av播播在线观看一区| 国产激情久久老熟女| 2021少妇久久久久久久久久久| 久久久久久人妻| 在线天堂最新版资源| 少妇的逼好多水| 麻豆乱淫一区二区| 99国产综合亚洲精品| 一边亲一边摸免费视频| 纯流量卡能插随身wifi吗| 99国产综合亚洲精品| 有码 亚洲区| 蜜桃国产av成人99| 久久女婷五月综合色啪小说| 男女边吃奶边做爰视频| 久久久久久久国产电影| 母亲3免费完整高清在线观看 | 中文精品一卡2卡3卡4更新| 天堂8中文在线网| 亚洲精品日韩在线中文字幕| 国产成人精品婷婷| 中文乱码字字幕精品一区二区三区| 乱人伦中国视频| 人妻一区二区av| 男女啪啪激烈高潮av片| 午夜久久久在线观看| 国产黄色视频一区二区在线观看| 亚洲第一av免费看| 99视频精品全部免费 在线| 亚洲欧美成人综合另类久久久| 日本猛色少妇xxxxx猛交久久| 高清不卡的av网站| 丝袜在线中文字幕| 久久久欧美国产精品| 永久免费av网站大全| av免费在线看不卡| 99久久综合免费| 中文字幕免费在线视频6| 久久久久国产网址| 国产熟女午夜一区二区三区| 成人漫画全彩无遮挡| 一边摸一边做爽爽视频免费| 秋霞伦理黄片| 永久网站在线| 十分钟在线观看高清视频www| 大片免费播放器 马上看| 亚洲精品国产av蜜桃| 国产探花极品一区二区| 亚洲精品美女久久久久99蜜臀 | 99热6这里只有精品| 一级爰片在线观看| 欧美日韩综合久久久久久| 亚洲国产av影院在线观看| 王馨瑶露胸无遮挡在线观看| av黄色大香蕉| 色网站视频免费| 精品一区二区三区四区五区乱码 | 欧美日韩国产mv在线观看视频| 久久久久久人妻| 成人免费观看视频高清| 久久精品久久久久久噜噜老黄| 久久97久久精品| 国产一区有黄有色的免费视频| 免费高清在线观看视频在线观看| 欧美国产精品va在线观看不卡| 深夜精品福利| 少妇的逼好多水| 女人被躁到高潮嗷嗷叫费观| 一级,二级,三级黄色视频| 国产av国产精品国产| 精品一区二区三卡| 亚洲美女视频黄频| 韩国高清视频一区二区三区| 肉色欧美久久久久久久蜜桃| 亚洲精品av麻豆狂野| 国产日韩欧美在线精品| 欧美bdsm另类| 亚洲五月色婷婷综合| 久久国产精品大桥未久av| 国产亚洲午夜精品一区二区久久| 黑人巨大精品欧美一区二区蜜桃 | 久久青草综合色| 亚洲精品国产av蜜桃| 男的添女的下面高潮视频| 久久国产亚洲av麻豆专区| h视频一区二区三区| 国产一区二区在线观看日韩| 黄网站色视频无遮挡免费观看| 99视频精品全部免费 在线| 精品第一国产精品| 精品人妻在线不人妻| 久久这里只有精品19| 少妇熟女欧美另类| 国产免费一级a男人的天堂| 国产av国产精品国产| 成年女人在线观看亚洲视频| 两性夫妻黄色片 | 免费av不卡在线播放| av一本久久久久| 亚洲av在线观看美女高潮| 1024视频免费在线观看| 一本久久精品| 亚洲久久久国产精品| 国产一区亚洲一区在线观看| av线在线观看网站| 亚洲天堂av无毛| 曰老女人黄片| 日韩熟女老妇一区二区性免费视频| 久久久精品94久久精品| 国产精品久久久av美女十八| 国产一区二区在线观看av| 黄网站色视频无遮挡免费观看| 免费观看av网站的网址| 黄色怎么调成土黄色| 日本欧美视频一区| 美女国产视频在线观看| 侵犯人妻中文字幕一二三四区| 国产不卡av网站在线观看| 满18在线观看网站| 亚洲国产精品国产精品| 熟女av电影| 只有这里有精品99| 国产男人的电影天堂91| 五月天丁香电影| 秋霞伦理黄片| 蜜臀久久99精品久久宅男| 日本av手机在线免费观看| 国产精品蜜桃在线观看| 久久久国产精品麻豆| 日本vs欧美在线观看视频| 国内精品宾馆在线| 午夜福利乱码中文字幕| 汤姆久久久久久久影院中文字幕| 天堂8中文在线网| 精品99又大又爽又粗少妇毛片| 插逼视频在线观看| 亚洲内射少妇av| 蜜桃国产av成人99| 大香蕉久久成人网| 99热全是精品| 国产女主播在线喷水免费视频网站| 精品少妇内射三级| www.av在线官网国产| 黑人猛操日本美女一级片| 免费黄频网站在线观看国产| 极品人妻少妇av视频| 国产精品久久久久久精品古装| 边亲边吃奶的免费视频| 国产精品嫩草影院av在线观看| 欧美最新免费一区二区三区| 国产成人a∨麻豆精品| 亚洲国产精品999| 大码成人一级视频| 久久久久久久亚洲中文字幕| 亚洲一级一片aⅴ在线观看| 中文字幕亚洲精品专区| 欧美人与性动交α欧美精品济南到 | 丝袜人妻中文字幕| 国产 精品1| 美女脱内裤让男人舔精品视频| 国产av国产精品国产| 亚洲欧美成人精品一区二区| 亚洲精品国产色婷婷电影| 国产免费又黄又爽又色| 国产麻豆69| 晚上一个人看的免费电影| 欧美亚洲 丝袜 人妻 在线| 9热在线视频观看99| 熟女电影av网| 午夜福利影视在线免费观看| 国产av码专区亚洲av| 卡戴珊不雅视频在线播放| 欧美人与善性xxx| √禁漫天堂资源中文www| 国产亚洲av片在线观看秒播厂| 视频区图区小说| 超碰97精品在线观看| 91午夜精品亚洲一区二区三区| 中文字幕精品免费在线观看视频 | 久久久久网色| 男女下面插进去视频免费观看 | 这个男人来自地球电影免费观看 | 美国免费a级毛片| 日韩熟女老妇一区二区性免费视频| 欧美日韩视频高清一区二区三区二| 国产成人午夜福利电影在线观看| 国产亚洲最大av| 精品国产一区二区三区四区第35| 免费不卡的大黄色大毛片视频在线观看| 三上悠亚av全集在线观看| 狠狠婷婷综合久久久久久88av| 亚洲国产av影院在线观看| 蜜臀久久99精品久久宅男| 777米奇影视久久| 日本欧美国产在线视频| 18+在线观看网站| 极品人妻少妇av视频| 9191精品国产免费久久| 午夜老司机福利剧场| 成人午夜精彩视频在线观看| 国产一区二区三区av在线| 国产精品久久久久成人av| a级片在线免费高清观看视频| freevideosex欧美| 一区二区三区四区激情视频| 熟妇人妻不卡中文字幕| 一级片免费观看大全| 国产男人的电影天堂91| 夜夜爽夜夜爽视频| 午夜激情久久久久久久| 精品一品国产午夜福利视频| 97在线视频观看| 免费人成在线观看视频色| av又黄又爽大尺度在线免费看| 老女人水多毛片| 亚洲av.av天堂| 丝瓜视频免费看黄片| 婷婷色综合大香蕉| 中文字幕亚洲精品专区| 搡女人真爽免费视频火全软件| 女的被弄到高潮叫床怎么办| 日韩伦理黄色片| 久久久久久久大尺度免费视频| 制服诱惑二区| 晚上一个人看的免费电影| 一区在线观看完整版| 婷婷成人精品国产| 国产欧美日韩综合在线一区二区| 99久久人妻综合| 少妇人妻 视频| 久久久久久久国产电影| 亚洲伊人久久精品综合| 丝袜喷水一区| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久av美女十八| 九色成人免费人妻av| 亚洲丝袜综合中文字幕| 久久久精品94久久精品| 尾随美女入室| 黄色怎么调成土黄色| 国产探花极品一区二区| 久久韩国三级中文字幕| 插逼视频在线观看| 国产老妇伦熟女老妇高清| 在线天堂中文资源库| 黄色怎么调成土黄色| 欧美日本中文国产一区发布| 久久韩国三级中文字幕| 国产熟女欧美一区二区| www日本在线高清视频| 亚洲一码二码三码区别大吗| 狂野欧美激情性bbbbbb| 蜜臀久久99精品久久宅男| 丰满饥渴人妻一区二区三| 中文欧美无线码| 亚洲精品456在线播放app| 观看av在线不卡| 欧美成人精品欧美一级黄| 日本色播在线视频| 99久久人妻综合| 桃花免费在线播放| 日韩电影二区| 91久久精品国产一区二区三区| 少妇的逼水好多| 在线观看国产h片| 国产成人一区二区在线| 91在线精品国自产拍蜜月| 大香蕉久久网| av国产久精品久网站免费入址| 亚洲av电影在线进入| 大片免费播放器 马上看| 美女脱内裤让男人舔精品视频| 我要看黄色一级片免费的| 午夜福利在线观看免费完整高清在| 热99国产精品久久久久久7| 亚洲一码二码三码区别大吗| 亚洲五月色婷婷综合| av播播在线观看一区| 久热久热在线精品观看| 亚洲国产精品专区欧美| 国产xxxxx性猛交| 少妇的逼好多水| 亚洲人成网站在线观看播放| 天美传媒精品一区二区| 久久久国产精品麻豆| 一级毛片我不卡| 最近最新中文字幕大全免费视频 | 亚洲av日韩在线播放| 水蜜桃什么品种好| 如日韩欧美国产精品一区二区三区| 欧美日韩亚洲高清精品| 99久久中文字幕三级久久日本| 又黄又爽又刺激的免费视频.| 国产亚洲欧美精品永久| 尾随美女入室| 精品少妇内射三级| 寂寞人妻少妇视频99o| 欧美xxxx性猛交bbbb| 中文欧美无线码| 成年av动漫网址| 最近最新中文字幕免费大全7| 男人添女人高潮全过程视频| 亚洲av电影在线观看一区二区三区| 日本黄色日本黄色录像| 午夜老司机福利剧场| 久久99一区二区三区| 亚洲五月色婷婷综合| 亚洲天堂av无毛| 亚洲av在线观看美女高潮| 亚洲精品色激情综合| 国产av一区二区精品久久| 国产综合精华液| 久久久久久久久久久免费av| 女性被躁到高潮视频| 亚洲av中文av极速乱| 免费看av在线观看网站| 日本wwww免费看| 日本欧美视频一区| 久久热在线av| 国产男女超爽视频在线观看| 在线 av 中文字幕| 中文欧美无线码| 黑人巨大精品欧美一区二区蜜桃 | 国产毛片在线视频| 亚洲婷婷狠狠爱综合网| 亚洲av欧美aⅴ国产| 国产精品一二三区在线看| 男男h啪啪无遮挡| 亚洲av电影在线进入| 热99国产精品久久久久久7| 亚洲av中文av极速乱| 欧美 日韩 精品 国产| av片东京热男人的天堂| 考比视频在线观看| 久久 成人 亚洲| 成年美女黄网站色视频大全免费| 黄色视频在线播放观看不卡| 两性夫妻黄色片 | 午夜av观看不卡| 十八禁高潮呻吟视频| 九色成人免费人妻av| 欧美人与善性xxx| 99九九在线精品视频| 国产免费一级a男人的天堂| 建设人人有责人人尽责人人享有的| 日韩精品有码人妻一区| 男女无遮挡免费网站观看| 黄片播放在线免费| 国产熟女欧美一区二区| 韩国av在线不卡| 人人妻人人澡人人看| 亚洲国产毛片av蜜桃av| 中国国产av一级| 人人妻人人澡人人爽人人夜夜| 中文乱码字字幕精品一区二区三区| 国产精品偷伦视频观看了| 丝瓜视频免费看黄片| 日韩av免费高清视频| 欧美+日韩+精品| 王馨瑶露胸无遮挡在线观看| 日韩一本色道免费dvd| 久久久久久人人人人人| 国产成人精品婷婷| 欧美日本中文国产一区发布| 热99国产精品久久久久久7| av在线观看视频网站免费| 内地一区二区视频在线| 色哟哟·www| 精品人妻熟女毛片av久久网站| 亚洲欧美清纯卡通| 日产精品乱码卡一卡2卡三| av线在线观看网站| 激情五月婷婷亚洲| 亚洲成色77777| 水蜜桃什么品种好| 国产国语露脸激情在线看| 精品视频人人做人人爽| 亚洲,一卡二卡三卡| 亚洲国产精品一区二区三区在线| 婷婷成人精品国产| 五月玫瑰六月丁香| 深夜精品福利| 2021少妇久久久久久久久久久| 一区二区三区精品91| av线在线观看网站| 狂野欧美激情性xxxx在线观看| 成人亚洲精品一区在线观看| 欧美xxxx性猛交bbbb| 国产精品国产三级专区第一集| 久久精品久久久久久久性| 波多野结衣一区麻豆| 啦啦啦视频在线资源免费观看| 久久久久久人妻| 国产精品秋霞免费鲁丝片| 九九爱精品视频在线观看| 日韩中字成人| 亚洲国产欧美日韩在线播放| 蜜臀久久99精品久久宅男| 日韩制服骚丝袜av| 天堂俺去俺来也www色官网| freevideosex欧美| 三级国产精品片| 汤姆久久久久久久影院中文字幕| 日韩一区二区三区影片| av播播在线观看一区| 美女视频免费永久观看网站| 波多野结衣一区麻豆| 久久久久久久国产电影| 国产精品秋霞免费鲁丝片| 久久精品国产a三级三级三级| 久久久久久伊人网av| 午夜影院在线不卡| 男女边吃奶边做爰视频| a级片在线免费高清观看视频| 午夜免费鲁丝| 国产在线免费精品| av国产久精品久网站免费入址| 肉色欧美久久久久久久蜜桃| 五月天丁香电影| 一级,二级,三级黄色视频| 国产乱来视频区| 婷婷色综合大香蕉| 内地一区二区视频在线| 成人国语在线视频| 在线观看免费视频网站a站| 91精品国产国语对白视频| 天天躁夜夜躁狠狠久久av| 乱人伦中国视频| 成人漫画全彩无遮挡| 精品熟女少妇av免费看| 女的被弄到高潮叫床怎么办| 欧美国产精品一级二级三级| 欧美最新免费一区二区三区| 国产精品一二三区在线看| 日韩三级伦理在线观看| 在线天堂中文资源库| 国产精品99久久99久久久不卡 | 全区人妻精品视频| 国产成人欧美| 日本猛色少妇xxxxx猛交久久| 亚洲成色77777| 国产精品一二三区在线看| 精品午夜福利在线看| 精品酒店卫生间| 99热6这里只有精品| 18禁裸乳无遮挡动漫免费视频| 国产探花极品一区二区| 午夜激情av网站| 美女国产视频在线观看| 91精品三级在线观看| a 毛片基地| 亚洲伊人久久精品综合| 亚洲成av片中文字幕在线观看 | 欧美日韩亚洲高清精品| 一级毛片电影观看| 国产伦理片在线播放av一区| 大香蕉97超碰在线| 久久久久精品久久久久真实原创| 国产不卡av网站在线观看| 国产亚洲精品久久久com| 99香蕉大伊视频| 国产成人午夜福利电影在线观看| 我要看黄色一级片免费的| 国产一区二区在线观看日韩| 国产成人午夜福利电影在线观看| 久久国内精品自在自线图片| 日韩伦理黄色片| 宅男免费午夜| 久久久久网色| 精品一区二区免费观看| 亚洲成av片中文字幕在线观看 | 久久久久久人人人人人| 精品一区在线观看国产| 永久网站在线| 黑人巨大精品欧美一区二区蜜桃 | 熟女电影av网| 中文字幕制服av| 狂野欧美激情性xxxx在线观看| 亚洲欧美一区二区三区黑人 | 免费高清在线观看日韩| 在线观看www视频免费| 国产激情久久老熟女| 在线天堂最新版资源| 久久精品久久久久久久性|