首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 437 毫秒
1.
Net primary production (NPP) supplies matter, energy, and services to facilitate the sustainable development of human society and ecosystem. The response mechanism of NPP to land use and climate changes is essential for food security and biodiversity conservation but lacks a comprehensive understanding, especially in arid and semi‐arid regions. To this end, taking the middle‐reaches of the Heihe River Basin (MHRB) as an example, we uncovered the NPP responses to land use and climate changes by integrating multisource data (e.g., MOD17A3 NPP, land use, temperature, and precipitation) and multiple methods. The results showed that (a) land use intensity (LUI) increased, and climate warming and wetting promoted NPP. From 2000 to 2014, the LUI, temperature, and precipitation of MHRB increased by 1.46, 0.58°C, and 15.76 mm, respectively, resulting in an increase of 14.62 gC/m2 in annual average NPP. (b) The conversion of low‐yield cropland to forest and grassland increased NPP. Although the widespread conversion of unused land and grassland to cropland boosted both LUI and NPP, it was not conducive to ecosystem sustainability and stability due to huge water consumption and human‐appropriated NPP. Urban sprawl occupied cropland, forest, and grassland and reduced NPP. (c) Increase in temperature and precipitation generally improved NPP. The temperature decreasing <1.2°C also promoted the NPP of hardy vegetation due to the simultaneous precipitation increase. However, warming‐induced water stress compromised the NPP in arid sparse grassland and deserts. Cropland had greater NPP and NPP increase than natural vegetation due to the irrigation, fertilizers, and other artificial inputs it received. The decrease in both temperature and precipitation generally reduced NPP, but the NPP in the well‐protection or less‐disturbance areas still increased slightly.  相似文献   

2.
Species occurrence databases and climate databases were used to examine differences in patterns of species experienced climate across latitude for wide‐ranging rodents in the central‐eastern and western North America. The accumulated data were used to address three questions: (1) Do rodent species in the central and eastern region of North America select habitat at range edges to remain closer to climate conditions at the range core? (2) Is there a trend toward species having greater experienced climate variation consistent with smaller effects of orbitally controlled climate oscillations in the south vs. north? (3) How do species experienced temperature, precipitation, and elevation means and variation in this region compare to rodent species in the adjacent but more heterogeneous western North America? Results showed that central‐eastern North American species occur in as wide a range of environmental conditions as available throughout their ranges. These patterns are different from previous findings for rodents in the adjacent western USA and highlight major differences in current structure of species experienced environmental means and variation over latitude for species in spatially heterogeneous, mountainous areas vs. those that occupy flatter lands. The differences are likely important for determining differential response to climate changes.  相似文献   

3.
4.
Land‐use change and climate change are driving a global biodiversity crisis. Yet, how species' responses to climate change are correlated with their responses to land‐use change is poorly understood. Here, we assess the linkages between climate and land‐use change on birds in Neotropical forest and agriculture. Across > 300 species, we show that affiliation with drier climates is associated with an ability to persist in and colonise agriculture. Further, species shift their habitat use along a precipitation gradient: species prefer forest in drier regions, but use agriculture more in wetter zones. Finally, forest‐dependent species that avoid agriculture are most likely to experience decreases in habitable range size if current drying trends in the Neotropics continue as predicted. This linkage suggests a synergy between the primary drivers of biodiversity loss. Because they favour the same species, climate and land‐use change will likely homogenise biodiversity more severely than otherwise anticipated.  相似文献   

5.
Ancient trees are considered one of the most important habitats for biodiversity in Europe and North America. They support exceptional numbers of specialized species, including a range of rare and endangered wood‐living insects. In this study, we use a dataset of 105 sites spanning a climatic gradient along the oak range of Norway and Sweden to investigate the importance of temperature and precipitation on beetle species richness in ancient, hollow oak trees. We expected that increased summer temperature would positively influence all wood‐living beetle species whereas precipitation would be less important with a negligible or negative impact. Surprisingly, only oak‐specialist beetles with a northern distribution increased in species richness with temperature. Few specialist beetles and no generalist beetles responded to the rise of 4°C in summer as covered by our climatic gradient. The negative effect of precipitation affected more specialist species than did temperature, whereas the generalists remained unaffected. In summary, we suggest that increased summer temperature is likely to benefit a few specialist beetles within this dead wood community, but a larger number of specialists are likely to decline due to increased precipitation. In addition, generalist species will remain unaffected. To minimize adverse impacts of climate change on this important community, long‐term management plans for ancient trees are important.  相似文献   

6.
Although species distribution modelling (SDM) is widely accepted among the scientific community and is increasingly used in ecology, conservation biology and biogeography, methodological limitations generate potential problems for its application in macroecology. Using amphibian species richness in North and South America, we compare species richness patterns derived from SDM maps and ‘expert’ maps to evaluate if: 1) richness patterns derived from SDM are biased toward climate‐based explanations for diversity when compared to expert maps, since SDM methods are typically based on climatic variables; and 2) SDM is a reliable tool for generating richness maps in hyperrich regions where point occurrence data are limited for many species. We found that although three widely used SDM methods overestimated amphibian species richness in grid cells when compared to expert richness maps in both North and South America due to systematic overestimation of range sizes, diversity gradients were reasonably robust at broad scales. Further, climatic variables statistically explained patterns of richness at similar levels among the different richness sources, although climatic relationships were stronger in the much better known North America than in South America. We conclude that in the face of the high deforestation rates coupled with incomplete data on species distributions, especially in the tropics, SDM represents a useful macroecological tool for investigating broad‐scale richness patterns and the dynamics between species richness and climate.  相似文献   

7.
Drivers of biodiversity at macroscales have long been of interest in ecology, and climate and topography are now considered to be major drivers. Because humans have transformed most of the Earth's land surface, land use may play a significant role as a driver of biodiversity at a macroscale. Here we disentangle the relationships among climate, topography, land use, available energy (measured by the normalized difference vegetation index [NDVI]), and species richness of Japanese forest birds. Species richness was better explained at 40‐ and 80‐km resolutions than at 5‐, 10‐ and 20‐km resolutions; it was explained by climate, topography, and land use, and the effects of land use were fully incorporated into those of climate and topography. As temperature increased and elevation decreased, natural forest area decreased, and this decrease intensified in warm lowland areas. With the loss of natural forest, species richness decreased below a certain threshold. As temperature increased and elevation decreased, species richness and NDVI increased slightly or were unchanged in cool highland areas and decreased in warm lowland areas. Species richness increased linearly with the increase in NDVI. Most effects of climate/topography on species richness in warm lowland areas were shared by those of land use, suggesting that the decrease in species richness in warm lowland areas has been caused by loss of natural forest. Therefore, it is suggested that climate and topography determined land use intensity, which in turn, drove species richness through the depletion of available energy. Increasing temperature and decreasing elevation leads to both benefits (increase in potential available energy) and costs (depletion of energy by human land‐use change) for forest birds. These costs seem to override benefits in warm lowland areas.  相似文献   

8.
Climate and land‐use change jointly affect the future of biodiversity. Yet, biodiversity scenarios have so far concentrated on climatic effects because forecasts of land use are rarely available at appropriate spatial and thematic scales. Agent‐based models (ABMs) represent a potentially powerful but little explored tool for establishing thematically and spatially fine‐grained land‐use scenarios. Here, we use an ABM parameterized for 1,329 agents, mostly farmers, in a Central European model region, and simulate the changes to land‐use patterns resulting from their response to three scenarios of changing socio‐economic conditions and three scenarios of climate change until the mid of the century. Subsequently, we use species distribution models to, first, analyse relationships between the realized niches of 832 plant species and climatic gradients or land‐use types, respectively, and, second, to project consequent changes in potential regional ranges of these species as triggered by changes in both the altered land‐use patterns and the changing climate. We find that both drivers determine the realized niches of the studied plants, with land use having a stronger effect than any single climatic variable in the model. Nevertheless, the plants' future distributions appear much more responsive to climate than to land‐use changes because alternative future socio‐economic backgrounds have only modest impact on land‐use decisions in the model region. However, relative effects of climate and land‐use changes on biodiversity may differ drastically in other regions, especially where landscapes are still dominated by natural or semi‐natural habitat. We conclude that agent‐based modelling of land use is able to provide scenarios at scales relevant to individual species distribution and suggest that coupling ABMs with models of species' range change should be intensified to provide more realistic biodiversity forecasts.  相似文献   

9.
We assess the importance of anthropogenic land‐use, altered productivity, and species invasions for observed productivity–richness relationships in California. To this end, we model net primary productivity (NPP) c. 1750 AD and at present (1982–1999) and map native and exotic vascular plant richness for 230 subecoregions. NPP has increased up to 105% in semi‐arid areas and decreased up to 48% in coastal urbanized areas. Exotic invasions have increased local species diversity up to 15%. Human activities have reinforced historical gradients in species richness but reduced the spatial heterogeneity of NPP. Structural equation modelling suggests that, prior to European settlement, NPP and richness were primarily controlled by precipitation and other abiotic variables, with NPP mediating richness. Abiotic variables remain the strongest predictors of present NPP and richness, but intermodel comparisons indicate a significant anthropogenic impact upon statewide distributions of NPP and richness. Exotic and native species each positively correlate to NPP after controlling for other variables, which may help explain recent reports of positively associated native and exotic richness.  相似文献   

10.
Climate and evolutionary factors (e.g. diversification, time‐for‐speciation, niche conservatism) are both thought to be major drivers of species richness in regional assemblages. However, few studies have simultaneously investigated the relative effects of climate and evolutionary factors on species richness across a broad geographical extent. Here, we assess their relative effects on species richness of angiosperm trees across North America. Species richness of angiosperm trees in 1175 regional assemblages were related to climate and phylogenetic structure using a structural equation modeling (SEM) approach. Climate was quantified based on the mean temperature of the coldest month and mean annual precipitation. Evolutionary factors (time‐for‐speciation vs diversification) were inferred from phylogeny‐based measures of mean root distance, phylogenetic species variability, and net relatedness index. We found that at the continental scale, species richness is correlated with temperature and precipitation with approximately similar strength. In the SEM with net relatedness index and phylogenetic species variability and with all the 1175 quadrats, the total direct effect size of phylogenetic structure on species richness is greater than the total direct effect size of climate on species richness by a factor of 3.7. The specific patterns of phylogenetic structure (i.e. greater phylogenetic distances in more species rich regions) are consistent with the idea that time and niche conservatism drive richness patterns in North American angiosperm trees. We conclude that angiosperm tree species richness in regional assemblages in North America is more strongly related to patterns of phylogenetic relatedness than to climatic variation. The results of the present study support the idea that climatic and evolutionary explanations for richness patterns are not in conflict, and that evolutionary processes explain both the relationship between climate and richness and substantial variation in richness that is independent of climate.  相似文献   

11.
Elevational patterns of species richness, local abundance and assemblage structure of rainforest birds of north‐eastern Australia were explored using data from extensive standardized surveys throughout the region. Eighty‐two species of birds were recorded with strong turnover in assemblage structure across the elevational gradient and high levels of regional endemism in the uplands. Both species richness and bird abundance exhibited a humped‐shaped pattern with elevation with the highest values being between 600 and 800 m elevation. While much of the variability in species richness could be explained by the species–area relationship, analyses of net primary productivity (NPP) and total daily energy consumption of the bird community (energy use) suggest that ecosystem energy flow or constraints may be a significant determinant of species richness. Species richness is positively correlated with local bird abundance which itself is closely related to total energy use of the bird community. We suggest the hypothesis that lower NPP limits bird abundance and energy use in the uplands (>500 m) and that low bird energy use and species richness in the lowlands is limited by a seasonal bottleneck in available primary productivity and/or a species pool previously truncated by an extinction filter imposed by the almost complete disappearance of rainforest in the lowlands during the glacial maxima. We suggest that some of the previously predicted impacts of global warming on biodiversity in the uplands may be partially ameliorated by increases in NPP because of increasing temperatures. However, these relationships are complex and require further data specifically in regard to direct estimates of primary productivity and detailed estimates of energy flow within the assemblage.  相似文献   

12.
Environmental gradients (EG) related to climate, topography and vegetation are among the most important drivers of broad scale patterns of species richness. However, these different EG do not necessarily drive species richness in similar ways, potentially presenting synergistic associations when driving species richness. Understanding the synergism among EG allows us to address key questions arising from the effects of global climate and land use changes on biodiversity. Herein, we use variation partitioning (also know as commonality analysis) to disentangle unique and shared contributions of different EG in explaining species richness of Neotropical vertebrates. We use three broad sets of predictors to represent the environmental variability in (i) climate (annual mean temperature, temperature annual range, annual precipitation and precipitation range), (ii) topography (mean elevation, range and coefficient of variation of elevation), and (iii) vegetation (land cover diversity, standard deviation and range of forest canopy height). The shared contribution between two types of EG is used to quantify synergistic processes operating among EG, offering new perspectives on the causal relationships driving species richness. To account for spatially structured processes, we use Spatial EigenVector Mapping models. We perform analyses across groups with distinct dispersal abilities (amphibians, non-volant mammals, bats and birds) and discuss the influence of vagility on the partitioning results. Our findings indicate that broad scale patterns of vertebrate richness are mainly affected by the synergism between climate and vegetation, followed by the unique contribution of climate. Climatic factors were relatively more important in explaining species richness of good dispersers. Most of the variation in vegetation that explains vertebrate richness is climatically structured, supporting the productivity hypothesis. Further, the weak synergism between topography and vegetation urges caution when using topographic complexity as a surrogate of habitat (vegetation) heterogeneity.  相似文献   

13.
Bird species richness is mediated by local, regional, and historical factors, for example, competition, environmental heterogeneity, contemporary, and historical climate. Here, we related bird species richness with phylogenetic relatedness of bird assemblages, plant species richness, topography, contemporary climate, and glacial‐interglacial climate change to investigate the relative importance of these factors. This study was conducted in Inner Mongolia, an arid and semiarid region with diverse vegetation types and strong species richness gradients. The following associated variables were included as follows: phylogenetic relatedness of bird assemblages (Net Relatedness Index, NRI), plant species richness, altitudinal range, contemporary climate (mean annual temperature and precipitation, MAT and MAP), and contemporary‐Last Glacial Maximum (LGM) change in climate (change in MAT and change in MAP). Ordinary least squares linear, simultaneous autoregressive linear, and Random Forest models were used to assess the associations between these variables and bird species richness across this region. We found that bird species richness was correlated negatively with NRI and positively with plant species richness and altitudinal range, with no significant correlations with contemporary climate and glacial–interglacial climate change. The six best combinations of variables ranked by Random Forest models consistently included NRI, plant species richness, and contemporary‐LGM change in MAT. Our results suggest important roles of local ecological factors in shaping the distribution of bird species richness across this semiarid region. Our findings highlight the potential importance of these local ecological factors, for example, environmental heterogeneity, habitat filtering, and biotic interactions, in biodiversity maintenance.  相似文献   

14.
Macroclimatic niche properties derived from species distribution ranges are fundamental for projections of climate change impacts on biodiversity. However, it has been recognized that changes in regional or local distribution patterns also depend on interactions with land use. The reliability and transferability of large scale geographic predictions to small scale plant performance need to be tested experimentally. Thus, we asked how grassland plant species pairs with different macroclimatic niche properties respond to increased spring temperature and decrease summer precipitation in three different land‐use types. An experiment was carried out in the framework of the German Biodiversity Exploratories simulating climate change in 45 experimental plots in three geographical regions (Schorfheide‐Chorin, Hainich‐Dün, Schwäbische Alb) and three grassland management types (meadow, pasture, mown pasture). We planted six plant species as phytometers, each two of them representing congeneric species with contrasting macroclimatic niches and recorded plant survival and growth over 1 year. To quantify the species macroclimatic niches with respect to drought tolerance, the species’ distribution ranges were mapped and combined with global climate data. The simulated climate change had a general negative effect on plant survival and plant growth, irrespective of the macroclimatic niche characteristics of the species. Against expectation, species with ranges extending into drier regions did not generally perform better under drier conditions. Growth performance and survival was best in mown pastures, representing a quite intensive type of land use in all study regions. Species with higher macroclimatic drought tolerance were generally characterized by lower growth rates and higher survival rates in land‐use types with regular mowing regimes, probably because of reduced competition in the growing season. In conclusion, plant species with similar climatic niche characteristics cannot be expected to respond consistently over different regions owing to complex interactions of climate change with land use practices.  相似文献   

15.
We model future changes in land biogeochemistry and biogeography across East Africa. East Africa is one of few tropical regions where general circulation model (GCM) future climate projections exhibit a robust response of strong future warming and general annual‐mean rainfall increases. Eighteen future climate projections from nine GCMs participating in the Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment were used as input to the LPJ dynamic global vegetation model (DGVM), which predicted vegetation patterns and carbon storage in agreement with satellite observations and forest inventory data under the present‐day climate. All simulations showed future increases in tropical woody vegetation over the region at the expense of grasslands. Regional increases in net primary productivity (NPP) (18–36%) and total carbon storage (3–13%) by 2080–2099 compared with the present‐day were common to all simulations. Despite decreases in soil carbon after 2050, seven out of nine simulations continued to show an annual net land carbon sink in the final decades of the 21st century because vegetation biomass continued to increase. The seasonal cycles of rainfall and soil moisture show future increases in wet season rainfall across the GCMs with generally little change in dry season rainfall. Based on the simulated present‐day climate and its future trends, the GCMs can be grouped into four broad categories. Overall, our model results suggest that East Africa, a populous and economically poor region, is likely to experience some ecosystem service benefits through increased precipitation, river runoff and fresh water availability. Resulting enhancements in NPP may lead to improved crop yields in some areas. Our results stand in partial contradiction to other studies that suggest possible negative consequences for agriculture, biodiversity and other ecosystem services caused by temperature increases.  相似文献   

16.
Aim Land use and climate are two major components of global environmental change but our understanding of their simultaneous and interactive effects upon biodiversity is still limited. Here, we investigated the relationship between the species richness of neophytes, i.e. non‐native vascular plants introduced after 1500 AD, and environmental covariates to draw implications for future dynamics under land‐use and climate change. Location Switzerland, Central Europe. Methods The distribution of vascular plants was derived from a systematic national grid of 1 km2 quadrates (n = 456; Swiss Biodiversity Monitoring programme) including 1761 species, 122 of which were neophytes. Generalized linear models (GLMs) were used to correlate neophyte species richness with environmental covariates. The impact of land‐use and climate change was thereafter evaluated by projections for the years 2020 and 2050 using scenarios of moderate and strong changes for climate warming (IPCC) and urban sprawl (NRP 54). Results Mean annual temperature and the amount of urban areas explained neophyte species richness best, with a high predictive power of the corresponding model (cross‐validated D2 = 0.816). Climate warming had a stronger impact on the potential increase in the mean neophyte species richness (up to 191% increase by 2050) than ongoing urban sprawl (up to 10% increase) independently from variable interactions and model extrapolations to non‐analogue environments. Main conclusions In contrast to other vascular plants, the prediction of neophyte species richness at the landscape scale in Switzerland requires few variables only, and regions of highest species richness of the two groups do not coincide. The neophyte species richness is basically driven by climatic (temperature) conditions, and urban areas additionally modulate small‐scale differences upon this coarse‐scale pattern. According to the projections climate warming will contribute to the future increase in neophyte species richness much more than ongoing urbanization, but the gain in new neophyte species will be highest in urban regions.  相似文献   

17.
Aim Broad‐scale spatial patterns of species richness are very strongly correlated with climatic variables. If there is a causal link, i.e. if climate directly or indirectly determines patterns of richness, then when the climatic variables change, richness should change in the manner that spatial correlations between richness and climate would predict. The present study tests this prediction using seasonal changes in climatic variables and bird richness. Location We used a grid of equal area quadrats (37 000 km2) covering North and Central America as far south as Nicaragua. Methods Summer and winter bird distribution data were drawn from monographs and field guides. Climatic data came from published sources. We also used remotely sensed NDVI (normalized difference vegetation index — a measure of greenness). Results Bird species richness changes temporally (between summer and winter) in a manner that is close to, but statistically distinguishable from, the change one would predict from models relating the spatial variation in richness at a single time to climatic variables. If one further takes into account the seasonal changes in NDVI and within‐season variability of temperature and precipitation, then winter and summer richness follow congruent, statistically indistinguishable patterns. Main conclusions Our results are consistent with the hypothesis that climatic variables (temperature and precipitation) and vegetation cover directly or indirectly influence patterns of bird species richness.  相似文献   

18.
Species–area relationships (SARs) provide an avenue to model patterns of species richness and have recently been shown to vary substantially across regions of different climate, vegetation, and land cover. Given that a large proportion of the globe has been converted to agriculture, and considering the large variety in agricultural management practices, a key question is whether global SARs vary across gradients of agricultural intensity. We developed SARs for mammals that account for geographic variation in biomes, land cover and a range of land‐use intensity indicators representing inputs (e.g. fertilizer, irrigation), outputs (e.g. yields) and system‐level measures of intensity (e.g. human appropriation of net primary productivity – HANPP). We systematically compared the resulting SARs in terms of their predictive ability. Our global SAR with a universal slope was significantly improved by the inclusion of any one of the three variable types: biomes, land cover, and land‐use intensity. The latter, in the form of human appropriation of net primary productivity (HANPP), performed as well as biomes and land‐cover in predicting species richness. Other land‐use intensity indicators had a lower predictive ability. Our main finding that land‐use intensity performs as well as biomes and land cover in predicting species richness emphasizes that human factors are on a par with environmental factors in predicting global patterns of biodiversity. While our broad‐scale study cannot establish causality, human activity is known to drive species richness at a local scale, and our findings suggest that this may hold true at a global scale. The ability of land‐use intensity to explain variation in SARs at a global scale had not previously been assessed. Our study suggests that the inclusion of land‐use intensity in SAR models allows us to better predict and understand species richness patterns.  相似文献   

19.
Aim In this study, I determine the relationships between net primary productivity (NPP), human population density, species richness and land use. I also examine the implications of human settlement patterns for species conservation. Location Australia. Methods I document the associations between NPP, human population density and the species richness of birds, butterflies and mammals using correlations and spatial regressions. I also assess changes in land‐use with NPP and population density, focussing particularly on protected areas. An initial exploration into the implications of the NPP‐population density relationship for regional conservation strategies is provided. Results Human population density increases with NPP suggesting that available energy may be a key driving force of human settlement patterns. The species richness of each taxonomic group and geographically restricted species also increases with NPP leading to substantial overlap between species diversity and populated regions. The percentage of land designated as minimal use decreases considerably with increasing human population density and NPP, while intensive agriculture is confined entirely to areas of high NPP. There are strong negative relationships between the size of Australia's National Parks and human population density and NPP. Small parks are often surrounded by relatively dense settlements, but have high average NPP, while large parks are mostly isolated and characterized by low productivity. There are no areas in the highest quartile of NPP that also occur in the most sparsely populated regions, presenting challenges for conservation strategies wanting to protect productive areas under the least threat of human development. Main conclusions Human population density and species richness respond similarly to variation in NPP, leading to spatial congruence between human settlements and productive, species rich regions. Planning strategies are required that minimize the potential threat posed by human development to diverse ecosystems and maximize the underlying productivity of protected areas. Reducing the level of threat may require stabilizing the size of the human population, while capturing larger areas of relatively high productivity in the conservation reserve system would lead to greater protection of local diversity.  相似文献   

20.
Multiple environmental factors are known to shape species distributions at the global scale, including climate and topography, but understanding current extents of occurrence and biodiversity patterns requires considering anthropogenic factors as well. Numerous studies have explored the relationship between contemporary human activities and different biodiversity metrics, but the influence of past activities, such as land‐use, remains poorly understood despite being one of the oldest human impacts. Here we evaluate the role of past land‐use modifications in the current distribution and conservation status of mammals worldwide using spatial data characterizing human land use from ca BC 6000 to ca AD 2000. First, we applied a clustering method that revealed three generalized past human land‐use trajectories that represent low‐, recently‐ and steadily‐used areas widely represented across the globe. Second, we fitted boosted regression trees to predict total and threatened mammalian richness, globally and within trajectory‐clusters, testing the role of environmental factors and multiple human land‐use metrics reflecting: total used area at different time spans, rates of land‐use change, and the occurrence of remarkable land‐use shifts. Environmental factors were identified as the main correlates of current mammalian richness, but several proposed metrics of past land‐use were also relevant predictors. Overall, these results highlight the likely existence of a land‐use legacy in some regions of the world that has influenced the distribution of extant mammals, particularly of those currently classified as threatened. Even if we cannot change that legacy, our results show that we need to account for past human impacts to understand present biodiversity patterns and, arguably, to guide future actions.  相似文献   

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

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