首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 343 毫秒
1.
Terrestrial ecosystems are an important part of Earth systems, and they are undergoing remarkable changes in response to global warming. This study investigates the response of the terrestrial vegetation distribution and carbon fluxes to global warming by using the new dynamic global vegetation model in the second version of the Chinese Academy of Sciences (CAS) Earth System Model (CAS-ESM2). We conducted two sets of simulations, a present-day simulation and a future simulation, which were forced by the present-day climate during 1981–2000 and the future climate during 2081–2100, respectively, as derived from RCP8.5 outputs in CMIP5. CO2 concentration is kept constant in all simulations to isolate CO2-fertilization effects. The results show an overall increase in vegetation coverage in response to global warming, which is the net result of the greening in the mid-high latitudes and the browning in the tropics. The results also show an enhancement in carbon fluxes in response to global warming, including gross primary productivity, net primary productivity, and autotrophic respiration. We found that the changes in vegetation coverage were significantly correlated with changes in surface air temperature, reflecting the dominant role of temperature, while the changes in carbon fluxes were caused by the combined effects of leaf area index, temperature, and precipitation. This study applies the CAS-ESM2 to investigate the response of terrestrial ecosystems to climate warming. Even though the interpretation of the results is limited by isolating CO2-fertilization effects, this application is still beneficial for adding to our understanding of vegetation processes and to further improve upon model parameterizations.  相似文献   

2.
Akihiko Ito 《Climatic change》2017,142(1-2):53-66
Geoengineering such as solar radiation management (SRM) can be an emergent option to avoid devastating climatic warming, but its ramifications are barely understood. The perturbation of the Earth’s energy balance, atmospheric dynamics, and hydrological cycling may exert unexpected influences on natural and human systems. In this study, I evaluate the impacts of SRM deployment on terrestrial ecosystem functions using a process-based ecosystem model (the Vegetation Integrative Simulator for Trace gases, VISIT) driven by the climate projections by multiple climate models. In the SRM-oriented climate projections, massive injection of sulphate aerosols into the stratosphere lead to increased scattering of solar radiation and delayed anthropogenic climate warming. The VISIT simulations show that canopy light absorption and gross primary production are enhanced in subtropics in spite of the slight decrease of total incident solar radiation. The retarded temperature rise during the deployment period leads to lower respiration, and consequently, an additional net terrestrial ecosystem carbon uptake by about 20%. After the SRM termination, however, along with the temperature rise, this carbon is released rapidly to the atmosphere. As a result of altered precipitation and radiation budget, simulated runoff discharge is suppressed mainly in the tropics. These SRM-induced influences on terrestrial ecosystems occurr heterogeneously over the land surface and differed among the ecosystem functions. These responses of terrestrial functions should be taken into account when discussing the costs and benefits of geoengineering.  相似文献   

3.
Terrestrial carbon fluxes are an important factor in regulating concentrations of atmospheric carbon dioxide (CO2). In this study, we use a coupled climate model with interactive biogeochemistry to benchmark the simulation of net primary productivity (NPP) and its response to elevated atmospheric CO2. Short-term field experiments such as Free-Air Carbon Dioxide Enrichment (FACE) studies have examined this phenomenon but it is difficult to infer trends from only a few years of field data. Here, we employ the University of Victoria's Earth System Climate Model (UVic ESCM) version 2.8 to compare simulated changes in NPP due to an elevated atmospheric CO2 concentration of 550 ppm to observed increases in NPP of 23% ±2% from four temperate forest FACE studies between 1997 and 2002. We further compare two scenarios: elevated CO2 with climate change, and elevated CO2 without climate change, the latter being consistent with FACE methodology. In the climate change scenario global terrestrial and forest-only NPP increased by 24.5% and 27.9%, respectively, while these increases were 21.0% and 17.2%, respectively, in the latitude band most representative of the location of the FACE studies. In the scenario without climate change, terrestrial and forest-only NPP increased instead by 28.3% and 30.6%, respectively, while these increases were 24.3% and 14.4%, respectively, in the FACE latitudes. This suggests that the model may underestimate temperate forest NPP increases when compared to results from temperate forest FACE studies and highlights the need for both increased experimental study of other forest biomes and further model development.  相似文献   

4.
The A.M.Obukhov Institute of Atmospheric Physics,Russian Academy of Sciences (IAP RAS) climate model (CM) of intermediate complexity is extended by a spatially explicit terrestrial carbon cycle module.Numerical experiments with the IAP RAS CM are performed forced by the reconstructions of anthropogenic and natural forcings for the 16th to the 20th centuries and by combined SRES (Special Report on Emission Scenarios) A2-LUH (Land Use Harmonization) anthropogenic scenarios for the 21st century.Hereby,the impact of uncertainty in land-use scenarios on results of simulations with a coupled climate-carbon cycle model is tested.The simulations of the model realistically reproduced historical changes in carbon cycle characteristics.In the IAP RAS CM,climate warming reproduced in the 20th and 21st centuries enhanced terrestrial net primary production but terrestrial carbon uptake was suppressed due to an overcompensating increase in soil respiration.Around year 2100,the simulations the model forced by different land use scenarios diverged markedly,by about 70 Pg (C) in terms of biomass and soil carbon stock but they differed only by about 10 ppmv in terms of atmospheric carbon dioxide content.  相似文献   

5.
Increasing concentrations of atmospheric CO2 influence climate, terrestrial biosphere productivity and ecosystem carbon storage through its radiative, physiological and fertilization effects. In this paper, we quantify these effects for a doubling of CO2 using a low resolution configuration of the coupled model NCAR CCSM4. In contrast to previous coupled climate-carbon modeling studies, we focus on the near-equilibrium response of the terrestrial carbon cycle. For a doubling of CO2, the radiative effect on the physical climate system causes global mean surface air temperature to increase by 2.14 K, whereas the physiological and fertilization on the land biosphere effects cause a warming of 0.22 K, suggesting that these later effects increase global warming by about 10 % as found in many recent studies. The CO2-fertilization leads to total ecosystem carbon gain of 371 Gt-C (28 %) while the radiative effect causes a loss of 131 Gt-C (~10 %) indicating that climate warming damps the fertilization-induced carbon uptake over land. Our model-based estimate for the maximum potential terrestrial carbon uptake resulting from a doubling of atmospheric CO2 concentration (285–570 ppm) is only 242 Gt-C. This highlights the limited storage capacity of the terrestrial carbon reservoir. We also find that the terrestrial carbon storage sensitivity to changes in CO2 and temperature have been estimated to be lower in previous transient simulations because of lags in the climate-carbon system. Our model simulations indicate that the time scale of terrestrial carbon cycle response is greater than 500 years for CO2-fertilization and about 200 years for temperature perturbations. We also find that dynamic changes in vegetation amplify the terrestrial carbon storage sensitivity relative to a static vegetation case: because of changes in tree cover, changes in total ecosystem carbon for CO2-direct and climate effects are amplified by 88 and 72 %, respectively, in simulations with dynamic vegetation when compared to static vegetation simulations.  相似文献   

6.
Probabilistic climate change projections using neural networks   总被引:5,自引:0,他引:5  
Anticipated future warming of the climate system increases the need for accurate climate projections. A central problem are the large uncertainties associated with these model projections, and that uncertainty estimates are often based on expert judgment rather than objective quantitative methods. Further, important climate model parameters are still given as poorly constrained ranges that are partly inconsistent with the observed warming during the industrial period. Here we present a neural network based climate model substitute that increases the efficiency of large climate model ensembles by at least an order of magnitude. Using the observed surface warming over the industrial period and estimates of global ocean heat uptake as constraints for the ensemble, this method estimates ranges for climate sensitivity and radiative forcing that are consistent with observations. In particular, negative values for the uncertain indirect aerosol forcing exceeding –1.2 Wm–2 can be excluded with high confidence. A parameterization to account for the uncertainty in the future carbon cycle is introduced, derived separately from a carbon cycle model. This allows us to quantify the effect of the feedback between oceanic and terrestrial carbon uptake and global warming on global temperature projections. Finally, probability density functions for the surface warming until year 2100 for two illustrative emission scenarios are calculated, taking into account uncertainties in the carbon cycle, radiative forcing, climate sensitivity, model parameters and the observed temperature records. We find that warming exceeds the surface warming range projected by IPCC for almost half of the ensemble members. Projection uncertainties are only consistent with IPCC if a model-derived upper limit of about 5 K is assumed for climate sensitivity.  相似文献   

7.
Global warming simulations are performed with a coupled climate model of reduced complexity to investigate global warming–marine carbon cycle feedbacks. The model is forced by emissions of CO2 and other greenhouse agents from scenarios recently developed by the Intergovernmental Panel on Climate Change and by CO2 stabilization profiles. The uptake of atmospheric CO2 by the ocean is reduced between 7 to 10% by year 2100 compared to simulations without global warming. The reduction is of similar size in the Southern Ocean and in low‐latitude regions (32.5°S‐32.5°N) until 2100, whereas low‐latitude regions dominate on longer time scales. In the North Atlantic the CO2 uptake is enhanced, unless the Atlantic thermohaline circulation completely collapses. At high latitudes, biologically mediated changes enhance ocean CO2 uptake, whereas in low‐latitude regions the situation is reversed. Different implementations of the marine biosphere yield a range of 5 to 16% for the total reduction in oceanic CO2 uptake until year 2100. Modeled oceanic O2 inventories are significantly reduced in global warming simulations. This suggests that the terrestrial carbon sink deduced from atmospheric O2/N2 observations is potentially overestimated if the oceanic loss of O2 to the atmosphere is not considered.  相似文献   

8.
As "the third pole", the Tibetan Plateau (TP) is sensitive to climate forcing and has experienced rapid warming in recent decades. This study analyzes annual and seasonal near-surface air temperature changes on the TP in response to transient and stabilized 2.0°C/1.5°C global warming targets based on simulations of the Community Earth System Model (CESM). Elevation-dependent warming (EDW) with faster warming at higher elevations is predicted. A surface energy budget analysis is adopted to uncover the mechanisms responsible for the temperature changes. Our results indicate a clear amplified warming on the TP with positive EDW in 2.0°C/1.5°C warmer futures, especially in the cold season. Mean TP warming relative to the reference period (1961–90) is dominated by an enhanced downward longwave radiation flux, while the variations in surface albedo shape the detailed pattern of EDW. For the same global warming level, the temperature changes under transient scenarios are ~0.2°C higher than those under stabilized scenarios, and the characteristics of EDW are broadly similar for both scenarios. These differences can be primarily attributed to the combined effects of differential downward longwave radiation, cloud radiative forcing, and surface sensible and latent heat fluxes. These findings contribute to a more detailed understanding of regional climate on the TP in response to the long-term climate goals of the Paris Agreement and highlight the differences between transient and stabilized warming scenarios.  相似文献   

9.
This study aims to demonstrate the potential of a process-based regional ecosystem model, LPJ-GUESS, driven by climate scenarios generated by a regional climate model system (RCM) to generate predictions useful for assessing effects of climatic and CO2 change on the key ecosystem services of carbon uptake and storage. Scenarios compatible with the A2 and B2 greenhouse gas emission scenarios of the Special Report on Emission Scenarios (SRES) and with boundary conditions from two general circulation models (GCMs) – HadAM3H and ECHAM4/OPYC3 – were used in simulations to explore changes in tree species distributions, vegetation structure, productivity and ecosystem carbon stocks for the late 21st Century, thus accommodating a proportion of the GCM-based and emissions-based uncertainty in future climate development. The simulations represented in this study were of the potential natural vegetation ignoring direct anthropogenic effects. Results suggest that shifts in climatic zones may lead to changes in species distribution and community composition among seven major tree species of natural Swedish forests. All four climate scenarios were associated with an extension of the boreal forest treeline with respect to altitude and latitude. In the boreal and boreo-nemoral zones, the dominance of Norway spruce and to a lesser extent Scots pine was reduced in favour of deciduous broadleaved tree species. The model also predicted substantial increases in vegetation net primary productivity (NPP), especially in central Sweden. Expansion of forest cover and increased local biomass enhanced the net carbon sink over central and northern Sweden, despite increased carbon release through decomposition processes in the soil. In southern Sweden, reduced growing season soil moisture levels counterbalanced the positive effects of a longer growing season and increased carbon supply on NPP, with the result that many areas were converted from a sink to a source of carbon by the late 21st century. The economy-oriented A2 emission scenario would lead to higher NPP and stronger carbon sinks according to the simulations than the environment-oriented B2 scenario.  相似文献   

10.
We analyzed interannual variability (IAV) of precipitation and air temperature over a 40-year period (1969–2008) for 11 sites along a precipitation gradient on the Tibetan Plateau. The observed IAV for both precipitation and air temperature decreases with increasing mean annual precipitation. Using Biome-BGC, a process-based ecosystem model, we simulated net primary production (NPP) along this gradient and find that the IAV of NPP is positively correlated to the IAV of precipitation and temperature. Following projected climate change scenarios for the Tibetan Plateau, our simulations suggest that with increasing IAV of precipitation and temperature, the IAV of NPP will also increase and that climate thresholds exist that, if surpassed, lead to ecosystem die-off. The impacts of these changes on ecosystem processes and climate-vegetation feedbacks on the rapidly warming Tibetan Plateau are potentially quite significant.  相似文献   

11.
Modeling potential global redistribution of terrestrial vegetation frequently is based on bioclimatic classifications which relate static regional vegetation zones (biomes) to a set of static climate parameters. The equilibrium character of the relationships limits our confidence in their application to scenarios of rapidly changing climate. Such assessments could be improved if vegetation migration and succession would be incorporated as response variables in model simulations. We developed the model MOVE (Migration Of VEgetation), to simulate the geographical implications of different rates of plant extirpation and in-migration. We used the model to study the potential impact on terrestrial carbon stocks of climate shifts hypothesized from a doubling of atmospheric greenhouse gas concentration. The model indicates that the terrestrial vegetation and soil could release carbon; the amount of this carbon pulse depends on the rate of migration relative to the rate of climate change. New temperate and boreal biomes, not found on the landscape today, increase rapidly in area during the first 100 years of simulated response to climate change. Their presence for several centuries and their gradual disappearance after the climate ceases to change adds uncertainty in calculating future terrestrial carbon fluxes.  相似文献   

12.
13.
The projected temperature and precipitationchange under different emissions scenarios using Coupled Model Intercomparison Project Phase 5 models over the northwestern arid regions of China(NWAC) were analyzed using the ensemble of three high-resolution dynamical downscaling simulations: the simulation of the Regional Climate Model version 4.0(Reg CM4) forced by the Beijing Climate Center Climate System Model version 1.1(BCC_CSM1.1); the Hadley Centre Global Environmental Model version 3 regional climate model(Had GEM3-RA) forced by the Atmosphere-Ocean coupled Had GEM version 2(Had GEM2-AO); and the Weather Research and Forecasting(WRF) model forced by the Norwegian community Earth System Model(Nor ESM1-M). Model validation indicated that the multimodel simulations reproduce the spatial and temporal distribution of temperature and precipitation well. The temperature is projected to increase over NWAC under both the 4.5 and 8.5 Representative Concentration Pathways scenarios(RCP4.5 and RCP8.5, respectively) in the middle of the 21 st century, but the warming trend is larger under the RCP8.5 scenario. Precipitation shows a significant increasing trend in spring and winter under both RCP4.5 and RCP8.5; but in summer, precipitation is projected to decrease in the Tarim Basin and Junggar Basin. The regional averaged temperature and precipitation show increasing trends in the future over NWAC; meanwhile, the large variability of the winter mean temperature and precipitation may induce more extreme cold events and intense snowfall events in these regions in the future.  相似文献   

14.
针对《巴黎协定》提出的温控目标,利用耦合模式比较计划第五阶段(CMIP5)模式在RCP2.6、RCP4.5和RCP8.5情景下的模拟结果,初步分析了全球升温情景下陆地生态系统净初级生产力(NPP)相对于参考时段(1986—2005年)的变化,重点分析了1.5℃和2℃升温时NPP相对于参考时段的变化量,并探讨了大气CO2浓度、气温、降水和辐射的变化及其对NPP变化的影响。CMIP5基于各典型浓度路径模拟的全球陆地生态系统NPP均呈增加趋势,且NPP增加量与升温幅度成正比。在相同的升温幅度下,基于各典型浓度路径模拟的各环境因子和NPP的变化量较为一致。陆地生态系统NPP总量增加主要由大气CO2浓度上升驱动,其他环境因子的影响相对较弱。中国东南部、非洲中部、美国东南部和亚马孙雨林西部地区NPP增加最明显。NPP变化量的空间格局主要由大气CO2浓度增加和升温控制,降水和辐射的影响相对较小。具体而言,大气CO2浓度上升对中低纬度的NPP变化贡献最大,对北方高纬度地区NPP变化贡献较小。温度上升有利于促进北方高纬度地区和青藏高原地区NPP,但对中低纬度地区的NPP有较强的抑制作用。鉴于既有典型浓度路径和地球系统模型的限制,本文对未来升温情景下陆地生态系统NPP的预估仍存在较大的不确定性,需要在未来的研究中进一步改进。  相似文献   

15.
Afforestation has been deployed as a mitigation strategy for global warming due to its substantial carbon sequestration, which is partly counterbalanced with its biogeophysical effects through modifying the fluxes of energy, water, and momentum at the land surface. To assess the potential biophysical effects of afforestation, a set of extreme experiments in an Earth system model of intermediate complexity, the McGill Paleoclimate Model-2 (MPM-2), is designed. Model results show that latitudinal afforestation not only has a local warming effect but also induces global and remote warming over regions beyond the forcing originating areas. Precipitation increases in the northern hemisphere and decreases in southern hemisphere in response to afforestation. The local surface warming over the forcing originating areas in northern hemisphere is driven by decreases in surface albedo and increases in precipitation. The remote surface warming in southern hemisphere is induced by decreases in surface albedo and precipitation. The results suggest that the potential impact of afforestation on regional and global climate depended critically on the location of the forest expansion. That is, afforestation in 0°–15°N leaves a relatively minor impact on global and regional temperature; afforestation in 45°–60°N results in a significant global warming, while afforestation in 30°–45°N results in a prominent regional warming. In addition, the afforestation leads to a decrease in annual mean meridional oceanic heat transport with a maximum decrease in forest expansion of 30°–45°N. These results can help to compare afforestation effects and find areas where afforestation mitigates climate change most effectively combined with its carbon drawdown effects.  相似文献   

16.
Summary Using a high resolution regional climate model we perform multiple January simulations of the impact of land cover change over western Australia. We focus on the potential of reforestation to ameliorate the projected warming over western Australia under two emission scenarios (A2, B2) for 2050 and 2100. Our simulations include the structural and physiological responses of the biosphere to changes in climate and changes in carbon dioxide. We find that reforestation has the potential to reduce the warming caused by the enhanced greenhouse effect by as much as 30% under the A2 and B2 scenarios by 2050 but the cooling effect declines to 10% by 2100 as CO2-induced warming intensifies. The cooling effect of reforestation over western Australia is caused primarily by the increase in leaf area index that leads to a corresponding increase in the latent heat flux. This cooling effect is localized and there were no simulated changes in temperature over regions remote from land cover change. We also show that the more extreme emission scenario (A2) appears to lead to a more intense response in photosynthesis by 2100. Overall, our results are not encouraging in terms of the potential to offset future warming by large scale reforestation. However, at regional scales the impact of land cover change is reasonably large relative to the impact of increasing carbon dioxide (up to 2050) suggesting that future projections of the Australian climate would benefit from the inclusion of projections of future land cover change. We suggest that this would add realism and regional detail to future projections and perhaps aid detection and attribution studies.  相似文献   

17.
基于部门间影响模式比较计划(ISI-MIP, Inter-Sectoral Impact Model Intercomparison Project)对CMIP5中5个气候(地球)系统模式模拟结果的降尺度数据,利用多模式集合预估了气候变化情景下21世纪环北极地区植被生长季与活动积温变化。研究发现:1)多模式集合模拟能够基本再现观测的初、终霜日及无霜期长度与>10°C积温的空间分布特征以及1979~2004年各指标变化趋势的空间分布特征,但其对气候变化年际变率的模拟能力较弱;2)至21世纪末,终霜日最多将提前60 d,初霜日将推迟20~40 d,无霜期延长幅度最高可达100 d,积温将增加1000~1200°C。其中RCP8.5情景下,各指标变幅最大,RCP2.6情景下变幅最小;3)各指标变幅呈现出较大的空间差异,亚欧大陆中西部的变幅普遍较大,随着气候变暖,>10°C积温增加幅度表现出明显的纬度地带性,南部增幅较大,北部增幅较小。  相似文献   

18.
This study explores natural and anthropogenic influences on the climate system, with an emphasis on the biogeophysical and biogeochemical effects of historical land cover change. The biogeophysical effect of land cover change is first subjected to a detailed sensitivity analysis in the context of the UVic Earth System Climate Model, a global climate model of intermediate complexity. Results show a global cooling in the range of –0.06 to –0.22 °C, though this effect is not found to be detectable in observed temperature trends. We then include the effects of natural forcings (volcanic aerosols, solar insolation variability and orbital changes) and other anthropogenic forcings (greenhouse gases and sulfate aerosols). Transient model runs from the year 1700 to 2000 are presented for each forcing individually as well as for combinations of forcings. We find that the UVic Model reproduces well the global temperature data when all forcings are included. These transient experiments are repeated using a dynamic vegetation model coupled interactively to the UVic Model. We find that dynamic vegetation acts as a positive feedback in the climate system for both the all-forcings and land cover change only model runs. Finally, the biogeochemical effect of land cover change is explored using a dynamically coupled inorganic ocean and terrestrial carbon cycle model. The carbon emissions from land cover change are found to enhance global temperatures by an amount that exceeds the biogeophysical cooling. The net effect of historical land cover change over this period is to increase global temperature by 0.15 °C.  相似文献   

19.
Although atmospheric greenhouse gas concentrations continuously increased, there was relatively little change in global-averaged surface temperatures from 1998 to 2013, which is known as atmospheric warming slowdown. For further understanding the mechanism involved, we explored the energy redistribution between the atmosphere and ocean in different latitudes and depths by using data analysis as well as simulations of a coupled atmosphere–ocean box model. The results revealed that, compared with observational changes of ocean heat content (OHC) associated with rapid warming, the OHC changes related to warming slowdown are relatively larger in multiple ocean basins, particularly in the deeper layer of the Atlantic. The coupled box model also showed that there is a larger increasing trend of OHC under the warming slowdown scenario than the rapid warming scenario. Particularly, during the warming slowdown period, the heat storage in the deeper ocean increases faster than the ocean heat uptake in the surface ocean. The simulations indicated that the warming patterns under the two scenarios are accompanied by distinct outgoing longwave radiation and atmospheric meridional heat transport, as well as other related processes, thus leading to different characteristics of ocean heat uptake. Due to the global energy balance, we suggest this slowdown has a tight relationship with the accelerated heat transport into the global ocean.  相似文献   

20.
Releases of halocarbons into the atmosphere over the last 50 years are among the factors that have contributed to changes in the Earth’s climate since pre-industrial times. Their individual and collective potential to contribute directly to surface climate change is usually gauged through calculation of their radiative efficiency, radiative forcing, and/or Global Warming Potential (GWP). For those halocarbons that contain chlorine and bromine, indirect effects on temperature via ozone layer depletion represent another way in which these gases affect climate. Further, halocarbons can also affect the temperature in the stratosphere. In this paper, we use a narrow-band radiative transfer model together with a range of climate models to examine the role of these gases on atmospheric temperatures in the stratosphere and troposphere. We evaluate in detail the halocarbon contributions to temperature changes at the tropical tropopause, and find that they have contributed a significant warming of ~0.4 K over the last 50 years, dominating the effect of the other well-mixed greenhouse gases at these levels. The fact that observed tropical temperatures have not warmed strongly suggests that other mechanisms may be countering this effect. In a climate model this warming of the tropopause layer is found to lead to a 6% smaller climate sensitivity for halocarbons on a globally averaged basis, compared to that for carbon dioxide changes. Using recent observations together with scenarios we also assess their past and predicted future direct and indirect roles on the evolution of surface temperature. We find that the indirect effect of stratospheric ozone depletion could have offset up to approximately half of the predicted past increases in surface temperature that would otherwise have occurred as a result of the direct effect of halocarbons. However, as ozone will likely recover in the next few decades, a slightly faster rate of warming should be expected from the net effect of halocarbons, and we find that together halocarbons could bring forward next century’s expected warming by ~20 years if future emissions projections are realized. In both the troposphere and stratosphere CFC-12 contributes most to the past temperature changes and the emissions projection considered suggest that HFC-134a could contribute most of the warming over the coming century.  相似文献   

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

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