首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
  总被引:1,自引:1,他引:1  
Aim  To produce a spatial clustering of Europe on the basis of species occurrence data for the land mammal fauna.
Location  Europe defined by the following boundaries: 11°W, 32°E, 71°N, 35°N.
Methods  Presence/absence records of mammal species collected by the Societas Europaea Mammalogica with a resolution of 50 × 50 km were used in the analysis. After pre-processing, the data provide information on 124 species in 2183 grid cells. The data were clustered using the k -means and probabilistic expectation maximization (EM) clustering algorithms. The resulting geographical pattern of clusters was compared against climate variables and against an environmental stratification of Europe based on climate, geomorphology and soil characteristics (EnS).
Results  The mammalian presence/absence data divide naturally into clusters, which are highly connected spatially and most strongly determined by the small mammals with the highest grid cell incidence. The clusters reflect major physiographic and environmental features and differ significantly in the values of basic climate variables. The geographical pattern is a fair match for the EnS stratification and is robust between non-overlapping subsets of the data, such as trophic groups.
Main conclusions  The pattern of clusters is regarded as reflecting the spatial expression of biologically distinct, metacommunity-like entities influenced by deterministic forces ultimately related to the physical environment. Small mammals give the most spatially coherent clusters of any subgroup, while large mammals show stronger relationships to climate variables. The spatial pattern is mainly due to small mammals with high grid cell incidence and is robust to noise from other subsets. The results support the use of spatially resolved environmental reconstructions based on fossil mammal data, especially when based on species with the highest incidence.  相似文献   

2.
    
Aim  Models of the potential distributions of invading species have to deal with a number of issues. The key one is the high likelihood that the absence of an invading species in an area is a false absence because it may not have invaded that area yet, or that it may not have been detected. This paper develops an approach for screening pseudo-absences in a way that is logical and defensible.
Innovation  The step-wise approach involves: (1) screening environmental variables to identify those most likely to indicate conditions where the species cannot invade; (2) identifying and selecting the most likely limiting variables; (3) using these to define the limits of its invasion potential; and (4) selecting points outside these limits as true absence records for input into species distribution models.
This approach was adopted and used for the study of three prominent Hakea species in South Africa. Models with and without the false absence records were compared. Two rainfall variables and the mean minimum temperature of the coldest month were the strongest predictors of potential distributions. Models which excluded false absences predicted that more of the potential distribution would have a high invasion potential than those which included them.
Main conclusions  The approach of applying a priori knowledge can be useful in refining the potential distribution of a species by excluding pseudo-absence records which are likely to be due to the species not having invaded an area yet or being undetected. The differences between the potential distributions predicted by the different models convey more information than making a single prediction, albeit a consensus model. The robustness of this approach depends strongly on an adequate knowledge of the ecology, invasion history and current distribution of that species.  相似文献   

3.
    
Aim Our aims were to test: (1) the extent to which vascular plant associations are related in space to mammalian associations, and (2) whether the plant associations are more closely related than the mammalian associations to climate and to a published environmental stratification of Europe. Location Europe, as defined by the following boundaries: 11° W, 32° E, 71° N and 35° N. Methods The analysis is based on presence/absence records of mammal species and plant species with a resolution of 50 km × 50 km. The similarity of the overall spatial structure was tested using a partial Mantel test while controlling for the effect of geographical proximity. To further identify the main spatial components in the datasets, we used k‐means clustering and principal components analysis. The resulting geographical patterns were compared with one another, with climate variables and with the environmental stratification of Europe. Results The clustering of the plant data forms coherent areas that can be interpreted as reflections of floristic regions that are controlled to a large extent by climate and topography. In terms of the correlation between distance matrices, the relationship between plants and mammals is relatively strong. The relationships between mammals and climate, and between plants and climate, are more complex but always statistically significant. There is no evidence that the plant clusters are more closely related than the mammalian clusters to climate, although plant clusters are closer to environmental data than to climate. Main conclusions The clustering patterns of mammals and plants form groups that agree with one another in their spatial extent. The forcing of floristic patterns into coherent entities appears mainly to be caused by climatic variables (temperature, temperature range and rainfall), mediated by elevation differences. The formation of individual plant clusters is also related to species numbers and to local and regional floristic differences. The close correlation between the floristic and faunal patterns suggests that the mammal and plant distributions are controlled by the same environmental variables, although the extent to which the mammals are controlled directly by climate or through the influence of vegetation requires more detailed study.  相似文献   

4.
ABSTRACT

The katipō is an endemic New Zealand spider that was previously common in the sand dunes at New Brighton. At sites on Banks Peninsula, katipō were detected under dried seaweed on the strandline 70% of the time. However, we detected no katipō among strandlines at New Brighton after 382 sampling visits. Incorporating these results into binomial and iterative Bayesian sampling models, it appeared highly unlikely that katipō still existed at New Brighton given so many non-detection events. However, when re-visiting the site, katipō were observed in the dunes at two locations, although they were still not found on the strandline. This specific habitat may be avoided at New Brighton due to high exposure to the prevalent strong easterly winds that occur at this site. The results emphasise that sampling models that use non-detection to indicate the likelihood of species absence can be highly specific to the sampling method used.  相似文献   

5.
6.
A simple sequential test based on the succession of ‘zero's or uninfested sampling units is proposed for use in verifying zero-infestation in pest control surveys. The method is applicable efficiently to both finite and infinite populations. Problems concerning its efficiency on application are discussed.  相似文献   

7.
    
The koala (Phascolarctos cinereus) occurs in the eucalypt forests of eastern and southern Australia and is currently threatened by habitat fragmentation, climate change, sexually transmitted diseases, and low genetic variability throughout most of its range. Using data collected during the Great Koala Count (a 1-day citizen science project in the state of South Australia), we developed generalized linear mixed-effects models to predict habitat suitability across South Australia accounting for potential errors associated with the dataset. We derived spatial environmental predictors for vegetation (based on dominant species of Eucalyptus or other vegetation), topographic water features, rain, elevation, and temperature range. We also included predictors accounting for human disturbance based on transport infrastructure (sealed and unsealed roads). We generated random pseudo-absences to account for the high prevalence bias typical of citizen-collected data. We accounted for biased sampling effort along sealed and unsealed roads by including an offset for distance to transport infrastructures. The model with the highest statistical support (wAICc ∼ 1) included all variables except rain, which was highly correlated with elevation. The same model also explained the highest deviance (61.6%), resulted in high R2(m) (76.4) and R2(c) (81.0), and had a good performance according to Cohen''s κ (0.46). Cross-validation error was low (∼ 0.1). Temperature range, elevation, and rain were the best predictors of koala occurrence. Our models predict high habitat suitability in Kangaroo Island, along the Mount Lofty Ranges, and at the tips of the Eyre, Yorke and Fleurieu Peninsulas. In the highest-density region (5576 km2) of the Adelaide–Mount Lofty Ranges, a density–suitability relationship predicts a population of 113,704 (95% confidence interval: 27,685–199,723; average density = 5.0–35.8 km−2). We demonstrate the power of citizen science data for predicting species'' distributions provided that the statistical approaches applied account for some uncertainties and potential biases. A future improvement to citizen science surveys to provide better data on search effort is that smartphone apps could be activated at the start of the search. The results of our models provide preliminary ranges of habitat suitability and population size for a species for which previous data have been difficult or impossible to gather otherwise.  相似文献   

8.
Continental‐scale assessments of 21st century global impacts of climate change on biodiversity have forecasted range contractions for many species. These coarse resolution studies are, however, of limited relevance for projecting risks to biodiversity in mountain systems, where pronounced microclimatic variation could allow species to persist locally, and are ill‐suited for assessment of species‐specific threat in particular regions. Here, we assess the impacts of climate change on 2632 plant species across all major European mountain ranges, using high‐resolution (ca. 100 m) species samples and data expressing four future climate scenarios. Projected habitat loss is greater for species distributed at higher elevations; depending on the climate scenario, we find 36–55% of alpine species, 31–51% of subalpine species and 19–46% of montane species lose more than 80% of their suitable habitat by 2070–2100. While our high‐resolution analyses consistently indicate marked levels of threat to cold‐adapted mountain florae across Europe, they also reveal unequal distribution of this threat across the various mountain ranges. Impacts on florae from regions projected to undergo increased warming accompanied by decreased precipitation, such as the Pyrenees and the Eastern Austrian Alps, will likely be greater than on florae in regions where the increase in temperature is less pronounced and rainfall increases concomitantly, such as in the Norwegian Scandes and the Scottish Highlands. This suggests that change in precipitation, not only warming, plays an important role in determining the potential impacts of climate change on vegetation.  相似文献   

9.
    
Aim The imperfect detection of species may lead to erroneous conclusions about species–environment relationships. Accuracy in species detection usually requires temporal replication at sampling sites, a time‐consuming and costly monitoring scheme. Here, we applied a lower‐cost alternative based on a double‐sampling approach to incorporate the reliability of species detection into regression‐based species distribution modelling. Location Doñana National Park (south‐western Spain). Methods Using species‐specific monthly detection probabilities, we estimated the detection reliability as the probability of having detected the species given the species‐specific survey time. Such reliability estimates were used to account explicitly for data uncertainty by weighting each absence. We illustrated how this novel framework can be used to evaluate four competing hypotheses as to what constitutes primary environmental control of amphibian distribution: breeding habitat, aestivating habitat, spatial distribution of surrounding habitats and/or major ecosystems zonation. The study was conducted on six pond‐breeding amphibian species during a 4‐year period. Results Non‐detections should not be considered equivalent to real absences, as their reliability varied considerably. The occurrence of Hyla meridionalis and Triturus pygmaeus was related to a particular major ecosystem of the study area, where suitable habitat for these species seemed to be widely available. Characteristics of the breeding habitat (area and hydroperiod) were of high importance for the occurrence of Pelobates cultripes and Pleurodeles waltl. Terrestrial characteristics were the most important predictors of the occurrence of Discoglossus galganoi and Lissotriton boscai, along with spatial distribution of breeding habitats for the last species. Main conclusions We did not find a single best supported hypothesis valid for all species, which stresses the importance of multiscale and multifactor approaches. More importantly, this study shows that estimating the reliability of non‐detection records, an exercise that had been previously seen as a naïve goal in species distribution modelling, is feasible and could be promoted in future studies, at least in comparable systems.  相似文献   

10.
SUMMARY 1. The prediction of species distributions is of primary importance in ecology and conservation biology. Statistical models play an important role in this regard; however, researchers have little guidance when choosing between competing methodologies because few comparative studies have been conducted. 2. We provide a comprehensive comparison of traditional and alternative techniques for predicting species distributions using logistic regression analysis, linear discriminant analysis, classification trees and artificial neural networks to model: (1) the presence/absence of 27 fish species as a function of habitat conditions in 286 temperate lakes located in south‐central Ontario, Canada and (2) simulated data sets exhibiting deterministic, linear and non‐linear species response curves. 3. Detailed evaluation of model predictive power showed that approaches produced species models that differed in overall correct classification, specificity (i.e. ability to correctly predict species absence) and sensitivity (i.e. ability to correctly predict speciespresence) and in terms of which of the study lakes they correctly classified. Onaverage, neural networks outperformed the other modelling approaches, although all approaches predicted species presence/absence with moderate to excellent success. 4. Based on simulated non‐linear data, classification trees and neural networks greatly outperformed traditional approaches, whereas all approaches exhibited similar correct classification rates when modelling simulated linear data. 5. Detailed evaluation of model explanatory insight showed that the relative importance of the habitat variables in the species models varied among the approaches, where habitat variable importance was similar among approaches for some species and very different for others. 6. In general, differences in predictive power (both correct classification rate and identity of the lakes correctly classified) among the approaches corresponded with differences in habitat variable importance, suggesting that non‐linear modelling approaches (i.e. classification trees and neural networks) are better able to capture and model complex, non‐linear patterns found in ecological data. The results from the comparisons using simulated data further support this notion. 7. By employing parallel modelling approaches with the same set of data and focusing on comparing multiple metrics of predictive performance, researchers can begin to choose predictive models that not only provide the greatest predictive power, but also best fit the proposed application.  相似文献   

11.
Data on the occurrence of species are widely used to inform the design of reserve networks. These data contain commission errors (when a species is mistakenly thought to be present) and omission errors (when a species is mistakenly thought to be absent), and the rates of the two types of error are inversely related. Point locality data can minimize commission errors, but those obtained from museum collections are generally sparse, suffer from substantial spatial bias and contain large omission errors. Geographic ranges generate large commission errors because they assume homogenous species distributions. Predicted distribution data make explicit inferences on species occurrence and their commission and omission errors depend on model structure, on the omission of variables that determine species distribution and on data resolution. Omission errors lead to identifying networks of areas for conservation action that are smaller than required and centred on known species occurrences, thus affecting the comprehensiveness, representativeness and efficiency of selected areas. Commission errors lead to selecting areas not relevant to conservation, thus affecting the representativeness and adequacy of reserve networks. Conservation plans should include an estimation of commission and omission errors in underlying species data and explicitly use this information to influence conservation planning outcomes.  相似文献   

12.
    

Aim

There is a wealth of information on species occurrences in biodiversity data banks, albeit presence‐only, biased and scarce at fine resolutions. Moreover, fine‐resolution species maps are required in biodiversity conservation. New techniques for dealing with this kind of data have been reported to perform well. These fine‐resolution maps would be more robust if they could explain data at coarser resolutions at which species distributions are well represented. We present a new methodology for testing this hypothesis and apply it to invasive alien species (IAS).

Location

Catalonia, Spain.

Methods

We used species presence records from the Biodiversity data bank of Catalonia to model the distribution of ten IAS which, according to some recent studies, achieve their maximum distribution in the study area. To overcome problems inherent with the data, we prepared different correction treatments: three for dealing with bias and five for autocorrelation. We used the MaxEnt algorithm to generate models at 1‐km resolution for each species and treatment. Acceptable models were upscaled to 10 km and validated against independent 10 km occurrence data.

Results

Of a total of 150 models, 20 gave acceptable results at 1‐km resolution and 12 passed the cross‐scale validation test. No apparent pattern emerged, which could serve as a guide on modelling. Only four species gave models that also explained the distribution at the coarser scale.

Main conclusions

Although some techniques may apparently deliver good distribution maps for species with scarce and biased data, they need to be taken with caution. When good independent data at a coarser scale are available, cross‐scale validation can help to produce more reliable and robust maps. When no independent data are available for validation, however, new data gathering field surveys may be the only option if reliable fine‐scale resolution maps are needed.  相似文献   

13.
    
Conservation planning assessments based on species atlas data are known to select planning units containing ecotones because these areas are relatively species‐rich. However, this richness is often dependent on the presence of adjoining core habitat, so populations within these ecotones might not be viable. This suggests that atlas data may also fail to distinguish between planning units that are highly transformed by agriculture or urbanization with those from neighbouring untransformed units. These highly transformed units could also be identified as priority sites, based solely on the presence of species that require adjoining habitat patches to persist. This potential problem was investigated using bird and mammal atlas data from Swaziland and a landcover map and found that: (i) there was no correlation between planning unit species richness and proportion of natural landcover for both taxa; (ii) the priority areas that were identified for both birds and mammals were no less transformed than if the units had been chosen at random and (iii) an approach that aimed to meet conservation targets and minimize transformation levels failed to identify more viable priority areas. This third result probably arose because 4.8% of the bird species and 22% of the mammal species were recorded in only one planning unit, reducing the opportunity to choose between units when aiming to represent each species. Therefore, it is suggested that using species lists to design protected area networks at a fine spatial scale may not conserve species effectively unless population viability data are explicitly included in the analysis.  相似文献   

14.
  总被引:1,自引:0,他引:1  
Aim Nowadays, large amounts of species distribution data and software for implementing different species distribution modelling methods are freely available through the internet. As a result, methodological works that analyse the relative performance of modelling techniques, as well as those that study which species characteristics affect their performance, are necessary. We discuss three important topics that must be kept in mind when modelling species distributions, namely (i) the distinction between potential and realized distribution, (ii) the effect of the relative occurrence area of the species on the results of the evaluation of model performance, and (iii) the general inaccuracy of the predictions of the realized distribution provided by species distribution modelling methods. Location Unspecific. Methods Using some recent papers as a basis, we illustrate the three issues mentioned above and discuss the negative implications of neglecting them. Results Considering a potential‐realized distribution gradient, different modelling methods may be arranged along this gradient according to their ability to model any concept. Complex techniques may be more suitable to model the realized distribution than simple ones, which may be more appropriate to estimate the potential distribution. Comparisons among techniques must consider this scenario. The relative occurrence area of the species conditions the results of the evaluation scores, implying that models of rare species will unavoidably yield higher discrimination values. Moreover, discrimination values that are usually reported in the literature may imply considerable over or underestimations of the distribution of the species. Main conclusions It is extremely important to establish a solid conceptual and methodological framework on which the emergent field of species distribution modelling can stand and develop.  相似文献   

15.
1. Evaluating the distribution of species richness where biodiversity is high but has been insufficiently sampled is not an easy task. Species distribution modelling has become a useful approach for predicting their ranges, based on the relationships between species records and environmental variables. Overlapping predictions of individual distributions could be a useful strategy for obtaining estimates of species richness and composition in a region, but these estimates should be evaluated using a proper validation process, which compares the predicted richness values and composition with accurate data from independent sources. 2. In this study, we propose a simple approach to estimate model performance for several distributional predictions generated simultaneously. This approach is particularly suitable when species distribution modelling techniques that require only presence data are used. 3. The individual distributions for the 370 known amphibian species of Mexico were predicted using maxent to model data on their known presence (66,113 presence-only records). Distributions were subsequently overlapped to obtain a prediction of species richness. Accuracy was assessed by comparing the overall species richness values predicted for the region with observed and predicted values from 118 well-surveyed sites, each with an area of c. 100 km(2), which were identified using species accumulation curves and nonparametric estimators. 4. The derived models revealed a remarkable heterogeneity of species richness across the country, provided information about species composition per site and allowed us to obtain a measure of the spatial distribution of prediction errors. Examining the magnitude and location of model inaccuracies, as well as separately assessing errors of both commission and omission, highlights the inaccuracy of the predictions of species distribution models and the need to provide measures of uncertainty along with the model results. 5. The combination of a species distribution modelling method like maxent and species richness estimators offers a useful tool for identifying when the overall pattern provided by all model predictions might be representing the geographical patterns of species richness and composition, regardless of the particular quality or accuracy of the predictions for each individual species.  相似文献   

16.
邓浩  纪力强 《生物多样性》2008,16(1):96-102
本研究设计并实现了一个基于地理信息系统(GIS)的仅用物种已知分布点数据预测物种潜在分布地的PSDS系统.在这一系统中,通过层次聚类算法对物种已知分布点数据进行处理,减少了异常值对预测结果的影响,从而解决了环境包络模型预测结果过于乐观的问题,在物种已知分布数据较少时也能取得较好的结果.该系统实现了数据加载与导出、图层浏览与显示、生态因子分析与分布地预测、结果展示等功能,操作方便,简单易用.本文以白冠长尾雉(Syrmaticus reevesii)为例,根据4个省的少量已知分布点数据对其在国内的潜在分布地进行了预测,获得了较理想的结果,可为该物种的保护提供科学依据.  相似文献   

17.
    
This pilot study of the rare Pagoda Rock Daisy (Leucochrysum graminifolium) in the western Blue Mountains of New South Wales (Australia) proposes a simple survey method combining timed meander and grid‐cell survey design to improve the survey effort required for monitoring of species growing in remote and/or inaccessible field locations. Where Pagoda Rock Daisies were known to be present, detection time was both rapid and effective (mean of 4.9 min for each 1 ha grid). Notably, the total survey effort remained constant for all grids, even though Pagoda Rock Daisies were unevenly distributed in the landscape (approximately 17 min/ha). Ultimately, the time required to traverse the landscape was deemed to be the primary limiting factor affecting survey effort. The application of this method is not restricted to challenging locations such as cliff edges; this method could be scaled according to the landscape or organism under investigation, providing a rapid method for surveying and monitoring rare, introduced or other plants from a site‐based scale to a broader geographic area.  相似文献   

18.
    
Aim In addition to the traditionally recognized Last Glacial Maximum (LGM, 21 ka) refuge areas in the Mediterranean region, more northerly LGM distributions for temperate and boreal taxa in central and eastern Europe are increasingly being discussed based on palaeoecological and phylogeographical evidence. Our aim was to investigate the potential refuge locations using species distribution modelling to estimate the geographical distribution of suitable climatic conditions for selected rodent species during the LGM. Location Eurasia. Methods Presence/absence data for seven rodent species with range limits corresponding to the limits of temperate or boreal forest or arctic tundra were used in the analysis. We developed predictive distribution models based on the species present‐day European distributions and validated these against their present‐day Siberian ranges. The models with the best predictors of the species distributions across Siberia were projected onto LGM climate simulations to assess the distribution of climatically suitable areas. Results The best distribution models provided good predictions of the present‐day Siberian ranges of the study species. Their LGM projections showed that areas with a suitable LGM climate for the three temperate species (Apodemus flavicollis, Apodemus sylvaticus and Microtus arvalis) were largely restricted to the traditionally recognized southern refuge areas, i.e. mainly in the Mediterranean region, but also southernmost France and southern parts of the Russian Plain. In contrast, suitable climatic conditions for the two boreal species (Clethrionomys glareous and Microtus agrestis) were predicted as far north as southern England and across southern parts of central and eastern Europe eastwards into the Russian Plain. For the two arctic species (Lemmus lemmus and Microtus oeconomus), suitable climate was predicted from the Atlantic coast eastward across central Europe and into Russia. Main conclusions Our results support the idea of more northerly refuge areas in Europe, indicating that boreal species would have found suitable living conditions over much of southern central and eastern Europe and the Russian Plain. Temperate species would have primarily found suitable conditions in the traditional southern refuge areas, but interestingly also in much of the southern Russian Plain.  相似文献   

19.
For a species to be able to respond to environmental change, it must either succeed in following its optimal environmental conditions or in persisting under suboptimal conditions, but we know very little about what controls these capacities. We parameterized species distribution models (SDMs) for 135 plant species from the Algerian steppes. We interpreted low false‐positive rates as reflecting a high capacity to follow optimal environmental conditions and high false‐negative rates as a high capacity to persist under suboptimal environmental conditions. We also measured functional traits in the field and built a unique plant trait database for the North‐African steppe. For both perennial and annual species, we explored how these two capacities can be explained by species traits and whether relevant trait values reflect species strategies or biases in SDMs. We found low false‐positive rates in species with small seeds, flowers attracting specialist pollinators, and specialized distributions (among annuals and perennials), low root:shoot ratios, wide root‐systems, and large leaves (perennials only) (R2 = .52–58). We found high false‐negative rates in species with marginal environmental distribution (among annuals and perennials), small seeds, relatively deep roots, and specialized distributions (annuals) or large leaves, wide root‐systems, and monocarpic life cycle (perennials) (R2 = .38 for annuals and 0.65 for perennials). Overall, relevant traits are rarely indicative of the possible biases of SDMs, but rather reflect the species' reproductive strategy, dispersal ability, stress tolerance, and pollination strategies. Our results suggest that wide undirected dispersal in annual species and efficient resource acquisition in perennial species favor both capacities, whereas short life spans in perennial species favor persistence in suboptimal environmental conditions and flowers attracting specialist pollinators in perennial and annual species favor following optimal environmental conditions. Species that neither follow nor persist will be at risk under future environmental change.  相似文献   

20.
    
  • 1 The effect of taxonomic level on the sensitivity of bioindicators has been widely investigated in aquatic ecosystems and, to a lesser extent, in terrestrial ecosystems. However, no studies have been conducted on the sensitivity of the different taxonomic levels of soil mites, especially Gamasina, to human activities.
  • 2 The present study aimed to assess the sensitivity of different taxonomic levels of soil Gamasina mites to anthropogenic disturbances in Europe and Argentina. We arranged the data from previous projects in a hierarchical system and conducted a study to identify the critical taxonomical levels that had the highest discriminative potential between sites (Europe and Argentina) or management types (forests, grasslands, fallows, succession, recultivation and agricultural sites).
  • 3 For the Gamasina community, geographical location was by far more important than the influence of any land use type. The analysis including only the European sites demonstrated that communities belonging to sites subjected to different land uses were also significantly different.
  • 4 The species data set provided a clearer separation of sites according to both the geographical and the land‐use gradients than the genus and family data sets. The genus and, to a lesser extent, the family approach may be sufficient to elucidate the influence of great geographical differences and also of certain land uses (e.g. grasslands from the forests and arable sites).
  • 5 Species presence/absence data provided valuable information in our analyses, although the use of quantitative data yielded a clearer separation of sites.
  相似文献   

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

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