首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
Background: Species-rich Nardus stricta grasslands are a priority habitat for conservation in Europe. They typically occur on siliceous substrates and less frequently are found on calcareous bedrock.

Aims: The present paper aimed to identify the environmental factors (i.e. bedrock type, topographic, and climatic factors) that are related with community diversity and to assess if differences in plant diversity between N. stricta communities on calcareous and siliceous bedrock occur. We hypothesised that Nardus grasslands on calcareous bedrock hosted a higher vascular plant diversity than those on siliceous bedrock.

Methods: Based on 579 vegetation surveys carried out in the south-western Alps, we assessed vascular plant diversity (species richness, Shannon diversity, and Pielou’s equitability index) of species-rich Nardus grasslands and compared it between N. stricta communities on calcareous and siliceous bedrock.

Results: Elevation was identified as the main factor related to species composition, while species diversity was mostly related to mean annual precipitation and bedrock type. Species richness, Shannon diversity, and Pielou’s equitability index were higher within the communities on calcareous rather than on siliceous bedrock and a total of 89 and 34 indicator species were detected, respectively.

Conclusions: Based on our results, we suggest to protect primarily, as a habitat of priority interest, N. stricta grasslands on calcareous substrates for the higher vascular plant diversity hosted.  相似文献   

2.
Numerous hybrid and polypoid species are found within the Triticeae. It has been suggested that the H subgenome of allopolyploid Elymus (wheatgrass) species originated from diploid Hordeum (barley) species, but the role of hybridization between polyploid Elymus and Hordeum has not been studied. It is not clear whether gene flow across polyploid Hordeum and Elymus species has occurred following polyploid speciation. Answering these questions will provide new insights into the formation of these polyploid species, and the potential role of gene flow among polyploid species during polyploid evolution. In order to address these questions, disrupted meiotic cDNA1 (DMC1) data from the allopolyploid StH Elymus are analyzed together with diploid and polyploid Hordeum species. Phylogenetic analysis revealed that the H copies of DMC1 sequence in some Elymus are very close to the H copies of DMC1 sequence in some polyploid Hordeum species, indicating either that the H genome in theses Elymus and polyploid Hordeum species originated from same diploid donor or that gene flow has occurred among them. Our analysis also suggested that the H genomes in Elymus species originated from limited gene pool, while H genomes in Hordeum polyploids have originated from broad gene pools. Nucleotide diversity (π) of the DMC1 sequences on H genome from polyploid species (π = 0.02083 in Elymus, π = 0.01680 in polyploid Hordeum) is higher than that in diploid Hordeum (π = 0.01488). The estimates of Tajima''s D were significantly departure from the equilibrium neutral model at this locus in diploid Hordeum species (P<0.05), suggesting an excess of rare variants in diploid species which may not contribute to the origination of polyploids. Nucleotide diversity (π) of the DMC1 sequences in Elymus polyploid species (π = 0.02083) is higher than that in polyploid Hordeum (π = 0.01680), suggesting that the degree of relationships between two parents of a polyploid might be a factor affecting nucleotide diversity in allopolyploids.  相似文献   

3.
4.

Background

Evaluating the influence of climate variability on enteric disease incidence may improve our ability to predict how climate change may affect these diseases.

Objectives

To examine the associations between regional climate variability and enteric disease incidence in New Zealand.

Methods

Associations between monthly climate and enteric diseases (campylobacteriosis, salmonellosis, cryptosporidiosis, giardiasis) were investigated using Seasonal Auto Regressive Integrated Moving Average (SARIMA) models.

Results

No climatic factors were significantly associated with campylobacteriosis and giardiasis, with similar predictive power for univariate and multivariate models. Cryptosporidiosis was positively associated with average temperature of the previous month (β =  0.130, SE =  0.060, p <0.01) and inversely related to the Southern Oscillation Index (SOI) two months previously (β =  −0.008, SE =  0.004, p <0.05). By contrast, salmonellosis was positively associated with temperature (β  = 0.110, SE = 0.020, p<0.001) of the current month and SOI of the current (β  = 0.005, SE = 0.002, p<0.050) and previous month (β  = 0.005, SE = 0.002, p<0.05). Forecasting accuracy of the multivariate models for cryptosporidiosis and salmonellosis were significantly higher.

Conclusions

Although spatial heterogeneity in the observed patterns could not be assessed, these results suggest that temporally lagged relationships between climate variables and national communicable disease incidence data can contribute to disease prediction models and early warning systems.  相似文献   

5.
Modeling and mapping possibilities of Shannon–Wiener, Simpson, and number of species (NS) indices were researched using geographic information systems (GIS) and remote sensing (RS) tools in Nallihan forest ecosystem of Turkey. The relationships between the indices and a number of independent variables such as topography, geology, soil, climate, normalized difference vegetation index (NDVI), and land cover were investigated to understand relationships between plant diversity and ecosystem. Georeferenced field data from the established 56 quadrats (50 × 20 m) were used to calculate the indices. Principle component analysis (PCA) and multiple regression were employed for data reduction and model development, respectively. Three diversity maps were produced using the developed models. Residual maps and logical interpretations in ecological point of view were used to test the validity of the models. Elevation and climatic factors formed the most important components that are effective determinants of plant species diversity, but geological formations, soil, land cover and land-use characteristics also influenced plant diversity. Considering the different responses of the models, Shannon–Wiener (SWI) and NS models were found suitable for rare cover types, while Simpson (SIMP) model might be appropriate for single dominant land covers in the study area.  相似文献   

6.
Climate-induced phenological shifts can influence population, evolutionary, and ecological dynamics, but our understanding of these phenomena is hampered by a lack of long-term demographic data. We use a multi-decade census of 5 salmonid species representing 14 life histories in a warming Alaskan stream to address the following key questions about climate change and phenology: How consistent are temporal patterns and drivers of phenology for similar species and alternative life histories? Are shifts in phenology associated with changes in phenotypic variation? How do phenological changes influence the availability of resource subsidies? For most salmonid species, life stages, and life histories, freshwater temperature influences migration timing – migration events are occurring earlier in time (mean = 1.7 days earlier per decade over the 3–5 decades), and the number of days over which migration events occur is decreasing (mean = 1.5 days per decade). Temporal trends in migration timing were not correlated with changes in intra-annual phenotypic variation, suggesting that these components of the phenotypic distribution have responded to environmental change independently. Despite commonalities across species and life histories, there was important biocomplexity in the form of disparate shifts in migration timing and variation in the environmental factors influencing migration timing for alternative life history strategies in the same population. Overall, adult populations have been stable during these phenotypic and environmental changes (λ ≈1.0), but the temporal availability of salmon as a resource in freshwater has decreased by nearly 30 days since 1971 due to changes in the median date of migration timing and decreases in intra-annual variation in migration timing. These novel observations advance our understanding of phenological change in response to climate warming, and indicate that climate change has influenced the ecology of salmon populations, which will have important consequences for the numerous species that depend on this resource.  相似文献   

7.
Upland calcareous grassland landscapes are typically comprised of a matrix of calcareous grassland, acid grassland and limestone heath plant communities. This matrix of habitats is produced by a combination of underlying geology, climate and management. These landscapes are typically managed through grazing, with management targeted to maintain particular plant communities in the calcareous grassland habitat, whilst patches of acid grassland and limestone heath are not targeted by conservation management. The biodiversity value of acid grassland and limestone heath patches within the calcareous grassland matrix are unknown. This study provides the first assessment of their biodiversity value by examining aspects of epigeal spider diversity supported by these non-target habitat patches in comparison to calcareous grassland. Spiders were sampled in each habitat from April to August 2014 using pitfall traps across three upland regions in Great Britain. Spider species assemblages were distinct between limestone heath and both grassland types. Distinction in species assemblages are likely due to differences in vegetation structure and microclimate, e.g., humidity, degree of shade. Each habitat type supported several rare species (e.g., Jacksonella falconeri, Agyneta subtilis) revealing the contribution to spider fauna. The distinct spider species assemblage and presence of rare species in limestone heath patches demonstrate their importance in the upland calcareous grassland matrix. This study highlights the value of monitoring biodiversity in non-target habitats within a habitat matrix alongside those that are actively targeted by management.  相似文献   

8.
The Neotropics harbors a high diversity of species and several hypotheses have been proposed to account for this pattern. However, while species of forested domains are frequently studied, less is known of species from open vegetation formations occupying, altogether, a larger area than the Amazon Forest. Here we evaluate the role of historical barriers and the riverine hypothesis in the speciation patterns of small mammals by analyzing an ancient rodent lineage (Thrichomys, Hystricomorpha). Phylogenetic and biogeographic analyses were carried out with mitochondrial and nuclear DNA markers to analyze the evolutionary relationships between Thrichomys lineages occurring in dry domains along both banks of the Rio São Francisco. This river is one of the longest of South America whose course and water flow have been modified by inland tectonic activities and climate changes. Molecular data showed a higher number of lineages than previously described. The T. inermis species complex with 2n = 26, FN = 48 was observed in both banks of the river showing a paraphyletic arrangement, suggesting that river crossing had occurred, from east to west. A similar pattern was also observed for the T. apereoides complex. Thrichomys speciation occurred in Late Miocene when the river followed a different course. The current geographic distribution of Thrichomys species and their phylogenetic relationships suggested the existence of frequent past connections between both banks in the middle section of the Rio São Francisco. The extensive palaeodune region found in this area has been identified as a centre of endemism of several vertebrate species and is likely to be a center of Thrichomys diversification.  相似文献   

9.
《PloS one》2014,9(9)
Understanding why host species differ so much in symbiont loads and how this depends on ecological host and symbiont traits is a major issue in the ecology of symbiosis. A first step in this inquiry is to know whether observed differences among host species are species-specific traits or more related with host-symbiont environmental conditions. Here we analysed the repeatability (R) of the intensity and the prevalence of feather mites to partition within- and among-host species variance components. We compiled the largest dataset so far available: 119 Paleartic passerine bird species, 75,944 individual birds, ca. 1.8 million mites, seven countries, 23 study years. Several analyses and approaches were made to estimate R and adjusted repeatability (Radj) after controlling for potential confounding factors (breeding period, weather, habitat, spatial autocorrelation and researcher identity). The prevalence of feather mites was moderately repeatable (R = 0.26–0.53; Radj = 0.32–0.57); smaller values were found for intensity (R = 0.19–0.30; Radj = 0.18–0.30). These moderate repeatabilities show that prevalence and intensity of feather mites differ among species, but also that the high variation within species leads to considerable overlap among bird species. Differences in the prevalence and intensity of feather mites within bird species were small among habitats, suggesting that local factors are playing a secondary role. However, effects of local climatic conditions were partially observed for intensity.  相似文献   

10.
In temperate mountains, fens have been reported as relict habitats subject to geographical fragmentation and broad climatic gradients, but few studies have analyzed the influence of these factors on plant diversity. Here we investigate the effect of isolation on the vegetation diversity of rich fens (Caricion davallianae) in the mountains of the Iberian Peninsula, the distribution limit of these habitats in south-western Europe. We used plot-based vegetation data from the Pyrenees and the Cantabrian mountain range to evaluate their regional species-pool, occurrence of specialists, beta-diversity and the effect of geo-climatic variables on their species-richness and species-composition. We found a lower ratio of rare specialists in the Pyrenees than in the Cantabrian range, but similar estimates in the species pools, total species-richness per plot and beta-diversity. The isolation of the two mountain regions resulted in different species assemblages best predicted by summer precipitation and bedrock types, showing region-based differences in the response of vegetation and plant specialists to the environment. The tighter correlation between local climate and diversity estimates in the Cantabrian range suggests relict character of rich fens in that region, where climatic conditions have restricted local distribution of formerly more widely distributed specialists. Although there is no relevant evidence of vegetation impoverishment in that region, historical isolation has probably resulted in the existence of fragmentary plant communities. We conclude that fen vegetation may experience long-time persistence in climatically sub-optimal mountain refugia, but related plant specialists may be sensitive to climatic changes and subject to the extinction of local populations.  相似文献   

11.
The seasonality of insect assemblages in Africa is poorly investigated. To provide information on the relationships between climate and insect assemblages in the Sudanian region, strongly affected by climate change, we studied Myrmeleontidae and Ascalaphidae assemblages (Insecta: Neuroptera) for 7 and 5 consecutive years respectively in southern Mali. To make the species inventory as exhaustive as possible, we performed weekly sampling by netting and light trapping. For both assemblages, results showed very similar patterns of variation in species diversity throughout the year. Adults of Myrmeleontidae and Ascalaphidae were active all year, and the species succession was influenced by a strong temporal segregation. Species diversity peaked at the end of the rainy season and surprisingly during the dry season. Principal component analysis of the climatic factors followed by co‐inertia analysis applied to two data sets, one comprising climatic factors and the other reporting presence/absence of species, showed a good association between the annual trend of climatic factors and the species diversity. But no well defined species grouping was clearly linked to a particular period of the year. This tight association between climate and species composition suggests that even small climate changes could modify significantly species assemblage characteristics.  相似文献   

12.
Campbell Island, an isolated island 600 km south of New Zealand mainland (52°S, 169°E) is oceanic (Conrad Index of Continentality  = −5) with small differences between mean summer and winter temperatures. Previous work established the unexpected result that a mean annual climate warming of c. 0.6°C since the 1940''s has not led to upward movement of the forest limit. Here we explore the relative importance of summer and winter climatic conditions on growth and age-class structure of the treeline forming species, Dracophyllum longifolium and Dracophyllum scoparium over the second half of the 20th century. The relationship between climate and growth and establishment were evaluated using standard dendroecological methods and local climate data from a meteorological station on the island. Growth and establishment were correlated against climate variables and further evaluated within hierarchical regression models to take into account the effect of plot level variables. Winter climatic conditions exerted a greater effect on growth and establishment than summer climatic conditions. Establishment is maximized under warm (mean winter temperatures >7 °C), dry winters (total winter precipitation <400 mm). Growth, on the other hand, is adversely affected by wide winter temperature ranges and increased rainfall. The contrasting effect of winter warmth on growth and establishment suggests that winter temperature affects growth and establishment through differing mechanisms. We propose that milder winters enhance survival of seedlings and, therefore, recruitment, but increases metabolic stress on established plants, resulting in lower growth rates. Future winter warming may therefore have complex effects on plant growth and establishment globally.  相似文献   

13.
Metrics that synthesize the complex effects of climate change are essential tools for mapping future threats to biodiversity and predicting which species are likely to adapt in place to new climatic conditions, disperse and establish in areas with newly suitable climate, or face the prospect of extirpation. The most commonly used of such metrics is the velocity of climate change, which estimates the speed at which species must migrate over the earth’s surface to maintain constant climatic conditions. However, “analog-based” velocities, which represent the actual distance to where analogous climates will be found in the future, may provide contrasting results to the more common form of velocity based on local climate gradients. Additionally, whereas climatic velocity reflects the exposure of organisms to climate change, resultant biotic effects are dependent on the sensitivity of individual species as reflected in part by their climatic niche width. This has motivated development of biotic velocity, a metric which uses data on projected species range shifts to estimate the velocity at which species must move to track their climatic niche. We calculated climatic and biotic velocity for the Western Hemisphere for 1961–2100, and applied the results to example ecological and conservation planning questions, to demonstrate the potential of such analog-based metrics to provide information on broad-scale patterns of exposure and sensitivity. Geographic patterns of biotic velocity for 2954 species of birds, mammals, and amphibians differed from climatic velocity in north temperate and boreal regions. However, both biotic and climatic velocities were greatest at low latitudes, implying that threats to equatorial species arise from both the future magnitude of climatic velocities and the narrow climatic tolerances of species in these regions, which currently experience low seasonal and interannual climatic variability. Biotic and climatic velocity, by approximating lower and upper bounds on migration rates, can inform conservation of species and locally-adapted populations, respectively, and in combination with backward velocity, a function of distance to a source of colonizers adapted to a site’s future climate, can facilitate conservation of diversity at multiple scales in the face of climate change.  相似文献   

14.
Tropical forest canopies house most of the globe''s diversity, yet little is known about global patterns and drivers of canopy diversity. Here, we present models of ant species density, using climate, abundance and habitat (i.e. canopy versus litter) as predictors. Ant species density is positively associated with temperature and precipitation, and negatively (or non-significantly) associated with two metrics of seasonality, precipitation seasonality and temperature range. Ant species density was significantly higher in canopy samples, but this difference disappeared once abundance was considered. Thus, apparent differences in species density between canopy and litter samples are probably owing to differences in abundance–diversity relationships, and not differences in climate–diversity relationships. Thus, it appears that canopy and litter ant assemblages share a common abundance–diversity relationship influenced by similar but not identical climatic drivers.  相似文献   

15.
Understanding the effects of climatic variability on marine mammals is challenging due to the complexity of ecological interactions. We used general linear models to analyze a 15-year database documenting marine mammal strandings (1994–2008; n = 1,193) and nine environmental parameters known to affect marine mammal survival, from regional (sea ice) to continental scales (North Atlantic Oscillation, NAO). Stranding events were more frequent during summer and fall than other seasons, and have increased since 1994. Poor ice conditions observed during the same period may have affected marine mammals either directly, by modulating the availability of habitat for feeding and breeding activities, or indirectly, through changes in water conditions and marine productivity (krill abundance). For most species (75%, n = 6 species), a low volume of ice was correlated with increasing frequency of stranding events (e.g. R2 adj = 0.59, hooded seal, Cystophora cristata). This likely led to an increase in seal mortality during the breeding period, but also to increase habitat availability for seasonal migratory cetaceans using ice-free areas during winter. We also detected a high frequency of stranding events for mysticete species (minke whale, Balaenoptera acutorostrata) and resident species (beluga, Delphinapterus leucas), correlated with low krill abundance since 1994. Positive NAO indices were positively correlated with high frequencies of stranding events for resident and seasonal migratory cetaceans, as well as rare species (R2 adj = 0.53, 0.81 and 0.34, respectively). This contrasts with seal mass stranding numbers, which were negatively correlated with a positive NAO index. In addition, an unusual multiple species mortality event (n = 114, 62% of total annual mortality) in 2008 was caused by a harmful algal bloom. Our findings provide an empirical baseline in understanding marine mammal survival when faced with climatic variability. This is a promising step in integrating stranding records to monitor the consequences of environmental changes in marine ecosystems over long time scales.  相似文献   

16.
The biogeographical relationships between far-separated populations, in particular, those in the mainland and islands, remain unclear for widespread species in eastern Asia where the current distribution of plants was greatly influenced by the Quaternary climate. Deciduous Oriental oak (Quercus variabilis) is one of the most widely distributed species in eastern Asia. In this study, leaf material of 528 Q. variabilis trees from 50 populations across the whole distribution (Mainland China, Korea Peninsular as well as Japan, Zhoushan and Taiwan Islands) was collected, and three cpDNA intergenic spacer fragments were sequenced using universal primers. A total of 26 haplotypes were detected, and it showed a weak phylogeographical structure in eastern Asia populations at species level, however, in the central-eastern region of Mainland China, the populations had more haplotypes than those in other regions, with a significant phylogeographical structure (N ST = 0.751> G ST = 0.690, P<0.05). Q. variabilis displayed high interpopulation and low intrapopulation genetic diversity across the distribution range. Both unimodal mismatch distribution and significant negative Fu’s FS indicated a demographic expansion of Q. variabilis populations in East Asia. A fossil calibrated phylogenetic tree showed a rapid speciation during Pleistocene, with a population augment occurred in Middle Pleistocene. Both diversity patterns and ecological niche modelling indicated there could be multiple glacial refugia and possible bottleneck or founder effects occurred in the southern Japan. We dated major spatial expansion of Q. variabilis population in eastern Asia to the last glacial cycle(s), a period with sea-level fluctuations and land bridges in East China Sea as possible dispersal corridors. This study showed that geographical heterogeneity combined with climate and sea-level changes have shaped the genetic structure of this wide-ranging tree species in East Asia.  相似文献   

17.
Despite that several studies have shown that data derived from species lists generated from distribution occurrence records in the Global Biodiversity Information Facility (GBIF) are not appropriate for those ecological and biogeographic studies that require high sampling completeness, because species lists derived from GBIF are generally very incomplete, Suissa et al. (2021) generated fern species lists based on data with GBIF for 100 km × 100 km grid cells across the world, and used the data to determine fern diversity hotspots and species richness–climate relationships. We conduct an evaluation on the completeness of fern species lists derived from GBIF at the grid–cell scale and at a larger spatial scale, and determine whether fern data derived from GBIF are appropriate for studies on the relations of species composition and richness with climatic variables. We show that species sampling completeness of GBIF is low (<40%) for most of the grid cells examined, and such low sampling completeness can substantially bias the investigation of geographic and ecological patterns of species diversity and the identification of diversity hotspots. We conclude that fern species lists derived from GBIF are generally very incomplete across a wide range of spatial scales, and are not appropriate for studies that require data derived from species lists in high completeness. We present a map showing global patterns of fern species diversity based on complete or nearly complete regional fern species lists.  相似文献   

18.
Abstract. Species composition patterns and vegetation-environment relationships were quantified for high-elevation rock outcrops of the Southern Appalachian Mountains, an infrequent and insular habitat in a forested landscape. Outcrops occur over a wide geographic range encompassing extensive variation in both geology and climate. Geographic-scale factors interact with site-scale factors to produce variation in vegetation among outcrops. Similarly, site-scale factors interact with micro-scale factors to produce variation in vegetation within outcrops. To provide a quantitatively-based classification of outcrop vegetation we used a TWINSPAN analysis of 154 100-m2 plots. We recognized nine communities that primarily correspond to different combinations of elevation, bedrock type, geography, and moisture. Within outcrops of a single bedrock type, vegetation composition of 100-m2 plots was consistently correlated with elevation and solar radiation, but relationships to soil nutrients varied with bedrock type. Both site-scale (100 m2) factors (e.g. elevation, slope, aspect, and bedrock type) and plot-scale (1-m2) microsite factors (e.g. soil depth, vegetation height, soil nutrients) were strongly correlated with species composition at the 1-m2 level. Environment can be used to predict composition more effectively for 100-m2 plots on a single bedrock type than either across bedrock types or at a 1-m2 scale. Composition-environment relationships resemble those described for outcrop systems from other regions with pronounced topographic relief more than they do those described for the nearby but flatter and lower-elevation outcrops of the Southeastern Piedmont. There is strong spatial autocorrelation in this community, perhaps owing to dispersal limitation. Consequently, a comprehensive conservation strategy must include reservation of both a range of geologic types and a range of geographic locations.  相似文献   

19.
We modelled the present and future sub-Saharan winter distributions of 64 trans-Saharan migrant passerines to predict the potential impacts of climate change. These predictions used the recent ensemble modelling developments and the latest IPCC climatic simulations to account for possible methodological uncertainties. Results suggest that 37 species would face a range reduction by 2100 (16 of these by more than 50%); however, the median range size variation is −13 per cent (from −97 to +980%) under a full dispersal hypothesis. Range centroids were predicted to shift by 500±373 km. Predicted changes in range size and location were spatially structured, with species that winter in southern and eastern Africa facing larger range contractions and shifts. Predicted changes in regional species richness for these long-distance migrants are increases just south of the Sahara and on the Arabian Peninsula and major decreases in southern and eastern Africa.  相似文献   

20.
Sinopodophyllum hexandrum is an important medicinal plant whose genetic diversity must be conserved because it is endangered. The Qinling Mts. are a S. hexandrum distribution area that has unique environmental features that highly affect the evolution of the species. To provide the reference data for evolutionary and conservation studies, the genetic diversity and population structure of S. hexandrum in its overall natural distribution areas in the Qinling Mts. were investigated through inter-simple sequence repeats analysis of 32 natural populations. The 11 selected primers generated a total of 135 polymorphic bands. S. hexandrum genetic diversity was low within populations (average He = 0.0621), but higher at the species level (He = 0.1434). Clear structure and high genetic differentiation among populations were detected by using the unweighted pair group method for arithmetic averages, principle coordinate analysis and Bayesian clustering. The clustering approaches supported a division of the 32 populations into three major groups, for which analysis of molecular variance confirmed significant variation (63.27%) among populations. The genetic differentiation may have been attributed to the limited gene flow (Nm = 0.3587) in the species. Isolation by distance among populations was determined by comparing genetic distance versus geographic distance by using the Mantel test. Result was insignificant (r = 0.212, P = 0.287) at 0.05, showing that their spatial pattern and geographic locations are not correlated. Given the low within-population genetic diversity, high differentiation among populations and the increasing anthropogenic pressure on the species, in situ conservation measures were recommended to preserve S. hexandrum in Qinling Mts., and other populations must be sampled to retain as much genetic diversity of the species to achieve ex situ preservation as a supplement to in situ conservation.  相似文献   

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

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