一种集成的地理空间建模工作流程可量化土地城市化与碳储量动态。建设用地扩张与碳储量损失及负向空间溢出效应相关,而生态规划情景则有助于保留高碳储量土地,并支持协同的低碳治理。
需要JoVE订阅才能观看此内容。 请登录或开始免费试用
一种集成的地理空间建模工作流程可量化土地城市化与碳储量动态。建设用地扩张与碳储量损失及负向空间溢出效应相关,而生态规划情景则有助于保留高碳储量土地,并支持协同的低碳治理。
本研究开发并应用了一种综合工作流程,结合遥感、土地利用模拟、生态系统服务评估和空间计量经济学方法,量化了中国西北地区新疆维吾尔自治区土地城市化与生态规划对碳储量动态的影响。利用2000–2020年的土地利用/土地覆被、不透水面、夜间灯光、植被、地形、社会经济以及规划约束数据,该工作流程识别了历史城市化模式,模拟了四种未来规划情景,采用InVEST碳模块估算碳储量,并通过空间计量经济模型评估了直接影响和空间溢出效应。结果表明,建设用地扩张主要通过占用耕地实现,并伴随着不透水面覆盖率和夜间灯光强度的增加。模拟结果显示,生态保护情景下碳储量保持最高,而自然发展情景导致的碳损失最大。林地被识别为最易受威胁的高碳土地利用类型,其碳密度为144.9 Mg C/ha,显著高于建设用地的29.4 Mg C/ha。空间杜宾模型结果表明,土地城市化对碳储量密度具有负向直接影响(−0.231)和间接溢出效应(−0.117),总效应为−0.348;而生态规划强度则表现出正向总效应(+0.245)。这些发现支持在行政边界之间开展协同生态规划,以促进区域低碳发展。
在全球气候变化和向低碳发展转型的背景下,快速城市化已成为区域土地覆盖格局和生物地球化学循环变化的主要驱动力1。随着人类活动在规模和强度上的持续扩大,发展中国家和地区的加速土地城市化不仅改变了土地利用格局,也重塑了生态系统的结构与功能。本质上,土地城市化涉及将高碳密度的自然植被转化为低碳密度的人工地表,导致区域碳储存能力大幅下降2。对于快速发展地区而言,如何在经济增长与生态保护空间维护之间实现平衡,已成为国土空间规划中的关键挑战。传统的土地利用研究主要关注物理扩张的模式,往往忽视了空间开发与碳循环过程之间的动态交互关系3。因此,亟需建立综合分析方法,量化土地城市化对碳储量动态的影响,为生态保护和低碳治理提供科学依据。
近年来,遥感技术和地理信息系统的发展显著提升了土地利用变化建模与生态系统服务评估的精度4,5。结合生态核算的多情景土地利用模拟已成为评估城市扩张环境影响的常用方法6。特别是,斑块生成式土地利用模拟模型(Patch-generating Land Use Simulation, PLUS)与生态系统服务综合评估模型(Integrated Valuation of Ecosystem Services and Tradeoffs, InVEST)已被广泛应用于量化规划约束对土地利用演变路径和碳储量结果的影响7,8,9,10,11,12。相关生态系统服务研究表明,土地覆被变化、城市发展政策以及生态保护措施会改变快速城市化区域的水源供给、碳固存、土壤保持功能及生态系统服务之间的权衡关系13,14,15,16,17。然而,大多数以往研究主要关注土地利用转换及其生态系统服务输出,较少关注社会经济发展压力如何在相邻行政单元之间引发具有空间依赖性的碳储量响应。
当前工作流程的提出主要基于三个关键的研究空白。首先,以往的 PLUS–InVEST 应用通常将研究区域视为封闭的物理系统,未明确评估土地城市化效应是否会溢出至相邻区域。其次,验证分析、不确定性分析和敏感性分析通常仅被简要报告,使得难以判断在不同土地利用模拟参数、碳库系数和空间权重设定条件下,情景排序是否保持稳定。第三,空间计量经济学研究常识别出开发强度与生态结果之间的关联,却未将这些关系与遥感观测、未来土地利用情景及碳储量核算相结合。因此,亟需一种集成化的工作流程,在可重复的分析框架内实现观测数据、情景模拟、碳储量核算与空间机制识别之间的有效衔接。
因此,本研究建立了一种可重复的分析工作流程,该流程整合了遥感技术、基于PLUS的情景模拟、InVEST碳储量评估以及空间计量经济分析。研究聚焦中国西北部的新疆维吾尔自治区,评估了2000年至2020年期间的历史土地城市化格局,探讨了不同规划情景对未来土地利用格局的影响,量化了其对碳储量分布的相应影响,并分析了土地城市化、生态保护规划强度、植被状况及社会经济因素对碳储量动态的直接影响和空间溢出效应。该工作流程的主要贡献在于将基于情景的预测与机制识别相结合,以支持生态红线划定、紧凑型增长规划、跨区域碳补偿以及低碳区域治理。
访问受限。请登录或开始试用以查看此内容。
本研究使用了公开可用的地理空间、遥感、规划以及汇总的社会经济数据集。未涉及人类受试者、动物实验、临床材料或可识别的个人信息,因此无需机构伦理审批。所有用于重现本工作流程所必需的资源,包括特定软件环境、计算软件包(含空间计量经济学和高分辨率可视化库)、数据集数字对象标识符(DOI)及来源网址,均在详尽的材料表中列出。
1. 研究区域与分析框架
分析在中国西北部的新疆维吾尔自治区进行,重点针对主要城市的划定城区范围(图1)开展具体分析,以准确捕捉局部扩张动态。该边界被统一用于栅格裁剪、土地利用模拟、碳储量评估以及行政单元的聚合。分析样本量(n)包括在三个离散时间点(2000年、2010年和2020年)上分析的105个县级行政单元,土地利用模拟中纳入了50次随机重复,以反映空间不确定性。
采用了一种整合栅格尺度与行政单元尺度分析的多尺度分析框架。在栅格尺度上,利用土地利用/土地覆被数据、夜间灯光强度、不透水面覆盖率、植被指数、地形变量以及生态规划约束条件,识别土地城市化模式,模拟未来土地利用情景,并评估碳储量空间分布。在行政单元尺度上,将社会经济、交通和环境变量纳入空间计量经济分析,以评估土地城市化和生态规划对碳储量动态的直接影响与溢出效应。
分析框架包含五个组成部分。首先,利用多时相遥感与土地利用数据,刻画土地城市化的时空格局,包括建设用地扩张、开发强度以及土地利用重构。其次,构建生态保护、耕地保护和低碳规划情景,以模拟未来的土地利用格局。第三,采用 InVEST 碳储量模块,量化在历史及未来土地利用条件下碳储量的分布及其变化。第四,开展全局与局部空间自相关分析,识别与碳储量变化相关的集聚模式及热点区域。第五,运用空间计量经济模型,评估土地城市化、生态规划强度以及环境与社会经济驱动因素对碳储量动态的影响程度、方向及溢出效应。
图1 展示了整体分析工作流程。该流程连接了研究区域界定、多源数据整合、土地城市化识别、生态规划情景构建、土地利用模拟、碳储量评估、空间自相关分析以及空间计量经济建模。多源空间数据作为输入层;土地城市化识别与历史土地利用分析构成模式识别层;生态规划情景与碳储量评估构成预测与影响评估层;空间自相关分析与计量经济分析则用于识别空间作用机制并支持政策解读。

图1:研究区域、数据整合与技术工作流程。(A)中国西北部新疆维吾尔自治区主要城市的多尺度空间背景及其精确的城市分析边界。(B)集成技术工作流程,整合多源空间数据、基于遥感的土地城市化识别、生态规划情景设计、PLUS土地利用模拟、InVEST碳储量评估、空间自相关分析、热点分析以及空间计量经济建模。ND = 自然发展;EP = 生态保护;CP = 耕地保护;LC = 低碳优化;OLS = 普通最小二乘法;SAR = 空间自回归模型;SEM = 空间误差模型;SDM = 空间杜宾模型。请点击此处查看该图的高清版本。
2. 数据来源与变量系统
整合了多源数据集以支持土地城市化识别、生态规划情景模拟、碳储量评估和空间计量经济分析。该数据集包括土地利用/土地覆被数据、不透水面比例、夜间灯光强度、归一化植被指数(NDVI)、地形变量、交通通达性指标、水文变量、人口密度、国内生产总值(GDP)密度、规划约束条件、气候变量以及碳库参数。表1汇总了数据集的类别、变量描述、单位、空间和时间分辨率、数据来源及分析应用。关键数据集分别获取了2000年、2010年和2020年的数据。根据分析需求,土地利用数据和主要遥感产品被统一重采样至30米至1公里不等的空间分辨率。
表1:数据来源与变量体系。 该表格列出了用于土地城市化识别、情景模拟、InVEST碳储量核算以及空间计量模型分析的数据集类别、变量、描述、单位、空间/时间分辨率、来源及分析用途。NDVI = 归一化植被指数;DEM = 数字高程模型;LULC = 土地利用/土地覆盖。请点击此处下载该表格。
所有空间图层均投影至统一的坐标系,按图1所示研究区域边界进行裁剪,并重采样或聚合至所需网格或行政单元尺度。在进行空间叠加分析前,已对缺失值进行筛查。PLUS 和 InVEST 分析所用的栅格数据集在像元级别上进行了对齐,而社会经济变量则聚合至行政单元以进行空间回归分析。该预处理流程确保了在时间序列分析、情景模拟、碳储量计算和计量经济建模中空间单元的一致性和时间点的可比性。
土地利用/土地覆盖数据被用作历史变化检测和未来土地利用模拟的主要输入数据。原始土地利用类别被重新分类为六类:耕地、林地、草地、水体、建设用地和未利用地。该分类体系在土地利用变化分析、情景构建以及InVEST碳储量计算中保持一致。建设用地比例、不透水表面覆盖率和夜间灯光强度被用作土地城市化的指标。不透水表面覆盖率和夜间灯光强度作为反映开发强度和人类活动的遥感代用指标。归一化植被指数(NDVI)被纳入后续空间计量经济分析中作为生态指标。
地形和区位变量包括高程、坡度、距主要道路的距离、道路密度以及距主要河流的 proximity。这些变量被用作土地利用模拟中的驱动因子,以及空间计量经济模型中的控制变量。它们共同表征了地形限制、交通可达性和水文连通性。
社会经济变量包括人口密度和GDP密度,并在数据可获得的情况下补充了道路网络密度。这些变量用于表征与碳储量变化相关的发展强度和人类活动压力。该变量框架整合了土地开发、社会经济压力和生态背景的指标,以反映影响碳储量动态的多种因素。
利用生态保护红线、自然保护区、水体缓冲区、坡度限制区及其他环境敏感区域,构建了生态规划强度指数。该指数既可作为未来土地利用模拟中的约束图层,也可作为空间计量经济分析中的解释变量。
在 InVEST 碳储量评估中,为每一类土地利用类型汇编了四个碳库参数:地上生物量碳、地下生物量碳、土壤有机碳和死有机物碳。参数值来源于已发表的区域研究、InVEST 指南文件以及当地土地覆盖特征,并被映射到统一的土地利用分类系统中。碳储量是指这四个碳库中所含碳的总量,而碳储量密度则是指单位面积上的碳储量。
变量被分为三类。第一类包括土地城市化指标,如建设用地比例、不透水面覆盖度和夜间灯光强度。第二类包括生态与规划变量,如归一化植被指数(NDVI)、生态规划强度、高程、坡度和水文邻近度。第三类包括社会经济驱动因素,如人口密度、GDP密度和道路密度。这些变量用于评估土地城市化、生态规划与碳储量动态之间的关系。
3. 基于遥感的城镇化土地识别
通过三个互补的指标识别土地城市化的时空格局:建成区扩张、不透水地表覆盖和夜间灯光强度。建成区扩张信息从连续的土地利用图中提取,用于界定城市开发从核心建成区向周边地区扩展的实际范围和演进过程。计算不透水地表覆盖以量化网格尺度上的开发强度,同时对夜间灯光强度栅格数据进行处理和归一化,以表征人类活动强度和功能集聚程度。随后将这些空间图层整合,构建历史土地城市化数据集。综合上述指标用于识别持续高强度开发区域、过渡带以及相对稳定的区域,并评估土地实体开发与功能城市化之间的空间对应关系。
构建了土地利用转移矩阵,以量化耕地、林地、草地、水域、建设用地和未利用地之间土地利用转换的强度和方向。特别关注了生态用地和农业用地向建设用地的转换。该矩阵用于识别主要的转换路径以及对城市扩张贡献最大的主要来源土地类型。
通过景观格局分析评估与土地城市化相关的结构变化。所采用的指标包括斑块密度、边缘密度、景观形状指数以及与破碎化相关的指标。这些指标针对每个研究时段进行计算,用于量化与建设用地扩张相关的景观格局配置、空间连续性及破碎化程度的变化。
4. 生态规划情景设计与土地利用模拟
在四种规划情景下,采用斑块生成型土地利用模拟(Patch-generating Land Use Simulation, PLUS)模型模拟了未来的土地利用格局,这四种情景分别为:自然发展(ND)、生态保护(EP)、耕地保护(CP)和低碳优化(LC)。情景假设、土地转换规则、限制性土地类型及预期碳储量结果汇总于表2中。
表2:情景控制规则与转移约束。 该表定义了自然发展(ND)、生态保护(EP)、耕地保护(CP)和低碳优化(LC)情景下的土地转换概率与空间约束条件。注:情景缩写严格对应文中PLUS模型模拟与结果部分所使用的缩写。情景设置定义了模拟中所采用的基于政策导向的土地转换规则及受限土地类型;所有法定禁止开发区域在最终的空间图层中均保持不可转换。 请点击此处下载该表格。
将历史土地利用地图和空间驱动变量整合到PLUS模型中,利用元胞自动机(CA)模块估算基线土地转移概率。对于生态保护(EP)、耕地保护(CP)和低碳优化(LC)情景,将生态保护红线、永久基本农田、水体缓冲区及其他规划约束作为空间限制图层纳入模型,以根据预设的情景规则限制土地转换。随后确定各情景下2030年的目标土地需求量,并通过执行50次随机重复模拟运行CA模块,生成最终的土地利用预测结果。
采用 PLUS 模型及其多类型随机斑块生成的元胞自动机(CA)框架进行土地利用模拟。利用历史土地利用图以及环境和社会经济驱动因子变量,估算每类土地利用类型的扩张概率。随后,根据各规划情景的需求设定未来的土地需求量,并生成目标模拟时段对应的土地利用图。
在进行未来模拟之前,通过历史回溯验证模型性能。利用早期的土地利用图及其相关驱动变量,模拟后期实际观测到的土地利用图,并通过总体精度(OA)、Kappa系数和优势度指数(FoM)评估模拟结果与观测结果之间的吻合程度。验证在总体层面以及主要土地利用类型层面分别进行。从2010年至2020年的历史回溯模拟获得了93.4%的总体精度(OA)、0.89的Kappa系数和0.26的优势度指数(FoM),表明该模型在后续多情景模拟中具有高度可靠的空間预测能力。
为评估模拟结果的稳健性,开展了不确定性分析和敏感性分析。PLUS 模型的敏感性分析考察了主要土地利用类型在不同转移阻力设置和邻域权重参数下的影响。InVEST 模型的敏感性分析评估了不同土地利用类型中碳库系数变化的影响。空间计量经济学敏感性分析则比较了不同空间权重矩阵的设定方案。这些分析用于判断在不同参数设置下,情景排序以及主要城市化与生态规划效应的方向是否保持一致。具体而言,结果的稳定性得到了严格验证:当转移阻力参数和碳库系数在 ±15% 范围内变动时,情景排序以及土地城市化的负向空间溢出效应均保持不变。
5. 碳储量评估
采用 InVEST 碳模型评估碳储量。根据表3汇总的生物物理参数,为每种重新分类的土地利用类型分配了四个碳库(地上生物量、地下生物量、土壤有机碳和死有机物)的基准参数。随后将历史土地利用图(2000–2020年)以及PLUS模型模拟的未来土地利用栅格数据集导入该模型,并与相应的碳密度参数进行整合。接着运行模型,以估算区域总碳储量(Tg C)、网格级别的碳储量密度(Mg C/ha),以及历史和未来情景下的碳储量变化空间分布图(ΔC)。
重新分类的土地利用图与表3中报告的相应基准碳库参数值进行了关联。通过将这些基准参数值调整±15%,检验了碳储量估算对参数不确定性的敏感性,并给出了详细的敏感性分析结果。该模型用于计算在历史条件和未来规划情景下,每种土地利用类型的总碳储量、碳储量密度以及碳储量变化。
评估了三类输出结果。首先,计算了区域总碳储量及其时间变化趋势,以量化碳储量随时间变化的幅度和方向。其次,绘制了碳储量及碳储量变化的空间分布图,以识别碳保留和碳损失的区域。第三,比较了不同规划情景下的碳储量估算结果,以评估生态保护、耕地保护和低碳优化策略对碳储量保护的相对影响。
表3:不同土地利用类型的基线碳库参数。 本表列出了 InVEST 碳储量模块中使用的地上生物量碳、地下生物量碳、土壤有机碳、死有机物碳以及总碳密度。单位为 Mg C/ha。注:数值代表 InVEST 模型中使用的基线参数。总碳密度等于四个碳库之和。对这些基线值进行 ±15% 扰动的敏感性分析见表6。请点击此处下载该表格。
6. 空间自相关与空间计量经济学分析
通过空间自相关分析,确定碳储量及其变化是否表现出显著的空间依赖性。计算全局Moran's I指数,以评估研究区域内碳储量分布及碳储量变化的整体空间聚类程度。随后采用局部Moran's I指数识别局部空间关联模式,包括高-高、低-低、高-低和低-高聚类。通过热点分析识别碳损失集中区域和碳储量保持区域。
采用空间计量经济模型,分析土地城市化、生态规划、环境条件、社会经济因素与碳储量动态之间的关系。以碳储量密度或碳储量变化作为因变量,解释变量包括土地城市化指数、建设用地比例、生态规划强度、归一化植被指数(NDVI)、人口密度、GDP密度、道路密度、海拔、坡度、年平均降水量和年平均气温。
采用普通最小二乘法(OLS)回归作为基线模型。在估计空间自回归(SAR)、空间误差(SEM)和空间杜宾(SDM)模型之前,评估了残差的空间依赖性。对不同模型设定下的模型性能和系数估计结果进行了比较。
构建了一个行标准化的空间权重矩阵,以表示行政单位之间的邻接关系。基本设定基于空间邻接性,而在可获得的情况下,稳健性分析比较了基于距离和最近邻的空间权重矩阵的替代方案。
通过SDM计算直接效应、间接效应和总效应,以评估局部和空间溢出关系。将地形、可达性、植被、气候和社会经济变量作为控制变量纳入模型,以减少遗漏变量偏差。模型系数被解释为条件性空间关联,而非确定性的因果效应。对于p值超过常规显著性阈值的变量,解释为证据较弱或具有提示性,不视为具有统计学意义的稳健效应。
访问受限。请登录或开始试用以查看此内容。
基于遥感的土地城市化识别
根据多时期遥感识别结果,研究区在2000年至2020年间经历了显著加剧的土地城市化进程,伴随土地利用格局的重构和景观格局的重塑。总体而言,建设用地的扩张呈现出由核心城区向外围区域辐射的同心圆式扩展模式。人类活动强度和地表不透水层水平的同步上升表明,快速城市化不仅改变了土地利用的数量结构,也深刻重构了区域空间组织格局。
以图2A中乌鲁木齐代表性子区域为例,土地城市化识别结果清晰地展示了2000–2020年间建设用地的阶段性扩张轨迹。2000年,建设用地主要集中于城市中心区域,空间格局相对紧凑。至2010年,外围扩张显著加剧,在核心区周围形成过渡带。到2020年,建设用地已突破原有紧凑边界,向多个方向扩展,并在远离核心区的位置形成离散的扩张斑块,反映出典型的向外扩张模式和多中心扩散趋势。与这些扩张模式一致,图2B显示,至2020年,研究区域内不透水面覆盖度呈现出显著的由中心向外围递减的梯度分布。最高值...
访问受限。请登录或开始试用以查看此内容。
本研究提出了一种综合分析工作流程,该流程结合了基于遥感的土地城市化识别、PLUS多情景土地利用模拟、基于InVEST的碳储量评估以及空间计量经济建模,以分析中国西北地区新疆维吾尔自治区土地城市化、生态规划与碳储量动态之间的关系。该工作流程并非将这些组成部分视为独立的技术操作,而是将观测到的土地转化、未来规划方案、碳库变化后果以及空间溢出机制整合于一个可重复的统一分析序列中。
时空分析表明,2000年至2020年间,研究区域的土地城市化进程遵循离心扩张的轨迹,从集中式的中心城区增长逐步转变为多方向、多中心的蔓延模式。这一模式与广泛的城市化文献一致,这些文献记录了快速工业化地区中跳跃式发展的现象,其中基础设施投资与经济集聚共同推动了城市开发边界的外扩。关键的是,主要的土地转化路径——耕地转为建设用地——证实了农业用地仍在承受城市侵占最沉重的压力,这一发现与中国及其他发展中经济体关于城乡结合部土地动态的研究结果相符。然而,与我国东部高度一体化的大型城市群不同,新疆地区的城市化受到绿洲边界和水资源可获得性的严格限制。因此,这种侵占对紧邻重要水文廊道的有限耕地和生态过渡带构成了不成比例的威胁...
访问受限。请登录或开始试用以查看此内容。
作者无任何利益冲突需要披露。
作者未获得此项工作的特定资助。
访问受限。请登录或开始试用以查看此内容。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 自定义空间计量经济学脚本与输出数据集 | 作者 | Zenodo 仓库:https://zenodo.org/records/21159171 | 用于重现本研究中描述的空间计量经济学分析及支持性计算工作流程。 |
| 高分辨率人口密度网格 | WorldPop | https://www.worldpop.org/ | 作为空间计量经济学分析的社会经济输入数据。 |
| InVEST(生态系统服务与权衡综合评估) | 自然资本项目,斯坦福大学 | 版本 3.14.0;https://naturalcapitalproject.stanford.edu/ | 用于估算历史和模拟土地利用情景下的碳储量及碳储量变化。 |
| Matplotlib & 枕头(Python 库) | Python 社区 | https://matplotlib.org/;https://python-pillow.org/ | 用于生成出版级图表并处理高分辨率(600 dpi)图形输出。 |
| 多时相土地利用/土地覆盖数据集 | 中国科学院资源与环境科学数据中心 | https://www.resdc.cn/ | 作为土地利用变化分析和PLUS模拟的主要输入数据集。 |
| PLUS(斑块生成型土地利用模拟)模型 | 高性能空间计算智能实验室 (HPSCIL) | https://github.com/HPSCIL/Patch-generating_Land_Use_Simulation_Model | 用于模拟多种生态规划情景下的未来土地利用模式。 |
| PySAL(Python 空间分析库) | PySAL 开发者 | https://pysal.org/ | 用于进行空间自相关分析和空间计量经济学建模。 |
| Python 编程环境 | Python 软件基金会 | 版本 3.9+;https://www.python.org/ | 用于数据预处理、空间分析、统计计算和工作流程实施。 |
| 地形(DEM)与遥感数据 | 美国地质调查局(USGS) | https://earthexplorer.usgs.gov/ | 用于推导土地城市化分析和模型参数化所需的地形变量及遥感输入数据。 |
访问受限。请登录或开始试用以查看此内容。