黄河流域碳排放量及影响因素的时空异质性研究.pdf

收稿日期：基金项目：国家社会科学基金重点项目（）作者简介：孙才志（），男，山东烟台人，教授，博士，主要从事水资源经济与海洋经济地理研究工作：【“黄河流域生态保护和高质量发展”专题】黄河流域碳排放量及影响因素的时空异质性研究孙才志，曾庆雨（教育部人文社科重点研究基地 辽宁师范大学 海洋经济与可持续发展研究中心，辽宁 大连；辽宁省海洋经济高质量发展高校协同创新中心，辽宁 大连；辽宁师范大学 地理科学学院，辽宁 大连）摘 要：低碳城市化对黄河流域生态保护和高质量发展起关键作用。利用莫兰指数测度黄河流域地级市碳排放量空间关联集聚性，并利用多尺度地理加权回归（）模型研究碳排放量影响因素时空异质性，结果表明：）年黄河流域碳排放量上升，呈东北部高、西南部低的空间分布格局，空间集聚特征显著，高高值集聚区位于鄂尔多斯市、榆林市、吕梁市、运城市和三门峡市，低低值集聚区位于西宁市和黄南藏族自治州等；）利用 模型研究碳排放量影响因素空间异质性优于 模型；）有显著空间异质性，、人口规模和土地城市化率对碳排放量影响为正效应，人口密度对碳排放量为负效应。关键词：碳排放量；模型；莫兰指数；黄河流域中图分类号：；文献标志码：引用格式：孙才志，曾庆雨黄河流域碳排放量及影响因素的时空异质性研究人民黄河，（）：，（，；，；，）：（），：），“”，），：；黄河流域是我国重要的能源、化工、原材料和基础工业基地，也是我国碳排放重点地区和生态治理重点区域。年 月，习近平总书记在中央财经委员会第九次会议上指出，要把碳达峰、碳中和纳入生态文明建设整体布局。年 月 日，生态环境部等四部门联合印发黄河流域生态环境保护规划，提出黄河流域 年前实现碳排放达峰，生态系统质量和稳定性全面提升。碳达峰目标和生态文明建设互为支撑、相互促进。一方面，生态文明建设是实现碳达峰目标的基础，高水平生态文明建设能促进经济绿色低碳发展、优化资源配置，推动高质量发展；另一方面，发展低碳经济是推动流域生态文明建设的重要路径。水生态文明建设作为生态文明发展的重要环节，是推动黄河流域经济社会高质量发展和实现碳达峰的重要影响因素。因此，研究黄河流域碳排放量影响因素时空异质性，对推动黄河流域的生态文明特别是水生态文明建设、实现黄河流域生态保护和高质量发展具有重要现实意义。当前对碳排放量影响因素的研究主要采用对数平第 卷第 期 人 民 黄 河 ，年 月 ，均迪 氏 指 数（）分 解 法、恒 等 式、模型等方法，研究范围包括工业和农业等部门、流域、省域、城市群、市域和县域等不同尺度。有学者采用经典地理加权回归（）模型和时空地理加权（）模型，分别对碳排放量及其影响因素时空异质性进行研究。通过梳理文献发现，目前对碳排放量影响因素的研究大多侧重于全局平均作用强度的估计，而 模型和 模型所有自变量带宽（带宽是衡量影响因素空间异质性的指标，带宽越小空间异质性越强）相同，无法反映各影响因素空间尺度差异，而多尺度地理加权回归（）模型则允许不同自变量拥有不同带宽，可真实有效地反映影响因素的空间异质性。鉴于此，本文以黄河流域地级市为研究对象，利用 模型研究其碳排放量影响因素的时空异质性，以期为黄河流域生态保护和高质量发展提供参考。研究区概况黄河全长 ，流经青海、四川、甘肃、宁夏、内蒙古、山西、陕西、河南、山东 个省（区），涉及 个地级行政单位，流域面积约为 万。黄河流域横跨我国西部、中部和东部，流域内城市经济发展水平差异较大。黄河流域 个省（区）年总人口亿，占全国总人口的；年地区生产总值 万亿元，约占全国的。黄河流域是我国煤炭主要生产和供应基地，煤炭产量约占全国煤炭总产量的。化石燃料燃烧是碳排放的主要来源，因此黄河流域面临巨大的碳减排压力。数据来源与研究方法 变量选取与数据来源 等研究表明，地级市尺度的、人口和城市化对碳排放量影响较大，因此本研究选取、人均、年末总人口、人口密度和土地城市化率 个指标作为自变量，研究其对碳排放量的影响。研究时段为 年，研究对象为黄河流域地级市，流域外部区域不在本研究范围内。为保证数据一致性，本研究采用的所有数据均从分辨率为 的卫星遥感影像按黄河流域边界范围内提取得到。用于提取城市建设用地面积的夜间灯光数据来自美国国家环境信息中心（：），用于计算城市行政区土地总面积的空间数据来自中国科学院资源环境科学数据中心（：），人口空间分布数据来自（：），排放空间分布数据来自人为二氧化碳开源数据清单（），空间分布数据来自全球高分辨率实际国内生产总值栅格数据。对处理后的自变量进行多重共线性检验，各指标的方差膨胀系数（）均小于，不存在显著共线性。研究方法 土地城市化率计算根据土地城市化的内涵，本研究用城市建设用地面积与城市行政区土地总面积的比例表征土地城市化率。利用 和 夜间灯光数据集，提取城市建设用地面积。处理步骤：）数据集（年）相互校准、年内融合和年际校正；）数据集（年）背景噪声和异常值去除、年度数据合成和年际校正；）根据 与 数据集在 年数据重叠的特点，建立二者之间的回归关系，得到校正后的 数据集（年）；）采用伪不变特征法提取城市建设用地面积。全局空间自相关性分析全局莫兰指数（）根据要素位置和要素属性值来度量要素的空间相关联程度，计算公式为 （）（）（）（）式中：为全局空间自相关指数，为地级市总数，和分别为地级市 和地级市 的样本值，为样本值的平均值，为空间权重。局部空间自相关性分析采用 的 工具，识别具有高值或低值要素的空间聚类情况，公式为 （）式中：和分别为地级市 和地级市 的碳排放量标准化值，为地级市 的局部空间自相关指数。和 模型 模型表达式为（，）（，）（）式中：为地级市 的碳排放量，（，）为地级市 的空间地理位置坐标，（，）为地级市 的截距，（，）为地级市 的第 个影响因素的回归系数，为地级市 的第 个影响因素的值，为误差项，为影响因素总数。人 民 黄 河 年第 期 模型表达式为（，）（，）（）式中：（，）为地级市 的第 个影响因素的局部回归系数。可以看出，模型的每个回归系数都是基于局部回归得到的，带宽具备特异性，这是 模型与 模型最大的不同。结果与分析 黄河流域碳排放量的时空格局及空间关联差异分析 黄河流域碳排放量时空格局黄河流域地级市平均碳排放量由 年的 万 增至 年的 万（见图）。年之前增速较快，之后增速变慢。碳排放量高值区主要位于上游与中游交错地带的鄂尔多斯市等地级市，整体呈东高西低的空间分布格局（见图）。图 年黄河流域地级市平均碳排放量图 黄河流域典型年碳排放量及变化速度空间分布 碳排放量空间自相关性分析采用全局莫兰指数表征黄河流域碳排放量空间自相关性（其计算结果见表），全局莫兰指数均为正值，且通过显著性水平为 的显著性检验，表明黄河流域地级市碳排放量存在显著正向空间自相关特征，即碳排放量水平接近的区域在空间上集中分布。采用 软件中的 工具进行局部自相关性分析，结果见图。高高值区主要分布在黄河上游与中游交错地带（鄂尔多斯市、榆林市和吕梁市）、中游与下游的交错地带（运城市、三门峡市），低低值区主要位于青海省西宁市和黄南藏族自治州。表 全局莫兰指数 计算结果年份全局莫兰指数 碳排放量标准化值显著性水平年份全局莫兰指数 碳排放量标准化值显著性水平图 黄河流域典型年碳排放量局部空间自相关性 模型对比与尺度分析采用 法对黄河流域碳排放量及影响因素进行相关性分析，结果见图（、分别表示通过显著性水平为、的显著性检验）。碳排放量与、年末总人口、土地城市化率的相关性均通过显著性水平为 的显著性检验，成显著正相关关系，与碳排放量的相关系数最大。人口密度与碳排放量的相关性在 年和 年分别通过显著性水平为和 的显著性检验。图 黄河流域典型年碳排放量与影响因素的相关关系人 民 黄 河 年第 期由空间自相关分析结果可知碳排放量有显著空间集聚特征，说明碳排放量与影响因素之间不满足普通最小二乘法要求，即区域之间相互独立的假设不成立。因此，需引入空间差异分析模型对经典线性模型进行改进。采用 模型和 模型对碳排放量及其影响因素异质性进行研究，结果见表。模型的（修正的赤池信息准则值）和（残差平方和）均比 模型的小，说明 模型优于 模型；模型的（拟合优度）和调整 均大于 模型的，说明 模型拟合效果优于 模型。因此，模型可降低模型残差的空间自相关性，是研究黄河流域地级市碳排放量影响因素空间异质性的较优模型。表 模型和 模型对比年份模型调整 模型和 模型变量带宽对比见表。模型变量带宽数值较大，不具有空间异质性。模型的 带宽最小，说明其具有显著空间异质性；土地城市化率的带宽从 年开始减小，说明其空间异质性逐步显现；常数项表示区位因素对碳排放量的影响程度，其带宽 年较小，表明在此期间区位因素对黄河流域地级市碳排放量的影响具有空间异质性；人均、年末总人口和人口密度的带宽较大，为全局尺度，说明这些因素对碳排放量的影响不具有空间异质性。表 模型和 模型变量带宽对比年份 模型（所有自变量）模型常数项人均年末总人口人口密度土地城市化率地级市拟合优度 反映各影响因素对黄河流域地级市碳排放量的解释程度。黄河流域 模型地级市 空间分布见图，可知黄河流域各地级市 均在 以上；年的地级市 大于，说明各影响因素对黄河流域地级市碳排放量空间分布特征具有较强解释程度；黄河中下游地区的地级市 相对较小，说明中下游地区地级市碳排放量还存在其他影响因素。图 黄河流域 模型地级市拟合优度 空间分布 基于 模型的黄河流域碳排放量影响因素异质性分析利用 模型计算 黄河流域碳排放量各影响因素回归系数年平均值，结果见图。由 相关性分析可知，碳排放量与人均 没有显著相关性，因此碳排放量影响因素异质性分析不包括人均。按回归系数大小排序：人口规模土地城市化率人口密度。、人口规模和土地城市化率为正效应，人口密度为负效应。图 基于 模型的黄河流域碳排放量影响因素回归系数人 民 黄 河 年第 期 黄河流域碳排放量与 回归系数空间分布见图。的回归系数随时间减小，即正效应强度逐年减弱，说明黄河流域各地级市在发展经济的同时，建设绿色低碳城市，可持续发展意识加强。正效应强度自东南向西北增强，说明黄河流域上游地区城市经济发展对碳排放量增加的贡献较大，这与张华明等的观点一致。原因是上游地区经济相对落后，碳减排能力提升难度较大；下游地区经济相对发达，生态保护和发展水平较高，碳减排能力较强。图 黄河流域典型年碳排放量与 回归系数空间分布 黄河流域碳排放量与年末总人口回归系数空间分布见图。人口规模正向影响碳排放量，说明城市人口规模扩大，居民能源消耗总量增长，会增大城市碳排放量。回归系数高值区自东向西移动至中部后向南收缩。年后，正效应最强的区域位于黄河流域南部的西安市等地级市，说明以西安为核心的经济发展水平较高的地区，由经济发展和人口激增带来的碳排放量高的问题较为突出。图 黄河流域典型年碳排放量与年末总人口回归系数空间分布 由于人口密度对碳排放量的影响为负效应，因此人口密度与碳排放量的回归系数绝对值越大，负效应越强。选择具有显著相关关系的 年、年和 年的回归系数进行空间分布分析（见图）。可以看出，人口密度负效应强度从西南向东北部增强，较强地区位于黄河流域下游以及中上游北部鄂尔多斯市及周边地级市，说明经济发展水平较高的城市规划合理、公共资源集聚，人口集聚有利于碳减排。图 黄河流域典型年碳排放量与人口密度回归系数空间分布 黄河流域典型年碳排放量与土地城市化率回归系数空间分布见图。由于 年和 年碳排放量与土地城市化率回归结果不显著，因此选择具有显著性的 年回归系数进行分析。土地城市化率的回归系数平均值从 年的 上升至 年的，其中经历 年和 年 个上升期。正效应强度自东向西逐渐增强，说明黄河上游地区土地城市化率的提高使碳排放量增大。该空间分布格局是以下因素共同作用的结果：）上游地区土地城市化基础薄弱，随着西部大开发战略的实施，需要更多的城市建设用地以发展工业企业，而工业企业以煤炭为主要燃料，导致碳排放量增加；）城市建设用地扩张的主要来源是耕地（由 年和 年土地利用分类图可知），而耕地整体表现为碳汇，耕地转化为城市建设用地导致陆地生态系统碳汇能力下降；）黄河上游地区碳减排能力低、提升难度大；）黄河下游地区城市发展较为完善，城市建设用地扩张空间有限，土地城市化率的影响较弱。人 民 黄 河 年第 期图 黄河流域典型年碳排放量与土地城市化率回归系数空间分布 结论在我国 年前实现碳达峰目标与黄河流域生态保护和高质量发展重大国家战略实施背景下，本文利用莫兰指数对黄河流域地级市碳排放量与影响因素的空间自相关特征进行分析。运用 模型，从经济、人口、土地城市化率等方面对碳排放量影响因素时空异质性进行研究，结论如下。）黄河流域地级市碳排放量呈上升趋势，具有显著空间集聚特征。高高值集聚区主要位于黄河上中游交错地带的鄂尔多斯市、榆林市、吕梁市，中下游交错地带的运城市和三门峡市；低低值集聚区位于黄河上游青海省西宁市和黄南藏族自治州等。）模型的和值比 模型的小，和调整 值比 模型的大，说明利用该模型研究黄河流域地级市碳排放量影响因素时空异质性优于 模型。由 模型得出：具有显著的空间异质性，土地城市化率的空间异质性从 年开始显现，人口规模和人口密度在地级市尺度上空间异质性较弱。）按回归系数大小排序：人口规模土地城市化率人口密度。、人口规模和土地城市化率呈显著正效应，人口密度呈显著负效应。空间分布上，的正效应强度从东南向西北逐渐增强，人口规模的正效应强度自东向西移动至中部后向南收缩，人口密度的负效应强度自西南向东北递增，土地城市化率的正效应强度自东南向西北递增。参考文献：任保平，张倩黄河流域高质量发展的战略设计及其支撑体系构建改革，（）：新华网习近平主持召开中央财经委员会第九次会议：？生态环境部 黄河流域生态环境保护规划：张国兴，王涵基于 模型的黄河流域中心城市水生态文明建设评价生态经济，（）：，王瑛，何艳芬中国省域二氧化碳排放的时空格局及影响因素世界地理研究，（）：王长建，汪菲，张虹鸥新疆能源消费碳排放过程及其影响因素：基于扩展的 恒等式生态学报，（）：公维凤，范振月，王传会，等黄河流域碳排放区域差异、成因及脱钩分析人民黄河，（）：袁伟彦，方柳莉，罗明中国工业碳排放驱动因素及其脱钩效应：基于时变参数 生产函数的分解和测算资源科学，（）：李天驹，贾磊，张晓慧中国农地碳排放时空分异及其驱动因素生态经济，（）：，韩梦瑶，刘卫东，谢漪甜，等中国省域碳排放的区域差异及脱钩趋势演变资源科学，（）：张永凯，田雨黄河流域城市群碳排放与经济增长脱钩状态及驱动因素人民黄河，（）：李云燕，张硕，中国城市碳排放强度时空演变与影响因素的时空异质性中国环境科学，（）：张赫，黄雅哲，王睿，等中国县域碳排放脱钩关系及其时空特征演变资源科学，（）：刘贤赵，高长春，张勇，等中国省域碳强度空间依赖格局及其影响因素的空间异质性研究 地理科学，（）：，：，：，（）：，：，（）：，（）：，人 民 黄 河 年第 期 ，：，：，（）：，：，：，：，（）：，：，（）：，（），（）：张华明，元鹏飞，朱治双黄河流域碳排放脱钩效应及减排路径资源科学，（）：武翠芳，赵慧，陈少炜黄河流域数字普惠金融与碳减排能力耦合协调分析（）：张杰，张剑勇，马蓝黄河流域生态保护与高质量发展耦合协调关系实证分析（）：，（）：，（），（）：【责任编辑 张华兴】（上接第 页）夏军，彭少明，王超，等气候变化对黄河水资源的影响及其适应性管理人民黄河，（）：，杨大文，张树磊，徐翔宇基于水热耦合平衡方程的黄河流域径流变化归因分析中国科学：技术科学，（）：黄建平，张国龙，于海鹏，等黄河流域近 年气候变化的时空特征水利学报，（）：刘昌明，王恺文，王冠，等 年黄河流域径流及其影响因素的变化分析人民黄河，（）：，周帅，王义民，畅建霞，等黄河流域干旱时空演变的空间格局研究水利学报，（）：王有恒，谭丹，韩兰英，等黄河流域气候变化研究综述中国沙漠，（）：黄婷婷，林青霞，吴志勇，等黄河流域干旱时空特征及其与 的关联性分析人民黄河，（）：高秉丽，巩杰，李焱，等基于 的黄河流域多尺度干湿特征分析干旱区研究，（）：张卓群，冯冬发，侯宇恒基于 函数的黄河流域干旱特征研究干旱区资源与环境，（）：，：，（）：，：，：，：宋松柏 函数在水文多变量分析计算中的问题人民黄河，（）：，【责任编辑 张智民】人 民 黄 河 年第 期