首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 484 毫秒
1.
To predict changes in species' distributions due to climate change we must understand populations at the poleward edge of species' ranges. Ecologists generally expect range shifts under climate change caused by the expansion of edge populations as peripheral conditions increasingly resemble the range core. We tested whether peripheral populations of two contrasting butterflies, a small-bodied specialist (Erynnis propertius) and a large-bodied generalist (Papilio zelicaon), respond favorably to warmer conditions. Performance of populations related to climate was evaluated in seven peripheral populations spanning 1.2 degrees latitude (160 km) using: (1) population density surveys, an indirect measure of site suitability; and (2) organismal fitness in translocation experiments. There was evidence that population density increased with temperature for P. zelicaon whose population density declined with latitude in 1 of 3 sample years. On the other hand, E. propertius showed a positive relationship of population density with latitude, apparently unrelated to climate or measured habitat variables. Translocation experiments showed increased larval production at increased temperatures for both species, and in P. zelicaon, larval production also increased under drier conditions. These findings suggest that both species may increase at their range edge with warming but the preference for core-like conditions may be stronger in P. zelicaon. Further, populations of E. propertius at the range boundary may be large enough to act as sources of colonists for range expansions, but range expansion in this species may be prevented by a lack of available host plants further north. In total, the species appear to respond differently to climate and other factors that vary latitudinally, factors that will likely affect poleward expansion.  相似文献   

2.
Gene flow may influence the formation of species range limits, and yet little is known about the patterns of gene flow with respect to environmental gradients or proximity to range limits. With rapid environmental change, it is especially important to understand patterns of gene flow to inform conservation efforts. Here we investigate the species range of the selfing, annual plant, Mimulus laciniatus, in the California Sierra Nevada. We assessed genetic variation, gene flow, and population abundance across the entire elevation‐based climate range. Contrary to expectations, within‐population plant density increased towards both climate limits. Mean genetic diversity of edge populations was equivalent to central populations; however, all edge populations exhibited less genetic diversity than neighbouring interior populations. Genetic differentiation was fairly consistent and moderate among all populations, and no directional signals of contemporary gene flow were detected between central and peripheral elevations. Elevation‐driven gene flow (isolation by environment), but not isolation by distance, was found across the species range. These findings were the same towards high‐ and low‐elevation range limits and were inconsistent with two common centre‐edge hypotheses invoked for the formation of species range limits: (i) decreasing habitat quality and population size; (ii) swamping gene flow from large, central populations. This pattern demonstrates that climate, but not centre‐edge dynamics, is an important range‐wide factor structuring M. laciniatus populations. To our knowledge, this is the first empirical study to relate environmental patterns of gene flow to range limits hypotheses. Similar investigations across a wide variety of taxa and life histories are needed.  相似文献   

3.
A central question in ecology is how individual fitness interacts with the spatial variation in population density and habitat characteristics across species' ranges. We used fluctuating asymmetry (FA) as a measure of developmental stability (DS) in individuals of Tyrannus forficatus to estimate the suitability of sites of varying abundance and position within the species' range. FA in the inner-tail feathers of males and females is not spatially correlated across the species' range. FA in males increases towards the centre of the range and is not correlated with abundance. FA in females is not correlated with position in the range or abundance. Our results suggest that optimal sites are found throughout the range of the species, whereas suboptimal sites are mainly found towards the centre of the range. Additionally, our results suggest that abundance may not reflect the suitability of sites across species' ranges.  相似文献   

4.
Many predictions of how climate change will impact biodiversity have focused on range shifts using species‐wide climate tolerances, an approach that ignores the demographic mechanisms that enable species to attain broad geographic distributions. But these mechanisms matter, as responses to climate change could fundamentally differ depending on the contributions of life‐history plasticity vs. local adaptation to species‐wide climate tolerances. In particular, if local adaptation to climate is strong, populations across a species’ range—not only those at the trailing range edge—could decline sharply with global climate change. Indeed, faster rates of climate change in many high latitude regions could combine with local adaptation to generate sharper declines well away from trailing edges. Combining 15 years of demographic data from field populations across North America with growth chamber warming experiments, we show that growth and survival in a widespread tundra plant show compensatory responses to warming throughout the species’ latitudinal range, buffering overall performance across a range of temperatures. However, populations also differ in their temperature responses, consistent with adaptation to local climate, especially growing season temperature. In particular, warming begins to negatively impact plant growth at cooler temperatures for plants from colder, northern populations than for those from warmer, southern populations, both in the field and in growth chambers. Furthermore, the individuals and maternal families with the fastest growth also have the lowest water use efficiency at all temperatures, suggesting that a trade‐off between growth and water use efficiency could further constrain responses to forecasted warming and drying. Taken together, these results suggest that populations throughout species’ ranges could be at risk of decline with continued climate change, and that the focus on trailing edge populations risks overlooking the largest potential impacts of climate change on species’ abundance and distribution.  相似文献   

5.
Range shifts due to annual variation in temperature are more tractable than range shifts linked to decadal to century long temperature changes due to climate change, providing natural experiments to determine the mechanisms responsible for driving long‐term distributional shifts. In this study we couple physiologically grounded mechanistic models with biogeographic surveys in 2 years with high levels of annual temperature variation to disentangle the drivers of a historical range shift driven by climate change. The distribution of the barnacle Semibalanus balanoides has shifted 350 km poleward in the past half century along the east coast of the United States. Recruits were present throughout the historical range following the 2015 reproductive season, when temperatures were similar to those in the past century, and absent following the 2016 reproductive season when temperatures were warmer than they have been since 1870, the earliest date for temperature records. Our dispersal dependent mechanistic models of reproductive success were highly accurate and predicted patterns of reproduction success documented in field surveys throughout the historical range in 2015 and 2016. Our mechanistic models of reproductive success not only predicted recruitment dynamics near the range edge but also predicted interior range fragmentation in a number of years between 1870 and 2016. All recruits monitored within the historical range following the 2015 colonization died before 2016 suggesting juvenile survival was likely the primary driver of the historical range retraction. However, if 2016 is indicative of future temperatures mechanisms of range limitation will shift and reproductive failure will lead to further range retraction in the future. Mechanistic models are necessary for accurately predicting the effects of climate change on ranges of species.  相似文献   

6.

Aim

The aim was to assess the sensitivity of butterfly population dynamics to variation in weather conditions across their geographical ranges, relative to sensitivity to density dependence, and determine whether sensitivity is greater towards latitudinal range margins.

Location

Europe.

Time period

1980–2014.

Major taxa studied

Butterflies.

Methods

We use long‐term (35 years) butterfly monitoring data from > 900 sites, ranging from Finland to Spain, grouping sites into 2° latitudinal bands. For 12 univoltine butterfly species with sufficient data from at least four bands, we construct population growth rate models that include density dependence, temperature and precipitation during distinct life‐cycle periods, defined to accommodate regional variation in phenology. We use partial R2 values as indicators of butterfly population dynamics' sensitivity to weather and density dependence, and assess how these vary with latitudinal position within a species' distribution.

Results

Population growth rates appear uniformly sensitive to density dependence across species' geographical distributions, and sensitivity to density dependence is typically greater than sensitivity to weather. Sensitivity to weather is greatest towards range edges, with symmetry in northern and southern parts of the range. This pattern is not driven by variation in the magnitude of weather variability across the range, topographic heterogeneity, latitudinal range extent or phylogeny. Significant weather variables in population growth rate models appear evenly distributed across the life cycle and across temperature and precipitation, with substantial intraspecific variation across the geographical ranges in the associations between population dynamics and specific weather variables.

Main conclusions

Range‐edge populations appear more sensitive to changes in weather than those nearer the centre of species' distributions, but density dependence does not exhibit this pattern. Precipitation is as important as temperature in driving butterfly population dynamics. Intraspecific variation in the form and strength of sensitivity to weather suggests that there may be important geographical variation in populations' responses to climate change.  相似文献   

7.
A frequent assumption in ecology is that biotic interactions are more important than abiotic factors in determining lower elevational range limits (i.e., the “warm edge” of a species distribution). However, for species with narrow environmental tolerances, theory suggests the presence of a strong environmental gradient can lead to persistence, even in the presence of competition. The relative importance of biotic and abiotic factors is rarely considered together, although understanding when one exerts a dominant influence on controlling range limits may be crucial to predicting extinction risk under future climate conditions. We sampled multiple transects spanning the elevational range limit of Plethodon shenandoah and site and climate covariates were recorded. A two‐species conditional occupancy model, accommodating heterogeneity in detection probability, was used to relate variation in occupancy with environmental and habitat conditions. Regional climate data were combined with datalogger observations to estimate the cloud base heights and to project future climate change impacts on cloud elevations across the survey area. By simultaneously accounting for species’ interactions and habitat variables, we find that elevation, not competition, is strongly correlated with the lower elevation range boundary, which had been presumed to be restricted mainly as a result of competitive interactions with a congener. Because the lower elevational range limit is sensitive to climate variables, projected climate change across its high‐elevation habitats will directly affect the species’ distribution. Testing assumptions of factors that set species range limits should use models which accommodate detection biases.  相似文献   

8.
Species often harbour large amounts of phenotypic variation in ecologically important traits, and some of this variation is genetically based. Understanding how this genetic variation is spatially structured can help to understand species' ecological tolerances and range limits. We modelled the climate envelopes of Arabidopsis thaliana genotypes, ranging from early- to late-flowering, as a function of several climatic variables. We found that genotypes with contrasting alleles at individual flowering time loci differed significantly in potential range size and niche breadth. We also found that later flowering genotypes had more restricted range potentials and narrower niche breadths than earlier flowering genotypes, indicating that local selection on flowering can constrain or enhance the ability of populations to colonise other areas. Our study demonstrates how climate envelope models that incorporate ecologically important genetic variation can provide insights into the macroecology of a species, which is important to understand its responses to changing environments.  相似文献   

9.
The incidence of butterflies on British islands and their geographical (latitudinal) ranges are regressed on ecological and life history variables. The objective has been to investigate the contribution of individual variables and to incorporate information on phylogenetic links. The findings confirm the close relationship of species' incidence on islands with their geographical ranges on mainland Britain and that of species' geographical ranges with ecological variables, particularly migration capacity, hostplant type (variety) and breeding-habitat range. The results for island incidence considering phylogenetic links are virtually identical to those disregarding them. For geographical range, the results are similar. The key variable in each case is dispersal, scored in either one of two different ways. However, hostplant type takes precedence over breeding-habitat range when phylogenetic links are considered. Species categorized for upper and lower quartiles for geographical range form isolated clusters in the first two axes of a principal components analysis on a set of seven ecological variables. This result suggests the combined influence of a number of ecological variables on range size. Species with wide geographical ranges tend to have high migration indices, a wide variety of hostplants and ubiquitous hostplants, long flight periods and are often multi-brooded; those with narrow geographical ranges tend to have low migration indices, sparse and limited hostplant resources, short flight periods and are often univoltine. A number of life history variables are found to correlate significantly with geographical range, but account for only small amounts of variation. The lack of any association between range size and population abundance may well reflect the difficulty of obtaining adequate measures for abundance. However, we caution against expecting a strong correlation between range size and abundance.  相似文献   

10.
Climate change is causing widespread geographical range shifts, which likely reflects different processes at leading and trailing range margins. Progressive warming is thought to relax thermal barriers at poleward range margins, enabling colonization of novel areas, but imposes increasingly unsuitable thermal conditions at equatorward margins, leading to range losses from those areas. Few tests of this process during recent climate change have been possible, but understanding determinants of species’ range limits will improve predictions of their geographical responses to climate change and variation in extinction risk. Here, we examine the relationship between poleward and equatorward range margin dynamics with respect to temperature‐related geographical limits observed for 34 breeding passerine species in North America between 1984–1988 and 2002–2006. We find that species’ equatorward range margins were closer to their upper realized thermal niche limits and proximity to those limits predicts equatorward population extinction risk through time. Conversely, the difference between breeding bird species’ poleward range margin temperatures and the coolest temperatures they tolerate elsewhere in their ranges was substantial and remained consistent through time: range expansion at species’ poleward range margins is unlikely to directly reflect lowered thermal barriers to colonization. The process of range expansion may reflect more complex factors operating across broader areas of species’ ranges. The latitudinal extent of breeding bird ranges is decreasing through time. Disparate responses observed at poleward versus equatorward margins arise due to differences in range margin placement within the realized thermal niche and suggest that climate‐induced geographical shift at equatorward range limits more strongly reflect abiotic conditions than at their poleward range limits. This further suggests that observed geographic responses to date may fail to demonstrate the true cost of climate change on the poleward portion of species’ distributions. Poleward range margins for North American breeding passerines are not presently in equilibrium with realized thermal limits.  相似文献   

11.
Poleward range expansions are widespread responses to recent climate change and are crucial for the future persistence of many species. However, evolutionary change in traits such as colonization history and habitat preference may also be necessary to track environmental change across a fragmented landscape. Understanding the likelihood and speed of such adaptive change is important in determining the rate of species extinction with ongoing climate change. We conducted an amplified fragment length polymorphism (AFLP)‐based genome scan across the recently expanded UK range of the Brown Argus butterfly, Aricia agestis, and used outlier‐based (DFDIST and BayeScan) and association‐based (Isolation‐By‐Adaptation) statistical approaches to identify signatures of evolutionary change associated with range expansion and habitat use. We present evidence for (i) limited effects of range expansion on population genetic structure and (ii) strong signatures of selection at approximately 5% AFLP loci associated with both the poleward range expansion of A. agestis and differences in habitat use across long‐established and recently colonized sites. Patterns of allele frequency variation at these candidate loci suggest that adaptation to new habitats at the range margin has involved selection on genetic variation in habitat use found across the long‐established part of the range. Our results suggest that evolutionary change is likely to affect species’ responses to climate change and that genetic variation in ecological traits across species’ distributions should be maximized to facilitate range shifts across a fragmented landscape, particularly in species that show strong associations with particular habitats.  相似文献   

12.
Predicted future climate change will alter species' distributions as they attempt to track the most suitable 'climate window'. Climate envelope models indicate the direction of likely range changes but do not incorporate population dynamics, therefore observed responses may differ greatly from these projections. We use simulation modelling to explore the consequences of a period of environmental change for a species structured across an environmental gradient. Results indicate that a species' range may lag behind its climate envelope and demonstrate that the rate of movement of a range can accelerate during a period of climate change. We conclude that the inclusion of both population dynamics and spatial environmental variability is vital to develop models that can both predict, and be used to manage, the impact of changing climate on species' biogeography.  相似文献   

13.
Aim It is often assumed that species reach their highest densities in the centre of their ranges and decline in abundance toward the edges of the range. Implicit in this notion, which we call the abundant centre hypothesis, is the assumption that the edges of the range are more stressful to organisms and are more likely to show responses to climate change. However, an earlier review and empirical study of patterns of abundance across the range of intertidal invertebrates show little support for the abundant centre hypothesis and further demonstrated that few studies have examined patterns in either abundance or stress across species ranges. In part this gap is due to the logistical difficulties of sampling species across large geographical ranges. Here we use intertidal invertebrates, which have relatively simple linear latitudinal ranges, and heat‐shock proteins, which have been shown to be an integrative measure of organismal stress, to test the hypothesis that species are more stressed at the edges of their range. We use complementary data on population density to test the relationship between stress proteins and overall species density across the species’ range. Location Our sampling programme covered the southern half of the large geographical ranges of two intertidal invertebrates on the Pacific Coast of North America. Sites were spread between northern Baja California, Mexico and Vancouver Island, Canada, a range of c. 22 degrees of latitude. Method We sampled levels of heat‐shock protein 70 (Hsp70) in eight to 12 individuals from each of 20 sites for the intertidal mussel Mytilus californianus and 11 sites for the intertidal snail Nucella ostrina, spread throughout the southern half of their geographical ranges. The relationships between levels of Hsp70 in individuals from a site and (1) latitude of the site, (2) the site's position in the species’ range and (3) average population density were determined. Results No significant relationship was found in either species between levels of Hsp70 and latitude, position in the range or population density. Complex patterns that did emerge may be explained by nonlinear gradients in environmental conditions along the Pacific coast. Specifically, we observed peak values of Hsp70 for both species in northern Oregon, where intertidal zones are disproportionately exposed to daytime emersion (exposure to air) in the summer months of collections. A second peak for M. californianus was found south of Point Conception, California, which is marked by dramatic shifts toward warmer sea temperatures and decreased wave exposure. Main conclusions Patterns that emerged were not predicted by simple models based on the abundant centre hypothesis. However, they are consistent with more complex pictures of heat stress, organismal condition and abundance along a latitudinal gradient that have been demonstrated in recent studies. We suggest that latitudinal complexity, species‐specific differences and local effects must be considered before generalizing the relationship between environmental stress, abundance, range limits and responses of ranges to climate change.  相似文献   

14.
Aim Bees are the most important pollinators of flowering plants and essential ecological keystone species contributing to the integrity of most terrestrial ecosystems. Here, we examine the potential impact of climate change on bees’ geographic range in a global biodiversity hotspot. Location South Africa with a focus on the Cape Floristic Region (CFR) diversity hotspot. Methods  Geographic ranges of 12 South African bee species representing dominant distribution types were studied, and the climate change impacts upon bees were examined with A2 and B2 climate scenarios of HadCM3 model, using MaxEnt for species distribution modelling. Results The predicted levels of climate change‐induced impacts on species ranges varied from little shifts and range expansion of 5–50% for two species to substantial range contractions between 32% and 99% in another six species. Four species show considerable range shifts. Bees of the winter rainfall area in the west of South Africa generally have smaller range sizes than in the summer rainfall area and generally show eastward range contractions toward the dry interior. Bee species prevalent in summer rainfall regions show a tendency for a south‐easterly shift in geographic range. Main conclusions The bee fauna of the CFR is identified as the most vulnerable to climate change due to the high level of endemism, the small range sizes and the island‐like isolation of the Mediterranean‐type climate region at the SW tip of Africa. For monitoring climate change impact on bees, we suggest to establish observatories in the coastal plains of the west coast that are predicted to be worst affected and areas where persistence of populations is most likely. Likely impacts of climate change on life history traits of bees (phenology, sociality, bee‐host plant synchronization) are discussed but require further investigation.  相似文献   

15.
Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well‐suited for studying local adaptation. Here, we analyzed genome‐wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.  相似文献   

16.
Ongoing changes in global climate are altering ecological conditions for many species. The consequences of such changes are typically most evident at the edge of a species’ geographical distribution, where differences in growth or population dynamics may result in range expansions or contractions. Understanding population responses to different climatic drivers along wide latitudinal and altitudinal gradients is necessary in order to gain a better understanding of plant responses to ongoing increases in global temperature and drought severity. We selected Scots pine (Pinus sylvestris L.) as a model species to explore growth responses to climatic variability (seasonal temperature and precipitation) over the last century through dendrochronological methods. We developed linear models based on age, climate and previous growth to forecast growth trends up to year 2100 using climatic predictions. Populations were located at the treeline across a latitudinal gradient covering the northern, central and southernmost populations and across an altitudinal gradient at the southern edge of the distribution (treeline, medium and lower elevations). Radial growth was maximal at medium altitude and treeline of the southernmost populations. Temperature was the main factor controlling growth variability along the gradients, although the timing and strength of climatic variables affecting growth shifted with latitude and altitude. Predictive models forecast a general increase in Scots pine growth at treeline across the latitudinal distribution, with southern populations increasing growth up to year 2050, when it stabilizes. The highest responsiveness appeared at central latitude, and moderate growth increase is projected at the northern limit. Contrastingly, the model forecasted growth declines at lowland‐southern populations, suggesting an upslope range displacement over the coming decades. Our results give insight into the geographical responses of tree species to climate change and demonstrate the importance of incorporating biogeographical variability into predictive models for an accurate prediction of species dynamics as climate changes.  相似文献   

17.
Global warming is now recognized as the dominant threat to biodiversity because even protected populations and habitats are susceptible. Nonetheless, current criteria for evaluating species' relative endangerment remain purely ecological, and the accepted conservation strategies of habitat preservation and population management assume that species can mount ecological responses if afforded protection. The insidious threat from climate change is that it will attenuate or preclude ecological responses by species that are physiologically constrained; yet, quantitative, objective criteria for assessing relative susceptibility of diverse taxa to warming-induced stress are wanting. We explored the utility of using interspecies physiological variation for this purpose by relating species' physiological phenotypes to landscape patterns of ecological and genetic exchange. Using a salamander model system in which ecological, genetic and physiological diversity are well characterized, we found strong quantitative relationships of basal metabolic rates (BMRs) to both macroecological and phylogeographic patterns, with decreasing BMR leading to dispersal limitation (small contemporary ranges with marked phylogeographic structure). Measures of intrinsic physiological tolerance, which vary systematically with macroecological and phylogeographic patterns, afford objective criteria for assessing endangerment across a wide range of species and should be incorporated into conservation assessment criteria that currently rely exclusively upon ecological predictors.  相似文献   

18.
Understanding the boundaries of species' ranges and the variations in population dynamics from the centre to margin of a species' range is critical. This study simulated spatial-temporal patterns of birth and death rates and migration across a species' range in different seasons. Our results demonstrated the importance of dispersal and migration in altering birth and death rates, balancing source and sink habitats, and governing expansion or contraction of species' ranges in changing environments. We also showed that the multiple equilibria of metapopulations across a species' range could be easily broken following climatic changes or physical disturbances either local or regional. Although we refer to our models as describing the population dynamics across whole species' range, they should also apply to small-scale habitats (metapopulations) in which species abundance follows a humped pattern or to any ecosystem or landscape where strong central-marginal (C-M) environmental gradients exist. Conservation of both central and marginal populations would therefore be equally important considerations in making management decisions.  相似文献   

19.
Although some organisms have moved to higher elevations and latitudes in response to recent climate change, there is little consensus regarding the capacity of different species to track rapid climate change via range shifts. Understanding species' abilities to shift ranges has important implications for assessing extinction risk and predicting future community structure. At an expanding front, colonization rates are determined jointly by rates of reproduction and dispersal. In addition, establishment of viable populations requires that individuals find suitable resources in novel habitats. Thus, species with greater dispersal ability, reproductive rate and ecological generalization should be more likely to expand into new regions under climate change. Here, we assess current evidence for the relationship between leading-edge range shifts and species' traits. We found expected relationships for several datasets, including diet breadth in North American Passeriformes and egg-laying habitat in British Odonata. However, models generally had low explanatory power. Thus, even statistically and biologically meaningful relationships are unlikely to be of predictive utility for conservation and management. Trait-based range shift forecasts face several challenges, including quantifying relevant natural history variation across large numbers of species and coupling these data with extrinsic factors such as habitat fragmentation and availability.  相似文献   

20.
Understanding and predicting how species' distributions will shift as climate changes are central questions in ecology today. The late Quaternary of North America represents a natural experiment in which we can evaluate how species responded during the expansion and contraction of the glaciers. Here, we ask whether species' range shifts differ because of taxonomic affinity, life-history traits, body size or topographic heterogeneity and whether the species survived the megafaunal extinction. There was no difference in range shifts between victims and survivors of the megafaunal extinction. In general, the change in the size of a species' range is not well correlated with any of the ecological or life-history traits evaluated. However, there are significant relationships between some variables and the movements of the centroids of ranges. Differences in the distances shifted exist among orders, although this is probably a result of body size differences as larger bodied species show larger shifts. Although there are a few exceptions, the distance that species shifted their range was weakly correlated with life-history traits. Finally, species in more topographically heterogeneous areas show smaller shifts than species in less-diverse areas. Overall, these results indicate that when trying to predict species range shifts in the future, body size, lifespan and the topographic relief of the landscape should be taken into account.  相似文献   

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

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