首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到19条相似文献,搜索用时 140 毫秒
1.
中国南北过渡带生态系统碳储量时空变化及动态模拟   总被引:3,自引:0,他引:3  
山地是全球变化的敏感地带,对生态安全与发展具有重要作用,山地生态系统服务变化和生态环境承载力是地理学与生态学的研究热点。以中国南北过渡带的主体秦巴山地为研究对象,采用CA-Markov模型与InVEST模型模拟和预测(2000—2040年)不同土地利用情景下秦巴山地生态系统碳储量变化,运用热点分析(Getis-Ord Gi*)探讨秦巴山地生态系统碳储量的空间分布差异。结果表明:(1)2000—2040年,研究区土地利用/土地覆被变化主要是耕地、林地、草地和建设用地。(2)2000—2020年,碳储量增加1.12×107 t;2020—2040年自然增长情景下,碳储量损失剧烈,减少50.24×107 t;生态保护情景下,碳损失幅度明显变弱,减少29.52×107 t,说明采取生态环境保护政策,能够有效控制碳储量减少。(3)土地利用/土地覆被与生态系统碳储量的变化呈现显著的一致性,土地利用数量变化决定了生态系统碳储量的质量和空间分布格局。(4)随着海拔抬升,碳储量呈现出“先增后减”的趋势;随着坡度升高,碳储量呈现出“W”型变化趋势。(5)热点分析结果显示,2000—2020年间,碳储量热点区和冷点区零散分布在研究区内;2040年自然增长情景下,碳储量冷热点分布范围有逐渐变大的趋势;2040年生态保护情景较2020年,秦巴山地生态系统碳储量的冷热点分布范围整体变化不大。  相似文献   

2.
科学划定“三线”(城镇开发边界、永久基本农田、生态保护红线)有利于合理安排国土空间资源,识别“三线冲突”强度影响因素,可加深对国土空间规划的认识。以省级空间规划试点区——江西省为例,通过“三线冲突强度指数”识别不同区域“三线冲突”差异,利用空间自相关分析探究空间特征,并应用地理加权回归模型分析不同冲突类型的影响因素。结果表明:研究区内,城镇—生态、城镇—农业、生态—农业、城镇—生态—农业冲突面积分别为4305.791 km2、2892.894 km2、1845.341 km2、159.444 km2,各类型冲突存在明显的空间聚集性;固定资产投资强度、建设用地适宜区面积占比和地形位指数分别为城镇—生态、城镇—农业和生态—农业冲突强度的主要影响因素,城镇化率和距河流距离对城镇—生态—农业冲突强度具有影响作用。研究结果可为今后合理布局“三线”、实现国土空间结构优化提供参考依据。  相似文献   

3.
文章在评述生态系统服务对全球变化的脆弱性研究概况的基础上,认为在生态系统服务退化和不可持续的状态下,改变土地利用变化对降低生态系统服务脆弱性具有重要意义。论文引入了生态系统服务对土地利用变化脆弱性评估方法,分析了生态系统服务对土地利用类型转变、土地利用类型渐变、土地管理方式变化的脆弱性响应机理,以江苏省环太湖地区碳储量为例进行实证研究。结果表明:①在1980—1990年、1990—2000年、2000—2005年3个时间段,农田、林地的显著减少和建设用地的急剧增加致使生态系统碳储量对土地利用强度变化表现得愈加脆弱,土地利用强度指数每增加一个单位,碳储量分别减少0.90 Tg、1.68 Tg、1.69 Tg;②2005—2015年间,通过对农田、林地、建设用地的规划控制,可以降低生态系统碳储量的脆弱性,土地利用强度指数每增加一个单位,碳储量分别由规划控制前减少0.50 Tg变为增加0.62Tg。据此得出结论,合理调整土地利用结构、加强规划管理可以降低生态系统碳储量对土地利用变化的脆弱性。最后,针对研究中存在的一些不足之处,提出了进一步研究的建议。  相似文献   

4.
基于土地利用变化的陕西省植被碳汇提质增效优先区识别   总被引:2,自引:0,他引:2  
在“双碳”目标背景下,陆地植被生态系统碳汇是实现碳中和目标的重要方式。为有效识别植被碳汇服务功能提质增效的优先区,利用InVEST模型定量评估陕西省植被碳储量时空演变特征及分布格局,分析土地利用/覆被类型变化对碳储量变化的影响,研究林草生态建设碳汇增长空间差异,确定林草生态建设提质增效对象区域。结果表明:(1)陕西省土地利用类型主要以耕地、林地、草地为主,土地利用类型转移变化也主要发生在三者之间;(2)1980― 2020年陕西省生态系统碳储量总体增加91.88×106 t,增幅3.16%,呈现出“总体上南高北低、局部地区明显过高或过低”的地带性分布特征;(3)退耕还林(草)工程对碳汇能力提升效果明显,存在全局空间相关性,表现为一定的空间趋同集聚现象;(4)陕北地区为生态保护修复工程极优先区和优先区,陕南地区为中等优先区,关中地区为一般优先区。研究基于不同区县生态系统碳汇年均增长率的差异,确定生态治理优先区域,可为实现生态修复工程主导模式的分区管理以及碳汇能力提质增效提供参考。  相似文献   

5.
陆地生态系统碳储量是表征生态系统碳储存服务的重要指标,与土地利用变化之间存在着密切的关系。退耕还林还草工程使区域土地利用格局发生巨大变化,并对生态系统碳储存服务造成了较大的影响。为了能简单快速的评估退耕还林还草工程对陆地生态系统碳储存服务所带来的影响,以位于黄土高原丘陵沟壑区的子长县为例,运用InVEST模型评估了退耕还林还草工程对陆地生态系统碳储量的影响,进一步耦合InVEST模型和FLUS模型,并设置四种不同的退耕还林还草实施情景,预测子长县2037年陆地生态系统碳储量变化和碳汇产生的经济价值。研究发现:(1)子长县退耕还林还草工程实施效果显著,17年间共有31627.98 hm2耕地退耕为林地和草地,境内的林草覆盖率由2000年的53.26%增长至2017年的64.20%;(2)退耕还林还草工程的实施显著提升了子长县陆地生态系统碳储存服务,碳储量由2000年的39.19×106 t增长至2017年42.34×106 t,增加量集中在工程实施主要阶段(2000—2008年);(3)未来子长县若继续实施退耕还林还草工程,其生态系统碳储存服务会得到进一步提升,且会获得一定的碳汇经济价值。预计到2037年子长县在退耕还林还草工程实施A、B、C、D四种情景下的陆地生态系统碳储量将分别达到:43.78×106 t、44.10×106 t、44.32×106 t和44.54×106 t,并将由此获得碳汇经济价值净收益分别为1627.88万美元、1979.89万美元、2231.39万美元和2471.67万美元。耦合InVEST-FLUS模型,不但能利用InVEST模型简单快速的对陆地生态系统碳储量进行评估,而且还能基于FLUS模型对未来土地利用变化情景下的陆地生态系统碳储量和碳汇经济价值做出测算。  相似文献   

6.
1985—2015年陕西黄土台塬表层土壤有机碳空间分布   总被引:1,自引:1,他引:0  
中尺度范围土壤有机碳的长期动态变化状态对大气温室气体的浓度、当地环境的生态效应至关重要。论文研究了陕西黄土台塬近30 a表层(0~20 cm)土壤有机碳密度的空间分布特征,并探讨了气候、地形、土壤类型、土地利用与土壤有机碳密度和储量的关系。结果表明:1)近30 a黄土台塬土壤有机碳密度和储量呈增加趋势,且增加幅度愈来愈明显,其中,1985—2000年有机碳密度和储量的增量分别为0.20 kg/m2和1.46×106 t,2000—2015年分别为0.75 kg/m2和10.87×106 t。2)1985—2015年有机碳密度随气温和降水量的增加而增加,随高程和坡度的增加呈现先增加后减少的趋势,在高程679~779 m、坡度10°~20°范围达到最大值。3)大多数土壤类型有机碳密度和储量随时间延长呈增长趋势,其中,黑垆土有机碳密度增加最大(1.59 kg/m2),黄绵土储量增量最多(5.64×106 t);不同土地利用类型有机碳密度和储量随时间延长也呈增加趋势,林草地有机碳密度增加量明显大于耕地。4)研究区表层土壤有机碳密度与气温、土壤类型、土地利用类型在P<0.001水平上相关性显著,与降水在P<0.05水平上显著相关。  相似文献   

7.
基于1990~2020年四平市土地利用类型数据,运用GeoSOS-FLUS模型,设定自然发展、耕地保护和生态优先三种情景,模拟不同情景下2030年四平市土地利用空间格局,同时结合InVEST模型定量分析研究区1990~2020年碳储量的时空分异特征,并探讨不同情境下土地利用变化对碳储量的影响,评估未来碳储量的潜力.结果表明:1990~2020年四平市耕地和林地分别减少了951.5 5km2和357.54km2,且以1990~2000年间的降幅最大.草地和建设用地呈增加趋势,分别增加了702.97km2和587.64km2.2030年在生态优先情景下,林地呈扩张态势,耕地有少量增加,在耕地保护情景下,耕地数量得到有效保障,而林地和草地有不同程度缩减.建设用地在三种情景下都呈现扩张的趋势,在自然发展情景下增长幅度最大.1990~2020年,四平市陆地生态系统的总碳储量及平均碳密度呈连续减少的势态,以1990~2000年的降幅最大,主要原因是该时段内土地利用变化以耕地的减少和建设用地的增加为主.四平市碳储...  相似文献   

8.
京津冀城市群土地利用变化对热环境的影响研究   总被引:5,自引:0,他引:5  
杨浩  王子羿  王婧  孟娜 《自然资源学报》2018,33(11):1912-1925
基于2000和2015年两期土地利用数据,运用WRF模式(Weather Research and Forecasting Model)模拟与数理统计方法,设计4组城市气温对比方案,定量模拟京津冀城市群土地利用时空格局变化对区域热环境的影响。结果表明: 1)2000—2015年城市群建设用地增加2 840 km2,94.37%由耕地转换而来; 2)北京土地利用变化对廊坊、天津和唐山气温变化的影响最为明显,天津的土地利用变化对廊坊气温变化的影响最为明显; 3)分季节来看,北京与天津的土地利用变化对城市群季均气温的影响排序为:秋季>夏季和春季>冬季; 4)其他地区土地利用变化使北京市年均气温增加0.07 ℃,且夏季降温作用明显,而该变化对天津市年均气温作用不明显,但秋季与冬季的增温区域形成了京津唐热环境连片发展的空间分布格局; 5)城市群土地利用变化对热环境影响表现出明显的“遥响应”特性,因而有必要从大尺度、泛区域的角度认识城市扩展对热环境的影响。  相似文献   

9.
海岸带国土空间开发适宜性评价及功能空间布局优化是推进陆海统筹战略目标的科学基础。以海南岛为例,基于环境损益分析思想构建海岸带“潜力—限制”评价指标体系,综合运用多要素乘积及空间叠置等方法,对海南岛海岸带地区的国土空间开发建设适宜性进行集成评价,并依据功能空间判别条件识别“港—工—城”区、滨海旅游区、农林牧渔区以及自然保护区四类海岸带功能空间。研究表明:海南岛海岸带国土空间开发适宜区总体呈现“西岸多,东岸少”的空间分布格局,适宜、基本适宜、基本不适宜、不适宜建设区面积依次为1785.32 km2、1724.57 km2、812.86 km2、4051.67 km2,适宜区与不适宜区重叠交错,给海岸带国土空间开发带来巨大压力。海岸带地区四类功能空间面积依次为1813.39 km2、2138.26 km2、1831.62 km2和2591.15 km2,分别占海岸带面积的21.65%、25.53%、21.87%和30.94%,海岸带的功能空间呈现“一环”“两极”和“多点”的空间格局。依据归纳出的三类典型功能空间冲突,对功能空间进行了调整优化,增强了各类功能空间的集聚性和连通性,但仍存在海域与陆域的相同功能空间不连片,海陆经济联系不够通畅,相互支撑不足等问题,国土空间开发还存在较多的陆海矛盾和冲突。研究结果指明了各类适宜性等级区域的开发或保护模式,丰富了适宜性的应用条件,也是对国土空间开发适宜性评价的适当拔高。  相似文献   

10.
基于1 km×1 km格网,论文综合运用景观指数、分散度、核密度、空间耦合关联等方法,研究河西走廊山丹县1998—2015年乡村聚落景观格局分异及时空演变特征,探讨绿洲乡村聚落与水土资源的空间耦合关系,结果表明:1)乡村聚落用地总面积扩张较快、空间分散化趋势明显,斑块数和斑块密度增大、平均斑块面积减小,乡村聚落景观指数在1998—2008年间相比2008—2015年间变化剧烈,后期逐渐趋于稳定;2)3个年份乡村聚落核密度空间分布格局基本一致,主干道路、渠系等线状地物沿线密集,外围区域稀疏,并以线状地物为中心向外呈阶梯状递减;3)1998—2008年乡村聚落分散度呈增加趋势,2008年后趋于稳定,乡村聚落空间分布由相对集中走向零星分散,聚落“散、空、乱”现象普遍;4)乡村聚落与耕地和渠系空间分布关系密切,存在显著的“耕地指向性”和“渠系指向性”;K指数1998年以“人多地少型”为主,2008和2015年转向“人少地多型”和“人地协调型”;L指数总体表现为“渠系偏多型”和“聚落-渠系协调型”。水土资源是约束和限制绿洲人口和聚落存在和发展的主要因素,人-水-土系统存在着高度的相互依存性和共生性。  相似文献   

11.
土地利用变化影响着浅层地下水水质,以滇池流域为研究对象,综合运用遥感影像解译、马尔科夫转移矩阵和冗余分析,对近20年(2002~2020)滇池流域土地利用和浅层地下水质变化进行分析,揭示长时间尺度下土地利用变化对浅层地下水水质的影响.结果表明:2002和2020年滇池流域土地利用类型以草地、林地和耕地为主,分别占总面积20.91%和17.43%、43.21%和37.99%、22.11%和17.08%,2002~2020年间耕地向建筑用地和林地、林地向草地和耕地、草地向林地和建筑用地转移概率分别为22.59%和20.72%、13.16%和10.49%、26.30%和15.65%,耕地面积减少146km2,建筑用地面积增加279km2.2002~2020年滇池流域浅层地下水化学类型由HCO3-·SO42--Mg2+型转变为HCO3-·SO42--Ca<...  相似文献   

12.
徐自为  张智杰 《环境科学研究》2018,31(11):1909-1917
为全面认识干旱区不同土地利用类型时空变化对区域生态系统碳储量的影响,以地处塔克拉玛干沙漠边缘生态脆弱区的新疆尉犁县为研究对象,基于详细的土地利用变更调查数据(2010-2016年),利用ArcGIS平台和InVEST模型,分析生态系统碳储量对土地利用变化尤其是二级土地利用类型变化的响应.结果表明:①研究区内,无论是区域平均碳密度还是灌木林地、其他林地、其他草地等主要土地利用类型的碳密度均较低,而面积较少的有林地和天然牧草地碳密度相对较高,因此对这些土地利用类型应着重加强保护.②2010-2016年新疆尉犁县碳存储量净减少24.23×104 t,这主要是由于其他草地、其他园地和果园等土地利用类型被开垦为水浇地,或被建设用地、交通用地占用所导致,而同时研究区内有林地和水浇地面积增加带来了碳储量的提高.③从空间变化看,受不同区域土地利用变化方式的影响,碳储量变化特征也有显著差异,但总体上变化敏感区域集中在塔里木河周边县、乡镇及兵团所在地等人类活动聚集区,这些区域平均碳密度较高,土地利用变化也更为剧烈.④尽管由于开垦行为带来耕地面积增加,从而使得研究区耕地总碳储量增加242.77×104 t,但由于塔里木河沿岸碳密度较高的耕地被建设用地占用,新增耕地多来源于土壤碳储量较低的其他草地等土地利用类型,导致耕地平均碳密度有所下降.研究显示,建设用地占用耕地、林地、草地等地类是尉犁县碳储量减少的重要原因,而林业建设能够带来碳储量增加.因此建议:一方面,推进林、草地建设提高区域生态系统固碳能力;另一方面,重点保护塔里木河周边碳密度较高区域,严控耕地开垦或建设用地占用,同时加强耕地保护,防止通过补充碳密度较低的耕地来弥补碳密度较高区域耕地的流失.   相似文献   

13.
Land use change is a main driver of carbon storage in terrestrial ecosystems. Based on land use data, research results related to carbon densities in vegetation and soil as well as government policies related to development in different regions of China, this paper optimized land use structure in 2020 for different regions with the goal of increasing terrestrial ecosystem carbon storage. We defined seven types of land use: (1) cultivated land, (2) garden land, (3) woodland, (4) pasture land, (5) other agricultural land, (6) urbanized land, and (7) a mixture of other land which we call mixed land which included open water, swamps, glaciers and other land as defined below. We found: (1) For most eastern regions, woodland has the highest carbon (C) densities while C densities of pasture land and cultivated land did not differ widely. Both have C densities higher than urbanized land while urbanized land has higher carbon densities than the areas placed in the mixed land type. (2) Under an optimized land use structure projected for 2020, the area of cultivated land will decrease compared with 2005 for most regions. The areas of garden land, pasture land and other agricultural land are much smaller compared with the mixed land use type, and the changes there are not obvious and their contributions to increased carbon storage are not significant. The area of woodland will increase the most obviously and it will contribute the most to increased carbon storage. The increasing urbanization of land and the decreasing trend of other land types make it difficult to change carbon storage patterns since the Chinese economy is expanding rapidly. (3) The optimized land use structure presented here will have effects on the entire country though with regional differences. Some inland regions will always have a larger potential to increase carbon storage than other areas because the potentialities in some coastal regions are limited by social and economic development.  相似文献   

14.
露天矿区开发活动主要是在景观层次上进行,其直接导致土地利用/覆被变化(LandUse/Cover Change,LUCC)。本文以黑岱沟大型露天煤矿1987年、1990年、2000年及2010年4期遥感影像的土地利用数据为基础,基于GIS空间分析功能,以破碎度、分离度和优势度景观格局指数和景观类型脆弱度为指标,分析景观格局现状。研究显示:开发初期,耕地和草地虽有所减少,但所占面积保持在80%以上;生态恢复开始后,土地类型向林地和草地转移的面积逐渐增多。景观水平上分析显示,景观上破碎化程度加大,煤矿建设用地斑块面积逐年增大。从景观指数看,耕地、草地、未利用地的破碎度、分离度指数变化呈上升或不稳定状态,建设用地的破碎度和分离度指数有所减小;优势景观类型有所改变,建设用地对于区域景观的控制作用逐渐加强。  相似文献   

15.
通过对太湖流域重污染区1999年和2007年的2期遥感影像数据解译结果的分析,获取了太湖流域重污染区主要土地利用类型的信息,分析了8年来研究区内土地利用与覆被变化趋势,并对区域内土地利用/覆被变化的驱动因子进行了探讨。结果表明,从1999~2007年间,太湖流域重污染区内耕地面积减少,而林地、建设用地、园地、水域总体呈增加的趋势,其中减少的耕地主要转化为建设用地。在区域结构中,影响土地利用结构变化的驱动力因素主要是经济的发展、城市化进程推进、人口增长以及政策因素这4个方面,而研究区内经济的快速发展和人口的增长成为土地利用/覆被变化的关键因素。  相似文献   

16.
科学地进行城镇发展空间与农业生产空间的划分,是解决城市无序发展问题和粮食安全问题的重要途径。本文以江西省临川区为例,借助ArcGIS、GeoDa等软件,对临川区耕地综合质量及其自相关关系进行评价,以此确定永久基本农田初步划定成果,采用多准则判断CA模型对城市扩张进行模拟,通过两方面结果的对比与调整,结合城市发展方向、发展模式和耕地国家利用等别,进行永久基本农田和城市开发边界的划定。结果表明:永久基本农田划定面积为593.46 km2,初步划定结果主要分布于上顿渡镇、崇岗镇、河埠乡和荣山镇等耕地综合质量较高,且多为高—高聚集的区域。经过城市扩张模拟结果可知,临川区城市建设用地面积较为显著增加,主要增加面积分布于展坪乡、上顿渡镇以及钟岭街道。城市发展方向和城市扩张模型分析结果显示,临川区城市建设用地主要向西北、东南、西南方向扩张,扩张类型主要为蔓延式扩展。根据耕地国家利用等别将初步划定永久基本农田区分为可调入区与可调出区,从而进行永久基本农田区与城市开发边界的优化与协调,形成永久基本农田区与城市开发边界的最终成果。研究结果可为保护耕地资源、避免城市的无序扩张提供参考,也可为今后的划定工作提供借鉴。  相似文献   

17.
为了有效评估城市群碳储量变化,以天山北坡城市群为研究对象,运用PLUS模型和InVEST模型,动态评估2000~2020年及2030年不同情景下土地利用变化及碳储量变化特征.结果表明,2000~2020年天山北坡城市群碳储量呈现持续增加趋势,且碳储量变化与土地利用变化密切相关,主要表现为2000~2010年林地面积的减少导致其碳储量减少约266×106t,2010~2020年草地面积的增加使其碳储量增加约69.14×106t.2030年自然发展情景、生态保护情景和经济快速发展情景下碳储量预测值分别为8875.88×106t、8895.58×106t和8841.58×106t;经济快速发展情景下碳储量最低,生态保护情景下碳储量最高.土地利用是影响碳储量空间变化分布的第一主导因素,贡献率接近于90%,土地利用强度与碳储量协调性分析与两者双变量空间自相关分析进一步验证了这一结论.土地利用变化在一定程度上能够对碳储量产生积极影响,对于本研究区而言,生态保护发展情景可能更符合未来城市发展模式,研究结果能够为土地利用规划提供参考.  相似文献   

18.
大清河阜平流域生态系统结构变化的服务价值响应研究   总被引:4,自引:0,他引:4  
生态系统通过物质循环、能量流动和信息传递为人类提供必要的生态服务功能,其结构变化可引起生态系统过程和服务功能的改变。研究区域的生态服务功能对生态系统结构变化的响应规律对高效配置环境资源和合理制定环境政策具有重要意义。论文以大清河阜平流域为研究对象,利用1985—2015年每隔5 a的遥感影像,采用监督分类法及价值量法,研究了该区域生态系统结构、服务功能近30 a的变化规律趋势,分析了研究区服务功能对生态系统结构变化的响应规律。结果表明:1)1985—2015年研究区生态系统结构发生明显变化,建设用地、林地、耕地、水域生态系统呈正变化,变化幅度分别为139.57%、90.66%、68.79%、6.55%;未利用地、草地生态系统呈负变化,变化幅度分别为-83.48%、-54.92%;其中林地面积的增加主要是由草地和未利用地转变而来,耕地面积增加主要由未利用地构成,建设用地的扩张主要占用耕地和未利用地。2)研究区服务价值变化趋势以2000年为界,前15 a呈增加趋势,后15 a呈减少趋势。2000年总服务价值最高,为58.08亿元;不同时期林地服务价值对研究区总服务价值贡献率最大,其次是耕地和草地;单项服务价值中,保持土壤服务价值贡献率最大,其次是维持生物多样性、气体调节和气候调节服务价值。3)研究区服务价值的增加得益于草地→林地的服务价值流向和未利用地→林地的价值流向,服务价值的损失来自于耕地转变为建设用地的服务价值损失和林地转变为耕地的服务价值损失。  相似文献   

19.
耦合InVEST与GeoSOS-FLUS模型的桂林市碳储量可持续发展研究   总被引:3,自引:0,他引:3  
为了量化桂林市碳储量并快速评估分级保护措施对区域碳储功能的影响,耦合InVEST模型碳储存模块和GeoSOS-FLUS模型,并基于土地利用数据和不同情景未来土地预测结果,对2000~2040年桂林市域范围六区十一县市内的碳储时空特征进行分析.结果表明:桂林市2000年、2010年和2020年的总碳储量分别为554.02×106t,553.58×106t,550.21×106t,呈现“逐年下降”的变化态势.同时,受人类活动和土地利用类型变化的影响,桂林市域各区县的碳储水平存在较大的时空差异,碳储量整体表现为“西北、西南及东部较高,东北、东南及中部较低”的空间分布特征.将桂林市碳储量高值区确定为碳储资源的优先保护区域,与自然变化情景相比,资源保护情景下桂林市林地得到有效保护,建设用地规模扩大受到限制.采取资源保护措施后,桂林市2040年总碳储量达到552.16×106t,较2020年增加了1.95×106t,中低密度碳储区所占比例明显下降,区域固碳能力大大增强.该研究结果可为桂林市国家可持续发展示范城市建设提供指导,也可为碳储资源精准保护和土地利用管理决策提供科学参考.  相似文献   

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

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