首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 93 毫秒
1.
太行山区植被NPP时空变化特征及其驱动力分析   总被引:2,自引:2,他引:2  
本文基于2000—2014年MODIS NPP数据,结合同期土地利用变化、气温、降水和DEM数据,运用趋势分析法、相关系数法及分区统计法等方法,研究了太行山区2000—2014年植被NPP时空变化特征,分析了气温、降水等气候因素和人为因素对植被NPP变化的影响,为太行山区植被资源管理及生态环境调控提供参考。研究结果表明:(1)太行山区植被NPP多年平均值为284.0 g(C)·m~(-2)·a~(-1),耕地、林地和草地的NPP均值分别为302.5 g(C)·m~(-2)·a~(-1)、258.1 g(C)·m~(-2)·a~(-1)、286.5 g(C)·m~(-2)·a~(-1)。(2)2000—2014年太行山区植被NPP整体呈上升趋势,但大部分植被NPP变化未达到显著水平;16.17%的植被NPP显著或极显著升高,主要分布在太行山区西侧;0.88%的植被NPP显著或极显著降低,零散分布在研究区内。(3)不同植被类型NPP变化速率为草地耕地林地。(4)基于区域平均计算,太行山区植被NPP与降水显著正相关(P0.05),与气温负相关(P0.05)。基于像元计算,植被NPP与降水显著或极显著正相关区面积比例为23.82%,主要分布在太行山区北段,几乎没有显著负相关区;植被NPP与气温显著或极显著负相关区面积比例为8.42%,主要分布在太行山区西侧,显著或极显著正相关区面积比例为0.81%,主要分布在太行山区最北端。(5)研究期内气候因子对植被NPP的升高整体上表现为促进作用,而人为因素主要表现为抑制作用。太行山区生态环境保护仍应以减少人为干扰为主。  相似文献   

2.
吉林省长时间序列陆地生态系统NPP时空变化特征   总被引:1,自引:0,他引:1  
植被净初级生产力(NPP)是陆地生态系统碳循环过程中的重要组成部分,也是评价碳源/汇的重要环节.基于2001—2019年MODIS地表反射率数据(MOD09A1)、气温和日照时间资料,利用光能利用率模型估算了吉林省陆地生态系统NPP,采用趋势分析、空间变异性分析等方法,分析了吉林省陆地生态系统NPP的时空变化特征及其与...  相似文献   

3.
对川西植被净初级生产力进行估算并分析了其时空格局变化及驱动因子与植被净初级生产力的变化关系,为深入认识川西植被生产力状况提供科学依据。在原CASA模型的基础上,通过区域实际蒸散量与区域潜在蒸散量对水分胁迫影响系数进行了改进,进而估算2000—2015年期间川西地区植被净初级生产力;运用逐像元趋势法分析了川西地区植被净初级生产力在研究期间内其空间变化情况;同时结合气象因子、土地利用变化、植被类型、地形因素、人类活动与植被净初级生产力进行了相关性分析。结果表明:川西地区植被空间差异明显,其NPP多年平均值为199 gC/(m2·a),在2000—2015年期间,大面积区域植被NPP呈显著上升,汶川、泸定、金川、康定局部区域受自然灾害及人类活动等因素NPP呈下降趋势;降水、气温、太阳辐射等气象因子对植被空间格局分布产生一定影响,不同土地利用及植被类型的NPP差异较大;海拔与研究区NPP相关性非常显著(R2=0.896,p<0.001);人类活动对汶川、泸定等局部地区负干扰明显。  相似文献   

4.
为了探明湖北省植被NPP空间分布特征及其驱动机制,利用MOD17A3H NPP时序数据、土地覆盖类型产品MCD12Q1、气象数据和SRTM DEM地形数据等,采用线性趋势分析、相关系数等方法,分析了2004-2015年湖北省植被NPP时空分布特征,探讨了不同土地利用类型、不同地形条件下植被NPP分布及对气候因素的响应特征。结果显示:(1)12年间,湖北省植被NPP总体呈上升趋势;植被NPP的空间分布大体呈东、西部山区高,中部略低的格局,全省植被NPP值大部介于400~600 g C/(m2·a);(2)全省植被NPP与年均气温呈现正相关性的面积占全省总面积的79.28%,与年降水量呈现正相关性的面积占全省总面积的59.07%;(3)气温是影响林地、草地植被NPP变化的主要驱动因子,耕地主要受降水量影响;(4)植被NPP与年均气温和年降水量的相关性随高程、坡度、坡向的不同表现出差异性,植被NPP变化是气温、降水和地形条件综合作用的结果。  相似文献   

5.
[目的]研究长江流域陆地植被生态系统初级净生产力(NPP)的时空演变特征,为政府部门建立和调整生态功能恢复项目提供参考。[方法]以长江流域为研究区,基于2000—2019年的降水和气温数据,采用周广胜—张新时模型(ZGS)和Thornthwaite Memorial (TM)模型估算NPP,并进一步利用皮尔逊相关分析、一元线性回归分析、Mann-Kendall显著性检验等,对比分析长江流域陆地植被生态系统NPP时空演变特征。[结果](1)基于上述两个模型模拟得到的长江流域NPP时空演变趋势基本一致,相关系数R为0.982,呈现显著正相关关系;(2)2000—2019年长江流域陆地植被生态系统实际NPP与潜在NPP均呈上升趋势,上升速率分别为6.85,2.74 g/(m~2·a)。(3)长江流域实际NPP和潜在NPP在空间上呈东南高西北低的分布格局,低值区域主要分布在草地生态系统;高值区域大部分分布在森林生态系统和农田生态系统。(4)2000—2019年长江流域实际NPP与潜在NPP呈上升趋势的面积分别占研究区总面积的80.65%和84.81%,主要分布在云南、青海、西藏、四川北部及浙江...  相似文献   

6.
李肖  袁金国    孟丹 《水土保持研究》2018,25(6):109-114,120
在GIS软件的支持下,利用2005—2014年MOD17A3H与MCD12Q1数据集,结合河北省2005—2012年的142个气象站点资料,利用一元线性回归法分析了河北省10年间NPP的时空演变特征,利用相关系数法分析了NPP与气候因子的相关关系。结果表明:河北省年NPP主要集中在200~400 g/(m2·a),其中NPP高值区分布在燕山北部区域和太行山东侧一带。2005—2014年河北省年NPP均值与最大值均在波动中微弱上升,年NPP均值在2012年最大为337.63 g/(m2·a)。2005—2014年NPP增长率大于10%的地区分布在冀西北区域和东南边界,其面积占比为43.03%。河北省5种植被覆盖类型按NPP年均值从大到小排列依次为:阔叶林、灌丛、草地、针叶林、农用地。河北省NPP与气温主要呈负相关,与降水主要呈正相关。  相似文献   

7.
土地利用/覆被变化是生态系统能量流动与物质循环的主要驱动力之一,土地利用/覆被变化最直观的表现就是植被覆盖度的变化,研究其对植被净初级生产力的影响,对区域土地利用的管理具有重要意义.研究利用MO-DIS反射率及NPP数据,分析了2008-2014年黑河流域植被覆盖度变化及其对NPP的影响.结果表明:(1) 2008-2014年,黑河流域年均NPP呈现上升的趋势,流域NPP总体呈现南多北少,上游地区的祁连山地带、中游地区的张掖、酒泉、临泽等地NPP值普遍较高,除沿河道外,下游NPP值普遍较低,并且在下游的北部地区NPP值达到了最低.(2)相比2008年、2014年黑河流域Ⅰ级(盖度小于20%)、Ⅳ级(盖度介于40%~60%)和Ⅴ级(盖度大于60%)植被覆盖等级所占比例有所上升,而Ⅱ级(盖度介于20% ~30%)和Ⅲ级(盖度介于30%~40%)比例下降,影响该流域植被生长的主导性因子不同,上游地区主要受降水影响,下游受气温影响较大.(3)植被覆盖度对NPP的产量具有很好的指示作用,二者的决定系数高达0.90,NPP产量随植被覆盖度的增大而增加.(4) 2008-2014年,黑河流域NPP总量由6 052.82 GgC/a增加为6 624.54 GgC/a,增加幅度为9.45%,其中,面积减少17 086 km2的低植被覆盖区NPP总量减少仅为25.65 GgC/a,而高植被覆盖度区增加最多,从2 861.96 GgC/a增加为3 206.86 GgC/a,增加量为344.90 GgC/a,较高植被覆盖度区的NPP总量增加次多,增加了251.84 GgC/a,增幅为12.57%.  相似文献   

8.
植被净初级生产力(NPP)可以直接反映植被在自然环境中的生产能力,利用遥感影像、气温降水数据结合简单差值法、趋势分析法以及线性回归分析对呼伦贝尔沙地河流湿地NPP进行了时空变化研究与驱动力分析.结果表明:时间上,呼伦贝尔沙地河流湿地NPP呈现波动式增长,NPP增长面积可达整个湿地面积的92.51%,年增长率为0~2 gC/(m2·a);空间上,河流湿地NPP值呈现东高西低格局.而15年间NPP增加最快的湿地为处于西部的克鲁伦河湿地,增长面积可达99.5%;气候变化与人类活动是导致NPP变化的两个重要因素,NPP增长与降水量呈显著正相关(R=0.703),与气温呈显著负相关(R=-0.559).经济产业结构的变化使得人类活动对生态环境的压力减小,促使湿地NPP的增加.同时,国家政策以及治沙工程同样影响着NPP的变化.  相似文献   

9.
2001—2010年中国农田生态系统NPP的时空演变   总被引:3,自引:0,他引:3  
植被净初级生产力(Net Primary Productivity,NPP)是地表碳循环的重要组成部分,而农业NPP代表了农田生态系统通过光合作用可固定大气中CO2的能力,决定了农田土壤可获得的有机碳含量.测算中国农田生态系统产生的NPP并分析其时空变异规律,探索其影响因子,对于了解全球碳循环,预测未来生态变化趋势具有重要意义.基于MOD17A3数据集,利用一元线性回归趋势线法和相关系数法定量分析了农田生态系统NPP的时空演变特征和气候因子对农田NPP的影响.结果表明,2001-2010年中国农田NPP平均值变化范围为C 0.21 ~ 17.24 Mg hm-2a-1,平均值为4.12 Mg hm-2 a-1.从时间变化来看,NPP年际变化呈现先增加后下降的趋势.从空间分布来看,表现为南方高,北方低,其中,甘新区、黄土高原区、东北区、内蒙古及长城沿线区和黄淮海区北部NPP值较低,<3 Mg hm-2 a-1;四川盆地、山东丘陵、长江中下游滨海和沿江平原、云南、贵州、海南等地NPP值较高,>5 Mghm-2 a-1.从变化趋势来看,黄淮海农业区、黄土高原区、新疆的绿洲农业区和长江中下游农业区南部上升趋势明显,东北农业区、四川盆地和长江中下游农业区下降趋势明显.2001-2010年,中国农田中有22%的像元NPP平均值与降水呈显著相关,7%的像元与气温呈显著相关,反映了农田生态系统受人类活动影响较大.  相似文献   

10.

气候变化和人类活动对植被NPP( net primary productivity,净初级生产力)的驱动作用是全球气候变化背景下的研究热点,并且在不同的时空尺度上尚不能达成共识。中国西北寒旱牧区植被生态系统脆弱,对气候变化和人类活动的响应十分敏感,该研究以青海省乌兰县作为代表性研究区,采用GIMMS-NDVI and MOD-NDVI数据融合构建了长时序归一化植被指数(normalized difference vegetation index,NDVI)数据集,并结合CASA(carnegie-ames-stanford approach)模型获取了研究区1982—2020年植被生长季NPP,利用Sen+MK趋势分析方法探究了研究区植被生长季NPP的时空演变特征,同时采用构建的ADE+Sen量化归因方法对多种气候要素和人类活动的驱动作用进行了定量分析。结果表明,研究区植被生长季NPP多年均值为(205.9 ± 11.5)g/(m2·a)(以C计),年际变化无显著趋势,不同植被类型的生长季NPP年际波动过程与全域生长季NPP基本一致。在空间上,植被生长季NPP自西向东逐渐增加,年际变化趋势具有明显的空间异质性,且整体以退化为主,平均变化率为−0.151 g/(m2·a2),其中表现出严重退化和轻度退化的面积占比分别达到了31.7%和29.5%。气候变化主导的植被面积占比达85.2%,其贡献值的绝对平均值为1.025 g/(m2·a2),约是人类活动贡献值绝对平均值的2倍,太阳辐射、降水量、平均气温和平均风速均是影响植被生长季NPP动态的主要气候因素。该研究表明研究区植被整体表现为退化状态,气候变化是导致该现象的主导驱动因素。研究结果可为中国西北寒旱牧区植被生态系统的可持续利用和保护提供参考。

  相似文献   

11.
植被净初级生产力对评价全球变化背景下植被生长状况、陆地生态系统质量、资源环境监测具有重要意义。基于2000—2015年的MODIS NPP产品MOD17A3数据,结合DEM数据、气象数据、土地利用数据,运用趋势分析、相关性分析、地理探测器等模型方法,探讨长江流域片植被NPP的时空特征,揭示了各驱动因子的贡献率。结果表明:(1)16 a间植被NPP均值在478.4~547.4 gC/(m2·a),平均值为516.5 gC/(m2·a),流域内NPP整体表现为缓慢上升趋势。(2)流域内NPP空间分布格局为自东南向西北减少,16 a间大部分地区NPP值基本不变,云南省迪庆州、贵州省毕节市西部增长最明显,下降区域多分布于贵州省东部。(3)研究区内植被NPP与气温、降水均呈正相关; NPP随海拔增加呈先增加后下降的趋势,2 000~3 500 m海拔范围内植被NPP值最高,且集中于横断山区; 坡度小于15°的区域对整个研究区NPP值贡献最大; 林地面积最大且NPP均值最大,对研究区的植被NPP值贡献最大。(4)各因子对NPP的解释力排序为海拔>气温>降水>土地利用类型>坡度,单因子海拔对NPP的解释力最强,海拔与气温双因子交互作用对NPP的解释力最强。研究结果可为长江流域生态修复及可持续发展提供数据支持。  相似文献   

12.
利用宁夏回族自治区20个气象站自建站以来到2012年共52 a的降水数据资料,运用趋势线、滑动平均等数理分析方法对宁夏52 a以来的降水变化特征进行了分析,并应用降水Z指数和正负20%的降水距平指标对宁夏的旱涝灾害特征进行了分析。最后基于信息预测理论的对称性及可公度法,对宁夏未来的旱涝灾害进行了趋势判断。结果表明:(1)宁夏52 a年降水量以8.57 mm/10 a的速度递减,从20世纪90年代起干旱化趋势明显。(2)宁夏旱涝灾害交替发生,旱涝出现的频率为32.7%和28.7%,降水正常年份仅占38.5%。(3)宁夏旱涝灾害的发生具有显著对称性周期特征,未来发生旱灾的可能年份是2014和2015年,发生涝灾的可能年份是2017年。  相似文献   

13.
植被NPP是判定生态系统碳源/碳汇及调节生态过程的主要因子,基于CASA模型估算NPP,探究时空尺度NPP的变化及其对气候变化的响应状况,可了解泾河流域植被恢复状况并为流域生态环境改善提供科学参考及建议。以泾河流域为研究区域,基于2000年、2009年、2018年MODIS NDVI数据、气象数据与植被分布数据等,运用CASA模型(Carnegie-Ames-Stanford approach,CASA)对2000年、2009年及2018年流域NPP进行了估算,探究流域内NPP时空动态及其对气候变化的响应,并分析了各植被类型下NPP分布规律。结果表明:(1)时间尺度上,2000年、2009年及2018年泾河流域NPP均值分别为521.81 gC/(m2·a),664.77 gC/(m2·a),719.78 gC/(m2·a),年际变化呈增长趋势; 月际变化曲线呈单峰型,4—8月呈较强上升趋势,8月后逐渐下降; 各季节的NPP均值由高到低依次为夏季>秋季>春季>冬季。(2)空间尺度上,NPP分布存在一定的地域差异性,水平方向呈“南高北低”的特点; 垂直方向上,NPP值随海拔高度的升高呈先下降后上升的趋势。(3)不同植被类型下NPP均值存在明显差异,其中常绿阔叶林年均NPP值最高,为1 544.50 gC/(m2·a)。(4)气候变化背景下,NPP主要受气温与降水的影响,且降水为主导因素。研究结果表明泾河流域NPP呈增加趋势,即植被覆盖情况在逐步改善; 流域北部植被覆盖状况仍有待改善,建议加大退耕还林政策实施力度,加大果树、茶树等防护型林地的种植; 且该流域NPP对降水的响应强于气温,故可加大植树造林恢复植被的力度,也可修建水库和水利工程,退田还湖,并加大湿地保护,从而保证空气湿度,增加降水,改善植被覆盖,实现人与自然和谐共生。  相似文献   

14.
根据宁夏回族自治区中卫市兴仁气象站1961-2008年逐月气温、降水数据,采用Miami和Thornthwaite Memorial模型模拟了中卫市环香山地区气候生产潜力及其变化。结果表明:(1) 环香山地区气候向"暖干型"方向变化,年平均气温呈明显上升趋势,降水量变化波动大,呈逐年递减的趋势;(2) 温度生产潜力逐年增加,降水生产潜力和气候生产潜力呈波动中减少的趋势;(3) 气候生产潜力与年平均降水量具有着极其显著的线性相关关系,与年平均气温没有明显的相关性,温度生产潜力分别是降水生产潜力和气候生产潜力的2.79和2.35倍,说明降水条件是制约环香山地区气候生产潜力的关键因子;(4) 对未来气候进行情景模拟分析,结果表明环香山地区气候主要以"暖干型"气候为主,这对当地农业生产极为不利,需大力发展节水农业技术以保证农业生产。  相似文献   

15.
植被碳利用效率(CUE)是评估陆地碳循环的重要指标,探讨其动态特征对气候变化的响应对于陆地生态系统碳循环研究具有重要的指示意义。基于MOD17数据计算中国植被CUE,辅以地统计学理论,利用趋势分析、变异系数及相关性分析等方法,研究了2000—2013年中国植被CUE的时空变化特征,并结合气象要素数据剖析植被CUE对气候变化的响应。结果表明:全国植被年CUE在研究年限内总体上表现为轻微的增长趋势,具体表现为在2000—2007年以0.000 6的变化率呈现波动下降的趋势,而2007—2013年植被CUE呈现显著上升的变化趋势,该变化特征可归因于2007年气温及降水格局的改变;植被CUE空间分布具有明显的空间异质性,大体呈现西部高东部低的状态。CUE增加较明显的区域主要分布在内蒙古呼伦贝尔地区,青藏高原大部分地区,东部沿海地区以及台湾岛。而CUE呈减少趋势的区域分布范围较广,其中减少趋势较为明显的地区主要包括东北平原及华北华中的大部分地区,另江南地区及新疆部分地区也有零星分布。草地及森林区域植被CUE波动变化较小,表明该部分地区生态系统处于良性循环。不同植被类型的CUE均值表现为:草地(0.21) > 农田(0.14) > 森林(0.09) > 灌丛(0.06)。全国大多数地区植被CUE与降水呈正相关,而与气温则呈负相关,而这种相关性会随着区域气候格局及植被类型的变化而变化。总体上,全国植被CUE的增加主要归因于降水量的增加,而气温的升高则造成植被CUE的降低。  相似文献   

16.
文琦  丁彩霞 《水土保持通报》2013,(5):285-289,308
利用宁夏回族自治区盐池县1981-2010年气候、社会经济数据,采用倾向率法、一元线性回归方程拟合等,深入分析了盐池县气候变化与人类活动对生态环境的影响。结果表明,盐池县正处于气候暖干化时期,1981-2010年,气温以0.22℃/10a的速率上升;降水量以1.44mm/10a的速率在波动中下降;蒸发量以9.85mm/a的速率下降。不同时期人类活动对生态环境的影响不同,20世纪80,90年代的滥挖、滥牧、滥采造成该地区生态环境恶化。2001年起,退耕还林(草)、禁牧等对盐池县生态环境产生了积极的影响,扭转了20世纪90年代生态环境恶化的局面,并开始表现出一定抵御自然灾害的能力。研究认为,人类活动是盐池县近30a生态环境演变的主导因素。因此,只有合理开发、利用自然资源,推进人与生态环境和谐共处,才能保障盐池县经济社会的可持续发展。  相似文献   

17.
利用BIOME—BGC模型估算北京山区刺槐林的净初级生产力(NPP),并分析了NPP对不同未来气候情景的响应。结果表明:模型模拟出的NPP总体上低于样地实测值,年际变动率为11.2%,并表现出波浪形年际变化规律;降水量与NPP呈现极显著正相关关系,是NPP年际变化的主要气候因子;CO2的浓度加倍、降水增加对NPP积累起到促进作用,温度增加不利于NPP的积累;降水增加对NPP积累的正向效应不敌温度增加对NPP积累的负向作用,CO2浓度加倍、降水增加对NPP积累的协同促进作用大于温度增加对NPP积累的抑制作用。  相似文献   

18.
近47年天山山区自然植被净初级生产力对气候变化的响应   总被引:3,自引:0,他引:3  
利用天山山区10个对各类天然草场气候和植被类型具有较好代表性的气象台站1961—2007年的历史气候资料,使用线性趋势、Modet小波和Mann—Kendall检测等方法在对年平均气温、降水量变化趋势、变化特征分析的基础上,分别采用Miami模型、Thornthwaite Memorial模型、Chikugo模型、朱志辉模型和周广胜模型对该地区的自然植被净初级生产力(NPP)及其对气候变化的响应进行了计算和对比分析。结果表明:近47a,天山山区的气候总体呈较明显的“暖湿化”趋势,并于1970年发生了由“冷干”向“暖湿”的突变;尽管基于各模型的NPP在数值上有一定差异,但其在时空领域对气候变化响应的趋势是一致的,即:47a来,在气候“暖湿化”背景下,天山山区的NPP总体呈较明显的增大趋势,并于1970年发生了突变性的增大;与年平均气温和年降水量一样,基于各模型的NPP也具有3~22a的不同时间尺度的周期变化;在各类草场中,以草甸草原和山地草甸的NNP最高,高寒草甸次之,荒漠草原再次,高寒草原NPP最低。最后建立了基于各模型的NPP估算值之间的相互换算关系。  相似文献   

19.
夏冰  马鹏宇  徐聪  张磊 《水土保持研究》2023,30(2):256-266,284
为探究黄河流域植被净初级生产力(Net Primary Productivity, NPP)变化及对极端天气指数的响应情况,利用MODIS NPP数据和极端气候数据,辅以斜率法和偏相关分析法分析了2000—2019年黄河流域植被NPP时空动态及其对极端降水指数和极端温度指数的响应情况。结果表明:(1)近20 a黄河流域植被NPP呈从北向南增加的趋势,大面积上表现为增加趋势。(2) NPP与极端降水指标以显著和极显著正相关为主。其中,除最长连续湿润天数呈显著正相关和极显著正相关的像元占比为16.5%,其他几种指标均在25%以上。(3)极端降水事件的水量和强度呈显著增加趋势,极端温度事件中与偏冷相关指标总体呈下降趋势,与变暖有关的事件呈明显上升趋势,其变化存在显著空间异质性,年际差异大。(4)日最高气温的最大值、日最低气温的最大值、冷夜日数、冷昼日数、冰冻天数和霜冻天与植被NPP以负相关为主,日最高气温的最小值、日最低气温的最小值、暖夜日数、暖昼日数、气温日较差和暖期与植被NPP以正相关为主。近20 a黄河流域植被NPP时空变化存在显著地域性差异,其中在干旱和半干旱地区极端降水增多在一定程...  相似文献   

20.
为理解中国植被GPP的变化规律,基于改进后的EC-LUE模型,分析了中国1982—2016年植被GPP的空间格局、变化特征及驱动因子。结果表明:(1)中国植被GPP年变化范围介于0~3 051.08 g C/(m2·a),总体呈现出南高北低的空间分布特征。高值区主要分布在华南、东南沿海地区,低值区分布在青藏高原、新疆、内蒙古西部等地区;(2)近35 a中国植被GPP整体呈增加趋势,增加速率为2.43 g C/(m2·a)(p<0.05),特别是黄土高原地区植被GPP增加趋势明显,但在东北及东南沿海的部分地区呈显著减少的趋势;(3)中国植被GPP变化与气温和降水呈正相关关系,气温和降水强驱动的面积占比8.1%,气温为主要驱动因子的面积占比25.3%,降水为主要驱动因子的面积占比15.1%,气温和降水弱驱动面积占比8.9%。虽然近年来中国植被固碳能力呈变好趋势,但在未来气候变化存在不确定性的背景下,生态保护工作任重道远。  相似文献   

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

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