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

    Determinants of species assemblages of insect pests in alpine forest ecosystems of western China

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

    Chun-Jing Wang,Rong Wang,Chun-Mei Yu,Yongcuo Pubu,Wan-Gui Sun,Xiao-Feng Dang,Qiang-Feng Li and Ji-Zhong Wan*

    Abstract

    Keywords: Alpine forest ecosystems, Community completeness, Dark diversity, Insect pest, Species co-occurrence,Species diversity

    Introduction

    The International Plant Protection Convention (IPPC)defined pests as species harmful to cultivated plants(agriculture, horticulture, and forestry), natural flora and plant products ranging from wood products and dried plant materials for handicrafts to fresh fruits and vegetables for consumption or processing and grains for food or feed (https://www.ippc.int). Insect pests are the surface feeders that damage foliage by chewing, piercing, or sucking and subsurface feeders, which can damage natural resources and social development across multiple ecosystems and regions worldwide (Clark and Kenna 2010; Savary et al. 2012; Venette 2017; Prather and Laws 2018; Cornelissen et al. 2019; Fei et al. 2019). Agricultural crops and forest resources are attacked by various insect pests, which pose a significant threat to global food security (Donatelli et al. 2017; Prather and Laws 2018; Fei et al. 2019). For example, insect pests (e.g.,Coleoptera and Hymenoptera) can damage crops and food production (L?vei and Sunderland 1996; Cebolla et al.2018).Furthermore,insect pests can threaten forest ecosystems jointly in many cases (Ji et al. 2011; Guo et al. 2019). Wood production declines because of insect pest damage (Sallé et al. 2014; Tiberi et al. 2016). Thus,the utilization efficiency of forest resources is reduced under a growing population. Hence, we should explore the mechanism of insect pest expansion in forest ecosystems.

    Alpine forest ecosystems are more sensitive to environmental changes than other forest ecosystems in temperate and tropical areas (Chapin III and K?rner 1994;Grabherr 2009; Minckley et al. 2012; Boix-Fayos et al.2020). Rising temperatures can affect insect physiology and associated fitness and shift populations in space and time, thereby altering community interactions and composition (Shah et al. 2020). Johnson et al. (2010) showed that forest insect outbreaks are generally more common in alpine regions due to global warming. In the context of global ecological change, plant health may be among the most critical determinants of ecological processes influenced by precipitation and temperature variability in alpine forest ecosystems (Pautasso et al. 2012; Milosavljevi? et al. 2016; Boix-Fayos et al. 2020; Grünig et al.2020; Lehmann et al. 2020). Numerous experimental and field studies (e.g., Frago et al. 2012; Grünig et al.2020; Lehmann et al. 2020) have shown that plant-pest interactions might vary sharply with global environmental changes. Furthermore, insect pests always occur together in similar areas (Guo et al. 2019; Grünig et al.2020). However, these previous studies (e.g., Logan et al.2003; Johnson et al. 2010; Guo et al. 2019; Wang et al.2021) have only concentrated on the effects of environmental changes on insect pests and have not considered their effects on pest species assemblages in alpine forest ecosystems. Modeling insect pest species assemblages can determine the risk of invasion (Worner and Gevrey 2006). Therefore, the relationships between patterns of species assemblages and environmental conditions in alpine forest ecosystems must be modeled to provide a reference for the prevention and control of insect pests in environment-sensitive areas. Here, we explored the determinants of insect pest species assemblages in alpine forest ecosystems to establish effective biological control systems.

    Species diversity is a critical index for insect species assemblages. As suggested by P?rtel et al. (2013), four indices (i.e., observed diversity, dark diversity, community completeness, and species pool) should be used to quantify species diversity. The concept of dark diversity could be crucial for developing effective tools for biological risk (Ronk et al. 2017). Dark diversity implies the absence of suitable species from a species pool despite the appropriate biogeographic history and current ecological and environmental conditions for their presence(P?rtel et al. 2013; Ronk et al. 2017). In consideration of dark diversity, community completeness expresses how much of the species pool is realized within a local community. More complete communities should exert greater biotic resistance to biological invasion because they are more saturated and have greater resource depletion. To assess the potential of biological risk, we should predict the dark diversity and community completeness of insect pests (P?rtel et al. 2013; Ronk et al. 2017; Venette 2017). Biological control systems could be improved for potential insect pest risks based on the concepts of dark diversity and community completeness. The appropriate selection of environmental determinants of species diversity can improve the robustness of mechanistic insect pest forecast modeling, which is key to improving prevention and control (Magarey and Isard 2017).

    Environmental changes (e.g., climatic changes) can drive the diversity of insect pests and shape their species assemblages in terrestrial ecosystems (Venette 2017;Cornelissen et al. 2019; Guo et al. 2019). Numerous studies have shown that abiotic factors (i.e., climate and human influence) are the main drivers of the species diversity and distribution of insect pests (Savary et al.2012; Venette 2017; Hill et al. 2019; Jactel et al. 2019).The main effects of climate change drivers on forest insect performance and damage (e.g., insect survival and voltinism; Jactel et al. 2019; Grünig et al. 2020). Warmer temperatures can lead to a higher number of generations per year and higher survival, and severe droughts can weaken the ability of tree species to resist insect attack during (Baeten et al. 2019; Guo et al. 2019; Jactel et al.2019; Grünig et al. 2020). In particular, strong human influences can increase the dispersal of insect pests worldwide (Cavers 2015). Thus, some insect pests with dark diversity may occur and then transfer the dark diversity of insect pest species to the observed diversity.Furthermore, species assemblages should form due to the positive interactions between different pest species in ecosystems (Cazelles et al. 2016; Nunez-Mir et al.2017; D’ Amen et al. 2018; Freilich et al. 2018; Grünig et al. 2020). Positive biotic interactions can shape spatial patterns of species aggregation (D’ Amen et al. 2018).Hence, positive biotic interactions may be the determinants of species assemblages of insect pests in alpine forest ecosystems. Plant-pest interactions also play an important role in species assemblages (Frago et al. 2012;Grünig et al. 2020). Host species diversity ranges are related to species abundance distributions of insect pests(Guo et al. 2019; Rasche and Taylor 2019; Delbac et al.2020). For example, tree species are likely to regulate the invasion of insect pests as their hosts (Guo et al. 2019).Hence, we sought to determine whether pest-pest and plant-pest interactions (i.e., biotic factors) coupled with environmental changes (i.e., abiotic factors) may be the determinants of species assemblages of insect pests in alpine forest ecosystems based on four indices (i.e., observeddiversity,darkdiversity,community completeness, and species pool). If we address this scientific issue, we could develop an effective indicator for insect pest risk assessment in alpine forest ecosystems based on relationships between abiotic factors and species diversity (Worner and Gevrey 2006).

    Here, we proposed the following four hypotheses: 1)Hypothesis 1 (H1): insect pest species interactions may be positive in alpine forest ecosystems; 2) Hypothesis 2(H2): dark diversity of insect pests may occur in alpine forest ecosystems; 3) Hypothesis 3 (H3): climatic factors and human influence may be the determinants of insect pest species diversity in alpine forest ecosystems; and 4)Hypothesis 4 (H4): the effects of abiotic factors on dark diversity and species pools may be larger than that on observed diversity. To test these four hypotheses, a field investigation was conducted for insect pests in the forest ecosystems of Qinghai Province, China, a region of the Qinghai-Tibet Plateau. We used multisite and multispecies measures of co-occurrence patterns to assess species assemblages in alpine forest ecosystems (Veech 2013;Griffith et al. 2016; Freilich et al. 2018). These measures were used to assess the interaction among different species (Veech 2013; Griffith et al. 2016). Co-occurrence methods are increasingly utilized in ecology to infer species interaction networks (Veech 2013; Griffith et al.2016; Freilich et al. 2018). Understanding species cooccurrence patterns can help determine the causes and consequences of the species evolution, histories, coexistence mechanisms, competition, and environments that lead to the community structure and assemblages(Veech 2013; Freilich et al. 2018).The approach of Lewis et al. (2016) was applied to quantify insect pest species diversity based on four indices of species diversity: observed diversity, dark diversity, community completeness, and species pool.

    Materials and methods

    Species data

    Qinghai Province, China, is a region of the Qinghai-Tibet Plateau from 103°04′–116°33′ E, 31°39′–39°19′ N.Elevations range from 1650 to 6860 m, and the average elevation is greater than 3000 m. Qinghai Province includes representative alpine forest ecosystems in western China (Hou 2001). Our field investigation was conducted for 3 years in a network of sites based on vegetation density across the entire province. The field investigation sites are shown in Fig. 1. This field investigation was managed by the Qinghai Forestry and Grassland Administration for 3 years. These plots can cover most forests across most counties of Qinghai Province.In each county of Qinghai Province, we conducted fieldwork for insect pests in alpine forest ecosystems based on the Vegetation Atlas of China 1:1,000,000 from robust field investigations (Hou 2001; Zhang et al. 2007).Our fieldwork covered most of the vegetation recorded in the Vegetation Atlas of China at 1:1,000,000 in Qinghai Province. The plot size was approximately 3 acres.Specifically, the sizes and numbers of the plots were determined according to the ranges of host species in one specific forest ecosystem (Wang et al. 2021).

    A given pest species can invade different parts of its plant hosts in alpine forest ecosystems. We identified the damage caused by a given insect pest based on the percentage of leaf loss, levels of stunting, and proportion of seedling death in host species in the field and conducted insect pest identification based on the knowledge of biocontrol and pest science experts. The experts assessed whether a specific insect species was a pest based on their knowledge and previous information from the first and second forestry pest survey projects in Qinghai Province, China. If the experts indicated that a given species was a pest, it was labeled as such. However, if insufficient evidence was available, a species was identified as a “potential insect pest”. The insect pest species studied could invade alpine woody plants (e.g.,Picea likiangensis, Populus cathayana, Picea crassifolia,Populus davidiana, and Hippophae rhamnoides) as harmful herbivores with high fecundity in forest ecosystems. Their invasion could lead to high percentages of leaf loss, stunting,and seedling death of host species.

    Fig.1 Species diversity(i.e.,observed diversity,dark diversity, community completeness,and species pool)of insect pests based on study sites in the forest ecosystems of Qinghai Province,China.The color of study sites(from blue to red) represents increasing levels of observed diversity(a),dark diversity(b),community completeness (c),and species pool(d) in the figures.The color of the figures from light to dark representedincreasing levels of human influence

    In each plot, we investigated insect pests in more than 30 host plant individuals twice in the seasons of plant growth in one given year from 2014 to 2016. We recoded the scientific names of insect pest species, the parts of plant hosts invaded by insect pests, and the number of insect pest individuals in each plot (Wang et al. 2021). In total, 172 insect pest species were studied. These 172 insect pest species could invade five parts of their plant hosts (i.e., leaves, roots, seeds and fruits,trunks, and twigs) in forest ecosystems and pose a large threat to forest ecosystems in Qinghai Province, China.

    Abiotic factors

    Previous studies have shown that abiotic factors, namely,climate factors and human influence, can determine insect pest species assemblages (Savary et al. 2012; Hill et al. 2019; Jactel et al. 2019). Data on climate factors were downloaded with a spatial resolution of 0.5 arcminutes from the CHELSA (Climatologies at high resolution for the Earth’s land surface areas) database(https://chelsa-climate.org/; Karger et al. 2017). We used 19 bioclimatic variables to assess the potential determinants of insect pest species assemblages in alpine forest ecosystems as follows: Bio1=annual mean temperature;Bio2=mean diurnal range; Bio3=isothermality; Bio4=temperature seasonality; Bio5=max temperature of warmest month; Bio6=min temperature of coldest month; Bio7=temperature annual range; Bio8=mean temperature of wettest quarter;Bio9=mean temperature of driest quarter; Bio10=mean temperature of warmest quarter; Bio11=mean temperature of coldest quarter;Bio12=annual precipitation; Bio13=precipitation of wettest month; Bio14=precipitation of driest month;Bio15=precipitation seasonality; Bio16=precipitation of wettest quarter; Bio17=precipitation of driest quarter;Bio18=precipitation of warmest quarter; and Bio19=precipitation of coldest quarter (Karger et al. 2017). We used principal component analysis (PCA) to extract the first two principal components (PCs) that explained the majority of 19 climatic variables. The first two PCs explained 65.93% of the climatic variation. Climatic PC1 was positively related to Bio1, Bio5, Bio8, and Bio10, and climatic PC2 was positively related to Bio6, Bio11, Bio12 Bio13, and Bio16.

    Sanderson et al. (2002) created a map of the Human Influence Index with a spatial resolution of 0.5 arcminutes based on human population pressure (population density), land use and infrastructure (built-up areas,night-time lighting, and land use/cover), as well as access (coastlines, roads, railroads, and navigable rivers;downloaded from https://sedac.ciesin.columbia.edu/data/set/wildareas-v2-human-footprint-geographic). The human footprint (HF) was strictly related to the Human Influence Index, with a high human footprint index indicating the intactness, naturalness, and function of natural communities. We used a 0.5-arcminute (~1 km)resolution for all abiotic data. We extracted the data on abiotic factors (i.e., climate factors and human influence)based on the longitude and latitude of each field investigation site.

    Species co-occurrence pattern

    We used studies by Veech (2013) and Griffith et al. (2016)to analyze species co-occurrence patterns of insect pests in alpine forest ecosystems. A presence-absence matrix was built for analyzing the species co-occurrence patterns of insect pest species based on study sites.We used a probabilistic model to test for statistically significant pairwise patterns of species co-occurrence based on the concept of effect sizes using sites × species and parts × species matrices of insect pest species (Veech 2013; Griffith et al. 2016).The effect size was defined as the absolute difference between observed and expected co-occurrence for insect pests in alpine forest ecosystems based on comparison with unbiased null model distributions (Veech 2013; Griffith et al.2016).The standardized range of effect size was from-1 to 1 and controlled for the number of sampling study sites(Veech 2013; Griffith et al. 2016). The effect size was used as an important index of interaction between two given species (Veech 2013; Cazelles et al. 2016; Freilich et al.2018). We used the number of association links (i.e., the number of one given species with the other species with positive species interaction) and average effect size(i.e.,the average values of effect sizes of association of one given species with the other species)to quantify the site-specific and part-specific interaction of the given species with the other species (Veech 2013; Freilich et al. 2018). Thus, the sitespecific and part-specific cooccurrence levels of insect pest species could be assessed by the interaction among different species based on the number of association links and average effect size (Veech 2013; Cazelles et al. 2016; Freilich et al.2018).

    Species diversity

    We used the approach of Lewis et al. (2016) to quantify insect pest species diversity based on four indices: observed diversity, dark diversity, community completeness, and species pool. We calculated the number of pest species (i.e., species richness) for observed diversity at the study sites. We quantified dark diversity based on the species co-occurrence patterns across the forest ecosystems of Qinghai Province, China. The assessment of species co-occurrence patterns was based on the probability of a species occurring in an investigation site using the Beals index (Lewis et al. 2016). According to the sites × species matrix for pest species, the Beals index was correlated with species frequency; therefore,each species was assigned a threshold value for inclusion in dark diversity (Lewis et al. 2016). These thresholds were derived from each species’ estimated probabilities of occurrence at sites where they were present, and they were set at the 5% quantile of those probabilities to omit outliers (Lewis et al. 2016). The details for dark diversity assessments were shared in the study by Lewis et al.(2016). We used the logarithmic ratio of observed and dark diversities to calculate the community completeness, and the species pool was the sum of the observed and dark diversities for a given study site, which indicated how much of the site-specific species pool was represented in local communities, i.e., the community completeness (P?rtel et al. 2013).

    Synthesis

    To test Hypothesis 1, we used ordinary least squares(OLS) regression modeling to explore the relationships between the number of association links and effect sizes for insect pests in forest ecosystems in Qinghai Province.We then compared the number of association links with positive, negative, and random effect sizes across 172 insect pest species. Following this, we applied Student’s ttest to check the differences in the number of association links and effect size among the five parts of the plant hosts that were invaded (i.e., leaves, roots, seeds and fruits, trunks, and twigs) in the ecosystems studied.To test Hypothesis 2, we compared four indices for pest species diversity (i.e., observed diversity, dark diversity,community completeness, and species pool) across different study sites. We used ANOVA tests to check the differences between the observed and dark diversities at these sites. To test Hypothesis 3, we used OLS regression modeling to explore the relationships between the same pest species diversity indices used to test Hypothesis 2 and abiotic variables, such as climate factors (PCs)and human influence. A boosted regression tree (BRT)was used to explore the contribution of abiotic factors(Bio1–19 and human influence) to the formation of insect pest species assemblages in alpine forest ecosystems.To test Hypothesis 4, we examined the results of OLS regression modeling to find the relationships between pest species diversity and abiotic factors. All analyses were conducted in R (https://www.r-project.org/) and JMP 11.0 (https://www.jmp.com/zh_cn/home.html).

    Results

    Fig.2 The number of association links and effect sizes for insect pest species based on species co-occurrence analysis across different parts of plant hosts(i.e.,leaves,roots,seeds and fruits,trunks,and twigs).There was a significant relationship between the number of association links and effect size in the forest ecosystems of Qinghai Province,China(R2=0.258;P <0.001).The line dash represents the 95%confidence intervals in regression modeling

    We found that the effect sizes of species cooccurrence among insect pests were positive for 839 pairs, negative for 152 pairs, and random for 1175 pairs, as shown in the Supporting Information. The number of association links ranged from 1 to 93 and the effect sizes ranged from -0.007 to 0.021 across different insect pest species(Fig. 2). The largest number of association links was 93,and the largest effect size was 0.021 (Fig. 2). There was a significant relationship between the number of association links and the effect size across all insect pest species (P <0.05; Fig. 2). We found that 86, 14, 46, and 60 pest species invaded the leaves, roots, trunks, and twigs,respectively (Fig. 2), and the average effect sizes were positive across the leaves, roots, seeds and fruits, trunks,and twigs (Fig. 3, details in the Supporting Information).We found that there were significant differences between seeds and fruit and the other parts (i.e., leaves and roots) for the number of association links (P <0.05;Table 1) and between twigs and trunks for the effect size(P <0.05; Fig. 3;Table 1).

    The observed diversity, dark diversity, community completeness, and species pool ranged from 1 to 52, 2 to 58, -4.060 to 1.553, and 4 to 63, respectively (Fig. 1).The areas with high species diversity, particularly dark diversity, community completeness, and species pool,were distributed in the eastern forest ecosystems of Qinghai Province (Fig. 1). Dark diversity was significantly higher than the observed diversity for insect pests in the Qinghai forest ecosystems(P <0.05; Fig.1).

    Fig.3 Effect sizes for insect pests across five parts (i.e.,leaves,roots,seeds and fruits,trunks,and twigs)of plant hosts in the forest ecosystems of Qinghai Province,China.The upper and lower bars represent the 95%confidence intervals of the contribution of environmental variables,and the lines represent the medium values of the contribution of environmental variables

    Table 1 Difference tests for the number of association links and effect sizes for insect pest species among five parts of plant hosts(i.e.,leaves, roots, seeds and fruits, trunks, and twigs). Difference represents the minus values of mean number of association links and effect sizes between two specific parts of plant hosts (Part1 minus Part2).The results of this table were shown based on the pairs of two specific parts of plant hosts

    There were significant positive relationships between climatic PCs (i.e., PC1 and PC2) and pest species diversity indices (P <0.05; Fig. 4). The effects of climatic PC1 on pest species diversity were higher than those of PC2(Fig. 4). High temperatures and precipitation could result in high pest species diversity, particularly in terms of dark diversity and species pools (Fig. 4). The effects of human influence on pest species diversity were significantly positive according to all four indices and significantly greater than the effects of climate factors (P <0.05; Fig. 4). High temperatures, precipitation, and human influence could lead to high community completeness (Fig. 4). The annual mean temperature and temperature seasonality contributed the most to the species diversity indices for insect pests in the alpine forest ecosystems studied (Table 2).

    The effects of human influence on pest species diversity were significantly larger than those of climate factors in Qinghai Province, China (P <0.05; Fig. 4). Furthermore, the effects of abiotic factors (i.e., climate factors and human influence) on the species pool were larger than the observed and dark species diversity (Fig. 4).However, the effects of abiotic factors on dark diversity were smaller than those on observed diversity (Fig. 4).All four indices of pest species diversity showed that the effect of abiotic factors on the species pool was the largest (Fig. 4).

    Discussion

    The determinants of insect pest species assemblages in alpine forest ecosystems were evaluated based on species interactions and the effects of abiotic factors on pest species diversity to improve the prevention and control of pest invasions. Species interactions shaped species assemblages in alpine forest ecosystems. Our results showed that the abiotic factors examined determined the levels of observed diversity, dark diversity, community completeness, and species pool of insect pests in the ecosystems studied. The results of our study, which provides a mechanism for shaping insect species assemblages for pest management, showed that high dark diversity and low community completeness among insect pests in alpine forest ecosystems were key factors for improving pest management in alpine forests. The insect pest species studied had great potential to generate large levels of leaf loss, stunting, and seedling death in their host species, thereby posing a threat to woody plant species in alpine forests. Thus, ecosystems with high dark diversity and low community completeness of insect pests present a high risk of damage. Therefore,high dark diversity and low community completeness of insect pests in alpine forest ecosystems can support pest management in alpine forests.

    Pest species interactions

    Fig.4 Relationships between abiotic factors(i.e.,climate factors and human influence)and species diversity(i.e.,observed diversity,dark diversity,community completeness,and species pool)in the forest ecosystems of Qinghai Province,China.Climatic PC1 was positively related to Bio1,Bio5,Bio8,and Bio10,and climatic PC2 was positively related to Bio6,Bio11,Bio12 Bio13,and Bio16,as described in the Materials and Methods section.All the relationships were significant(P <0.01)

    Our results described above and the significant relationship observed between the number of association links and effect size for insect pest species indicated that positive pest species interactions were detected widely in alpine forest ecosystems. Furthermore, the number of positive effect sizes of species co-occurrence was higher than the number of negative effect sizes. The high number of association links (i.e., 839 pairs) suggest that insect pest species could invade similar areas jointly.These pest species invaded the leaves, trunks, and twigs of plants. Species interactions were mainly the determinant of species assemblages of insect pests in alpine forest ecosystems. However, 1175 association pairs were random for species co-occurrence among insect pests, indicating that there were no associations among some specific insect pest species in alpine forest ecosystems. Most of the nonrandom positive and negative species associations could be attributed to environmental filtering and/or dispersal limitations (D’ Amen et al.2018;Freilich et al. 2018;Blanchet et al. 2020).Biological interpretations, including dispersal limitations, environmental niche differences, and biotic interactions, cannot contribute to insect pest species assemblages in alpine forest ecosystems (D’ Amen et al. 2018; Blanchet et al.2020). Therefore, our study could provide evidence thatspecies interactions among insect pests may not occur in some cases for forest ecosystems.

    Table 2 The contribution of abiotic factors(Bio1–19,and human influence) to observed diversity,dark diversity, community completeness,and species pool for insect pests based on Boosted Regression Tree (BRT). The codes of abiotic factors were shown in Materials and Methods

    Species co-occurrence patterns are an important indicator of species assemblages of insect pests across different spatial scales (Horner-Devine et al. 2007; Veech 2013; Cazelles et al. 2016; D’ Amen et al. 2018; Freilich et al. 2018). Species interactions are the basic driver of species assemblages in alpine areas (Cazelles et al. 2016;D’ Amen et al. 2018; Freilich et al. 2018). The biotic resistance hypothesis and fluctuating resource hypothesis could support our findings (Maron and Vilà, 2001;Müller-Sch?rer et al. 2004; Sax et al. 2007). Biotic interactions with native species can affect the establishment and spread of invasive plant species, and alpine forest ecosystems provide hosts that could promote the invasion of pest species (Nunez-Mir et al. 2017; Venette 2017; Guo et al. 2019; Grünig et al. 2020). Tree diversity,as an effective biotic factor, plays an important role in insect pest species assemblages in alpine regions (Guyot et al. 2019; Guo et al. 2019; Jactel et al. 2021). Higher tree species diversity may lead to higher resistance to insect pest invasion (Guyot et al. 2016, 2019; Jactel et al.2021). Several main mechanisms have been proposed to explain associational resistance effects (Guyot et al.2016, 2019; Jactel et al. 2019, 2021; Poeydebat et al.2021). In a diverse host community, host individuals are more likely isolated from one another due to the presence of nonhost species (Guo et al. 2019; Jactel et al.2021). Furthermore, low diversity resident communities with low biomass have a weak ability to resist the establishment of pest invaders in alpine forest ecosystems(Forister et al. 2012; Nunez-Mir et al. 2017; Fei et al.2019). Thus, negative relationships may occur between tree species diversity and insect pest species diversity.Furthermore, heterospecific neighborhoods favor a greater abundance of herbivore natural enemies (predators and parasitoids) because of the increased availability of habitats or alternative resources, which leads to reduced herbivore abundance and damage (Guo et al.2019; Guyot et al. 2019; Jactel et al. 2021). Hence, pest species community assemblages can invade alpine forest ecosystems widely due to their hosts and weak competition.

    The fluctuating resource hypothesis assumes that an invading species has access to available resources (e.g.,light, nutrients, and water for pests) and greater success in invading a community if it does not encounter intense competition for resources from resident species (Klapwijk et al. 2016; Pearson et al. 2018; Guo et al. 2019).Limited resource availability causes competition in various parts of plants as pests compete for resources(Sequeira and Mackauer 1992; Eyles et al. 2010). However, we found that the average effect sizes were positive across all five invaded parts of plant hosts (i.e., leaves,roots, seeds and fruits, trunks, and twigs) in the forest ecosystems of Qinghai Province, China, indicating that there were positive species interactions for insect pests across different parts of plant hosts in alpine forest ecosystems. The results on the number of association links and effect size indicated that there were significant differences in species interactions for insect pests across the different invaded parts of hosts for insect pests. Host resource availability could affect species interactions among insect pests in cultivated and forest ecosystems(Eyles et al. 2010; Klapwijk et al. 2016; Nunez-Mir et al.2017; Pearson et al. 2018; Delbac et al. 2020). Plant host parts (e.g., seeds and fruits) provide sufficient resources for insect pests in alpine forest ecosystems (Eyles et al.2010; Guo et al. 2019). Therefore, species cooccurrences of insect pests could occur widely in alpine forest ecosystems due to the availability of resources from their hosts. Thus, species assemblages of insect pests are determined by invaded parts of plant hosts in alpine forest ecosystems, which indicates that the invaded parts of plant hosts should be monitored depending on different assemblages of insect pests.

    Effects of abiotic factors on species diversity

    We found that the areas with high levels of dark diversity, community completeness, and species pools were distributed in the eastern forest ecosystems of Qinghai Province, China, and dark diversity was significantly higher than the observed diversity for insect pests in the forest ecosystems. This finding indicates a high potential invasion risk in alpine forest ecosystems (P?rtel et al.2013; Ronk et al. 2017). Dark diversity is the set of absent species from a study site but present in the surrounding region and potentially able to inhabit the same ecological conditions (Lewis et al. 2017; P?rtel et al.2013). High dark diversity may lead to insect pest outbreaks in alpine forest ecosystems (Lewis et al. 2017;Ronk et al. 2017), and this high potential risk may be caused by abiotic factors (i.e., climate factors and human influence; Gallardo et al. 2015; Kulhanek et al. 2011).We found that the effects of abiotic factors on the species pool were larger than the effects of observed and dark species diversity,indicating that insect pest invasion could occur widely in alpine forest ecosystems (P?rtel et al. 2013; Lewis et al. 2017; Ronk et al. 2017). It is necessary to determine the predictors of pest species diversity to prevent and control pest risk (Gurr et al. 2017;Lewis et al. 2017).

    Climatic factors, particularly annual mean temperature and temperature seasonality, are the main drivers of insect pest species assemblages in alpine forest ecosystems. Our analysis showed that high temperatures and precipitation levels could lead to high pest species diversity, particularly in terms of dark diversity and species pools. Numerous studies have shown that climatic factors are the determinants of observed insect diversity at large scales (Venette 2017; Grünig et al. 2020; Lehmann et al. 2020; Schneider et al. 2021; Wang et al. 2021).Similarly, our study provides evidence of the effects of climatic conditions on dark diversity and species pools for insect pests in alpine forest ecosystems. Temperature shifts can directly influence insect populations by modifying the rates of development, reproduction, and mortality, while reductions in host tree resistance caused by changes in precipitation can contribute to the growth of insect populations (Bale et al. 2002; Amarasekare and Savage 2012; Régnière et al. 2012; Schneider et al. 2021;Wang et al. 2021). Suboptimal tree growth conditions and reduced tree vigor can also lead to outbreaks of forest insect attacks due to changes in temperature and precipitation (Waring & O’Hara 2005). Therefore, climatic variability (i.e., shifts in temperature and precipitation) can result in high dark diversity and species pools.We found that the effects of climatic PC1 on the pest species diversity indices studied were higher than those of PC2, indicating that temperature changes should be monitored to prevent invasions of and otherwise control insect pests in alpine forest ecosystems. Insect pests can benefit from increasing minimum winter temperatures(Schneider et al. 2021), with a particular risk of invasions likely due to the effects of temperature changes on dark diversity.

    Our results show that the effects of human influence on pest species diversity were significantly positive, indicating that it was the main driver in alpine forest ecosystems. Disturbance can promote invasion success for some insects in native habitats, and human influence might be associated with both climatic niche emptiness and expansion (Eschtruth and Battles 2008; González-Moreno et al. 2015; Hill et al. 2017, 2019). Anthropogenic changes to the landscape create niche opportunities for insect pests (Cariveau and Winfree 2015; Hill et al. 2017). Small patches can lead to novel environments for insect pests within a broader landscape (Cariveau and Winfree 2015). The novel niche hypothesis suggests that increased environmental heterogeneity may favor the expansion of pest species and coexistence with their hosts and competitors via niche differences(Lu et al. 2007; Cariveau and Winfree 2015; Hill et al.2017). Hence, human influence could drive the high observed diversity of insect pests. Our study provides evidence that human influence leads to high levels of dark diversity and species pools of insect pests in alpine forest ecosystems due to novel niches under environmental changes. Human activities (e.g., agricultural development and urbanization) could increase the invasion risk of insect pests in alpine forest ecosystems by impacting dark diversity and species pools (Sanderson et al. 2002; Gallardo et al. 2015;Hill et al. 2019).

    The concept of insect pest community completeness could contribute to biodiversity and ecosystem conservation in alpine forest ecosystems (P?rtel et al.2013; Ronk et al. 2017). The community completeness of insect pests was low in the forest ecosystems of Qinghai Province, China, because of changes in temperature, precipitation, and human influence.Community completeness is an interface between observed and dark diversities within the species pool(P?rtel et al. 2013). Low levels of community completeness of insect pests could indicate increasing pest invasion risk in alpine forest ecosystems (P?rtel et al.2013; Ronk et al. 2017). Positive species interactions and suitable abiotic factors favored high dark diversity and species pools; therefore, alpine forest ecosystems were exposed to insect pest invasion (Lewis et al.2017; P?rtel et al. 2013; Ronk et al. 2017). A community completeness index could be an effective indicator for macroscale risk assessments of insect pests in alpine forest ecosystems (P?rtel et al. 2013; Ronk et al. 2017). High invasion risk could occur widely in areas with low community completeness for insect pests (P?rtel et al. 2013; Ronk et al. 2017). Our results indicate that the effect of long-term changes in temperature, precipitation, and human influence on insect pests should be monitored in alpine forest ecosystems.

    Conclusions

    Our study provides new evidence for the determinants of insect pest species assemblages in alpine forest ecosystems from the perspectives of pest species interactions and abiotic factors. Species interactions shape species assemblages due to the positive effect sizes of species co-occurrence between insect pests and the high number of association links. Positive species interactions determined the species assemblages of insect pests in alpine forest ecosystems. Furthermore, climate factors and human influence were the abiotic determinants of species assemblages because high temperatures, precipitation, and human influence increased the observed diversity, dark diversity, and species pools in alpine forest ecosystems. Community completeness could be an effective indicator for insect pest risk assessments in alpine forest ecosystems due to the positive relationships between human influence and community completeness.This study will be useful for supporting the prevention and control of insect pests and provides new information for species assemblages in alpine forest ecosystems.

    Supplementary Information

    The online version contains supplementary material available at https://doi.org/10.1186/s40663-021-00351-7.

    Additional file 1:.

    Acknowledgements

    Many thanks for the helpful comments from the editor and reviewers on the improvement of our manuscript. We thank the conservation staff from the government bodies, universities, and research institutes who have made the database on insect pest species in Qinghai Province, China, and have conducted the projects of the first and second forestry pest survey of Qinghai Province, China. We acknowledge Pablo A Marquet’s Lab for supporting our work.

    Authors’contributions

    Chun-Jing Wang, Rong Wang and Ji-Zhong Wan conceived the ideas and designed the study; Rong Wang and Chun-Mei Yu collected and managed the data; Chun-Mei Yu, Yongcuo Pubu, Xiao-Peng Dang, Wan-Gui Sun,and Qiang-Feng Liconducted the fieldwork;Chun-Jing Wang, Rong Wang,and Ji-Zhong Wan collaborated with the statistical analysis and interpretation of data; Chun-Jing Wang and Ji-Zhong Wan wrote the first version of the manuscript with substantial contribution from Rong Wang. All authors contributed to subsequent drafts and gave final approval for publication.

    Funding

    This work has been supported by the National Natural Science Foundation of China(Nos. 31800449 and 31800464), and the project of the third forestry pest survey of Qinghai Province, China.

    Availability of data and materials

    All environmental data used in this manuscript are publicly available and sources described in the text with the exception of field investigation data on insect pests, which are maintained by Forest Pest Control and Quarantine Station of Qinghai Province, China with compliance of confidentiality laws in China. All the resulting data are available in Supporting Information.

    Declarations

    Ethics approval and consent to participate

    All the authors have approved the manuscript and agreed with submission to your esteemed journal.

    Consent for publication

    Not applicable.

    Competing interests

    There are no conflicts of interest to declare.

    Author details

    1State Key Laboratory of Plateau Ecology and Agriculture, College of Agriculture and Animal Husbandry, Qinghai University, Xining 810016, China.2Forestry and Grassland Planning Institute of Qinghai Province, Xining 810016, China.3Forest Pest Control and Quarantine Station of Qinghai Province, Xining 810016, China.4Forestry Station of Yushu City, Yushu 815000, China.

    Received: 14 July 2021 Accepted: 14 October 2021

    日韩一区二区视频免费看| 男女床上黄色一级片免费看| 美女脱内裤让男人舔精品视频| 91国产中文字幕| 亚洲av福利一区| 国产1区2区3区精品| 国产片特级美女逼逼视频| 亚洲国产成人一精品久久久| 国产不卡av网站在线观看| 狠狠婷婷综合久久久久久88av| 两性夫妻黄色片| 中文字幕色久视频| 啦啦啦中文免费视频观看日本| 女性被躁到高潮视频| 国产成人欧美| 亚洲精品美女久久久久99蜜臀 | 亚洲av日韩在线播放| 亚洲欧洲国产日韩| 国产1区2区3区精品| 建设人人有责人人尽责人人享有的| 青春草视频在线免费观看| 亚洲成色77777| 男人操女人黄网站| 国产99久久九九免费精品| 亚洲熟女毛片儿| 精品一品国产午夜福利视频| 国产精品久久久久久人妻精品电影 | 亚洲国产欧美日韩在线播放| 深夜精品福利| 欧美日韩视频精品一区| 黑人欧美特级aaaaaa片| 免费看av在线观看网站| 午夜免费观看性视频| 美女脱内裤让男人舔精品视频| 日韩精品有码人妻一区| 综合色丁香网| 欧美亚洲日本最大视频资源| 国产欧美亚洲国产| 观看av在线不卡| 国产精品免费大片| 国产免费视频播放在线视频| 国产精品一二三区在线看| 中文字幕亚洲精品专区| 午夜免费男女啪啪视频观看| 亚洲精品aⅴ在线观看| 亚洲精品,欧美精品| 菩萨蛮人人尽说江南好唐韦庄| 最近中文字幕2019免费版| 欧美日韩亚洲综合一区二区三区_| 国产有黄有色有爽视频| 看免费成人av毛片| 日韩一区二区视频免费看| 精品亚洲乱码少妇综合久久| 啦啦啦啦在线视频资源| 国产亚洲一区二区精品| 一边摸一边抽搐一进一出视频| 黑人欧美特级aaaaaa片| 免费观看a级毛片全部| 亚洲,欧美,日韩| 最近手机中文字幕大全| 成人亚洲欧美一区二区av| 晚上一个人看的免费电影| 汤姆久久久久久久影院中文字幕| 王馨瑶露胸无遮挡在线观看| 日韩电影二区| 9191精品国产免费久久| 男的添女的下面高潮视频| 国产爽快片一区二区三区| av网站在线播放免费| 亚洲av日韩精品久久久久久密 | 不卡av一区二区三区| 观看av在线不卡| 赤兔流量卡办理| 日韩精品有码人妻一区| 在线亚洲精品国产二区图片欧美| 久久久久精品国产欧美久久久 | 在现免费观看毛片| 欧美国产精品一级二级三级| 99国产综合亚洲精品| 免费黄频网站在线观看国产| 一边摸一边做爽爽视频免费| 无限看片的www在线观看| 女人精品久久久久毛片| 日韩大码丰满熟妇| 十八禁网站网址无遮挡| 久久久久久人人人人人| 日本91视频免费播放| xxxhd国产人妻xxx| 国产精品无大码| 国产成人一区二区在线| 亚洲欧美精品综合一区二区三区| 国产精品人妻久久久影院| 欧美xxⅹ黑人| 免费在线观看视频国产中文字幕亚洲 | 久久人人97超碰香蕉20202| 亚洲视频免费观看视频| 青春草亚洲视频在线观看| 欧美日韩亚洲高清精品| 美女脱内裤让男人舔精品视频| 亚洲国产精品一区三区| 日本欧美视频一区| √禁漫天堂资源中文www| 国产乱人偷精品视频| 精品国产乱码久久久久久小说| 国产一卡二卡三卡精品 | 夫妻午夜视频| 日韩大片免费观看网站| av免费观看日本| 男女边摸边吃奶| 免费高清在线观看日韩| 丝瓜视频免费看黄片| 亚洲精品成人av观看孕妇| 亚洲精品av麻豆狂野| 日韩成人av中文字幕在线观看| 日韩精品有码人妻一区| 91国产中文字幕| 国产成人精品无人区| 国产探花极品一区二区| avwww免费| 欧美 日韩 精品 国产| 日韩欧美精品免费久久| 国产精品免费大片| 亚洲熟女精品中文字幕| 中文字幕人妻丝袜制服| 国产日韩欧美亚洲二区| 777久久人妻少妇嫩草av网站| 国产av国产精品国产| 国产精品欧美亚洲77777| 丝袜人妻中文字幕| 国产高清不卡午夜福利| 午夜福利,免费看| 在线观看免费视频网站a站| 精品视频人人做人人爽| 一级黄片播放器| av女优亚洲男人天堂| 天天躁日日躁夜夜躁夜夜| 精品少妇久久久久久888优播| 熟女少妇亚洲综合色aaa.| 亚洲av福利一区| 日韩av免费高清视频| 亚洲综合色网址| 国产成人91sexporn| 国产成人91sexporn| 一区二区三区乱码不卡18| 亚洲欧美日韩另类电影网站| 午夜福利,免费看| 亚洲第一区二区三区不卡| 国产一区亚洲一区在线观看| 亚洲四区av| 亚洲专区中文字幕在线 | 一本色道久久久久久精品综合| 国产亚洲最大av| 男女国产视频网站| 亚洲欧洲日产国产| 亚洲成人手机| 久久精品亚洲熟妇少妇任你| 一级黄片播放器| 老司机靠b影院| 婷婷成人精品国产| 九九爱精品视频在线观看| 美国免费a级毛片| 欧美日韩福利视频一区二区| 欧美在线一区亚洲| 成人手机av| 乱人伦中国视频| 国产精品久久久久久精品古装| 女人高潮潮喷娇喘18禁视频| 久久精品国产亚洲av涩爱| 视频在线观看一区二区三区| 男女高潮啪啪啪动态图| 9热在线视频观看99| 国产亚洲av片在线观看秒播厂| 大香蕉久久网| 别揉我奶头~嗯~啊~动态视频 | 国产在线视频一区二区| 欧美日韩视频精品一区| 一二三四中文在线观看免费高清| 黄色 视频免费看| 午夜免费男女啪啪视频观看| 国产黄色免费在线视频| 黄片播放在线免费| 啦啦啦啦在线视频资源| 午夜福利一区二区在线看| 国产免费福利视频在线观看| e午夜精品久久久久久久| 亚洲 国产 在线| 97超级碰碰碰精品色视频在线观看| 两个人视频免费观看高清| 久久久久久久久久久久大奶| 成人亚洲精品av一区二区| 在线十欧美十亚洲十日本专区| 国产精品久久久av美女十八| 欧美黑人欧美精品刺激| 亚洲人成电影观看| 老司机午夜福利在线观看视频| 夜夜躁狠狠躁天天躁| 欧美成人午夜精品| 久久国产乱子伦精品免费另类| 日韩成人在线观看一区二区三区| 精品久久久久久久毛片微露脸| 中文字幕人妻熟女乱码| 亚洲国产精品久久男人天堂| 激情在线观看视频在线高清| 色在线成人网| 国产黄a三级三级三级人| 亚洲成av人片免费观看| 天堂影院成人在线观看| 看黄色毛片网站| 波多野结衣高清无吗| 大陆偷拍与自拍| 国产欧美日韩综合在线一区二区| 窝窝影院91人妻| 国产乱人伦免费视频| 国产主播在线观看一区二区| 国产伦人伦偷精品视频| 精品一区二区三区四区五区乱码| 精品久久久精品久久久| 成人三级黄色视频| 国产成人精品久久二区二区91| 三级毛片av免费| 我的亚洲天堂| 成人国产一区最新在线观看| 国产成人精品在线电影| 免费在线观看黄色视频的| 黄片播放在线免费| ponron亚洲| 久久香蕉国产精品| 一区二区三区精品91| 国产麻豆69| 欧美另类亚洲清纯唯美| 男男h啪啪无遮挡| 国产三级在线视频| 丝袜美腿诱惑在线| 青草久久国产| 精品无人区乱码1区二区| 精品国产乱子伦一区二区三区| 一级毛片高清免费大全| 国产精品亚洲一级av第二区| 不卡一级毛片| 久久久国产欧美日韩av| 人人妻人人澡欧美一区二区 | 午夜福利免费观看在线| 国产av又大| 我的亚洲天堂| 免费在线观看视频国产中文字幕亚洲| 久久 成人 亚洲| 99香蕉大伊视频| 老司机福利观看| 精品国产乱子伦一区二区三区| 美女国产高潮福利片在线看| 国产av一区在线观看免费| 成人欧美大片| 91在线观看av| 欧美国产精品va在线观看不卡| 欧美激情久久久久久爽电影 | 久久天堂一区二区三区四区| 宅男免费午夜| 国产精品一区二区三区四区久久 | 老司机午夜福利在线观看视频| 国产高清videossex| 欧美日本视频| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 99香蕉大伊视频| 久久久久精品国产欧美久久久| 老熟妇乱子伦视频在线观看| 十八禁网站免费在线| 精品熟女少妇八av免费久了| 国产亚洲精品一区二区www| 午夜免费观看网址| 成熟少妇高潮喷水视频| 国产视频一区二区在线看| 婷婷六月久久综合丁香| 国产精华一区二区三区| 国产精品野战在线观看| 一a级毛片在线观看| 视频在线观看一区二区三区| 两个人看的免费小视频| 国产三级在线视频| cao死你这个sao货| 九色国产91popny在线| 此物有八面人人有两片| 亚洲人成电影观看| 亚洲午夜精品一区,二区,三区| 国产亚洲欧美在线一区二区| 日韩欧美国产在线观看| 午夜精品国产一区二区电影| 18禁国产床啪视频网站| 国产亚洲欧美98| 久久精品亚洲精品国产色婷小说| 99re在线观看精品视频| 免费看十八禁软件| 色综合婷婷激情| 日韩成人在线观看一区二区三区| 欧美成人一区二区免费高清观看 | √禁漫天堂资源中文www| 777久久人妻少妇嫩草av网站| 丰满的人妻完整版| 如日韩欧美国产精品一区二区三区| 欧美精品啪啪一区二区三区| 一二三四社区在线视频社区8| 精品国产一区二区久久| 欧美乱色亚洲激情| 亚洲成人国产一区在线观看| 国内精品久久久久久久电影| 色老头精品视频在线观看| 国产av一区二区精品久久| 高清毛片免费观看视频网站| 国产精品国产高清国产av| 两个人免费观看高清视频| 女性生殖器流出的白浆| av片东京热男人的天堂| 麻豆成人av在线观看| 中文字幕人妻丝袜一区二区| 这个男人来自地球电影免费观看| 久久久久久亚洲精品国产蜜桃av| 欧美日本中文国产一区发布| 精品久久久久久,| 亚洲精品美女久久久久99蜜臀| 欧美乱妇无乱码| 亚洲av电影不卡..在线观看| 黑人欧美特级aaaaaa片| www日本在线高清视频| 亚洲电影在线观看av| 999精品在线视频| 男人舔女人下体高潮全视频| 久久中文字幕一级| 久久久精品国产亚洲av高清涩受| 成人国语在线视频| 色播在线永久视频| 老汉色∧v一级毛片| 国产xxxxx性猛交| av视频免费观看在线观看| 黄色视频,在线免费观看| 极品教师在线免费播放| a级毛片在线看网站| 一级作爱视频免费观看| 精品日产1卡2卡| 精品国内亚洲2022精品成人| 真人做人爱边吃奶动态| 色老头精品视频在线观看| 变态另类成人亚洲欧美熟女 | 国产精品影院久久| 色综合亚洲欧美另类图片| 亚洲片人在线观看| 啦啦啦 在线观看视频| 久久天堂一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 在线观看午夜福利视频| 女人被狂操c到高潮| 久久精品国产99精品国产亚洲性色 | 在线av久久热| av福利片在线| 母亲3免费完整高清在线观看| 亚洲国产精品成人综合色| 一级片免费观看大全| 中出人妻视频一区二区| 丝袜美足系列| 久久久久国产精品人妻aⅴ院| 神马国产精品三级电影在线观看 | 午夜福利欧美成人| 久久久久久大精品| 欧美乱妇无乱码| 精品久久久久久成人av| av电影中文网址| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产精品999在线| 啦啦啦 在线观看视频| 1024视频免费在线观看| 校园春色视频在线观看| 9色porny在线观看| 精品不卡国产一区二区三区| 久久狼人影院| 老汉色av国产亚洲站长工具| 亚洲中文av在线| a在线观看视频网站| 桃色一区二区三区在线观看| 99精品在免费线老司机午夜| av免费在线观看网站| 日韩欧美国产一区二区入口| 亚洲第一青青草原| 女性生殖器流出的白浆| 欧美色欧美亚洲另类二区 | 妹子高潮喷水视频| 国产高清激情床上av| 久久人妻熟女aⅴ| 少妇的丰满在线观看| 真人做人爱边吃奶动态| 国产午夜福利久久久久久| 天堂√8在线中文| 精品欧美一区二区三区在线| 国产精品综合久久久久久久免费 | 99久久综合精品五月天人人| 日本免费一区二区三区高清不卡 | 一本综合久久免费| 欧美黑人欧美精品刺激| 亚洲国产中文字幕在线视频| 人人妻人人澡欧美一区二区 | 12—13女人毛片做爰片一| 久久人妻熟女aⅴ| 精品久久久精品久久久| 99热只有精品国产| 夜夜夜夜夜久久久久| 亚洲欧美日韩另类电影网站| 日韩精品免费视频一区二区三区| 久久久精品国产亚洲av高清涩受| 亚洲专区中文字幕在线| 高清毛片免费观看视频网站| 9191精品国产免费久久| 久久中文看片网| 欧美丝袜亚洲另类 | 啦啦啦观看免费观看视频高清 | 嫩草影院精品99| 亚洲成av人片免费观看| 老司机在亚洲福利影院| 精品国产美女av久久久久小说| 亚洲九九香蕉| √禁漫天堂资源中文www| 国内久久婷婷六月综合欲色啪| 岛国在线观看网站| 一级a爱片免费观看的视频| 亚洲黑人精品在线| 亚洲精品国产一区二区精华液| 亚洲第一电影网av| 亚洲全国av大片| 午夜亚洲福利在线播放| 精品第一国产精品| 欧美不卡视频在线免费观看 | 久久精品aⅴ一区二区三区四区| 亚洲七黄色美女视频| 欧美最黄视频在线播放免费| 最好的美女福利视频网| 大陆偷拍与自拍| 99精品久久久久人妻精品| 黄色a级毛片大全视频| 一边摸一边抽搐一进一出视频| 女生性感内裤真人,穿戴方法视频| 操出白浆在线播放| 热re99久久国产66热| 色老头精品视频在线观看| 一区二区日韩欧美中文字幕| av网站免费在线观看视频| 国产单亲对白刺激| 国产精品免费视频内射| avwww免费| 一区在线观看完整版| 色综合婷婷激情| 国产三级在线视频| 91国产中文字幕| 国内久久婷婷六月综合欲色啪| 国语自产精品视频在线第100页| 美女高潮喷水抽搐中文字幕| 淫妇啪啪啪对白视频| 人人妻人人澡欧美一区二区 | 一级,二级,三级黄色视频| 午夜久久久久精精品| 日韩成人在线观看一区二区三区| 99精品久久久久人妻精品| 久久婷婷成人综合色麻豆| 中国美女看黄片| 不卡av一区二区三区| cao死你这个sao货| 超碰成人久久| 欧美激情高清一区二区三区| 欧美日韩瑟瑟在线播放| 国内毛片毛片毛片毛片毛片| 免费一级毛片在线播放高清视频 | 成人18禁高潮啪啪吃奶动态图| 免费少妇av软件| 99久久精品国产亚洲精品| 亚洲欧美日韩无卡精品| 精品无人区乱码1区二区| 91成人精品电影| 久久久久国内视频| 亚洲成人免费电影在线观看| 国产成人影院久久av| 久久精品亚洲熟妇少妇任你| 免费搜索国产男女视频| 女性被躁到高潮视频| 精品久久久久久久毛片微露脸| 午夜久久久久精精品| 精品一品国产午夜福利视频| 人妻久久中文字幕网| 嫩草影视91久久| 久久久国产成人精品二区| 久久久久亚洲av毛片大全| 久久精品成人免费网站| 国产精品一区二区精品视频观看| 亚洲七黄色美女视频| 亚洲欧美精品综合一区二区三区| 国产精品久久久久久人妻精品电影| 亚洲中文日韩欧美视频| 叶爱在线成人免费视频播放| 麻豆一二三区av精品| 国产成年人精品一区二区| 久久国产亚洲av麻豆专区| 欧美黑人欧美精品刺激| 日韩欧美国产在线观看| 天堂动漫精品| 黄片播放在线免费| 午夜福利,免费看| 制服人妻中文乱码| 亚洲 国产 在线| av视频在线观看入口| 伦理电影免费视频| 亚洲专区中文字幕在线| 日韩成人在线观看一区二区三区| 一边摸一边抽搐一进一小说| 精品国内亚洲2022精品成人| 精品国产一区二区三区四区第35| 中文字幕色久视频| 悠悠久久av| 极品教师在线免费播放| 久久青草综合色| 午夜影院日韩av| 国内毛片毛片毛片毛片毛片| 多毛熟女@视频| 国产精品综合久久久久久久免费 | 中亚洲国语对白在线视频| 欧美最黄视频在线播放免费| 精品人妻1区二区| 久久久久久久午夜电影| 久久久精品欧美日韩精品| 首页视频小说图片口味搜索| 亚洲专区国产一区二区| 久久久国产精品麻豆| 精品久久久久久成人av| 每晚都被弄得嗷嗷叫到高潮| 香蕉丝袜av| 国产精品一区二区在线不卡| 此物有八面人人有两片| 色精品久久人妻99蜜桃| 色综合婷婷激情| 国产真人三级小视频在线观看| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 美国免费a级毛片| 欧美激情极品国产一区二区三区| 国产片内射在线| 久久国产精品男人的天堂亚洲| 欧美黄色淫秽网站| 97碰自拍视频| 又大又爽又粗| 不卡av一区二区三区| 久久中文字幕一级| 精品国产国语对白av| 国产欧美日韩一区二区三| 欧美丝袜亚洲另类 | 国产亚洲欧美在线一区二区| 精品欧美一区二区三区在线| 老汉色∧v一级毛片| 国内毛片毛片毛片毛片毛片| 国内精品久久久久精免费| 午夜福利欧美成人| 午夜老司机福利片| 久久香蕉精品热| 少妇裸体淫交视频免费看高清 | 99riav亚洲国产免费| 午夜两性在线视频| 午夜免费激情av| 激情在线观看视频在线高清| 免费看美女性在线毛片视频| 女同久久另类99精品国产91| 午夜福利视频1000在线观看 | 国产成人精品久久二区二区91| 国产亚洲精品一区二区www| 国产精品二区激情视频| 色播亚洲综合网| 久久久久久久精品吃奶| 在线播放国产精品三级| 久久狼人影院| 无人区码免费观看不卡| 国产亚洲精品久久久久久毛片| 欧美日韩亚洲综合一区二区三区_| 欧美国产日韩亚洲一区| 亚洲少妇的诱惑av| 国产精品 国内视频| 中文字幕高清在线视频| 国产私拍福利视频在线观看| 欧美丝袜亚洲另类 | 叶爱在线成人免费视频播放| 高清在线国产一区| 亚洲 国产 在线| 中亚洲国语对白在线视频| 男男h啪啪无遮挡| 久久狼人影院| 一卡2卡三卡四卡精品乱码亚洲| 久久婷婷人人爽人人干人人爱 | 一边摸一边抽搐一进一出视频| 日韩av在线大香蕉| 亚洲自偷自拍图片 自拍| 日日干狠狠操夜夜爽| 色婷婷久久久亚洲欧美| 欧美成人性av电影在线观看| 亚洲欧美精品综合一区二区三区| 国产三级黄色录像| 91在线观看av| 日韩成人在线观看一区二区三区| 美女扒开内裤让男人捅视频| 宅男免费午夜| 国产成人av激情在线播放| 国产午夜精品久久久久久| 欧美不卡视频在线免费观看 | 亚洲成人精品中文字幕电影| 久久精品国产99精品国产亚洲性色 | 一级毛片高清免费大全| 嫩草影院精品99| 在线十欧美十亚洲十日本专区| 久久精品影院6| 9热在线视频观看99| 视频在线观看一区二区三区| 女人高潮潮喷娇喘18禁视频| 欧美成人免费av一区二区三区| 变态另类丝袜制服| 在线观看舔阴道视频|