首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 46 毫秒
1.
Given species inventories of all sites in a planning area, integer programming or heuristic algorithms can prioritize sites in terms of the site's complementary value, that is, the ability of the site to complement (add unrepresented species to) other sites prioritized for conservation. The utility of these procedures is limited because distributions of species are typically available only as coarse atlases or range maps, whereas conservation planners need to prioritize relatively small sites. If such coarse‐resolution information can be used to identify small sites that efficiently represent species (i.e., downscaled), then such data can be useful for conservation planning. We develop and test a new type of surrogate for biodiversity, which we call downscaled complementarity. In this approach, complementarity values from large cells are downscaled to small cells, using statistical methods or simple map overlays. We illustrate our approach for birds in Spain by building models at coarse scale (50 × 50 km atlas of European birds, and global range maps of birds interpreted at the same 50 × 50 km grid size), using this model to predict complementary value for 10 × 10 km cells in Spain, and testing how well‐prioritized cells represented bird distributions in an independent bird atlas of those 10 × 10 km cells. Downscaled complementarity was about 63–77% as effective as having full knowledge of the 10‐km atlas data in its ability to improve on random selection of sites. Downscaled complementarity has relatively low data acquisition cost and meets representation goals well compared with other surrogates currently in use. Our study justifies additional tests to determine whether downscaled complementarity is an effective surrogate for other regions and taxa, and at spatial resolution finer than 10 × 10 km cells. Until such tests have been completed, we caution against assuming that any surrogate can reliably prioritize sites for species representation.  相似文献   

2.
3.
Regions of mediterranean-type climate represent most extra-tropical biodiversity hotspots, being both highly diverse and highly endangered. Though renowned for their plant richness, these regions’ insects constitute the bulk of their alpha diversities. Data on insect distribution and rarity are generally lacking for such regions, and are often considered unattainable. Intensive field inventories combined with statistical extrapolation methods can provide a sufficient understanding of alpha, beta, and gamma diversity components for application to conservation planning. These are essential to assessing the adequacy of a regional reserve network for the conservation of insect diversity. Here the beetle faunas of three protected areas spanning three major ecoregions in the California Floristic Province were inventoried and analyzed for species richness, complementarity, and uniqueness. These surveys produced collectively nearly 1,200 species from all three sites, estimated to be about 80% of their total faunas. Diversity was highly partitioned among sites, no one site containing more than 60% of the species. Dissimilarity was moderate to high for all comparisons, and all sites contained >40% unique species. Comparison of these results with those based on species of co-occurring plants reveal incongruent species richness but congruent similarities among sites. These results provide quantitative support to the perception that mediterranean insect faunas show high spatial variability. Along with online specimen level data on distribution, rarity, and seasonality, these results will help speed the incorporation of insect data into serious conservation planning.  相似文献   

4.
Here we report that prioritizing sites in order of rarity-weighted richness (RWR) is a simple, reliable way to identify sites that represent all species in the fewest number of sites (minimum set problem) or to identify sites that represent the largest number of species within a given number of sites (maximum coverage problem). We compared the number of species represented in sites prioritized by RWR to numbers of species represented in sites prioritized by the Zonation software package for 11 datasets in which the size of individual planning units (sites) ranged from <1 ha to 2,500 km2. On average, RWR solutions were more efficient than Zonation solutions. Integer programming remains the only guaranteed way find an optimal solution, and heuristic algorithms remain superior for conservation prioritizations that consider compactness and multiple near-optimal solutions in addition to species representation. But because RWR can be implemented easily and quickly in R or a spreadsheet, it is an attractive alternative to integer programming or heuristic algorithms in some conservation prioritization contexts.  相似文献   

5.
Indicator species groups are often used as surrogates for overall biodiversity in conservation planning because inventories of multiple taxa are rare, especially in the tropics where most biodiversity is found. At coarse spatial scales most studies show congruence in the distribution of species richness and of endemic and threatened species of different species groups. At finer spatial scale levels however, cross-taxon congruence patterns are much more ambiguous. In this study we investigated cross-taxon patterns in the distribution of species richness of trees, birds and bats across four tropical forest types in a ca. 100 × 35 km area in the Northern Sierra Madre region of Luzon Island, Philippines. A non-parametric species richness estimator (Chao1) was used to compensate for differential sample sizes, sample strategies and completeness of species richness assessments. We found positive but weak congruence in the distribution of all and endemic tree and bird and tree and bat species richness across the four forest types; strong positive congruence in the distribution of all and endemic bat and bird species richness and low or negative congruence in the distribution of globally threatened species between trees, birds and bats. We also found weak cross-taxon congruence in the complementarity of pairs of forest types in species richness between trees and birds and birds and bats but strong congruence in complementarity of forest pairs between trees and bats. This study provides further evidence that congruence in the distribution of different species groups is often ambiguous at fine to moderate spatial scales. Low or ambiguous cross-taxon congruence complicates the use of indicator species and species groups as a surrogate for biodiversity in general for local systematic conservation planning.  相似文献   

6.
We examined the spatial fidelity in the pattern in species richness, abundance and composition of ants, birds, mammals and reptiles in a Eucalyptus vegetation type in Australian tropical savanna woodland. We sampled 32 sites representing intact (uncleared) vegetation, and three treatments of different clearing age (2, 12 and 18 years). We investigated whether each fauna taxon could act a surrogate for pattern in the other, and whether they responded in parallel to time since clearing. Reptiles and mammals were combined into a single group. The correlation between taxa was low for abundance (the best being between rept/mamm and ant abundance R = 0.34) and moderate for species richness (varied from 0.36 to 0.44). Mantel tests identified a moderately strong relationship between bird and rept/mamm composition (R = 0.48). Recent clearing (2–12 years) depleted both abundance and species richness of the vertebrate fauna, whereas ant abundance spiked. Ant species richness was consistent over time since clearing. ANOSIM revealed that clearing age was a strong a priori predictor of birds and rept/mamm composition (Global R = 0.48, 0.65, respectively), whereas it was not for ants (Global R = 0.20). Pair-wise comparisons suggested that for birds, rept/mamm, there was consistent changes in composition from intact woodlands to different aged clearing. This was not the case for ants; most distinct variation in composition occurred between treatments of increasing in clearing age (2–18 years R = 0.69, 12–18 years R = 0.65). The results of this study do not support the use of ants as a surrogate of vertebrates at a local scale in a cleared and intact vegetation mosaic. We acknowledge that identification of universal surrogates and indicators are difficult to find, but funding and political expediency that demand rapid solutions to conservation planning necessitate continued investigation of the merits of using surrogate species.  相似文献   

7.
Increasingly large presence‐only survey datasets are becoming available for use in conservation assessments. Potentially, these records could be used to determine spatial patterns of plant species rarity and endemism. We test the integration of a large South Korean species record database with Rabinowitz rarity classes. Rabinowitz proposed seven classes of species rarity using three variables: geographic range, habitat specificity, and local population size. We estimated the range size and local abundance of 2,215 plant species from species occurrence records and habitat specificity as the number of landcover types each species’ records were found in. We classified each species into a rarity class or as common, compared species composition by class to national lists, and mapped the spatial pattern of species richness for each rarity class. Species were classed to narrow or wide geographic ranges using 315 km, the average from a range size index of all species (Dmax), based on maximum distance between observations. There were four classes each within the narrow and wide range groups, sorted using cutoffs of local abundance and habitat specificity. Nationally listed endangered species only appeared in the narrow‐range classes, while nationally listed endemic species appeared in almost all classes. Species richness in most rarity classes was high in northeastern South Korea especially for species with narrow ranges. Policy implications. Large presence‐only surveys may be able to estimate some classes of rarity better than others, but modification to include estimates of local abundance and habitat types, could greatly increase their utility. Application of the Rabinowitz rarity framework to such surveys can extend their utility beyond species distribution models and can identify areas that need further surveys and for conservation priority. Future studies should be aware of the subjectivity of the rarity classification and that regional scale implementations of the framework may differ.  相似文献   

8.
We provide total vascular plant species counts for three 1-ha plots in deciduous, semi-deciduous and evergreen forests in central Bolivia. Species richness ranged from 297 species and 22,360 individuals/ha in the dry deciduous forest to 382 species and 31,670 individuals/ha in the evergreen forest. Orchidaceae, Pteridophyta and Leguminosae were among the most species-rich major plant groups in each plot, and Peperomia (Piperaceae), Pleurothallis (Orchidaceae) and Tillandsia (Bromeliaceae), all epiphytes, were the most species-rich genera. This dominance of a few but very diverse and/or widespread taxa contrasted with the low compositional similarity between plots. In a neotropical context, these Central Bolivian forest plots are similar in total species richness to other dry deciduous and humid montane forests, but less rich than most Amazonian forests. Nevertheless, lianas, terrestrial herbs and especially epiphytes proved to be of equal or higher species richness than most other neotropical forest inventories from which data are available. We therefore highlight the importance of non-woody life-forms (especially epiphytes and terrestrial herbs) in Andean foothill forest ecosystems in terms of species richness and numbers of individuals, representing in some cases nearly 50% of the species and more than 75% of the individuals. These figures stress the need for an increased inventory effort on non-woody plant groups in order to accurately direct conservation actions. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.  相似文献   

9.
The emerging interest in the biological and conservation significance of locally rare species prompts a number of questions about their correspondence with other categories of biodiversity, especially global rarity. Here we present an analysis of the correspondence between the distributions of globally and locally rare plants. Using biological hotspots of rarity as our framework, we evaluate the extent to which conservation of globally rare plants will act as a surrogate for conservation of locally rare taxa. Subsequently, we aim to identify gaps between rarity hotspots and protected land to guide conservation planning. We compiled distribution data for globally and locally rare plants from botanically diverse Napa County, California into a geographic information system. We then generated richness maps highlighting hotspots of global and local rarity. Following this, we overlaid the distribution of these hotspots with the distribution of protected lands to identify conservation gaps. Based on occupancy of 1 km2 grid cells, we found that over half of Napa County is occupied by at least one globally or locally rare plant. Hotspots of global and local rarity occurred in a substantially smaller portion of the county. Of these hotspots, less than 5% were classified as multi-scale hotspots, i.e. they were hotspots of global and local rarity. Although, several hotspots corresponded with the 483 km2 of protected lands in Napa County, some of the richest areas did not. Thus, our results show that there are important conservation gaps in Napa County. Furthermore, if only hotspots of global rarity are preserved, only a subset of locally rare plants will be protected. Therefore, conservation of global, local, and multi-scale hotspots needs serious consideration if the goals are to protect a larger variety of biological attributes, prevent extinction, and limit extirpation in Napa County.  相似文献   

10.
We related pteridophytes versus tree species composition to identify surrogate measures of diversity, and complementarity of seven cloud forest fragments. Forest structure, and fern and tree composition were determined in 70 (2 × 50 m) transects. Fern density (10,150–25,080 individuals/ha) differed among sites. We recorded 83 fern species in the transects. Nonparametric richness estimators indicated that more sampling effort was needed to complete fern inventories (14 more species). However, ferns recorded outside of the transects increased richness to 103 species (six more species than predicted). Twenty-eight species were unique and rare due to special habitat requirements (Diplazium expansum, Hymenophyllum hirsutum, Melpomene leptostoma, Terpsichore asplenifolia), or were at a geographical distribution edge (Diplazium plantaginifolium, Lycopodium thyoides, Pecluma consimilis, Polypodium puberulum). Correlations between fern richness and tree richness and density were not significant, but were significant between fern richness and fern density, between epiphytic fern density and tree richness and density. Tree richness is not a good surrogate for fern diversity. Only three species were recorded in all fragments (Polypodium lepidotrichum, P. longepinnulatum, P. plebeium); thus fragments pteridophytes compositions are highly complementary, but more similar for ferns than for trees. A regional conservation approach which includes many small reserves needs to focus supplementarity on patterns of tree and fern species richness.  相似文献   

11.
We studied species richness and density of rattan palms in 50 plots of 10 × 100 m2 each between 250 m and 2420 m in eight sites in Lore Lindu National Park (LLNP), Central Sulawesi. Rattans were observed in all sample sites, representing three genera and 34 species. The elevational patterns for species richness and density were humped-shaped with maxima around 1000 m. Polynomial models of second order explained 59 and 32% of species richness and density with the factor elevation, respectively. A majority of rattan species (65%) overlapped between 1000 and 1100 m elevation, while a pronounced change in the rattan flora occurred above 1100 m. Commercially important rattan species (Calamus zollingeri, C. ornatus var. celebicus, Daemonorops macroptera) were not observed above 1250 m. The change of species assemblage was significantly related to elevation (56%), followed by geographical distance (47%) and precipitation (40%). Less than 10% of LLNP is lowland forests, much of which is threatened by agricultural intensification. In contrast, montane forests are well represented in the park and high elevation forests are not subject to agricultural conversion or intensive harvesting of rattan and other forest products.  相似文献   

12.
The use of indicator taxa as biodiversity surrogates has received widespread attention in conservation planning, but remains a highly contentious issue. Here we assess biodiversity surrogacy in the two most important biomes of tropical northeastern Brazil, Atlantic forest and Caatinga. We examine the extent to which species richness is correlated among taxonomic groups, and assess relationships between the richness of individual groups, and combinations of groups and total species richness. We introduce a new approach to dealing with autocorrelation between focal taxon richness and total species richness, using standardized data such that each taxon is given equal weight. Our Atlantic forest data covered seven taxa (bryophytes, pteridophytes, trees, ants, euglosine bees, birds, and mammals; total of 768 species) sampled from 12 sites; in Caatinga it was four taxa (trees, spiders, beetles and ants; total of 184 species) from 25 sites. Our results showed that: (1) in nearly all cases the species richnesses of individual taxa were significantly correlated with each other; (2) the species richnesses of most individual taxa were significantly correlated with total species richness in both biomes; (3) only two taxa were required for excellent (R 2 > 80%) surrogacy of total species richness in both biomes; and (4) the same two taxa (trees and ants) can provide reasonable (R 2 > 60%) surrogacy for total richness in these contrasting biomes. Our findings therefore suggest that the ‘shopping basket of taxa’ required for effective biodiversity surrogacy may not only need to be very small (two taxa), but may also be very limited in composition.  相似文献   

13.
Abstract Processes acting on different spatial and temporal scales may influence local species richness. Ant communities are usually described as interactive and therefore determined by local processes. In this paper we tested two hypotheses linked to the question of why there is local variation in arboreal ant species richness in the Brazilian savanna (‘cerrado’). The hypotheses are: (i) there is a positive relationship between ant species richness and tree species richness, used as a surrogate of heterogeneity; and (ii) there is a positive relationship between ant species richness and tree density, used as a surrogate of resource availability. Arboreal ants were sampled in two cerrado sites in Brazil using baited pitfall traps and manual sampling, in quadrats of 20 m × 50 m. Ant species richness in each quadrat was used as the response variable in regression tests, using tree species richness and tree density as explanatory variables. Ant species richness responded positively to tree species richness and density. Sampling site also influenced ant species richness, and the relationship between tree density and tree species richness was also positive and significant. Tree species richness may have influenced ant species richness through three processes: (i) increasing the variety of resources and allowing the existence of a higher number of specialist species; (ii) increasing the amount of resources to generalist species; and (iii) some other unmeasured factor may have influenced both ant and tree species richness. Tree density may also have influenced ant species richness through three processes: (i) increasing the amount of resources and allowing a higher ant species richness; (ii) changing habitat conditions and dominance hierarchies in ant communities; and (iii) increasing the area and causing a species–area pattern. Processes acting on larger scales, such as disturbance, altitude and evolutionary histories, as well as sampling effect may have caused the difference between sites.  相似文献   

14.
Abstract: Little quantitative information exists about the survey effort necessary to inventory temperate bat species assemblages. We used a bootstrap resampling algorithm to estimate the number of mist net surveys required to capture individuals from 9 species at both study area and site levels using data collected in a forested watershed in northwestern California, USA, during 1996–2000. The mean number of simulated surveys required to capture individual species varied with species' rarity and ranged from 1.5 to 44.9. We retrospectively evaluated strategies to reduce required survey effort by subsampling data from 1996 to 1998 and tested the strategies in the field during 1999 and 2000. Using data from 1996 to 1998, the mean number of simulated surveys required to capture 8 out of 9 species was 26.3, but a 95% probability of capture required >61 surveys. Inventory efficiency, defined as the cumulative proportion of species detected per survey effort, improved for both the study area and individual sites by conducting surveys later in summer. We realized further improvements in study area inventory efficiency by focusing on productive sites. We found that 3 surveys conducted between 1 July and 10 September at each of 4 productive sites in this 10-km2 study area resulted in the capture of 8 species annually. Quantitative estimation of the survey effort required to assess bat species occurrence improves the ability to plan and execute reliable, efficient inventories. Results from our study should be useful for planning inventories in nearby geographical areas and similar habitat types; further, the analytical methods we used to assess effort are broadly applicable to other survey methods and taxa.  相似文献   

15.
Species richness and abundance are central in biodiversity inventories and in measuring the structure of communities. Neglecting the assessment of sampling efficiency may lead to spurious estimates of species richness and conservation value. Our aim was to examine species richness, sampling effectiveness, species-abundance distribution (SAD) and rarity of a boreal, mire-dwelling crane fly (Diptera, Tipuloidea) assemblage in western Finland. 12 Malaise traps dispersed in 4 subplots and standardized sweep net samples were used to collect adult flies from the mire. A total of 23 species and 1,569 specimens were identified. In general all species richness estimators were highly correlated and indicated rather good sampling effort. Sample completeness, expressed as percentage of observed richness divided by estimated richness, was higher for mire-dwellers (mean 75 %) than for all species (mean 63 %). Crane fly assemblages of subplots and combined data fitted best with log-series SAD. Species spatial distribution was positively correlated with average abundance. In other words, the most abundant species occurred in the most of Malaise traps. Seven mire-dwelling species greatly outnumbered (94 % of the collected specimens) all other members in the assemblage, and only one observed species was rare by several definitions (local abundance, extent of occurrence in Finland and area of occupancy). Although the studied assemblage was characterized by commonness, five of the species have threatened status in Europe south of Finland. Separate species richness estimation of all species (vagrants and occasional species included) and focal species (here mire-dwellers) is supported if ecological information is available on the taxonomic group being studied.  相似文献   

16.
As time and money is limited, explicit, cost-effective, quick, and appropriate methods are needed to assist conservation planners and managers for making quick decisions. Butterflies promise to be a good model for rapid assessment and habitat monitoring studies because they are widespread, conspicuous, and easily recognizable and they are effective indicators of forest health. We conducted a rapid assessment of butterflies at five disturbance gradient sites that varied in elevation from 900 m a.s.l. to 3500 m a.s.l. for 20 days during March–April 2010 and recorded 79 butterfly species and 1504 individuals in the Tons valley in Western Himalayas. We were able to sample approximately 77% (123 species) of the estimated species richness on continuing the sampling until July 2010. Species richness at the study site is estimated to be 159 (95% CI: 145–210) species. Diversity was highest in heterogeneous habitats and decreased towards homogeneous habitats. Unique species were highly restricted to lowest disturbed sites. Using Pearson's correlation analysis, the strongest vegetative predictors of butterfly richness were plant species richness, canopy cover, and herb and shrub density. Butterfly species richness and abundance were highly correlated with altitude, temperature, relative humidity, fire signs, and livestock abundance. We also found positive cross-taxon correlation among butterflies, moths, and beetles across sites, indicating that butterflies can be used as surrogate or indicator taxa for insect conservation. Short sampling periods providing comprehensive estimates of species richness were reliable for identifying habitats and sites with the most conservation value in the Tons valley landscape.  相似文献   

17.
Amphibians are one of the most threatened vertebrate classes, yet at the same time new species are being described every year, demonstrating that the number of existing species is grossly underestimated. In groups such as amphibians, with high extinction rates and poorly known species boundaries, DNA barcoding is a tool that can rapidly assess genetic diversity and estimate species richness for prioritizing conservation decisions. However, reliable recovery of the 5′ region of the cytochrome c oxidase subunit 1 (COI) gene is critical for the ongoing effort to gather DNA barcodes for all amphibian species. Here, we provide new PCR conditions and tested new primers that increase the efficiency of barcode recovery in amphibians. We found that a low extension temperature for PCR cycles significantly improves the efficiency of amplification for all combinations of primers. Combining low PCR extension temperature and primers AnF1 + AnR1, we were able to recover COI sequences for 100% of the species analysed (N = 161), encompassing ~15% of the species known from Brazil (representing 77 genera and 23 families), which is an important improvement over previous studies. The preliminary assessment of species diversity suggested that number of species might be underestimated by about 25%. We conclude that DNA barcoding is an efficient, simple, and standardized protocol for identifying cryptic diversity in amphibians and advocate for its use in biodiversity inventories and across widespread populations within known species.  相似文献   

18.
Fifty-five macroinvertebrate metrics were tested for their response to pond condition in 41 ponds of northwest Spain to develop a preliminary multimetric index for ecological assessment of Mediterranean flatland ponds. Stressor specific response of individual attributes to eutrophication and habitat alteration was also investigated to identify differences in the responses of metrics to single stressors and elucidate how this might affect the performance of the final index. Several combinations were tested using discrimination efficiency (25th percentile of slightly impaired sites for metrics decreasing with perturbation and 75th percentile of slightly impaired sites for metrics increasing with perturbation) and Mann–Whitney U-test with Bonferroni adjustment (P < 0.001). The final index comprised five measures (generic richness of Chironominae, generic richness of Dytiscidae + Odonata + Tanypodinae, relative richness of Chironomidae, % Macropelopini and Shannon index) and discriminated between acceptable (good) and unacceptable (moderate) conditions with more than 86% efficiency. Moreover, all the five measures included in the final index showed unidirectional responses to eutrophication, decreasing as eutrophication increased. In contrast, the effect of habitat alteration was less clear, especially in ponds in best available conditions where a vegetation belt of shrubs and trees prevented growing of macrophitic vegetation on shores and consequently associated fauna. Interestingly, none of the functional groups (e.g. % predators and % collector–gatherers) were sensitive to degradation. Handling editor: K. Martens  相似文献   

19.
Identifying nonrandom species composition patterns predicted by assembly rules has been a central theme in community ecology. Few studies have investigated the prevalence of multiple drivers on species composition patterns in small mammal assemblages in the Old World. This study investigated seasonal changes in rodent and shrew diversity in eleven savannah vegetation types in South Africa. We tested whether species composition patterns are nonrandom with respect to predictions from Diamond's assembly rules, niche limitation hypothesis and nestedness hypothesis. Species richness estimators indicated that inventories for the rodents (80%) and shrews (100%) were relatively complete. Rodent (n = 11 species) diversity and shrew (n = 5 species) diversity were highest in summer and lowest in autumn. Rodent richness was highest in the Terminalia sericea bushveld and woodlands and lowest in the Drypetes arguta sand forest, whilst shrew richness was highest in the T. sericea bushveld and woodlands and lowest in the Acacia nilotica/Dichrostachys cinerea open shrub savannah. We found no support for the predictions of competition and nestedness hypotheses and suggest that this was probably due to the high seasonal and annual variability in rodent and shrew diversity.  相似文献   

20.
The arbuscular mycorrhizal status of fifteen mangroves and one mangrove associate was investigated from 27 sites of three inundation types namely, diurnal, usual springtide and summer springtide. Roots and rhizospheric soil samples were analysed for spore density, frequency of mycorrhizal colonization and some chemical characteristics of soil. Relative abundance, frequency and spore richness of AMF were assessed at each inundation type. All the plant species except Avicennia alba exhibited mycorrhizal colonization. The study demonstrated that mycorrhizal colonization and spore density were more influenced by host plant species than tidal inundation. Forty four AMF species belonging to six genera, namely Acaulospora, Entrophospora, Gigaspora, Glomus, Sclerocystis and Scutellospora, were recorded. Glomus mosseae exhibited highest frequency at all the inundation types; Glomus fistulosum, Sclerocystis coremioides and Glomus mosseae showed highest relative abundance at sites inundated by usual springtides, summer springtides and diurnal tides, respectively. Spore richness of AMF was of the order usual springtide > diurnal > summer springtide inundated sites. The mean spore richness was 3.27. Diurnally inundated sites had the lowest concentrations of salinity, available phosphorus, exchangeable potassium, sodium and magnesium. Statistical analyses indicated that mycorrhizal frequency and AMF spore richness were significantly negatively correlated to soil salinity. Spore richness was also significantly negatively correlated to available phosphorus. The soil parameters of the usual springtide inundated sites appeared to be favourable for the existence of maximum number of AMF. Glomus mosseae was the predominant species in terms of frequency in the soils of the Sundarbans.  相似文献   

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

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