关中-天水经济区生态系统固碳服务空间流动及格局优化

李婷，李晶,3，王彦泽，曾莉



关中-天水经济区生态系统固碳服务空间流动及格局优化

李婷1,2，李晶1,2,3，王彦泽1,2，曾莉1,2

（1陕西师范大学地理科学与旅游学院，西安710119；2地理学国家级实验教学示范中心（陕西师范大学），西安710119；3宝鸡文理学院/陕西省灾害监测与机理模拟重点实验室，陕西宝鸡721013）

量化区域生态系统固碳服务（以下简称“固碳服务”）的供需平衡状况，模拟固碳服务的空间流动，揭示区域固碳服务功能转变的空间规律，并据此给出区域固碳空间布局优化策略，为引导区域的低碳发展提供直观的科学参考。通过多源数据融合模型模拟人口空间密度，进而估算关中-天水经济区生态系统固碳服务需求量；使用CASA模型计算研究区生物碳库的固碳服务供给量，在此基础上结合遥感碳循环过程模型估算土壤碳库的固碳服务供给量；使用流动比率量化区域固碳服务的供需平衡状况，以空间可视化的方式梳理固碳服务从供给区到受益区的流动过程；基于贝叶斯原理计算环境变量对固碳服务供给的状态条件概率，将其可视化表达为条件概率像素图，通过熵减模型筛选关键变量，利用关键因子最优状态子集分布探讨固碳格局的不确定性，并给出固碳空间布局优化策略。（1）研究区固碳服务总体上供大于求，供需平衡空间差异明显，需求的高值区主要分布在关中平原的人口高度聚集区，供给的高值区则主要沿秦岭山脉和北山山系分布。（2）根据流动比率分布，研究区大致可分为三大碳源集聚中心：以天水市区为中心的碳源集聚中心（R>0.04），以彬县为中心的碳源集聚中心（R>0.04），以及以最高值西安市为中心的多级碳源集聚中心（R>0.20）。其空间流动则可分为四组供需对应区域：秦岭中段、东段及北山的永寿县等流向以西安市为中心的关中城市群，秦岭西段流向天水市，麟游县、旬邑县等流向彬县，铜川市、澄城县、华县等流向蒲城县。（3）根据条件概率和熵减度计算，研究区可取{DEM=3, PET=1}作为生物碳库的关键变量最优状态子集，该子集主要分布在宝鸡市南部秦岭山脉和天水市西南角，该区域中生物固碳状态为最优的概率可以达到54.36%；取{NPP=3, DEM=3}作为土壤碳库的关键变量最优状态子集，该子集主要分布在秦岭山脉沿线、天水市西南角和咸阳市东北角，这些区域中土壤固碳状态为最优的概率高达92.84%。与现有固碳格局对比得出，生物碳库的适宜优化区主要分布在天水市武山县和秦州区，土壤碳库的适宜优化区主要分布在秦岭中段各区县。研究区固碳服务总体上需求小于供给，区域内存在较明显的固碳服务空间流动，未来的固碳格局优化中以天水和秦岭中段作为固碳功能优化的主要区域，能够得到良好固碳格局优化成效的概率较大。

固碳；生态系统服务；空间流动；贝叶斯；关中-天水经济区

0 引言

【研究意义】生态系统固碳服务（以下简称“固碳服务”）通过捕获大气中的碳并把捕获到的碳固定起来，它可以抵消人类向大气中排放的一部分二氧化碳，从而起到调节气候的作用，是生态系统服务中调节服务的重要组成部分[1]。【前人研究进展】关于生态系统服务功能的受益人群、生态补偿等研究一直都广受研究者关注[2-3]，而生态服务供给与人类需求的相互关系及其空间平衡问题近年来逐渐引起国内外专家的重视[4-8]。服务供给与需求的空间匹配问题是梳理生态系统服务空间流动过程研究的主要障碍之一[9-10]。目前得到较多学者认同的空间匹配模型是服务路径属性网络（SPANs）模型[11-13]，但其研究主要还处于概念阶段[14-15]。从空间流动的角度来看，固碳服务供需的研究主要可分为典型生态系统与大气之间的固碳服务供给（以下简称“碳供给”）研究[16-17]和区域之间的固碳服务需求（以下简称“碳需求”）研究[18]。这些研究促进了人们对于碳供需的理解，但对于以人类活动和自然环境共同组成的复杂系统而言，单纯的生态碳循环研究或经济学碳供需研究不能全面反映区域中固碳服务的供需平衡[19-20]。近年来该研究领域逐渐突破单一视角的局限，开始从人类福祉的角度出发，开展区域固碳服务供需平衡研究[21]，一定程度上填补了区域碳平衡研究的空白，但对碳管理实践而言仍有一定局限性。【本研究切入点】作为影响地表植被覆盖和人们能源消费水平的重要因素，科学合理的土地利用格局配置无疑是实现有效的区域碳管理、引导区域低碳发展的重要工具[22]。国内外学者已经在土地利用的碳排放效应、城市化进程中的碳收支等领域做了大量的实践研究，但多侧重于土地类型面积的线性组合和土地利用预测，对区域的空间结构优化研究仍需继续探索。【拟解决的关键问题】本文以关中-天水经济区为研究区域，估算研究区碳需求与供给，量化固碳服务供需平衡状况，模拟区域内固碳服务的流动状况，分析环境变量对碳供给的影响，探讨区域中碳供给的不确定性，并给出空间布局优化策略，促进对固碳服务空间流动的客观认识，为区域碳管理提供科学直观的理论依据和数据支撑。

1 材料与方法

1.1 研究区域

关中-天水经济区（简称“关天经济区”）（图1）的经度范围是104°34′47″E—110°48′38″E，纬度范围是33°21′37″N—35°51′15″N，面积约8.01×104km2，年均气温6—13℃，年降水量500—800 mm，海拔200—3 700 m。根据地形地貌和气候差异可将研究区分为三大部分：秦岭山脉、关中平原和北山山系。

图1 研究区行政区划图

作为陕西省重点生态保护区和南水北调重要水源涵养基地，秦岭的大部分区域都属于限制开发区，生态基础优异，是关中地区的主要林地和天然碳库。其土壤类型以棕壤、黄棕壤、褐土为主，其间散布山地棕壤、暗棕壤、粗骨土和土等，主要植被类型包括常绿阔叶林、常绿针叶林、落叶阔叶林、落叶针叶林等。北山山系耸立于关中平原与黄土高原之间，主要由中生界和古生界砂页岩、砾岩组成，土壤类型以黄墡土、褐土为主，其间零星分布小面积红壤、粗骨土和棕壤，具有较高的植被覆盖率。渭河自西向东流经研究区，关中平原则是自宝鸡峡沿渭河向东延展至潼关的冲积平原，地势平坦开阔、土壤肥沃、水源充足，适宜农业发展，被称为“八百里秦川”。总的来说，秦岭山脉和北山山系的生态基础较好，而关中平原的人类社会经济活动密集，生态环境则相对较差。近年来，研究区经济发展迅速，西咸新区建设如火如荼，宝鸡等次核心城市群和产业集聚带已经初具规模。然而，日新月异的经济发展使得城市建设与生态环境建设之间的矛盾不可避免。因此，关天经济区是开展固碳服务平衡及空间流动、固碳格局优化等相关研究的典型区域。

1.2 数据来源

研究选用的数据源主要包括遥感数据、地理基础数据、气象数据、土地利用专题地图和Harmonized World Soil Database（HWSD）。其中，遥感数据为2014年250 m空间分辨率的MODIS数据产品。地理基础数据包括DEM、坡度、坡向、土壤类型等。气象数据来源于中国气象科学数据共享服务网，包括年平均降水、年平均气温、7月平均气温、地表太阳辐射等。HWSD为本文提供中国1∶100万土壤有机碳含量分布图。土壤有机碳实测数据为2015年和2016年8月分别于研究区内野外采样并通过重铬酸钾法测定所得。

1.3 固碳服务供需模型

1.3.1 碳需求模型 碳排放的测算方法主要有实测法、物料衡算法和排放系数法，其中排放系数法（Emission-Factor Approach）是IPCC提出的、适用于多尺度的碳排放估算方法，但单纯的碳排放测算难以体现碳排放的空间差异。为了更加直观地表现区域固碳服务需求的空间分布，本文基于能源统计中人均碳排放数据和人口密度分布计算区域碳源的空间分布。

式中，C为人类社会经济活动中的碳排放，即固碳服务需求；()为像元的人口空间密度；()为像元的人均碳排放;表示研究区的像元总数。

应用较广的人口数据空间化的主要模型有负指数模型及其改进模型、核密度估计模型、分区密度模型、多元回归模型、多因素融合模型和智能化模型[23]。为了综合考虑自然、经济因素对人口分布的影响程度，本文选用多因素融合模型，基于遥感影像解译结果和居民点分布，将区域分为城镇区域和农村区域；基于多种基础数据分别计算城镇人口密度系数和农村人口密度系数；结合市县人口统计数据和人口密度系数计算出区域人口密度分布状况。基于城镇面积的城镇人口密度系数模型为[24]：

式中，V为第个城市中第个栅格的人口密度系数，A为栅格中的城市用地面积，A为第个城市的面积，r为第个栅格中心至城市中心的距离，参数反映城市的不同发展阶段。而农村人口密度系数则是根据调研所获得的数据和公式（3）求D与D比值后，经单要素加权融合、多元逐步回归分析、归一化处理等后得出[25]。

式中，P、P、A、D、A、D分别表示各县总人口、城镇人口、农村居民点的总面积、农村居民点平均人口密度、远离居民点地区的总面积和平均人口密度。

1.3.2 碳供给模型 碳供给部分包括生物固碳和土壤固碳两部分。考虑到模型复杂度、估算精度、所需数据的可获得性[26]及研究区域实际情况，本研究选用CASA模型计算生物部分的碳库存量，选用由周涛等[27]改进的碳循环过程模型来反演土壤基础呼吸进而估算土壤部分的碳库存量。

CASA模型[28]通过光合有效辐射（）和实际光能利用率()的乘积来表示每个月的植被初级生产力的变化，计算公式如下：

式中，为植物净初级生产力，()表示像元在月份吸收的光合有效辐射（MJ·m-2），()表示像元在月份的实际光能利用率（gc·MJ-1）。根据不同植被地上生物量与地下生物量的经验比例值，最终计算区域生物碳库的碳供给量。

通过碳循环过程模型来反演土壤基础呼吸，建立土壤中CO2排放与土壤有机碳（SOC）之间的关系模型。该模型在原碳循环过程的基础上增加了水分因子，用年降水量和年潜在蒸散量结合温度敏感因子描述土壤水分对土壤呼吸的影响。其公式如下：

式中，A为土壤基础呼吸，为净初级生产力，为年平均温度，为温度敏感常数因子，为水分对土壤呼吸的限制因子。

由于HWSD中的土壤数据来自全国第二次土壤普查，与现有的实际数据可能存在一定的偏差。因此本文先计算了2000年的土壤基础呼吸，将其与HWSD土壤有机碳数据做回归分析，得出土壤基础呼吸与土壤有机碳含量的回归模型；将2014年的土壤基础呼吸代入回归模型，得出2014年的土壤有机碳含量。使用野外调研实测数据验证模拟精度，标准误差为4.9726，研究结果具有较好可信度。

1.4 固碳服务的空间流动模型

从生态学的角度，生态学家们关于“流”的概念最早可以追溯到将生命的代谢过程看作是能量、物质与周围环境不断发生交换的过程，即能量流和物质流[29-30]。生态系统中的“能量流”通常用“林德曼定律”来表示[31]，但关于碳的空间流动目前尚没有成熟的理论和方法。“场”在物理学中是指物体在空间中的分布情况，可以用来描述某种物理量在空间的分布和变化规律，用空间位置函数来表征。若物理量是标量，那么空间上每一点都对应着该物理量的一个确定数值，则称此空间为标量场（如电势场、温度场）；若物理量是矢量，那么空间每一点都存在着它的大小和方向，则此空间为矢量场。物理学中常见的有电场、速度场等，地理学中常见的有风场等。受此启发，固碳服务的空间分布也可以理解成一个场：由于环境变量空间分布的差异，区域的固碳服务势能也有所差异，在这一势能的驱使下，区域内部以及区域与环境之间会产生固碳服务在空间上的流动。

区域的碳源与碳汇受到来自自然生态和社会经济等多方面的影响，二者之间的平衡数量关系也比较复杂。在此引入流动比率（R）的概念，通过计算碳源和碳汇的比值，以表征同一区域中碳源与碳汇平衡的数量关系[19]。其概念模型为：

式中，R为像元的源汇流动比率，O为像元某时段的碳排放量，即通过碳源模型计算得到的碳排放量；I为同一时段像元的碳吸收量，即生物固碳量。R能够表达区域固碳服务的自给率，R越小，区域的固碳服务自给率越强，反之则越弱。当R<1，表明区域碳供给倾向大于碳需求倾向，即为供给型区域；反之，当R＞1，表明区域碳需求倾向大于碳供给倾向，即需求型区域。运用地理信息系统手段，对固碳服务的空间格局及其流动比率做空间表面分析，可绘制区域固碳服务空间流动的方向和趋势。

1.5 条件概率与因子状态可视化

生态系统固碳过程涉及大量环境因子，其格局优化需以各因子对固碳功能的影响程度为前提。本文基于贝叶斯信念网络（Bayesian Belief Networks，BBN），以概率统计为理论基础，借助python编程，模拟和计算各固碳因子的不确定性，并将其分布格局空间可视化（图2）。

首先，参考自然分级结果，综合考虑研究区实际情况与坡度等相关等级划分标准，将备选因子的栅格数据按照离散化标准（表1）做离散化处理。考虑到生物碳库和土壤碳库的生物循环和代谢过程存在较大区别，因此本文中二者的备选因子也有所差异。其中，生物碳库的备选因子有高程、7月平均气温、潜在蒸散发、降水、坡度、土壤类型、太阳辐射，土壤碳库的备选因子则包括高程、净初级生产力、降水、坡度、土壤类型、太阳辐射、年平均气温。

其次，基于贝叶斯网络原理，计算碳库中各变量状态的条件概率。设有栅格变量，记作A={A，…，A}，其中A为栅格变量A的不同状态。各变量不同状态的栅格数量记为SUM，研究区栅格总数记为，那么栅格变量A的不同状态概率P为:

又有事件B={B1，…，B}，事件A和B同时发生记作AB，满足此条件的栅格总数记为SUM 。那么事件A和B同时发生的联合概率（A B）为：

那么事件B发生的条件下，事件A发生的条件概率（A|B）为：

通过Python编程统计各因子对应的栅格数，计算备选变量两两配对的条件概率，生成条件概率表，并利用Python将其可视化为像素图。根据像素图，选取关键变量状态子集，并将关键变量状态组合空间可视化。

然后，筛选关键变量最优状态子集。其中，最优状态是指当固碳状态等级达到最优（状态等级=3）时所对应的最高概率备选变量状态，即最适宜于固碳的状态，可通过分析条件概率像素图得出；关键变量是指备选因子中对固碳结果影响最大的几个变量，本文根据实际情况将关键变量数量设置为2个，可通过计算熵减度筛选出。

表1 碳库备选因子及其数据状态离散化标准

图2 生态系统固碳服务功能关键变量筛选及其可视化

通过熵差法计算每个变量对固碳结果的熵减度：

其中，为生物（或土壤）固碳节点，为生物（土壤）碳库的备选因子；E()为生物（或土壤）固碳的熵，E()为生物（或土壤）固碳与备选因子联合发生的熵；()为生物（或土壤）固碳状态发生的概率，()为备选因子状态发生的概率，()为其联合概率。为熵减度，取值范围在0—1，越大则该因子对固碳结果的影响越大，反之越小则该因子对固碳结果的影响越小。

根据上述方法，可以挑选出关键变量的最优状态子集，在关键变量的最优状态能够同时发生的区域，获得高固碳效益的概率最高。

2 结果

2.1 固碳服务供给与需求

土地利用方式是人类社会活动作用于自然环境和气候变化的重要介质，不同的土地利用方式具有不同的生物碳固定能力和土壤碳固定能力。在关天经济区，草地、林地和耕地的生物固碳总量分别占全区的41.90%、28.54%和28.18%，其土壤固碳总量分别占全区的27.59%、36.41%和34.73%，是研究区的主要碳汇区域。而城市用地的单位面积固碳能力为269.47 g·m-2·a-1，单位面积平均碳排放能力却高达128.59 g·m-2·a-1，总碳排放占全区总量的49.35%，是研究区最主要的碳源区域。耕地的单位面积碳排放能力虽不强（4.81 g·m-2·a-1），但因面积广大而导致其总碳排放能力占全区的33.94%，成为研究区的主要碳源之一。由此可见，土地利用格局对于区域碳源汇平衡具有重要意义。

研究区固碳服务在整体上供大于求，但各区县的供需状况差异明显（图3）。分布在秦岭山脉上的凤县、太白县等，依靠秦岭良好的生态基础，拥有较高的固碳服务供给量（1.10×1010—1.33×1010gC），而其经济发展水平有限，所需要的固碳服务量相对较低（0.41×105—0.89×105gC）。天水市麦积区、周至县、宝鸡市陈仓区和洛南县虽然固碳服务需求量较高（8.23×105—10.80×105gC），但由于供给量（0.90×1010—1.42×1010gC）远远大于需求，因此固碳平衡仍为供大于求。而西安多个市区、咸阳市各区、宝鸡市金台区等的固碳服务供给量（≤5.57×108gC）明显低于其他市区，而需求量却较高（≥6.84×105gC）。

1：白水县；2：宝鸡市陈仓区；3：宝鸡市金台区；4：宝鸡市渭滨区；5：彬县；6：长武县；7：澄城县；8：淳化县；9：大荔县；10：丹凤县；11：凤县；12：凤翔县；13：扶风县；14：富平县；15：甘谷县；16：高陵县；17：韩城市；18：合阳县；19：华县；20：华阴市；21：户县；22：泾阳县；23：蓝田县；24：麟游县；25：礼泉县；26：陇县；27：洛南县；28：眉县；29：蒲城县；30：乾县；31：千阳县；32：秦安县；33：清水县；34：岐山县；35：三原县；36：商洛市商州区；37：太白县；38：天水市麦积区；39：天水市秦州区；40：铜川市王益区；41：铜川市耀州区；42：铜川市印台区；43：潼关县；44：渭南市临渭区；45：武功县；46：武山县；47：西安市灞桥区；48：西安市碑林区；49：西安市长安区；50：西安市莲湖区；51：西安市临潼区；52：西安市未央区；53：西安市新城区；54：西安市阎良区；55：西安市雁塔区；56：咸阳市秦都区；57：咸阳市渭城区；58：兴平市；59：旬邑县；60：杨陵区；61：宜君县；62：永寿县；63：张家川自治州；64：柞水县；65：周至县

2.2 固碳服务供需关系的空间格局

通过固碳服务流动比率（R）空间分布，可以更直观地了解区域的碳平衡状况（图4）。关中平原各市县的城区普遍表现为R≥1，分布特征为R以西安市、咸阳市和渭南市为高值中心向四周递减；而研究区大部分区域R<1，属于供给型区域，秦岭山脉和北山山系的R则明显偏低，部分森林核心区域R几乎为零。根据研究区固碳服务供需的格局分区（图5）可知，研究区可分为三个主要的碳源集聚中心，分别是以天水市（R>0.06）和彬县（R>0.06）为中心的两个碳源集聚中心以及以西安市（R>0.20）为中心的多级碳源集聚中心。

其中，西安市（碳源中心A，R>0.20）为一级碳源聚集区，杨陵区（碳源中心B，R>0.14）、咸阳市和高陵县（碳源中心C，R>0.10）为二级碳源聚集区，渭南市、武功县等（碳源中心D，R>0.04）为三级碳源聚集区。研究区各县级以上行政区的碳源汇流动比率总和分布情况主要表现为以关中平原为中心轴，南部流动比率的变化率略大于北部，尤其是眉县、周至县、武功县和杨陵区之间，流动比率呈现出骤降趋势（图6）。供给特征最强的区域（∑R≤200）主要为秦岭山脉和北山山系的各县，供给特征较明显的区县（400<∑R≤800）主要分布在研究区西侧的宝鸡市和天水市，而关中平原则分布着供给特征较弱的区县（∑R>800）。

固碳服务空间分布的不均衡致使需求型区域与供给型区域的流动比率相差上百倍，根据物质循环原理，需求型区域所消耗的碳将由其周边的供给型区域提供，这在一定程度上会引起固碳服务由供给型区到需求区的空间流动。研究区的固碳服务空间流动整体上表现为：在秦岭中段和东段，大量的固碳服务向以西安市为中心的关中城市群涌进，关中平原北侧也有固碳服务流入；在秦岭西段，固碳服务则主要流向天水市，且流量较小；北山的彬县呈现出一定的弱供给特征，麟游县、旬邑县等有小支固碳服务向彬县流入补给；东北侧蒲城县表现出较明显的弱供给特征，其固碳服务的补给源则是铜川市、澄城县、华县等（图7）。

图4 研究区生态系统固碳服务流动比率空间分布

2.3 固碳的条件概率与格局分析

2.3.1 条件状态子集筛选 将生物碳库和土壤碳库的固碳能力与其备选因子的条件概率可视化表达为像素图（图8、图9），图中将每个因子离散化为3个等级，横纵坐标是每个因子的不同状态，纵坐标是事件A，横坐标是事件B，两两状态所对应的色块表示B条件下A发生的条件概率。颜色由蓝到红表示概率从0到1。图中最右一列九宫格表示当某一地理位置的固碳功能处于不同状态时，对应位置的各种变量状态的组合情况。最下方一行九宫格则可以表示当某一地理位置的变量在不同状态时，对应位置的固碳功能的条件概率；其中当固碳能力状态为3时，其所对应的高概率备选变量状态集合为最优固碳状态的高概率条件状态子集（简称“最优状态子集”）。

令生物固碳功能处于高中低三种状态时，其对应的备选变量状态组合称为“备选变量状态子集Ⅰ、Ⅱ、Ⅲ”。那么，状态子集Ⅰ= {DEM=1, T_July=2, PET=2, Prec=2, Slope=1, Soil=2, Sol=2}。状态子集Ⅱ= {DEM=2, T_July=2, PET=2, Prec=2, Slope=1, Soil=2,Sol=2}。状态子集Ⅲ= {DEM=3, T_July=2, PET=2, Slope=2, Soil=2, Sol=2}。生物碳库最优状态子集= {DEM=3, T_July=3, PET=1, Prec=1, Slope=3, Soil=3, Sol=2}。其中，P(DEM=3)=0.4513, P(T_July=3)=0.3474, P(PET=1)=0.5169, P(Prec=1)=0.3475, P(Slope=3)= 0.5254, P(Soil=3)= 0.4939, P(Sol=2)=0.3362。

图5 研究区生态系统固碳服务平均流动比率空间分区

图6 研究区生态系统固碳服务流动比率总和空间分区

图7 研究区生态系统固碳服务的空间流动示意图

图8 生物碳库变量间的条件概率像素图

图9 土壤碳库变量间的条件概率像素图

同样地，令土壤固碳功能处于高中低三种状态时，其对应的备选变量状态组合称为“关键变量状态子集Ⅳ、Ⅴ、Ⅵ”。那么子集Ⅳ= {DEM=1, NPP=1, Prec=2, Slope=1, Soil=2, Sol=2, T=2}，子集Ⅴ= {NPP=2, Prec=2, Slope=1, Soil=2, Sol=2, T=2}，子集Ⅵ= {DEM=3, NPP=3, Slope=2, Soil=2, Sol=2, T=2}。土壤碳库最优状态子集= {DEM=3, NPP=3, Prec=1, Slope=3, Soil=3, Sol=2, T=3}。其中，P(DEM=3)= 0.4089，P(NPP=3)= 0.9013，P(Prec=1)=0.2736, P(Slope=3)= 0.4397, P(Soil=3)= 0.4118, P(Sol=2)=0.0994，P(T=3)= 0.3387。

2.3.2 状态子集的空间格局 将各变量状态子集分别做空间可视化表达可知，生物碳库子集Ⅰ主要分布在富平县东部韩城县中部以及三原县北部和宝鸡市渭滨区北部；子集Ⅱ主要分布在研究区北部及秦岭北向分支；子集Ⅲ主要分布在研究区北部本山山系。生物碳库最优状态子集则主要分布在研究区西南侧秦岭山脉沿线及秦岭北向支线的部分区域（图10）。

土壤碳库的各组变量子集同样也主要分布在秦岭山脉沿线及本山山系（图11）。子集Ⅳ分布区域由研究区的东北角经关中平原延伸至宝鸡市辖区。子集Ⅴ该子集的分布较前一子集而言更为零散，主要分布在研究区北部，秦岭的北向分支有零星分布。子集Ⅵ主要分布在旬邑县的北部及其与耀州区、淳化县交界处和凤县的北部。土壤碳库最优状态子集则主要分布在研究区西南侧秦岭山脉沿线及秦岭北向支线的部分区域。

2.4 固碳服务空间格局优化

2.4.1 关键变量的最优状态子集 经过熵减算法得出生物、土壤固碳与其备选因子的熵减度（表2），熵减度越高表示该因子对固碳能力的影响越大。

其中高程（DEM）对生物固碳影响最大，熵减度为0.0160，其次为潜在蒸散发（PET），熵减度为0.0082。坡度和降水对生物固碳的影响程度相当，熵减度同为0.0032，太阳辐射和7月平均温度的熵减度紧随其后，土壤类型对生物固碳的影响程度最不显著。而净初级生产力（NPP）对土壤固碳影响最大，熵减度为0.0795，其次为高程（DEM），熵减度为0.0106。气温（T）和太阳辐射（Sol）对土壤固碳的影响程度相当，熵减度同为0.0023，降水（Prec）和坡度（Slope）的熵减度紧随其后，土壤类型（Soil）对土壤固碳的影响程度最不显著。

图10 生物固碳功能关键变量状态子集空间分布图

图11 土壤碳固定功能关键变量状态子集空间分布图

表2 备选因子对固碳的敏感性（按熵减度从大到小排序）

取熵减度前二的因子作为碳库的关键因子，即生物碳库的关键因子为DEM和潜在蒸散发，土壤碳库的关键因子为NPP和DEM。

结合最优状态分析可知：生物碳库的关键变量最优状态子集={DEM=3, PET=1}，其中P(DEM=3)= 0.4513, P(PET =1)= 0.3475；土壤碳库的关键变量最优状态子集={NPP=3, DEM=3}，其中P(NPP=3)= 0.9013, P(DEM=3)= 0.4089。如果能够同时满足某一碳库的两个关键变量最优状态条件，那么该区域是该碳库关键变量的最优状态分布区；如果只满足其中一个，那么该区域是该碳库关键变量的次优状态分布区。

生物碳库关键变量的最优状态子集（图12）主要分布在宝鸡市南部秦岭山脉和天水市西南角，该区域中生物固碳状态为最优的概率可以达到54.36%。而生物碳库关键变量的次优状态子集则主要在天水市、秦岭山脉和北山山系，该区域中生物固碳状态为最优的概率为42.49%。其中分布于宝鸡的关键变量最优子集大多已是省级或国家级森林保护单位，受保护程度较高，但天水市西南侧的区域则有待加强。如果要进一步发展生物固碳服务的储备区域，则可以从关键变量次优状态子集的广大区域着手。

土壤碳库关键变量的最优状态子集（图13）主要分布在秦岭山脉沿线、天水市西南角和咸阳市东北角，这些区域中土壤固碳状态为最优的概率高达92.84%。而土壤碳库关键变量的次优状态子集则主要分布在天水市境内，该区域中生物固碳状态为最优的概率为19.32%。土壤固碳能力的强度很大程度上依赖于植被对土壤的有机质输入，因此保护土壤固碳能力的关键在于保护植被。

图12 研究区生物碳库关键变量最优与次优状态子集分布

图13 研究区土壤碳库关键变量最优与次优状态子集分布

2.4.2 固碳格局优化建议 低碳目标导向的固碳格局应以生物碳库和土壤碳库的关键变量最优状态分布区域作为维持该区域基本碳汇功能的碳汇核心保护区。在此基础上，以生物碳库和土壤碳库的关键变量次优状态分布区域作为继续提升研究区碳汇功能的优选区域，充分利用固碳关键变量最优状态分布区域，能够在维持现有固碳服务供需格局的基础上，进一步增强研究区的碳汇功能，使研究区的低碳经济发展更加有保障。

对比研究区2014年固碳估算结果与生物、土壤碳库关键变量最优状态子集的空间分布情况，得出研究区生物固碳和土壤固碳的适宜优化区域（图14）。这些区域满足固碳关键变量最优状态组合条件，但固碳能力未达到最优状态，可作为研究区低碳目标导向下的固碳格局优化的首选区域。其中，生物固碳的适宜优化区域明显大于土壤固碳适宜优化区。生物碳库的优化应将天水市武山县和秦州区划为重点碳汇格局优化区，加强绿化工程的建设。此外，麦积区、清水县、张家川回族自治县以及凤县和太白县需根据实际情况开展不同程度的生态修复工程。而土壤固碳格局的优化则需要从秦岭中段着手。由于土壤固碳能力在很大程度上依赖于土壤有机质的输入量，因此需要秦岭中段各市县在保护已有森林的同时，改善适宜优化区的植被覆盖，增加其土壤的有机质输入，从而提高适宜优化区的土壤固碳量。

图14 研究区生物固碳与土壤固碳的适宜优化区域

3 讨论

3.1 区域生态系统固碳服务的空间流动

目前关于生态系统服务评估的研究大多数依赖于静态指标，例如Bryan等基于3-PG-2模型预测全球气候变化和土地利用变化背景下的澳大利亚碳汇供给[32]和未来国土部门对于温室气体减排需求[33]。对于固碳服务在供给区与受益区之间的空间流动研究主要还在探索阶段。Serna-Chavez等[34]提出一个全球尺度的多种生态系统服务流分析框架，在气候调节服务部分中分析了碳封存与碳固定的供给热点（provisioning hotspots，p）、受益区（benefiting areas，b）以及流动区（flow area，f）及其相互关系，认为与b、p不重合的区域都是流动区，固碳服务的空间流动可能具有空间依赖性。但关于固碳服务空间流动的实质性研究成果还较少[14]。

本文受到“场”概念和“资源流动”相关研究的启发，以碳源汇流动比率的等值线为基础，借助“等压梯度力”的原理给出研究区碳源汇空间流动示意图，在碳源汇的空间流动上做了可取的尝试，为生态安全、生态补偿等研究领域提供可参考的思路。但是根据SPANs的概念模型，生态系统服务可以通过分布式空间模拟详细地描绘出生态服务的流动路径、流量大小和流动过程等[11-13]，这为本文的下一步研究提供了重要的思路。精确地可视化服务流还需要更系统地考虑固碳服务产生和被消费的过程机制，对服务流的载体、路径、流量、衰减和转化等做具体的描述和研究。以此为基础准确界定固碳服务的供给区和受益区，结合流动路径、流量等将生态系统服务流转化为经济价值，将为生态补偿的应用提供有力的科学依据。

3.2 基于贝叶斯与熵减度的空间格局优化

关于土地利用结构的低碳优化模式与策略的相关研究近年来有所发展，例如赵荣钦等[22]基于线性规划评估了南京市土地利用碳减排潜力，并给出低碳土地利用优化建议。但此类研究多注重各土地利用类型面积的分配上，关于土地利用空间布局上的优化讨论比较欠缺。本文基于贝叶斯原理，借助GIS手段计算关键变量对碳固定影响的不确定性，并将不确定性和关键因子分布做可视化分析，研究结果对于区域的空间格局优化具有一定的参考价值。

当然，该模型也存在一定的缺陷和不足之处：备选因子的选取、判断最优关键变量子集的阈值以及模型本身的结构优化问题依然有待探讨和完善。在未来的相关研究中，应综合考虑更多的因子（例如人文、社会经济发展等因素）以使研究结果更加贴合实际情况。关于空间格局优化，则可以结合未来气象预测、经济发展等对未来土地利用进行更加细致的预测。另外，受方法所限，文章关于不确定性分析还不够深刻，未来研究中应结合多情景分析探讨变化中的多个因子之间的不确定性关系。

3.3 不确定性分析

本文耦合了多个模型以模拟研究区生态系统固碳能力，从空间流动和不确定性的角度较客观地分析了固碳空间格局，但多个模型的耦合也给结果带来了一定的不确定性。

目前人口数据的空间化模型和方法趋于多样化，本文选用的多源数据融合的方法，能够从自然、人文等多方面反映人口的空间分布，对于后续的碳源空间布局分析具有较好的参考价值。但该模型也存在一定的不足，例如参数设置复杂，融合权重的确定较为主观等。目前应用较广的还有基于夜间灯光数据和NDVI及土地利用类型数据的人口数据空间化模型，但该方法更适合于中大尺度的研究，且存在灯光像元过饱和及溢出等问题[23]。随着移动基站等新型数据源的发展，基于“大数据”的人口数据空间化模型也日益引起相关领域的关注[35-36]，这可能是未来人口数据空间化的一大趋势。

在固碳供给估算部分，研究通过野外采样和实验，模拟和验证了土壤表层（30 cm）土壤有机碳含量，提高了土壤固碳估算数据的可靠性。然而受客观条件所限，研究缺乏对NPP估算的验证，这给估算过程带来了一定程度的不确定性。在未来的研究中，将更加注重研究的验证部分，增加净初级生产力的野外考察和验证实验。

此外，受先验知识和研究水平所限，本研究只考虑了关天经济区内部各区县的固碳服务流动，而没有讨论研究区边缘区县与研究区外部的固碳服务交流。在未来研究中，我们将把区域与外界的固碳服务交流纳入研究范畴，力争更加全面和系统地分析区域的固碳服务空间流动。同时，生态系统服务流的研究不应局限于单种生态系统服务的研究和探讨，多种生态系统服务流之间的时空耦合与权衡/协同关系是该领域研究未来发展的重要方向。

4 结论

本文结合多源数据融合模型、CASA模型、遥感碳循环过程模型，量化关天经济区固碳服务的供需平衡，引入流动比率，借助贝叶斯原理和熵模型等多种模型和方法，模拟固碳服务空间流动，探讨研究区固碳服务的格局优化对策。结果表明：研究区固碳服务需求小于供给，具有较明显的空间流动，生物碳库的适宜优化主要分布在天水市武山县和秦州区，土壤碳库的适宜优化区主要分布在秦岭中段各区县。本研究以空间可视化的方式梳理生态系统固碳服务从供给区到受益区的流动过程，揭示区域生态系统固碳服务功能转变的空间规律；使用贝叶斯原理和熵减模型，从概率的角度探讨研究区固碳格局优化的不确定性问题，为生态系统服务流和区域固碳格局优化提供了一种新的思路。

[1] MA B.. Washington DC: New Island, 2005.

[2] KONTOGIANNI A, LUCK G W, SKOURTOS M. Valuing ecosystem services on the basis of service-providing units: A potential approach to address the ‘endpoint problem’and improve stated preference methods., 2010, 69(7): 1479-1487.

[3] ZANDER K K, STRATON A. An economic assessment of the value of tropical river ecosystem services: Heterogeneous preferences among Aboriginal and non-Aboriginal Australians., 2010, 69(12): 2417-2426.

[4] KROLL F, MULLER F, HAASE D, FOHRER N. Rural–urban gradient analysis of ecosystem services supply and demand dynamics., 2012, 29(3): 521-535.

[5] 李双成, 刘金龙, 张才玉, 赵志强. 生态系统服务研究动态及地理学研究范式. 地理学报, 2011, 66(12): 1618-1630.

LI S C, LIU J L, ZHANG C Y, ZHAO Z Q. The research trends of ecosystem services and the paradigm in geography., 2011, 66(12): 1618-1630.(in Chinese)

[6] ANTON C, YOUNG J, HARRISON P A, MUSCHE M, BELA G, FELD C K, Harrington R, Haslett J R, Pataki G, Rounsevell M D. Research needs for incorporating the ecosystem service approach into EU biodiversity conservation policy., 2010, 19(10): 2979-2994.

[7] FISHER B, TURNER R K, MORLING P. Defining and classifying ecosystem services for decision making., 2009, 68(3): 643-653.

[8] BURKHARD B, KROLL F, NEDKOV S, MULLER F. Mapping ecosystem service supply, demand and budgets., 2012, 21: 17-29.

[9] HEIN L, VAN KOPPEN K, DE GROOT R S, VAN IERLAND E C. Spatial scales, stakeholders and the valuation of ecosystem services., 2006, 57(2): 209-228.

[10] BRAUMAN K A, DAILY G C, DUARTE T K E, MOONEY H A. The nature and value of ecosystem services: an overview highlighting hydrologic services.., 2007, 32: 67-98.

[11] PALOMO I, MARTIN-LOPEZ B, POTSCHIN M, HAINES-YOUNG R, MONTES C. National parks, buffer zones and surrounding lands: mapping ecosystem service flows., 2013, 4: 104-116.

[12] TURNER W R, BRANDON K, BROOKS T M, GASCON C, GIBBS H K, LAWRENCE K S, MITTERMEIER R A, SELIG E R. Global biodiversity conservation and the alleviation of poverty., 2012, 62(1): 85-92.

[13] WENDLAND K J, HONZAK M, PORTELA R, VITALE B, RUBINOFF S, RANDRIANARISOA J. Targeting and implementing payments for ecosystem services: Opportunities for bundling biodiversity conservation with carbon and water services in Madagascar., 2010, 69(11): 2093-2107.

[14] 肖玉, 谢高地, 鲁春霞, 徐洁. 基于供需关系的生态系统服务空间流动研究进展. 生态学报, 2016, 36(10): 3096-3102.

XIAO Y, XIE G D, LU C X, XU J. Involvement of ecosystem service flows in human wellbeing based on the relationship between supply and demand.,2016, 36(10): 3096-3102.(in Chinese)

[15] 刘慧敏, 范玉龙, 丁圣彦. 生态系统服务流研究进展. 应用生态学报, 2016, (07): 2161-2171.

LIU H M, FAN Y L, DING S Y. Research progress of ecosystem service flow., 2016, (07): 2161-2171. (in Chinese)

[16] 刘春雨. 省域生态系统碳源汇的时空演变及驱动机制——以甘肃省为例[D]. 兰州大学, 2015.

LIU C Y. The temporal-spatial changes and dynamic mechanism of carbon source/sink of provincial ecosystem—— a case of Gansu province[D]., 2015. (in Chinese)

[17] 于贵瑞, 王秋凤, 朱先进. 区域尺度陆地生态系统碳收支评估方法及其不确定性. 地理科学进展, 2011, 30(01): 103-113.

YU G R, WANG Q F, ZHU X J. Methods and uncertainties in evaluating the carbon budgets of regional terrestrial ecosystems., 2011, 30(1): 103-113. (in Chinese)

[18] 李博, 张文忠, 余建辉. 碳排放约束下的中国农业生产效率地区差异分解与影响因素. 经济地理, 2016, 36(09): 150-157.

LI B, ZHANG W Z, YU J H. Decomposition and influence factors of district difference of China agricultural production efficiency under the constraint of carbon emission., 2016, 36(09): 150-157. (in Chinese)

[19] 蒋金亮, 徐建刚, 吴文佳, 周亮, 孙东琪. 中国人-地碳源汇系统空间格局演变及其特征分析. 自然资源学报, 2014, 29(05): 757-768.

JIANG J L, XU J G, WU W J, ZHOU L, SUN D Q. Patterns and dynamics of China’s human-nature carbon source-sink system., 2014, 29(5): 757-768. (in Chinese)

[20] 李晶, 任志远. 基于GIS的陕北黄土高原土地生态系统固碳释氧价值评价. 中国农业科学, 2011, 44(14): 2943-2950.

LI J, REN Z Y. Research on the values of CO2fixation and O2release by landuse ecosystem in Loess Plateau in Northern Shaanxi Province., 2011, 44(14): 2943-2950. (in Chinese)

[21] 孙伟, 乌日汗. 长三角核心区碳收支平衡及其空间分异. 地理研究, 2012, 31(12): 2220-2228.

SUN W, WR H. Study on the balance of carbon budget and its spatial differentiation in Yangtze River Delta., 2012, 31(12): 2220-2228. (in Chinese)

[22] 赵荣钦, 黄贤金, 钟太洋, 揣小伟. 区域土地利用结构的碳效应评估及低碳优化. 农业工程学报, 2013, 29(17): 220-229.

ZHAO R Q, HUANG X J,ZHONG T Y,CHUAI X W.Carbon effect evaluation and low-carbon optimization of regional land use., 2013, 29(17): 220-229. (in Chinese)

[23] 董南, 杨小唤, 蔡红艳. 人口数据空间化研究进展. 地球信息科学学报, 2016, (10): 1295-1304.

DONG N, YANG X H, CAI H Y. Research progress and perspective on the spatialization of population data., 2016, 18(10): 1295-1304. (in Chinese)

[24] 田永中, 陈述彭, 岳天祥, 朱莉芬, 王英安, 范泽孟, 马胜男. 基于土地利用的中国人口密度模拟. 地理学报, 2004, 59(02): 283-292.

TIAN Y Z, CHEN S P, YUE T X, ZHU L F,WANG Y A, FAN Z M, MA S N. Simulation of Chinese population density based on land use., 2004, 59(2): 283-292. (in Chinese)

[25] 廖顺宝, 孙九林. 基于GIS的青藏高原人口统计数据空间化. 地理学报, 2003, 58(01): 25-33.

LIAO S B, SUN J L. GIS based spatialization of population census data in Qinghai-Tibet Plateau., 2003, 58(1): 25-33. (in Chinese)

[26] YUAN J, NIU Z, WANG C. Vegetation NPP distribution based on MODIS data and CASA model—A case study of northern Hebei Province., 2006, 16(4): 334-341.

[27] 周涛, 史培军, 罗巾英, 邵振艳. 基于遥感与碳循环过程模型估算土壤有机碳储量. 遥感学报, 2007, 11(1): 127-136.

ZHOU T, SHI P J, LUO J Y, SHAO Z Y. Estimation of Soil Organic Carbon Based on Remote Sensing and Process Model., 2007, 11(1): 127-136. (in Chinese)

[28] 李婷, 李晶, 杨欢. 基于遥感和碳循环过程模型的土壤固碳价值估算——以关中-天水经济区为例. 干旱区地理, 2016, 39(2): 451-459.

LI T, LI J, YANG H. Estimate value of soil organic carbon based on remote sensing and process model in Guanzhong-Tianshui Economic Region., 2016, 39(2): 451-459. (in Chinese)

[29] LINDEMAN R L. Experimental simulation of winter anaerobiosis in a senescent lake., 1942, 23(1): 1-13.

[30] LOTKAAJ.. USA:Williams andWilkins, 1925.

[31] 张新林, 赵媛. 基于空间视角的资源流动内涵与构成要素的再思考. 自然资源学报, 2016, 31(10): 1611-1623.

ZHANG X L, ZHAO Y. Rethinking of connotation and constituent elements of resources flow based on spatial perspective., 2016, 31(10): 1611-1623. (in Chinese)

[32] BRYAN B A, NOLAN M, HARWOOD T D, CONNOR J D, NAVARRO-GARCIA J, KING D, SUMMERS D M, NEWTH D, CAI Y, GRIGG N, HARMAN I, CROSSMAN N D, GRUNDY M J, FINNIGAN J J, FERRIER S, WILLIAMS K J, WILSON K A, LAW E A, HATFIELD-DODDS S. Supply of carbon sequestration and biodiversity services from Australia's agricultural land under global change., 2014, 28: 166-181.

[33] BRYAN B A, CROSSMAN N D, NOLAN M, LI J, NAVARRO J, CONNOR J D. Land use efficiency: anticipating future demand for land-sector greenhouse gas emissions abatement and managing trade-offs with agriculture, water, and biodiversity., 2015, 21(11): 4098-4114.

[34] SERNA-CHAVEZ H M, SCHULP C J E, VAN BODEGOM P M, BOUTEN W, VERBURG P H, DAVIDSON M D. A quantitative framework for assessing spatial flows of ecosystem services., 2014, 39: 24-33.

[35] 洪东升. 基于定位数据的人口分布特征研究[D]. 中国地质大学(北京), 2015.

HONG D S. Research on characteristics of population distribution based on positioning data[D]., 2015. (in Chinese)

[36] 毛夏, 徐蓉蓉, 李新硕, 王煜, 李程, 曾波, 何宇华, 刘锦泉. 深圳市人口分布的细网格动态特征. 地理学报, 2010, (4): 443-453.

MAO X, XU R R, LI X S, WANG Y, LI C, ZENG B, HE Y H, LIU J Q. Fine grid dynamic features of population distribution in Shenzhen., 2010, 65(4): 443-453. (in Chinese)

（责任编辑 李云霞）

The Spatial Flow and Pattern Optimization of Carbon Sequestration Ecosystem Service in Guanzhong-Tianshui Economical Region

LI Ting1,2, LI Jing1,2,3, WANG YanZe1,2, ZENG Li1,2

(1School of Geography and Tourism, Shaanxi Normal University, Xi’an 710119;2National Demonstration Center for Experimental Geography Education (Shaanxi Normal University), Xi’an 710119;3Baoji University of Arts and Sciences/Key Laboratory of Disaster Monitoring and Mechanism Simulation of Shaanxi Province, Baoji 721013, Shaanxi)

The aims of this paper were to quantify the balance of supply and demand situation of carbon sequestration service in regional ecosystem(hereinafter referred to as "carbon sequestration service"), simulate the spatial flow of carbon sequestration service, reveal the function change spatial rule of carbon sequestration service in regional ecosystem, and give the regional spatial layout optimization strategy which will provide an intuitive scientific reference for guiding the low carbon development of regions.Population spatial density was simulated with a multi-source data fusion model, and then the demand quantity of carbon sequestration service was estimated in Guanzhong-Tianshui Economic Region. The supply quantity of biologic carbon sequestration service in study area was calculated using CASA model, the soil carbon sequestration service supply was estimated using remote sensing model of carbon cycle. Based on these, the regional carbon balance was quantified with current ratio and the flow process of carbon sequestration service from generation to the use of land was elaborated with spatial visualization method. The conditional probability of environmental variables to the supply of carbon sequestration was calculated rely on Bayes principle, and the key variables were screened by entropy reduction model, the uncertainty of carbon sequestration pattern was discussed by using the distribution of key factor subsets with optimal state, and at last, the carbon sequestration space layout optimization strategy was given.(1) Carbon sequestration service in the study area overall oversupply, the spatial difference of balance between supply and demand is obvious, the high value areas of demand are mainly distributed in the Guanzhong plain high population areas, the high value areas of supply are mainly distributed in the areas along the Qinling Mountains and the Beishan mountains. (2) The spatial flow of carbon sequestration service in ecosystems is obviously different. According to the current ratio distribution, the study area can be divided into three carbon source concentration centers: Tianshui carbon source concentration center (R>0.04), Binxian carbon source concentration center (R>0.04), and Xi’an multi-level carbon source concentration center, with the highest value (R>0.20). The carbon sequestration service space flow in the study area can be divided into four groups: the flow from the Middle and East of Qinling Mountains, Yongshou County in Beishan Muontain to Guanzhong City Group, the flow from the west of Qinling Mountains to Tianshui city, the flow from Linyou county and Xunyi County to Binxian County, the flow from Tongchuan City, Chengcheng County, Hua County to Pucheng County. (3) According to the conditional probability and entropy calculation, taking set{DEM=3, PET=1} as a key variable optimal state subset of biology carbon pool, which is mainly distributed in Baoji City in southern Qinling Mountains and the southwest corner of Tianshui, where the optimal biological carbon fixation probability can reach 54.36%; taking set {NPP=3, DEM=3} as a key variable in the optimal state subset of soil carbon pool, which is mainly distributed in the Qinling Mountains, along the southwest corner of Tianshui and the northeast corner of the city of Xianyang, where the optimal probability is as high as 92.84%. The suitable areas of biology carbon pool mainly distributed in Tianshui and Wushan County, Qinzhou District and the suitable areas of soil carbon pool are mainly distributed in the middle part of Qinling Mountains district.In general, carbon sequestration service demand is less than the supply in study area with obvious spatial flow. Taking Tianshui City and the middle areas of Qinling Mountains as the main region to optimize the carbon sequestration function can get a higher probability of good carbon sequestration pattern.

carbon sequestration; ecosystem services; spatial flow; Bayes; Guanzhong-Tianshui Economic Region

2017-04-27；接受日期：2017-06-16

国家自然科学基金（41371020）、中央高校特别支持项目（GK2015020210）、陕西省教育厅重点实验室科研计划（12JS014）

李婷，E-mail：lite@snnu.edu.cn。通信作者李晶，E-mail：lijing@snnu.edu.cn