基于InVEST和CA-Markov模型的黄河流域碳储量时空变化研究.docxVIP

基于InVEST和CA-Markov模型的黄河流域碳储量时空变化研究 摘要：区域土地利用/覆被变化是导致生态系统碳储量变化的主要原因,预测未来土地利用/覆盖变化及其对碳储量的影响对区域陆地生态系统的认识具有重要意义。本研究基于黄河流域2005—2018年土地利用/覆被变化规律,运用CA-Markov模型分别预测了生态保护情景(EVC)和自然变化情景(NVC)下的土地利用/覆被空间格局,采用修正后的碳密度,运用InVEST模型评估黄河流域2005—2030年6期碳储量。结果表明:2005—2018年黄河流域林地、水域和建设用地面积持续增加,耕地、草地和未利用土地面积减少, 13 a间全流域碳储量减少28.734×10 陆地生态系统的碳储量在全球碳循环和气候变化中起着非常重要的作用黄河流域是我国重要的生态屏障,是连接青藏高原、黄土高原、华北平原的生态廊道,拥有多个国家公园和国家重点生态功能区,流域内林草面积占流域总面积的50%以上,是我国陆地生态系统重要的碳汇和碳储存区域。伴随黄河流域人口增加和城镇化进程,国土开发利用活动强烈而且频繁,为了保护流域生态环境,国家实施了退耕还林还草等一系列生态保护工程,流域内土地利用覆被发生了显著变化。王光谦等1 研究区概况与研究方法1.1 研究区概况黄河是我国第二大河流,位于96°～119°E、32°～42°N,发源于青藏高原巴颜喀拉山,在山东省垦利县注入渤海,流经青海、四川、甘肃、宁夏、内蒙古、山西、陕西、河南、山东9省区,干流全长5464 km,东西长约1900 km,南北宽约1100 km,落差4480 m,流域面积7.95×101.2 研究方法1.2.1 In VEST模型碳储量模块In VEST模型(Integrated Valuation of Ecosystem Services and Trade-offs)即“生态系统服务功能与权衡交易综合评价模型”,In VEST模型碳储量模块将生态系统的碳储量划分为4个基本碳库:地上生物碳(土壤以上所有存活的植物材料中的碳)、地下生物碳(存在于植物活根系统中的碳)、土壤碳(分布在有机土壤和矿质土壤中的有机碳)、死亡有机碳(凋落物、倒立或站立的己死亡树木中的碳)。根据土地利用的分类情况,分别对不同地类地上碳库(C基于各地类的碳密度和土地利用数据,流域内每种土地利用类型的碳储量计算公式为:式中:i为每种土地利用的平均碳密度,A流域内总碳储量为所有地类碳储量之和。1.2.2 碳密度及其修正碳密度是In VEST模型准确评估碳储量非常重要的输入参数。本研究不同地类碳密度数据来源于国家生态科学数据中心(/)并参考了部分文献由于本研究的碳密度数据来源于全国和一些地方的研究结果,并不是实际的测量结果,而碳密度值随气候、土壤性质和土地利用的不同而不同式中:C式中:K1.2.3 CA-Markov模型及情景设定1.2.3. 1 CA-Markov模型原理CA-Markov模型将CA和Markov模型结合起来,利用一个转移概率矩阵来模拟土地利用随时间的变化。Markov链是以马尔可夫随机过程为理论基础的预测方法式中:S但在马尔科夫链过程中没有充分考虑空间参数,无法识别土地利用的空间变异性式中:S为元胞有限、离散状态的集合,f为元胞状态的转换规则函数,N为每个元胞的邻域,t、t+1为两个不同的时刻。CA-Markov模型结合Markov模型长时间序列模拟预测的优势和CA模型在模拟空间变化的优势,可较好地预测和模拟LUCC在数量和空间的时空格局1.2.3. 2 CA-Markov模型预测情景设置及过程在本研究中,使用IDRISI 17.2软件运行CA-Markov模型,进而预测2030年土地利用格局,具体经过以下操作步骤:1)首先IDRISI 17.2软件中的Markov模块以2005年和2018年的土地利用图为初始年份和最后年份,得到转移概率矩阵和转移面积。其中转移概率矩阵反映了各土地利用转化为其他类型的概率,转移面积反映了下一阶段区域内各土地利用转化为其他土地利用的面积。2)利用多标准评价模块(Multi-criteria Evaluation,MCE)建立土地利用变化适宜性规则图集。以限制性因素和驱动因素作为自变量进行分析,得到每种地类的空间分布概率图。本研究根据未来对黄河流域生态环境的保护、生态工程的实施以及土地利用的管控方向,设定自然变化和生态保护两种可能的土地利用覆被变化情景,并分别设置约束条件。情景1:生态保护情景(EVC)。模拟到2030年,根据生态保护的要求,对林地、草地以及水域转出实施严格的控制,林地和草地不能转为其他地类,但其他地类可以转为