首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
Identification of the primary factors that influence the ecological distribution of species groups is important to managers of lowland‐mountain forests in northern Iran. The aim of this study was to identify main ecological species groups, describe the site conditions associated with these species groups and the relationships between environmental factors and the distribution of ecological species groups using multi‐variate analysis (Detrended correspondence analysis (DCA) and Canonical correspondence analysis (CCA)). For this purpose, 50 relevés (400 m2 each) were sampled using the Braun‐Blanquet method. Vegetation was classified into three ecological species groups using a modified two‐way indicator species analysis (TWINSPAN). In each relevé, environmental factors (topographic and soil variables) were measured and analysed using one‐way ANOVA and Pearson r statistics. Further, species diversity indices were determined for the identified ecological species groups. Our results show that the environmental factors, e.g. elevation, slope, slope aspect, soil texture, pH and organic matter, were the most important factors explaining the distribution of the three ecological species groups in the study area. The diversity of the ecological species groups decreased with elevation. The results provide an ecological basis for forest management and for developing strategies for forest conservation in the study area.  相似文献   

2.
Knowledge about relationships between specialization degree of species, i.e. the width of their realized niche and functional traits, may have important implications for the assessment of future population developments under environmental change. In this study, we used a recently introduced method to calculate ecological niche widths of plant species in mixed broad-leaved deciduous forests and to investigate the dependence between niche widths of plants and their functional traits and Ellenberg indicator values. The research is based on a dataset of 4556 phytosociological relevés of mixed broad-leaved deciduous forests in Slovenia. We calculated theta indices for 326 species, which ranks them along a continuous gradient of habitat specialization. For 272 species, we compiled 26 functional traits and Ellenberg indicator values. We found some significant correlations between theta indices of species and their functional traits and Ellenberg indicator values; habitat specialists thrive primarily on the highest altitudes, on colder, dry sites and achieve the age of first flowering later than generalists. They also have smaller seed diameter, lower leaf dry matter content, lower mean canopy height and bigger specific leaf area than generalists. Two species groups, chamaephytes and spring green species, are particularly characterized as specialist species. The added value of our work is in complementing the knowledge about the niche differentiating along different environmental gradients and species coexistence in mixed broad-leaved deciduous forests.  相似文献   

3.
Understanding how land‐use changes affect different facets of plant biodiversity in seminatural European grasslands is of particular importance for biodiversity conservation. As conclusions of previous experimental or synchronic observational studies did not converge toward a general agreement, assessing the recent trends in vegetation change in various grassland systems using a diachronic approach is needed. In this resurvey study, we investigated the recent changes in grassland vegetation of the French Jura Mountains, a region with a long tradition of pastoralism. We compared the floristic composition of 150 grassland plots recorded between 1990 and 2000 with new relevés made in 2012 on the same plots. We considered taxonomic, phylogenetic and functional diversity as well as ecological characteristics of the plant communities derived from ecological indicator values and life strategies of the species. PCA of the floristic composition revealed a significant general trend linked to the sampling year. Wilcoxon paired tests showed that contemporary communities were generally more dominated by grass species and presented a higher tolerance to defoliation, a higher pastoral value, and a higher nutrient indicator value. Comparisons revealed a decrease in phylogenetic and functional diversity. By contrast, local species richness has slightly increased. The intensity of change in species composition, measured by Hellinger distance between pairs of relevés, was dependent on neither the time lag between the two surveys, the author of the first relevé nor its location or elevation. The most important changes were observed in grasslands that previously presented low pastoral value, low grass cover, low tolerance to defoliation, and high proportion of stress‐tolerant species. This trend was likely linked to the intensification of grassland management reported in the region, with a parallel increase in mowing frequency, grazing pressure, and fertilization level. More restrictive specifications should be applied to agricultural practices to avoid overexploitation of mountain species‐rich grasslands and its negative consequences on their biodiversity and resilience.  相似文献   

4.
In this study a new field of action for outcome-oriented payment schemes (payments by results) is introduced. That is the promotion of extensive summer grazing activities of livestock. Combining this objective with methods of outcome-oriented rewarding means to tie the payment to vegetation criteria indicating regular grazing activities instead of regularly controlling the presence of livestock on the field. In comparison to action-oriented concepts, the main advantage of this approach is that only one inspection walk is required. The paper describes and evaluates the first attempt to determine suitable indicator species for such a scheme. We pre-tested the approach with existing vegetation data only. A database of relevés, plant species and agriculturally used grassland communities of the regional district Doberan, Mecklenburg-Western Pomerania (Germany) was utilised. Grazing activities can be identified by the existence of a specific number of vascular species indicating grazing activities on a particular grassland field. We utilised the indicator function of grassland communities for management conditions to resolve the deficiency of data base relevés for providing no data about the management regime. The pasture community Lolio perennis-Cynosuretum cristati was used as a reference system to pre-test the suitability of different combinations of potential indicator species for separating pasture from non-pasture grassland. A tentative quantitative selection criterion for pasture grassland was proposed. That is the number of indicator species a relevé/grassland field has to hold to be identified as pasture grassland. Database relevés are not sufficient for calibrating such a concept, though. It is necessary to adjust the indicator species list to other drivers like use characteristics, nutrient availability, and soil moisture level. Here additional fieldwork is needed. As a first evidence for the validity of the concept, the suitability of the indicator list was measured by testing the correct identification of pasture grassland in a long-term grassland experiment built by the treatments “grazing” and “mowing”. ANOVA tests on the precision of the indicator species showed significant results. In the discussion a modification of the inventory method of grassland fields was recommended. The precision of the developed concept could be further enhanced by introducing the species dominance as a second factor deciding about the correct identification of pasture use.  相似文献   

5.
6.
Question: Which are the gradients of floristic differentiation in Greek beech (Fagus sylvatica) forests? Which is the role of geographical and ecological factors in this differentiation? Location: Beech forests of the plant geographical regions Northeast, North Central and East Central Greece. Methods: A total of 1404 published and unpublished phytoso‐ciological relevés were used in the analyses. TWINSPAN and DCA were applied to classify and ordinate the relevés. Altitude, Indicator Values of relevés and their X and Y coordinates were used in a posteriori interpretation of the ordination axes. Kendall's correlation coefficients were calculated between DCA relevé scores and explanatory variables. Multiple linear regression was used to partition the variation explained by the first two DCA axes, between the geographical and the ecological variables. Results: Classification resulted in 14 vegetation units defined by species composition. Two types of gradients, ecological and geographical, were revealed by the DCA of all releves. The partition of the variation accounted for by the first and second DCA axis was attributed mainly to ecological and geographical variables, respectively. Conclusions: Beech forests of northeast and Central Greece show phytogeographical differences, while ecologically similar vegetation units occur in both regions. A west‐east gradient is revealed in Greek beech forest vegetation. The extent of the study area, its position along regional gradients and the comprehensiveness of the data set that is analysed determine the types of the gradients which can be revealed in a vegetation study.  相似文献   

7.
Question: Is it possible to mathematically classify relevés into vegetation types on the basis of their average indicator values, including the uncertainty of the classification? Location: The Netherlands. Method: A large relevé database was used to develop a method for predicting vegetation types based on indicator values. First, each relevé was classified into a phytosociological association on the basis of its species composition. Additionally, mean indicator values for moisture, nutrients and acidity were computed for each relevé. Thus, the position of each classified relevé was obtained in a three‐dimensional space of indicator values. Fitting the data to so called Gaussian Mixture Models yielded densities of associations as a function of indicator values. Finally, these density functions were used to predict the Bayesian occurrence probabilities of associations for known indicator values. Validation of predictions was performed by using a randomly chosen half of the database for the calibration of densities and the other half for the validation of predicted associations. Results and Conclusions: With indicator values, most reléves were classified correctly into vegetation types at the association level. This was shown using confusion matrices that relate (1) the number of relevés classified into associations based on species composition to (2) those based on indicator values. Misclassified relevés belonged to ecologically similar associations. The method seems very suitable for predictive vegetation models.  相似文献   

8.
Species Distribution Models are key in modern ecological studies. They employ information about species locations and environmental factors to generate statistical functions that predict the potential distribution of species on the basis of landscape suitability. Although these models are powerful and useful tools, often the required information about species distribution is lacking, and the only resources are pre-collected museum data. Phytosociological databases contain a myriad of relevés with precious information, but are often considered to be the exclusive ownership of vegetation scientists. Our study tested the efficiency of a phytosociological database in the building of Species Distribution Models, including spatial autocorrelation (SAC) as a predictor to evaluate its effects on model performance. Spatial autocorrelation (SAC) is a natural characteristic of species distribution that depends on exogenous and endogenous processes. The latter??s effects could be overestimated by a subjective sample choice. We chose Festuca riccerii, an Italian endemic species. We split the whole dataset (671 relevés) into a calibration (443 relevés) and testing set (228 relevés) and performed a GLM on these data to identify the main ecological factors that lead distribution in order to build a Species Distribution Model. The dataset??s efficiency was assessed by testing the predicting power of the calibrated model on the testing subset. The phytosociological database proved to be good for building model (AUC?=?0.821), providing a useful basis for fast and low cost ecological analysis, and could be used subsequently for more detailed analyses.  相似文献   

9.
Abstract

The Centaureo-Portenschlagiellion describes the cliff vegetation of the Southeastern Adriatic. This habitat type harbours many endemic taxa and has a high value for biodiversity conservation. Notwithstanding its importance, knowledge of the syntaxonomy of this alliance is still poor. This article aims at revising the synchorological, coenological and floristic relationships of the associations of the Centaureo-Portenschlagiellion. The revision is based on a data set of 103 relevés of Mediterranean xerothermic cliffs from Croatia, Bosnia i Herzegovina and Montenegro. The relevés were clustered by using the Flexible beta method. An indicator species analysis was used to identify the diagnostic taxa of the main clusters of relevés and non-metric multidimensional scaling ordination was undertaken to visualize the floristic relationships among them. Results revealed that only seven associations belong to the Centaureo-Portenschlagiellion, whereas the others already described in phytosociological literature were invalidly described or should be treated as syntaxonomic synonyms of the previous ones. Two associations (Inulo-Centaureetum cuspidatae and Portenschlagiello-Campanuletum portenschlagianae) were reduced to the rank of subassociations of the Moltkio-Inuletum verbascifoliae. Two other subassociations were described for the first time. Data provided with this revision may be considered as essential base-line information that should aid in evaluating the state of this vegetation type in the future.  相似文献   

10.
Question: Species optima or indicator values are frequently used to predict environmental variables from species composition. The present study focuses on the question whether predictions can be improved by using species environmental amplitudes instead of single values representing species optima. Location: Semi‐natural, deciduous hardwood forests of northwestern Germany. Methods: Based on a data set of 558 relevés, species responses (presence/absence) to pH were modelled with Huisman‐Olff‐Fresco (HOF) regression models. Species amplitudes were derived from response curves using three different methods. To predict the pH from vegetation, a maximum amplitude overlap method was applied. For comparison, predictions resulting from several established methods, i. e. maximum likelihood/present and absent species, maximum likelihood/present species only, mean weighted averages and mean Ellenberg indicator values were calculated. The predictive success (squared Pearson's r and root mean square error of prediction) was evaluated using an independent data set of 151 relevés. Results: Predictions based upon amplitudes defined by maximum Cohen's x probability threshold yield the best results of all amplitude definitions (R2= 0.75, RMSEP = 0.52). Provided there is an even distribution of the environmental variable, amplitudes defined by predicted probability exceeding prevalence are also suitable (R2= 0.76, RMSEP = 0.55). The prediction success is comparable to maximum likelihood (present species only) and – after rescaling – to mean weighted averages. Predicted values show a good linearity to observed pH values as opposed to a curvilinear relationship of mean Ellenberg indicator values. Transformation or rescaling of the predicted values is not required. Conclusions: Species amplitudes given by a minimum and maximum boundary for each species can be used to efficiently predict environmental variables from species composition. The predictive success is superior to mean Ellenberg indicator values and comparable to mean indicator values based on species weighted averages.  相似文献   

11.
Floristic differentiation of the oriental beech (Fagus orientalis Lipsky) forests in Turkey and Bulgaria was investigated and the role of geographical and topographical factors in this differentiation was assessed. After geographical and ecological stratification of the available 922 relevés, 288 remained. Classification, by applying cluster analysis, resulted in seven vegetation units defined by species composition which represent the geographical and ecological variation of Fagus orientalis forests. DCA ordination was applied to these units by passively projecting their chorological structure, as supplementary variables. For more detailed interpretation of vegetation types with similar geographic distribution patterns, PCA was applied by passively projecting the chorological elements, life-forms and topographical factors as supplementary variables. Seven vegetation units representing the geographical and ecological variety of Fagus orientalis forests were described. Four vegetation units represent the core area of Fagus orientalis distribution on the western and middle coast of the Black Sea region (Euxine region); the remaining three types represent the distribution in the eastern Black Sea region (Colchic region), the distribution in western and southern Anatolia under the influence of the Mediterranean climate and the distribution in the transitional zone from the Euxine region to the continental parts of Inner Anatolia, respectively. The four vegetation types in Euxine region reflect the decreasing effect of Black Sea towards Inner Anatolia, as well as altitudinal differences, except the forest type representing forests on calcareous sites. The other three vegetation units represent ravine, lowland to montane and altimontane forests in Euxine region. Fagus orientalis forests could be distinguished by their floristic composition, their chorological elements and life-forms spectra, which reflect a geographical and ecological gradients.  相似文献   

12.
The present study is focused on Potamogeton polygonifolius – the species considered to be endangered in Germany, the Czech Republic, Poland and Switzerland. Field studies were carried out in Lower Silesia and Lusatia in southwestern Poland in order to verify the present state of the marginal populations of P. polygonifolius, to determine the habitat preferences of the species in the study area, to compare the ecological data between marginal and central populations and finally to recognise the factors constraining the easternmost limit of the P. polygonifolius. In total, 18 populations were examined, and 11 relevés were collected. The relevés gathered were compared with 95 relevés made in northwestern Germany, and with 10 relevés collected in the German part of Lusatia. There were no significant differences between the ecology of populations of P. polygonifolius on the eastern edge of its range and those in the central part. P. polygonifolius was found growing both in natural habitats like unspoiled peat bogs, and in anthropogenic habitats such as exploited peat bogs and drainage ditches. In both types of habitats its populations were abundant and competent in terms of generative propagation. The lack of geographical barriers as well as biology and fitness of P. polygonifolius individuals from the study area suggest that marginal populations are possibly a part of a large metapopulation, covering the whole (Polish and German) area of Lusatia, however, this hypothesis was not tested in detail. The main factor shaping the easternmost limit of the general range of P. polygonifolius is most likely climate, which prevents P. polygonifolius from spreading farther eastwards.  相似文献   

13.
Abstract. A new method of indicator species analysis is suggested that can be used to demonstrate changes in the species composition of forests resulting from atmospheric deposition of acidifying substances and nitrogen. The method is not dependent on an exact re-localization of sites. 2162 published relevés of Fagus sylvatica and Quercus robur-Carpinus betulus forests from two study areas in Germany (S Niedersachsen and Nordrhein-Westfalen) were compiled and divided into relevés from before 1975 and more recent relevés. Weighted averages of Ellenberg indicator values for light (L), moisture (M), reaction (R), and nitrogen (N) were calculated and for each of the four data sets - mesic and dry sites in two areas - N was regressed on R to give a predicted value (Npred) for any given R. The differences between observed values Nobs and Npred were calculated separately for old and new relevés (Ndev-old and Ndev-new). The same analysis was applied to model data including 16 different scenarios of acidification and/or eutrophication. The model runs resulted in six different types of regression plots for Ndev and R. All models had in common that the regression lines for Ndev-new lay, completely or largely, above the zero line ofno change and above the lines for Ndev-old. Also, in the four data sets, Ndev-new was generally larger than Ndev-old, and particularly so in Nordrhein-Westfalen. The values of Ldev-old and Ldev-new indicate that the differences in Ndev are unlikely to be caused by altered forest management, but must result from atmospheric pollution.  相似文献   

14.
The actual state of neophyte invasion in Moldavia (Eastern Romania) is described in this paper on the basis of 11,055 phytosociological relevés. We analyzed the i) proportion of relevés with neophytes, ii) mean proportion of neophytes per relevés, and iii) mean coverage of neophytes per relevé for 36 EUNIS habitat types to identify general plant invasion patterns. The level of invasion differed considerably between habitats. The invasion of neophytes especially affected habitats strongly determined or influenced by man, such as anthropogenic woodlands, ruderal habitats, arable lands or trampled areas. Most natural habitats are either slightly invaded or entirely free of neophytes. Only riverine willow stands and wet tall-herb stands are relatively highly invaded. However, the absence of neophytes in some natural habitats less represented in the phytosociological dataset could be of artifactual nature. No significant relationship between the number of neophytes and non-neophytes was found in the analysis across different habitats. When the analyses were made within-habitats, both negative and positive relationships were found, which confirm that the relationship between alien and native species richness depends on the habitats. A total number of 105 neophyte species were recorded in the phytosociological relevés used in this study. Among these, 13 species that are currently considered invasive in Moldavia occur in at least 10 types of habitats.  相似文献   

15.
In the phytosociological literature, there are numerous different approaches to the designation of diagnostic species. Frequently, this results in discrepancies between the lists of diagnostic species published for one and the same community. We examined different approaches to determining diagnostic species using as an examplePicea abies forests within the broader context of all Central European forests. Diagnostic species of spruce forests were determined from a data set of 20,164 phytosociological relevés of forests from the Eastern Alps, Western Carpathians, and the Bohemian Massif, which included 3,569 relevés of spruce forests. Phi coefficient of association was used to measure species fidelity, and species with the highest fidelities were considered as diagnostic. Diagnostic species were determined in four ways, including (A) comparison of spruce forests among the three mountain ranges, (B) comparison between spruce forests and other forests, performed separately in each of the mountain ranges, (C) simultaneous comparison of spruce forests of each of the mountain ranges with spruce forests of the other two ranges and with the other forests of all ranges, (D) comparison of spruce forests with the other forests, using pooled data sets from the three mountain ranges. The sets of diagnostic species of spruce forests yielded in comparisons A and B were sharply different; the set resulting from comparison C was intermediate between the first two and comparison D resulted in similar diagnostic species as comparison B. In comparison A, spruce forests of the Eastern Alps had a number of diagnostic species, while the spruce forests of the other two mountain ranges had only few diagnostic species. In comparison B, by contrast, the number and quality of diagnostic species decreased from the Bohemian Massif to the Eastern Alps. This exercise points out that lists of diagnostic species published in phytosociological literature are dependent on the context, i.e. the underlying data sets and comparisons: some of these lists are useful for identification of vegetation units at a local scale, some others for distinguishing units within a narrowly delimited community type over a large area. The thoughtless application of published lists of diagnostic species outside of the context for which they were intended should therefore be avoided.  相似文献   

16.
Hédl  Radim 《Plant Ecology》2004,170(2):243-265
From 1941–;1944 nearly 30 phytosociological relevés were completed by F. K. Hartmann in the Rychlebské Mountains, a typical mountainous area in northeastern Czech Republic. Of the original plots still covered with adult grown beech (Fagus sylvatica) forest, 22 were resampled in 1998 and 1999. In order to describe the recent vegetation variability of the sites 57 relevés were recorded. Changes in vegetation were estimated using relative changes in species density and ordinations (PCA, RDA). Environmental changes were assessed using Ellenberg indicator values when no direct measurements were available. A decline in species diversity has been documented, particularly, many species occurring frequently in deciduous forests with nutrient and moisture well-supplied soils around neutral have decreased. In contrast, several light-demanding, acid- and soil desiccation-tolerant species have increased. Natural succession, quantified as forest age, contributed slightly to these changes. In Ellenberg indicator values, a decline in F (soil moisture), R (soil calcium) and N (ecosystem productivity), and an increase in L (understorey light) were shown. This is interpreted as the influence of modified forestry management and of airborne pollutants. Intensified logging caused the canopy to open and soil conditions to worsen. The latter is most likely also due to acid leaching of soil cations (Ca, K, Na). This caused a decline in soil productivity, thus the effect of nitrification could not be detected. The original relevés may have differed in size influencing the results. This revised version was published online in August 2006 with corrections to the Cover Date.  相似文献   

17.
This study provides a view of vegetation types dominated by Sesleria species in the Western Carpathians (Sesleria caerulea, S. heufleriana, S. tatrae, and S. uliginosa). We also took into account characteristics/traits such as plant life forms, autochthonous status, endemism, and ploidy level occurring within each relevé in our data set. Altitude and Ellenberg indicator values derived for each relevé were considered as well. Eight vegetation types/formations/habitats were recognized in a data set of 942 phytosociological relevés: (1) synanthropic vegetation, (2) rock fissures and screes, (3) alpine grasslands and heaths, (4) mesic grasslands and pastures, (5) springs and fens, (6) forests, (7) xerophilous shrubland, and (8) high-mountain scrubs (krummholz). Results corroborated and clearly emphasized that Sesleria caerulea has the widest ecological amplitude of all studied species because the species occurred among all studied vegetation types. Sesleria tatrae was present only in several vegetation types occurring from montane to alpine vegetation belts in the highest mountains. Sesleria heufleriana and S. uliginosa were recorded only in low-altitude areas. The difference between them lies mainly in the dampness of each locality. Sesleria heufleriana was frequently found in xerophilous communities, whereas S. uliginosa preferred humid habitats of springs and fens.  相似文献   

18.
An analysis of indicator values of four Saussurea DC. species, S. baicalensis (Adams) B. L. Rob., S. frolowii Ledeb., S. salicifolia (L.) DC., and S. schanginiana (Wydl.) Fisch. ex Serg., have been carried out by the method of indicator scales of soil moisture and nutrient availability and salinity in more than 1600 relevés. The size of their ecological ranges and the values of their competitiveness in phytocenoses are determined. As a result of a comparison of the locality of local populations within the ecological range and their occurrence, the vulnerability of S. baicalensis, S. frolowii, S. salicifolia, and S. schanginiana in the Kuznetsk Alatau is evaluated. The types of rarity for these species in the Kuznetsk Alatau are determined.  相似文献   

19.
The large, comprehensive vegetation database of Mecklenburg-Vorpommern/NE Germany with 51,328 relevés allowed us to study an entire regional flora of 133 non-native plants (NNP, immigration after 1492 AD) with regard to their preferences to all kinds of habitats and along different ecological gradients. For each relevé, we computed average Ellenberg indicator values (EIV) for temperature, light, moisture, reaction, nutrients and salt as well as plant strategy type weights. We partitioned the dataset into relevés with and without occurrences of NNP and compared them with respect to the relative frequencies of EIVs and strategy type weights. We identified deviations from random differences by testing against permuted indicator values. To account for bias in EIV between community types, NNP preferences were differentiated for 34 phytosociological classes. We tested significance of preferences for the group of NNP as a whole, as well as for single NNP species within the entire dataset, as well as differentiated by phytosociological classes and formations. NNP as a group prefer communities with high EIVs for temperature and nutrients and low EIVs for moisture. They avoid communities with low EIV for reaction and high EIV for salt. NNP prefer communities with high proportions of ruderal and low proportion of stress strategists. The differentiation by phytosociological classes reinforces the general trends for temperature, nutrients, moisture, R and S strategy types. Nevertheless, preferences of single species reveal that NNP are not a congruent group but show individualistic ecological preferences.  相似文献   

20.
We compared the coenological information of the only Italian population of Malcolmia littorea (L.) R. Br. with published phytosociological relevés, including ones of this species, throughout its European range. With the aim of highlighting the main climatic features influencing the distribution patterns of M. littorea, we integrated coenological data with some climatic variables and considered major drivers for plant distribution at the European scale. Finally, we analysed the population extent of M. littorea in Italy, in order to assess its conservation status at regional level. The DCA analysis, performed on a matrix 139 relevés × 183 species, separated the relevés according to their floristic composition, showing a geographic gradient from Portugal to Mediterranean coasts, until Italy; with Mediterranean relevés clearly separated from Atlantic ones as well. Along the beach-inland gradient, the analysis highlights that the species is typical of fixed dunes habitats (more inland, mainly stabilised dunes), although in the Atlantic it can also be found in mobile dunes. The analysis of climatic variables in relation to M. littorea distribution, suggested that the species is sensitive to low winter temperature and to summer drought. The only Italian population of M. littorea is subjected to many threats, due to its small dimensions (<1 ha), isolation from the rest of its distribution area, that ranges from Portugal to France (until the Camargue region) and intensive human disturbances. Using both field and remote sensing information, we showed a considerable decrease of the occupied surface in Italy, leading us to suggest that the IUCN threat category of M. littorea in Italy should be reassessed from endangered to critically endangered.  相似文献   

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

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