首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Freshwater biodiversity has shown to be highly vulnerable to climate warming, alpine cold stenotherm populations being especially at risk of getting extinct. This paper aims at identifying the environmental factors favouring cold stenotherm species in alpine ponds. This information is required to provide management recommendations for habitats restoration or creation, needed for the mitigation of the effects of climate warming on alpine freshwater biodiversity. Cold stenotherm species richness as well as total (i.e. stenotherm and eurytherm) richness were analyzed for aquatic plants, Coleoptera and Odonata in 26 subalpine and alpine ponds from Switzerland and were related to environmental factors ecologically relevant for pond biodiversity. Our results confirmed that the set of environmental variables governing pond biodiversity in alpine or subalpine ponds is specific to altitude. Altitude and macrophyte presence were important drivers of cold stenotherm and total species richness, whereas connectivity did not show any significant relation. Therefore, the management of pond biodiversity has to be ‘altitude-specific’. Nevertheless, cold stenotherm species from the investigated alpine ponds do not show some specific requirements if compared to the other species inhabiting these ponds. Therefore, both total and cold stenotherm species richness could be favoured by the same management measures.  相似文献   

2.
玛曲县植被覆被变化及其对环境要素的响应   总被引:5,自引:0,他引:5  
植被覆被变化是气象要素和人类活动综合作用的结果,能够反映区域内生态系统的演替趋势。玛曲高寒生态区作为黄河上游重要的水源涵养和补给区,具有维持区域生物多样性和生态安全,保障经济社会健康发展的重要作用,因此厘清该区域植被变化与气候及人类活动等环境要素的相互关系有助于为玛曲县生态治理与恢复提供科学参考。鉴于此,以2000-2015年MODIS/NDVI数据为基础,结合同期气象与人类活动数据,应用趋势分析法、相关分析以及通径分析等方法,分析了玛曲县植被NDVI的时空变化规律,并详细探讨了气象要素和人类活动对植被覆被变化的影响。结果表明:(1)2000-2015年玛曲县NDVI呈波动上升趋势,上升速率为0.01015/10a;各土地利用/覆被类型中,增加幅度由大到小依次为高山稀疏植被、湿地、沙化草甸、山地疏林地、高寒草甸、亚高山硬叶灌丛、亚高山阔叶灌丛和高寒草原;增加面积占相应地类总面积比例由大到小分别是高山稀疏植被(75.57%)、山地疏林地(71.45%)、沙化草甸(71.18%)、湿地(70.66%)、高寒草甸(68.15%)、亚高山硬叶灌丛(66.96%)、亚高山阔叶灌丛(66.24%)和高寒草原(66.05%);(2)气象要素中,气温与NDVI间具有显著正相关(P < 0.05),是影响植被覆被的决定性因子,利于植被的生长与发育;降水与NDVI间相关不显著(P > 0.05),对植被覆被的影响较小;(3)人类活动要素中,与放牧强度密切相关的大牲畜存栏数和羊存栏数是植被生长的主控因子,其中大牲畜存栏数呈显著的抑制作用(P < 0.05),羊存栏数具有较强的促进作用(P < 0.05);(4)通径分析发现,气温、大牲畜存栏数和羊存栏数的决定系数依次为0.3005,-0.0563和0.0128,说明气温对NDVI的综合作用强度最高、大牲畜存栏数次之,羊存栏数最低;此外,剩余通径系数为0.53。该数值较大,表明仍有部分对NDVI增加存在影响的环境要素未考虑到,需在今后的研究中给予关注。  相似文献   

3.
The transitions between ecosystems (ecotones) are often biodiversity hotspots, but we know little about the forces that shape them. Today, often sharp boundaries with low diversity are found between terrestrial and aquatic ecosystems. This has been attributed to environmental factors that hamper succession. However, ecosystem properties are often controlled by both bottom-up and top-down forces, but their relative importance in shaping riparian boundaries is not known. We hypothesize that (1) herbivores may enforce sharp transitions between terrestrial and aquatic ecosystems by inhibiting emergent vegetation expansion and reducing the width of the transition zone and (2) the vegetation expansion, diversity, and species turnover are related to abiotic factors in the absence of herbivores, but not in their presence. We tested these hypotheses in 50 paired grazed and ungrazed plots spread over ten wetlands, during two years. Excluding grazers increased vegetation expansion, cover, biomass, and species richness. In ungrazed plots, vegetation cover was negatively related to water depth, whereas plant species richness was negatively related to the vegetation N:P ratio. The presence of (mainly aquatic) herbivores overruled the effect of water depth on vegetation cover increase but did not interact with vegetation N:P ratio. Increased local extinction in the presence of herbivores explained the negative effect of herbivores on species richness, as local colonization rates were unaffected by grazing. We conclude that (aquatic) herbivores can strongly inhibit expansion of the riparian vegetation and reduce vegetation diversity over a range of environmental conditions. Consequently, herbivores enforce sharp boundaries between terrestrial and aquatic ecosystems.  相似文献   

4.
Fluctuations of local but connected populations may show correlation or synchrony whenever they experience significant dispersal or correlated environmental biotic and abiotic variability. Synchrony may be an important variable in multispecies systems, but its nature and implications have not been explicitly examined. Because the number of locally coexisting species (richness) affects the population variability of community members, we manipulated richness under different regimes of environmental fluctuation (EF). We predicted that the temporal synchrony of populations in a species should decline with increasing richness of the metacommunity they live in. Additionally, we predicted that specialist species that are sensitive to a specific environmental factor would show higher synchronization when EF increases. We thus created experimental communities with varied richness, EF, and species specialization to examine the synchronizing effects of these factors on three aquatic invertebrate species. We created four levels of richness and three levels of EF by manipulating the salinity of the culture media. Monocultures exhibited higher population synchrony than metacommunities of 2–4 species. Furthermore, we found that species responded differently to EF treatments: high EF enhanced population synchrony for the specialist and intermediate species, but not for the generalist species. Our findings emphasize that the magnitude of EF and species richness both contribute to determine population synchrony, and importantly, our results suggest that biotic diversity may actually stabilize metacommunities by disrupting synchrony.  相似文献   

5.
Fire is known to facilitate the invasion of many non-native plant species, but how invasion into burnt areas varies along environmental gradients is not well-understood. We used two pre-existing data sets to analyse patterns of invasion by non-native plant species into burnt areas along gradients of topography, soil and vegetation structure in Yosemite National Park, California, USA. A total of 46 non-native species (all herbaceous) were recorded in the two data sets. They occurred in all seven of the major plant formations in the park, but were least common in subalpine and upper montane conifer forests. There was no significant difference in species richness or cover of non-natives between burnt and unburnt areas for either data set, and environmental gradients had a stronger effect on patterns of non-native species distribution, abundance and species composition than burning. Cover and species richness of non-natives had significant positive correlations with slope (steepness) and herbaceous cover, while species richness had significant negative correlations with elevation, the number of years post-burn, and cover of woody vegetation. Non-native species comprised a relatively minor component of the vegetation in both burnt and unburnt areas in Yosemite (percentage species = 4%, mean cover < 6.0%), and those species that did occur in burnt areas tended not to persist over time. The results indicate that in many western montane ecosystems, fire alone will not necessarily result in increased rates of invasion into burnt areas. However, it would be premature to conclude that non-native species could not affect post-fire succession patterns in these systems. Short fire-return intervals and high fire severity coupled with increased propagule pressure from areas used heavily by humans could still lead to high rates of invasion, establishment and spread even in highly protected areas such as Yosemite.  相似文献   

6.
该研究通过对甘肃省兴隆山自然保护区进行野外实地调查和相关文献资料收集,依据海拔梯度将保护区划分为6个植被垂直带[草原带(1800~2000 m)、山地灌丛带(2000~2200 m)、亚高山针叶林带(2200~2900 m)、亚高山矮林带(2900~3000 m)、高山灌丛带(3000~3500 m)、高山草甸带(>3500 m)],在整理6个植被垂直带植物名录的基础上,以种子植物为研究对象,对各植被垂直带的种子植物丰富度、生活型、区系成分和系统发育结构进行分析,探讨植物多样性沿植被垂直带海拔升高的垂直变化规律,以揭示植物对环境的生态适应性,为山地植物多样性保护与开发利用提供理论依据。结果表明:(1)按最新的分类系统划分,兴隆山自然保护区种子植物隶属于87科387属889种,植物科、属、种丰富度随植被垂直带海拔升高呈单峰分布格局,在亚高山针叶林带达到峰值(81科304属661种);各植被垂直带间的Jaccard相似性系数呈中等不相似和极不相似水平,海拔越靠近的植被垂直带间相似性系数越高。(2)保护区内种子植物不同生活型的垂直变化格局存在差异,木本植物所占比例沿植被垂直带海拔升高呈先升后降的变化趋势,而草本植物呈相反的变化格局,且各植被垂直带中草本植物所占比例始终高于木本植物。(3)保护区植物在属水平上,热带成分所占比例随植被带海拔升高呈下降的变化趋势,而温带成分占比呈上升的变化趋势。(4)系统发育结构在中低海拔区域的亚高山针叶林带呈发散型,在高海拔区域(>2900 m)的3个植被带中呈聚集型,说明兴隆山自然保护区非随机分布格局在群落构建机制中发挥主要作用。  相似文献   

7.
Although there is much research showing a strong negative effect of habitat fragmentation and deterioration on the viability of different insect populations and on species richness, the effect of fragmentation is modified by other local and landscape factors. One of the most substantial gaps in knowledge is whether species are similar in their response to the same environmental factors and if their response mirrors response of the entire community. From the conservation point of view this knowledge is of primary importance in planning conservation actions, yet these studies are rare. In this paper we test the relative effects of habitat patch and landscape characteristics on butterflies inhabiting calcareous grasslands in southern Poland. Butterfly species richness and abundance were positively affected by patch size and wind shelter. In the case of species richness there was also a positive effect of plant species richness. Butterfly diversity was enhanced in wind sheltered patches, and commonness (non-rarity) enhanced by distance to buildings and by shorter vegetation. Multivariate analysis suggested differences in the responses of individual species to the examined environmental variables, with some species more responsive to patch size and shelter and others to sward height. The conservation of butterfly communities requires sensible and complex management to ensure high habitat diversity. The most important challenge for future studies on calcareous grasslands is to formulate a model of management that guarantees high species richness and conservation of each individual species.  相似文献   

8.
We fill a gap in understanding wetland vegetation diversity and relationship with environmental determinants in Bulgarian high mountains. A total of 615 phytosociological samples were taken from springs, mires, wet meadows and tall-forb habitats throughout Bulgaria, of which 234 relevés are from mire and spring vegetation above timberline. The vegetation was classified by TWINSPAN and the resulting vegetation types were reproduced by the formal definitions using the combination of Cocktail species groups based on phi-coefficient of joint co-occurrence of the species. Nine vegetation types of springs and fens have been clearly delimited above the timberline. All vegetation types include Balkan endemic species, the representation of which varies. Fens generally harbour more Balkan endemics than do springs, with the exception of species-poor high-altitude Drepanocladetum exannulati. The gradient structure of the vegetation was revealed by DCA and by CCA with forward selection of environmental factors. The major determinants of vegetation variation strongly differ above and below the timberline and likewise between springs and fens. The base-richness gradient controls the floristic variation of Bulgarian submontane fens, whereas the complete data set including both submontane and subalpine fens is governed by the altitude gradient from lowland and basin fens to subalpine fens rich in Balkan endemics. When focusing on sites above the timberline only, the first DCA axis separates fens from springs without organic matter. The major species turnover in springs follows the variation in water pH and mineral content in water, whereas fen vegetation variation is primarily controlled by succession gradient of peat accumulation. Altitude remains an important factor in all cases. Weak correlation between water pH and conductivity was found. This correlation was even statistically insignificant in fens above the timberline. Water pH is not influenced by mineral richness in Bulgarian high mountains, since it is buffered by decomposition of organic matter in fens. In springs, pH reaches maximum values due to strong aeration caused by water flow. The plant species richness decreases significantly with increasing altitude. The increase of species richness towards circumneutral pH, often found in mires, was not confirmed in Bulgarian high mountains. The correlation between species richness and pH was significant only when arctic-alpine species and allied European high-mountain species were considered separately. The richness of boreal species was independent on pH. Some of them had their optima shifted to more acidic fens as compared to regions below the timberline. Our results suggest that subalpine spring and fen vegetation should be analysed separately with respect to vegetation-environment correlations. Separate analysis of fens below and above timberline is quite appropriate.  相似文献   

9.
10.
In order to clarify how vegetation types change along the environmental gradients in a cool temperate to sub-alpine mountainous zone and the determinant factors that define plant species richness, we established 360 plots (each 4 × 10 m) within which the vegetation type, species richness, elevation, topographic position index (TPI), slope inclination, and ground light index (GLI) of the natural vegetation were surveyed. Mean elevation, TPI, slope inclination, and GLI differed across vegetation types. Tree species richness was negatively correlated with elevation, whereas fern and herb species richness were positively correlated. Tree species richness was greater in the upper slope area than the lower slope area, whereas fern and herb species richness were greater in the lower slope area. Ferns and trees species richness were smaller in the open canopy, whereas herb species richness was greater in the open canopy. Vegetation types were determined firstly by elevation and secondary by topographic configurations, such as topographic position, and slope inclination. Elevation and topography were the most important factors affecting plant richness, but the most influential variables differed among plant life-form groups. Moreover, the species richness responses to these environmental gradients greatly differed among ferns, herbs, and trees.  相似文献   

11.
Geographic variation in species richness has been explained by different theories such as energy, productivity, energy–water balance, habitat heterogeneity, and freezing tolerance. This study determines which of these theories best account for gradients of breeding bird richness in China. In addition, we develop a best-fit model to account for the relationship between breeding bird richness and environment in China. Breeding bird species richness in 207 localities (3271 km2 per locality on average) from across China was related to thirteen environmental variables after accounting for sampling area. The Akaike's information criterion (AIC) was used to evaluate model performance. We used Moran's I to determine the magnitude of spatial autocorrelation in model residuals, and used simultaneous autoregressive model to determine coefficients of determination and AIC of explanatory variables after accounting for residual spatial autocorrelation. Of all environmental variables examined, normalized difference vegetation index, a measure of plant productivity, is the best variable to explain the variance in breeding bird richness. We found that species richness of breeding birds at the scale examined is best predicted by a combination of plant productivity, elevation range, seasonal variation in potential evapotranspiration, and mean annual temperature. These variables explained 47.3% of the variance in breeding bird richness after accounting for sampling area; most of the explained variance in richness is attributable to the first two of the four variables.  相似文献   

12.
13.
Nutrient enrichment weakens the stabilizing effect of species richness   总被引:2,自引:0,他引:2  
With global freshwater biodiversity declining at an even faster rate than in the most disturbed terrestrial ecosystems, understanding the effects of changing environmental conditions on relationships between biodiversity and the variability of community and population processes in aquatic ecosystems is of significant interest. Evidence is accumulating that biodiversity loss results in more variable communities; however, the mechanisms underlying this effect have been the subject of considerable debate. We manipulated species richness and nutrients in outdoor aquatic microcosms composed of naturally occurring assemblages of zooplankton and benthic invertebrates to determine how the relationship between species richness and variability might change under different nutrient conditions. Temporal variability of populations and communities decreased with increasing species richness in low nutrient microcosms. In contrast, we found no relationship between species richness and either population or community variability in nutrient enriched microcosms. Of the different mechanisms we investigated (e.g. overyielding, statistical averaging, insurance effects, and the stabilizing effect of species richness on populations) the only one that was consistent with our results was that increases in species richness led to more stable community abundances through the stabilizing effect of species richness on the component populations. While we cannot conclusively determine the mechanism(s) by which species richness stabilized populations, our results suggest that more complete resource-use in the more species-rich low nutrient communities may have dampened population fluctuations.  相似文献   

14.
Assessing how environmental changes affect the distribution and dynamics of vegetation and animal populations is becoming increasingly important for terrestrial ecologists to enable better predictions of the effects of global warming, biodiversity reduction or habitat degradation. The ability to predict ecological responses has often been hampered by our rather limited understanding of trophic interactions. Indeed, it has proven difficult to discern direct and indirect effects of environmental change on animal populations owing to limited information about vegetation at large temporal and spatial scales. The rapidly increasing use of the Normalized Difference Vegetation Index (NDVI) in ecological studies has recently changed this situation. Here, we review the use of the NDVI in recent ecological studies and outline its possible key role in future research of environmental change in an ecosystem context.  相似文献   

15.
The niche-based argument that species are filtered from environments in which they cannot sustain viable populations is the basis of the Richness-Heterogeneity Relationship (RHR). However, the multi-dimensionality of niches suggests that the RHR may take different shapes along different environmental axes, with potential confounding effects if filtering along the axes is not equally strong. Here, we explore how different structural and landscape variables drive the RHR as the accumulative outcome of environmental preferences at the species-level while considering the intercorrelation between heterogeneity levels along three niche axes. We used occurrence data of avifauna from 226 sites situated along a grassland-to-woodland gradient in a Midwestern USA study area. In each site, we quantified horizontal (habitat cover type), vertical (vegetation height structure), and spatial (habitat configuration) heterogeneity and explored the shape of the observed RHR at the landscape scale, as well as the correlations among heterogeneity levels at different axes. We then fitted species distribution models to environmental variables from the three axes separately and compared the stacked probabilities of occurrences of all species to the observed species richness. We found that predictions of richness patterns improved when more than one heterogeneity axis was included in RHR models, and that habitat suitability along different axes is not equally strong. Furthermore, a unimodal RHR along the vegetation height axis, which the species distribution models revealed to be a weak predictor for most species, may arise through intercorrelation with heterogeneity along the two other axes, along which we recorded stronger signals of environmental preference at the species level. Our results emphasize the importance of selecting relevant niche axes in studies of species richness patterns because ultimately, these patterns reflect the various environmental preferences of individual species.  相似文献   

16.
The relative involvement of larval dietary tolerance to the leaf-litter toxic polyphenols in shaping population genetic structure of the subalpine mosquito Aedes rusticus was examined. This was compared with other parameters such as geographical range, type of vegetation surrounding the breeding site, and occurrence of annual larvicidal treatments. Population genetic structure was analysed at 10 presumed neutral polymorphic isoenzyme loci. Toxicological comparisons involved standard bioassays performed on larvae fed on toxic decomposed leaf litter. Significant overall genetic differentiation was observed among the 22 studied populations and within the five defined geographical groups. Analysis of molecular variance revealed an absence of relation between genetic and environmental parameters, genetic variance being essentially found within populations. This suggested that the larval dietary tolerance to the toxic leaf litter and the other studied parameters poorly influence population genetic structure. The local adaptation of subalpine mosquito populations to the surrounding vegetation thus appears as a labile trait. Such a dynamic adaptation is also suggested by the correlation between geographical and toxicological distances and the correlation between dietary tolerance to the leaf-litter toxic polyphenols and annual larvicidal treatments.  相似文献   

17.
Many ecological hypotheses have been widely used to explain species richness variation across the globe. We investigated lizard species richness patterns in China, and identified areas of high species richness. Furthermore, we tested hypotheses concerning the relationships between lizard richness and environmental variables. A large data including 30,902 records of point locality data for 151 lizard species occurring in China were retrieved from Herpetology museums of CIB/CAS and other museums through HerpNET, and published sources, and then predicted distributions maps were generated using ecological niche modeling. We overlaid all species prediction maps into a composite map to describe species richness patterns. A multiple regression analysis using eigenvector-based spatial filtering (SEVM) was performed to examine the best environmental predictors of species richness. Richness peaked mainly in southern China located in the Oriental realm. Our best multiple regression models explained a total of 80.1% variance of lizard richness (r2 = 0.801; F = 203.47; P < 0.001). Among related factors in shaping species richness distribution, the best environmental predictors of species richness were: frost-day frequency, elevation, vegetation, and wet-day frequency. Based on models selection, our results revealed that underlying mechanisms related to different ecological hypotheses might work together and best explain lizard richness in China. We are in an initial step to develop a large data set on species richness, and provide the necessary conservation implications from habitat loss. Additional studies that test species richness at different geographical scale are required to better understand the factors that may influence the species richness distribution in East Asia.  相似文献   

18.
Genetic markers provide potentially sensitive indicators of changes in environmental conditions because the genetic constitution of populations is normally altered well before populations become extinct. Genetic indicators in populations include overall genetic diversity, genetic changes in traits measured at the phenotypic level, and evolution at specific loci under selection. While overall genetic diversity has rarely been successfully related to environmental conditions, genetically based changes in traits have now been linked to the presence of toxins and both local and global temperature shifts. Candidate loci for monitoring stressors are emerging from information on how specific genes influence traits, and from screens of random loci across environmental gradients. Drosophila research suggests that chromosomal regions under recent intense selection can be identified from patterns of molecular variation and a high frequency of transposable element insertions. Allele frequency changes at candidate loci have been linked to pesticides, pollutants and climate change. Nevertheless, there are challenges in interpreting allele frequencies in populations, particularly when a large number of loci control a trait and when interactions between alleles influence trait expression. To meet these challenges, population samples should be collected for longitudinal studies, and experimental programmes should be undertaken to link variation at candidate genes to ecological processes.  相似文献   

19.
Wang D  Jiang Z  Shen Z  Wang H  Wang B  Shou W  Zheng H  Chu X  Shi J  Huang W 《PloS one》2011,6(10):e24900
Variations in the activities of Cytochrome P450s are one of the major factors responsible for inter-individual differences in drug clearance rates, which may cause serious toxicity or inefficacy of therapeutic drugs. Various mRNA level is one of the key factors for different activity of the major P450 genes. Although both genetic and environmental regulators of P450 gene expression have been widely investigated, few studies have evaluated the functional importance of cis- and trans-regulatory factors and environmental factors in the modulation of inter-individual expression variations of the P450 genes. In this study, we measured the mRNA levels of seven major P450 genes (CYP1A1, CYP1A2, CYP2C9, CYP2C19, CYP2D6, CYP3A4 and CYP3A5) in 96 liver biopsy samples from Chinese population. Both trans-acting (mRNA levels and non-synonymous SNPs of putative regulator genes) and cis-acting (gene copy number and functional SNPs) factors were investigated to identify the determinants of the expression variations of these seven P450 genes. We found that expression variations of most P450 genes, regulator genes and housekeeping genes were positively correlated at the mRNA level. After partial correlation analysis using ACTB and GAPDH expression to eliminate the effect of global regulators, a UPGMA (Unweighted Pair Group Method with Arithmetic Mean) tree was constructed to reveal the effects of specific regulation networks potentially masked by global regulators. Combined with the functional analysis of regulators, our results suggested that expression variation at the mRNA level was mediated by several factors in a gene-specific manner. Cis-acting genetic variants might play key roles in the expression variation of CYP2D6 and CYP3A5, environmental inducers might play key roles in CYP1A1 and CYP1A2 variation and global regulators might play key roles in CYP2C9 variation. In addition, the functions of regulators that play less important roles in controlling expression variation for each P450 gene were determined.  相似文献   

20.
Mexico has higher mammalian diversity than expected for its size and geographic position. High environmental hetero geneity throughout Mexico is hypothesized to promote high turnover rates (β‐diversity), thus contributing more to observed species richness and composition than within‐habitat (α) diversity. This is true if species are strongly associated with their environments, such that changes in environmental attributes will result in changes in species composition. Also, greater heterogeneity in an area will result in greater species richness. This hypothesis has been deemed false for bats, as their ability to fly would reduce opportunities for habitat specialization. If so, we would expect no significant relationships between 1) species composition and environmental variables, 2) species richness and environmental heterogeneity, 3) β‐diversity and environmental heterogeneity. We tested these predictions using 31 bat assemblages distributed across Mexico. Using variance partitioning we evaluated the relative contribution of vegetation, climate, elevation, horizontal heterogeneity (a variate including vegetation, climate, and elevational heterogeneity), spatial variation (lat‐long), and vertical hetero geneity (of vegetation strata) to variation in bat species composition and richness. Variation in vegetation explained 92% of the variation in species composition and was correlated with all other variables examined, indicating that bats respond directly to habitat composition and structure. Beta‐diversity and vegetational heterogeneity were significantly correlated. Bat species richness was significantly correlated with vertical, but not horizontal, heterogeneity. Nonetheless, neither horizontal nor vertical heterogeneity were random; both were related to latitude and to elevation. Variation in bat community composition and richness in Mexico were primarily explained by local landscape heterogeneity and environmental factors. Significant relationships between β‐diversity and environmental variation reveal differences in habitat specialization by bats, and explain their high diversity in Mexico. Understanding mechanisms acting along environmental or geographic gradients is as important for understanding spatial variation in community composition as studying mechanisms that operate at local scales.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号