首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到17条相似文献,搜索用时 218 毫秒
1.
研究西南地区陆地植被生态系统净初级生产力(NPP)的时空演变特征及其驱动力,对区域生态环境保护具有重要的现实意义.利用2000~2021年MODIS NPP、1999~2021年基于站点的气象数据和2000~2020年土地利用类型等数据,结合主成分分析、残差分析、Theil-Sen Median趋势分析和偏相关分析等方法,研究西南地区陆地植被生态系统NPP时空演变及其对驱动力的响应特征.结果表明,时间上,2000~2021年西南地区植被NPP呈波动上升趋势,速率为3.54g·(m2·a)-1.气候变化和人类活动影响下,农田、草地和森林生态系统NPP均呈上升趋势,但农田生态系统NPP的上升趋势最为显著.空间上,西南地区植被NPP呈上升趋势的面积占比为89.06%,显著上升和极显著上升的区域主要分布在广西南部、四川东部、重庆西部,以及云南和贵州交界处.气候变化和人类活动对西南地区植被生长具有双重影响,气候变化和人类活动影响下农田生态系统NPP呈上升趋势的面积占比均高于草地和森林生态系统.西南地区植被NPP与各气象因子的相关性呈明显地域差异.区域尺...  相似文献   

2.
利用GIMMS和MODIS两种遥感数据,分析了1982~2008年东北冻土区植被生长季平均NDVI的时空特征,并探讨不同类型冻土区和不同植被类型归一化差值植被指数(NDVI)对气候变化和CO2体积分数增加的响应.研究表明,不同冻土类型区植被生长季NDVI均值从大到小依次为:连续多年冻土区大片多年冻土区季节性冻土区岛状多年冻土区.东北冻土区不同植被类型生长季平均NDVI值由大到小依次为:森林灌丛沼泽农田草地,其中森林植被生长季平均NDVI值为0.61,草地为0.46.过去27年间,东北冻土区植被生长季平均NDVI年际变化曲线可分为3个变化阶段:①1982~1990年,小幅上升阶段;②1990~2000年,缓慢下降阶段;③2000~2008年,明显上升阶段.1982~2008年期间,连续多年冻土区及大片多年冻土区植被生长季平均NDVI值呈显著上升趋势(p0.05).对于不同植被类型而言,除森林植被NDVI呈显著上升趋势外(p0.05),其它植被类型NDVI值无显著变化趋势.过去27年间,东北冻土区年均气温显著升高,年降水量显著下降,CO2浓度显著升高.研究区全区平均NDVI与年平均气温呈显著正相关(p0.05),气温是影响东北冻土区生长季植被NDVI的主要气候因子.森林和沼泽湿地植被生长季平均NDVI与年平均气温呈显著正相关,与降水量呈显著负相关(p0.05);5种植被类型中仅森林植被受CO2浓度影响显著.年平均气温对不同植被类型的影响由高到低的顺序为:森林沼泽湿地灌丛农田草地;降水的影响为:森林沼泽湿地草地灌丛农田;CO2浓度的影响为:森林沼泽湿地草地农田灌丛.  相似文献   

3.
利用GIMMS和MODIS两种遥感数据,分析了1982~2008年东北冻土区植被生长季平均NDVI的时空特征.并探讨不同类型冻土区和不同植被类型归一化差值植被指数(NDVI)对气候变化和CO2体积分数增加的响应.研究表明.不同冻土类型区植被生长季NDVI均值从大到小依次为:连续多年冻土区>大片多年冻土区>季节性冻土区>岛状多年冻土区.东北冻土区不同植被类型生长季平均NDVI值由大到小依次为:森林>灌丛>沼泽>农田>草地,其中森林植被生长季平均NDVI值为0.61,草地为0.46,过去27年间,东北冻土区植被生长季平均NDVI年际变化曲线可分为3个变化阶段:①1982~1990年,小幅上升阶段;②1990~2000年,缓慢下降阶段;③2000~2008年.明显上升阶段.1982~2008年期间,连续多年冻土区及大片多年冻土区植被生长季平均NDVI值呈显著上升趋势((P<0.05)对于不同植被类型而言,除森林植被NDVI呈显著上升趋势外(P<0.05),其它植被类型NDVI值无显著变化趋势.过去27年间,东北冻土区年均气温显著升高,年降水量显著下降,CO2浓度显著升高.研究区全区平均NDVI与年平均气温呈显著正相关(P<0.05),气温是影响东北冻土区生长季植被NDVI的主要气候因子.森林和沼泽湿地植被生长季平均NDVI与年平均气温呈显著正相关,与降水量呈显著负相关(P<0.05);5种植被类型中仅森林植被受CO2浓度影响显著.年平均气温对不同植被类型的影响由高到低的顺序为:森林>沼泽湿地>灌丛>农田>草地;降水的影响为:森林>沼泽湿地>草地>灌丛>农田;CO2浓度的影响为:森林>沼泽湿地>草地>农田>灌丛.关健词:NDVI;气候变化;CO2浓度增加;东北冻土区  相似文献   

4.
长江流域植被净初级生产力对未来气候变化的响应   总被引:3,自引:0,他引:3  
苗茜  黄玫  李仁强 《自然资源学报》2010,25(8):1296-1305
研究基于气象观测和B2气候变化情景数据,利用大气-植被相互作用模型(AVIM2)模拟了1981—2000年和2010—2050年两个时段内植被NPP的空间分布格局及其时间变化趋势并分析了其时空变化与气温和降水量的关系。研究表明1981—2000年流域内植被NPP的空间分布大致呈现自西向东、自北向南递增的趋势。未来长江流域气温将整体增加,但各地增温幅度不同。流域降水量有增有减,主要增加区域位于长江源头和上游及中游的江北地区。未来在气温增加幅度较小而降水量增加的区域,如长江源头和上游的青海、西藏、川西及云南的部分地区的植被NPP将增加。在气温增幅较大而降水量减少或者降水量增加不多的区域如长江中游和下游的广大地区植被NPP将减少。从植被类型来看,长江流域大部分森林、郁闭灌丛和农作物的NPP在B2气候变化情景下将减少,每年减少量分别在0~4.5 gC.m-2、0~2 gC.m-2和0~2.5 gC.m-2之间。高寒草甸、草地和稀疏灌丛的NPP将增加,每年增长量介于0~2 gC.m-2之间。  相似文献   

5.
单振东  刘顿  骆汉  刘建伟  张丽梅  魏宇航 《环境科学》2023,44(11):6215-6225
承德接坝区位于农牧过渡区,对气候变化和人类活动极为敏感.以植被净初级生产力(NPP)作为评价指标,基于Thornthwaite Memorial模型计算潜在NPP和MODIS NPP遥感影像获取实际NPP数据,利用潜在NPP与实际NPP间的差值衡量人类活动作用下NPP的大小,运用Slope趋势和变异系数法分析实际NPP、潜在NPP和人类活动作用下NPP的变化趋势及稳定性分布,并采用相关系数法分析实际NPP与年降水量和年均气温间的相关性,最终量化气候变化和人类活动对该区域植被变化的影响.结果表明,潜在NPP自西北向东南递增,其变化趋势和稳定性均为自西向东递增.实际NPP与年降水量和年均气温呈正相关区域面积占比分别是99.87%和91.66%.该区域99.85%的植被得到改善且变化稳定,主要是由气候因素和人类活动共同主导(99.71%),而植被退化完全是由人为因素所导致(0.15%).  相似文献   

6.
基于MOD17数据,采用趋势分析和相关分析法,研究了2000~2019年粤港澳大湾区植被生产力和植被碳利用率(CUE)变化及其与气候变化的关系.结果表明,粤港澳大湾区总初级生产力(GPP)和净初级生产力(NPP)均值分别为1.80和0.89kg·C/m2,呈中部低四周高的空间格局,CUE均值为0.51,呈中部和东南部略高于四周的空间格局;GPP和NPP总体呈增加趋势,其年际变化率均值分别为0.01和0.001kg·C/(m2·a),GPP年增长率远高于NPP年增长率,不同植被类型GPP和NPP也呈增加趋势.植被CUE则呈逐年下降趋势,其年际变化率均值为-0.002a-1,由于气候因子变化对光合作用速率和呼吸作用速率产生的不同程度影响,从未来发展趋势看,68.22%的区域表明CUE仍呈下降趋势,植被碳固定能力减弱;不同植被类型CUE存在差异,其中农田CUE平均值最高,为(0.511±0.014),森林和草地CUE分别为(0.500±0.019)和(0.501±0.020),从变化趋势看,不同植被类型CUE呈极显著下降趋势(P < 0.01);GPP与气温、累积降雨和累积净太阳辐射量总体呈正相关关系,其所占比例分别为94.52%、53.36%和90.58%,NPP与气温、累积降雨和累积净太阳辐射量总体上也表现为正相关关系,其所占比例分别为86.86%、71.10%和85.97%.然而,CUE与气候因子的关系却有所不同.CUE与气温和累积降雨呈正相关关系,但与累积太阳辐射呈负相关关系,不同植被类型GPP、NPP和CUE与气温、累积降雨量呈正相关关系,GPP和NPP与累计净太阳辐射呈正相关关系,但CUE与累积太阳辐射却呈负相关关系.  相似文献   

7.
利用1982~2006年的NOAA AVHRR-GIMMS和MODIS 2种数据集的归一化植被指数(NDVI)数据对东北多年冻土区植被NDVI年际动态和空间差异进行分析,并结合气象数据和土地利用/覆被数据分析了植被NDVI对气候变化和土地利用/覆被变化的响应.研究表明,东北多年冻土区植被NDVI值较高,且空间差异明显;森林为该区主要植被类型,NDVI值较高,主要分布于大小兴安岭和伊春地区;草地集中分布于西南部, NDVI值相对较低.东北多年冻土区过去25a间植被生长的变化趋势为:伴随着气温的显著升高和降水量减少,植被NDVI显著下降.较气温而言,降水量是影响植被NDVI的主要因子(r = 0.77, P < 0.01).在气候变化和人类活动的双重作用下,东北多年冻土区植被NDVI在1982~2006年间表现为4个阶段:1982~1990年,植被NDVI虽有小幅波动,但整体上呈持续增加的趋势;1990~1993年,植被NDVI呈迅速下降趋势; 1993~1997年,植被NDVI呈现回升态势,表现出缓慢上升的趋势;1998~2006年,植被NDVI呈现总体下降趋势.不同植被类型表现出不同的NDVI年际变化规律,尤以草地NDVI值波动最大.植被NDVI变化空间异质性显著.气候变化和多年冻土退化影响了东北多年冻土区植被NDVI动态.年均气温升高和年降水量降低影响了植被的生长.从像元尺度来看,研究区植被NDVI与气温和降水均具有较显著的相关性.研究区土地利用/覆被变化的分析结果表明,不同的土地利用类型间的转变对植被NDVI的大小和空间分布产生了重要影响.  相似文献   

8.
张宏 《自然资源学报》2001,16(3):216-220
利用作者建立的塔里木盆地北部盐化草甸植被净初级生产力模型 ,模拟了极端干旱气候下盐化草甸植被净初级生产力对全球变化的响应。结果表明 ,在一定土壤质地条件下 ,盐化草甸植被净初级生产力随地下水埋深的增加而逐渐下降。地下水埋深越大 ,盐化草甸植被净初级生产力对地下水埋深变化的响应越敏感。全球变化造成的温度升高对盐化草甸植被净初级生产力的影响 ,也是依地下水埋深的不同而有所差异。地下水埋深较小时,盐化草甸植被净初级生产力对温度升高的反应较小 ;随着地下水埋深的加大,其响应程度明显增加。因而 ,地下水埋深越大 ,盐化草甸植被净初级生产力对全球变化的响应就越明显。  相似文献   

9.
基于空间尺度效应的西南地区植被NPP影响因子探测   总被引:2,自引:1,他引:1  
植被净初级生产力(NPP)是评价陆地生态系统质量的重要参数,研究植被NPP时空演变特征及其驱动力对区域生态环境保护和可持续发展具有重大意义.基于MODIS NPP数据、气象数据、 DEM数据、人口密度数据、 GDP数据和土地利用类型数据,采用一元线性回归分析、 R/S分析和地理探测器模型,分析西南地区及其六大地貌单元植被NPP时空演变特征及未来变化趋势,探究植被NPP空间分异的影响因子.结果表明,2000~2020年西南地区植被NPP整体呈极显著上升趋势.地貌单元中,除青藏高原南部外,其余地貌单元植被NPP均表现为改善态势,其中四川盆地和云贵高原表现为极显著改善.西南地区的植被NPP变化斜率整体呈现“东高西低”的分布格局.西南地区及各地貌单元植被NPP呈上升趋势的区域面积均大于呈下降趋势的区域面积,但未来植被NPP变化趋势均以下降为主.地理探测器结果表明,除云贵高原植被NPP空间分异主要受气温影响外,海拔是西南地区及各地貌单元植被NPP空间分异的主导因子.交互探测结果表明,影响因子之间的交互作用均表现为双因子增强或非线性增强,其中,海拔∩温度对西南地区植被NPP空间分异的解释力最大.地...  相似文献   

10.
采用M-K趋势检验与R/S分析法,探究了近70a全球全年及各季节的降水、气温、蒸发量的演变特征与空间分布格局.研究发现:(1)1948~2016年间,全球72.7%地区的降水呈现出非显著的上升或下降趋势,全球总降水量呈现上升趋势,趋势率为1.9mm/10a;全球气温呈现显著性、持续的上升趋势,趋势率为0.23℃/10a,1980年后气温上升速率变快;1980~2016年间全球蒸发量在大部分地区呈现上升趋势,总蒸发量的变化趋势率为8.2mm/10a;Hurst指数显示气温与蒸发变化的持续性明显大于降水量的变化;(2)降水量在北半球高纬地区多呈现出非常显著的增加趋势,在低纬度地区多呈现波动或下降趋势,且在12月~次年5月(DJF、MAM)的上升趋势的显著性普遍高于其他两个季节;气温在DJF和MAM的上升趋势相对微弱,中高纬地区呈现非显著性的变化趋势,亚洲中部部分地区在JJA呈现出显著性的下降趋势;蒸发量在沿海湿润地区的上升趋势显著,美洲北部在DJF呈现出显著下降的趋势,格陵兰岛、尼罗河流域在一年四季均为下降趋势;(3)各大洲的气温在1948~2016年均呈现出显著的上升趋势,其中北美洲的平均上升率最高;降水除非洲外均为上升趋势,非洲地区的降水量呈现出下降趋势,南美洲降水呈现出3个阶段的下降趋势;各大洲的蒸发量均呈现上升趋势,其中欧洲的平均上升率最高;除大洋洲外,各大洲降水量的上升速率低于蒸发量的上升速率.  相似文献   

11.
草原区作为影响我国陆地生态系统碳水循环和生态安全的重点区域,对气候变化和人类活动极为敏感.然而,有关气候变化和人类活动对该区域植被恢复相对贡献的认识尚存分歧.以生态系统净初级生产力(NPP)为评价指标,通过对比MODIS观测的实际NPP和基于Thornthwaite Memorial模型估算的潜在NPP的趋势差异,量化了2000~2020年气候变化和人类活动对我国典型草原区(北方温性草原区和青藏高原高寒草原区)植被恢复的相对作用.结果表明,研究区内93%的草地植被呈恢复趋势,NPP平均增加速率达(以C计)2.12 g·(m2·a)-1(P<0.01),其中,近一半植被恢复区受气候变化和人为活动共同控制,约36%和10%植被恢复区分别受气候变化和人类活动的独立控制;此外,不同草地类型气候变化主导植被恢复的面积占比差异大,主要表现为高寒草地明显高于温性草地,气候条件越干旱,气候主导面积占比越大.人类活动不是北方温性草原区和青藏高原高寒草原区植被恢复的主要原因,但在气候条件恶化地区,人类活动可降低甚至抵消气候变化对植被的负面影响.未来需加强长...  相似文献   

12.
朱士华  艳燕  邵华  李超凡 《自然资源学报》2017,32(11):1844-1856
中亚干旱区分布着世界80%以上的温带荒漠,受气候变化影响显著。论文首先收集实验观测数据验证了干旱区生态系统模型(AEM),然后运用AEM开展数值模拟实验量化研究了1980—2014年中亚净初级生产力(NPP)的时空格局,评估了不同环境因子(降水、温度、CO2)的相对贡献率及其交互效应。结果表明:过去35 a中亚干旱区年均NPP总量为1 125±129 Tg C(1 T=1012)或218±25 g C/m2。哈萨克斯坦北部地区年NPP值较高(349±39 g C/m2),而南疆地区年NPP值较低(123±45 g C/m2)。1980—2014年间,中亚NPP总体呈减少趋势 [-0.71 g C/(m2·a)],南疆极端干旱区的NPP降低最为显著 [-2.05 g C/(m2·a)]。相较于1980—1984年NPP均值,在1985—2014年中亚区域NPP总体降低了118 Tg(-10%)。其中CO2施肥效应促进NPP增加了99.7 Tg (+8%),气温升高的正效应促进NPP增加了35.4 Tg(+2%),而降水减少导致NPP降低了221 Tg(-18%)。研究区内9%的地区的NPP主要控制因子为温度,主要分布在天山和哈萨克斯坦北部等高纬高寒地区。降水主控区面积占整个研究区的69%,主要分布在荒漠平原特别是南疆等植被受水分限制的区域。CO2主控区占研究区面积的20%,主要分布在天山中山带森林区和低海拔地区等水热条件好的区域。研究表明新疆南部地区是中亚的关键生态脆弱区,其生态安全面临着气候变化的挑战,但21世纪的升温不大可能因刺激自养呼吸而对中亚区域NPP造成显著影响。  相似文献   

13.
全球变化背景下,青藏高原作为我国乃至全球气候变化的“天然实验室”,植被生态系统发生了深刻变化。引入重心模型等方法分析和探讨2000—2015年青藏高原植被NPP时空变化格局及其驱动机理,并定量区分NPP变化过程中气候变化和人类活动的相对作用。研究发现:(1)2000—2015年,青藏高原植被NPP年均值总体上呈现从东南向西北递减的趋势。在年际变化方面,近16年植被NPP呈现波动上升趋势,其中在2005年出现上升陡坡,并在2005—2015年表现为高位波动的态势。(2)青藏高原植被NPP增加区(变化率>10%)主要集中于三江源地区、横断山区北部、雅鲁藏布江中下游以及那曲地区的中东部,而植被NPP减小区(变化率<-10%)则主要分布于雅鲁藏布江上游和阿里高原。(3)近16年青藏高原植被NPP重心总体向西南方向移动,表明西南部植被NPP在增量和增速上大于东北部。(4)青藏高原植被NPP与气候因子相关性的地区差异显著,其中植被NPP与降水显著相关的区域主要位于青藏高原中部、青藏高原东南部及雅鲁藏布江流域中下游,而植被NPP与气温显著相关的区域主要位于藏南地区、横断山区北部、青藏高原中部和北部。(5)气候变化和人类活动在青藏高原植被NPP变化过程中的相对作用存在显著的时空差异性,在空间上呈现“四线—五区”的格局。研究成果能够为揭示青藏高原区域生态系统对全球变化的响应机制提供理论和方法支撑。  相似文献   

14.
Shifts in the carbon balance of high-latitude ecosystems could result from differential responses of vegetation and soil processes to changing moisture and temperature regimes and to a lengthening of the growing season. Although shrub expansion and northward movement of treeline should increase carbon inputs, the effects of these vegetation changes on net carbon exchange have not been evaluated. We selected low shrub, tall shrub, and forest tundra sites near treeline in northwestern Alaska, representing the major structural transitions expected in response to warming. In these sites, we measured aboveground net primary production (ANPP) and vegetation and soil carbon and nitrogen pools, and used these data to parameterize the Terrestrial Ecosystem Model. We simulated the response of carbon balance components to air temperature and precipitation trends during 1981–2000. In areas experiencing warmer and dryer conditions, Net Primary Production (NPP) decreased and heterotrophic respiration (R H ) increased, leading to a decrease in Net Ecosystem Production (NEP). In warmer and wetter conditions NPP increased, but the response was exceeded by an increase in R H ; therefore, NEP also decreased. Lastly, in colder and wetter regions, the increase in NPP exceeded a small decline in R H , leading to an increase in NEP. The net effect for the region was a slight gain in ecosystem carbon storage over the 20 year period. This research highlights the potential importance of spatial variability in ecosystem responses to climate change in assessing the response of carbon storage in northern Alaska over the last two decades.  相似文献   

15.
湖北省神农架林区是全国唯一以林区命名的行政单位,拥有全球中纬度地区唯一一块保存完好的原始森林,量化其森林生态系统NPP(net primary productivity,净初级生产力)对县域生态系统评估工作十分重要.基于CEVSA2(carbon exchange between vegetation,soil and atmosphere 2)模型模拟1981-2015年神农架林区森林生态系统NPP,并利用中国生态系统研究网络神农架站观测数据和野外调查数据进行验证,进而分析其NPP时空变化特征及其主要环境影响因子.结果表明:①1981-2015年神农架林区森林生态系统年均NPP为628.27 g/m2(以C计),空间分布表现为中部较低、东部以及周边较高,具有明显的空间异质性.②1981-2015年神农架林区年均NPP的增长速率为2.58 g/(m2·a)(R2=0.65,P < 0.001);运用Mann-Kendall突变检验法发现,1998年前后是NPP增长速率变化的突变点,虽然1998年前后两个时段NPP均呈上升趋势,但1999-2015年NPP增长率较1981-1998年下降了7.01%;从空间上来看,林区中部和北部NPP增长率[4~6 g/(m2·a)]相对较高,南部和东部部分地区NPP呈下降趋势,其变化速率在-3~0 g/(m2·a)之间.③神农架林区NPP年际变化与年均温、总辐射年际变化均呈正相关,与年降水量年际变化呈负相关,其中年均温年际变化对NPP年际变化的解释率最高,为43%(P < 0.01);在空间尺度上,林区森林生态系统约67.83%区域的NPP年际变化由年均温年际变化控制,主要分布在林区中部和东南部,可见年均温是该地区森林生态系统NPP的主要影响因素.   相似文献   

16.
以遥感和气象数据为主要数据源,运用改进后的CASA模型,估算了俄罗斯布里亚特共和国2000-2008年的植被NPP,并验证了模型的精度,分析了该地区植被NPP的时空变化规律及其与气候因子的相互关系。研究结果表明:时间上,植被NPP年际上呈现为在波动中上升,月份上表现为先升后降的趋势;空间上,植被NPP随经度的增加而增大,随纬度的增加而减小,由西南到东北逐渐递增的趋势;不同植被类型的NPP也不同,从大到小依此为:草地与沼泽林>森林>森林与草原>稀树草原>高山植被。其变化主要受气温和降水量变化的作用。改进后的CASA模型运用于布里亚特共和国植被NPP估算的精度较高。该研究对中国北方植被NPP估算和生态跨境研究具有借鉴意义。  相似文献   

17.
柴达木盆地植被生长时空变化特征及其对气候要素的响应   总被引:7,自引:0,他引:7  
以2001-2010 年MODIS NDVI 植被数据为基础,并结合该区同期气温和降水量数据,采用线性趋势分析、偏相关分析、Hurst 指数等数理分析方法,研究了近10 a 来柴达木盆地植被时空变化特征、影响因素及未来可能的变化趋势。结果表明:①2001-2010 年柴达木盆地植被生长季NDVI平均值(NDVI)呈显著增加趋势,线性增长率为0.041/10 a。柴达木盆地主要植被类型灌丛、草原、荒漠NDVI 的线性增长率分别为0.043/10 a、0.034/10 a、0.028/10 a;②柴达木盆地植被变化具有阶段性特征,整体呈“S”型增长,具有两次明显的植被快速增长期;③柴达木盆地植被改善面积占研究区总面积的12.43%,并主要分布在布尔汉布达山、祁漫塔格山、鄂拉山、柴达木山、宗务隆山的高山区,冬给措纳湖周边和绿洲核心区。植被退化面积仅占研究区总面积的0.35%,并主要分布在绿洲边缘区,特别是柴达木河绿洲和那陵格勒河绿洲的中下游地区;④植被NDVI 与1-5 月平均气温以及5-8 月可利用降水量呈显著正相关关系,植被生长与温度呈显著正相关的区域主要分布在绿洲核心区及高山区,占植被区面积的8.36%,植被生长与可利用降水量呈显著正相关的区域主要分布在低山区及山地河谷地带,占植被区面积的30.95%;⑤植被生长季延长和生长加速是柴达木盆地植被NDVI 增加的主要原因,气候的暖湿化是促使柴达木盆地植被改善的主要驱动力;⑥柴达木盆地植被改善具有强持续性,未来大部分区域植被将持续改善。  相似文献   

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

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