国土空间作为人类和社会经济活动的主要空间载体, 对于维持粮食安全、生态系统健康、土地资源优化配置、社会协调发展等起着举足轻重的作用。但在社会经济和工业化的飞速发展进程中, 建设空间肆意扩张、农业用地利用强度增大、景观破碎化程度增加等现象日益凸显, 致使区域生境质量的抗干扰能力下降, 进而对农业产出效率和居民生活质量造成一定影响, 国土空间生态质量的优劣逐渐成为制约人与社会和谐发展、资源可持续利用的主要因素^[1-3]。与此同时, 国家出台有关政策明确开展“五级三类”体系的国土空间规划工作, 并将国土空间生态修复提升为国家战略工程, 先后多次提及生态文明建设方针政策^[4-5]。在此背景下进行国土空间土地生态质量调查与综合评价研究, 能够精准地识别出国土空间土地生态优劣差异情况, 可为后续开展国土空间规划、生态修复等提供有针对性的科学参考依据, 对于促进生态-社会-经济和谐发展具有十分重要的现实意义。

20世纪40年代末美国环保主义者奥尔多· 利奥伯德率先提出“土地健康”概念^[6], 随着环境管理和资源保护意识的提高, 土地生态评价研究得到广泛关注^[7-8]。目前国内外关于土地质量的研究多集中于土地生态安全评价^[8-12]、生态风险评价^[13-14]、生态脆弱性评价^[15-18]、生态敏感性评价^[19]、生态系统健康评价^[20]、生态适宜性评价^[21]等研究角度, 运用压力-状态-响应模型^{[10, 22-23]}、改进TOPSIS(Technique for Order Preference by Similarity to Ideal Solution)模型^[24]、多目标综合评价(Multi-objective comprehensive evaluation, MCE)模型^{[8, 18]}等研究方法, 融合数理统计、GIS空间分析技术和景观生态学^[25-26]等学科知识, 从不同的研究尺度开展土地生态环境定量评价、空间差异以及内外在影响机制等的研究。这些现有研究多集中于土地本底条件和社会经济属性的综合评估, 研究视角较为单一, 缺乏以国土空间综合特性为出发点的土地生态质量评估实例研究。

国土空间生态质量的优劣受到自然和社会等多方面因素的制约, 因此需要构建多维度多层次的综合评估体系。多目标综合评价法能够围绕研究目的构建多准则综合评价体系, 定量测算出准确的评价结果, 原理易懂, 操作简便, 现已广泛运用于资源承载力评价、土地开发利用适宜评价等多个研究领域。本研究结合国土空间规划思想, 以饶河流域涉及的县市为研究区域, 着眼于该区域城镇空间、农业空间和生态空间自身的结构、功能和效应, 从自然基础属性、生态结构特征、生态胁迫压力与生态效益水平四个层面构建国土空间生态质量综合评价指标体系, 基于多种数据分析软件, 结合多目标综合评价法对饶河流域国土空间土地生态质量进行综合评估, 探究饶河流域土地生态质量现状及空间分异特征, 并定量诊断影响区域土地生态质量的障碍因素, 综合评价和障碍度诊断结果提出提升研究区生态质量的改进措施, 以期为研究区国土空间土地整治、空间规划、生态修复等工程项目的实施提供切实可行、科学合理的参考建议, 促进研究区人-自然-社会三者的和谐发展, 为进一步优化研究区国土空间生态安全格局打下良好基础。

1 研究区概况

饶河流域为鄱阳湖五大流域之一, 地处江西省东北部, 涉及鄱阳县、万年县、婺源县、弋阳县、德兴市、景德镇市等地区, 本研究以饶河流域所涉县市的127个乡镇为研究区域(图 1), 该区域国土总面积为1716064.12 hm², 截止2018年末, 研究区总人口483.85万人, 地区生产总值达1589.87亿元。境内主要地形为山地和丘陵, 约占70%, 整体地势东北部高, 西南部及中部地势较低, 土壤类型以红壤和水稻土为主, 各占比68.85%和23.22%。境内光照充足, 年均降水量1070.47 mm, 年均气温19.7℃, 属于亚热带季风性湿润气候区。国家高度重视长江经济带的生态环境保护, 为此江西省相继开展了山水林田湖草生态保护修复、森林湿地湖泊生态保护和修复、矿山地质环境修复治理等试点工程, 研究区各县市积极响应, 水土流失面积持续减少, 矿区水质得到明显改善。工业化和社会经济的飞速发展, 以及国土空间规划浪潮的来袭, 使得研究区国土空间面临着生态安全、资源利用效率以及人与自然和谐共生等多重考验。

2 数据来源与研究方法 2.1 数据来源与预处理

研究数据主要包括矢量数据、栅格数据以及统计数据三类。2018年土地利用现状矢量数据来自于六个县市自然资源局以及江西省土地整理中心提供, 按照土地利用代码表, 基于ArcGIS 10.2软件重新划分为农用地、建设用地、水域湿地、其他土地四种土地利用类型, 参考相关文献进一步划分出生产用地、生活用地和生态用地三生空间用地类型。栅格数据包括遥感影像和DEM数据, 均是从地理空间数据云(http://www.gscloud.cn/)下载得来, 运用ENVI 5.3软件进行图像预处理操作。统计数据包括社会经济数据、气象数据和土壤样点检测数据, 社会经济数据来源于2018年六个县市统计年鉴以及国民发展与社会经济统计公报；气象数据来源于江西省气象局提供的2018年气象数据, 基于ArcGIS 10.2软件进行空间插值得到气象基础数据；土壤样点数据来源于2018年江西省耕地质量等别数据, 共有528个土壤样点, 参照《农用地土壤污染风险管控标准》和内梅罗综合污染指数法, 基于ArcGIS 10.2软件对数据进行空间插值, 得到土壤污染、土壤质地等基础数据。

2.2 指标体系构建与权重确定 2.2.1 评价指标体系构建

评价指标的选取应遵循科学性、综合性、代表性、可获得性等原则, 本研究以农业空间、城镇空间和生态空间为研究视角, 着眼于这三类空间自身的结构、功能和效应, 从自然基础、生态结构、生态胁迫和生态效益四个方面选取25个指标, 构建研究区国土空间生态质量评价指标体系(表 1)。

2.2.2 评价单元及指标权重确定

运用ArcGIS 10.2的创建渔网工具将研究区行政边界划分为3 km×3 km格网作为评价单元, 共计1905个。土地利用、水文、土地退化/恢复等矢量因子数据, 通过Fragstats 4.2的景观格局指数计算以及ArcGIS 10.2的提取分析、转换、分区统计等操作得到相应数值。

指标赋权的方法有层次分析法、德尔菲法、熵权法等, 本研究结合层次分析法和德尔菲法来确定指标权重, 由于指标具有不同的性质和量级, 在确定权重之前首先参考极差归一化法对指标数据进行归一化处理, 公式详见参考文献^[8], 计算得到每个评价指标的权重值(表 1)。

2.3 国土空间生态质量综合评价

本研究采取多目标综合评价法, 对饶河流域国土空间生态质量进行综合评估, 公式如下：

(1)

式中, TEI为国土空间生态质量综合指数, y_ij为指标归一化后的数值, W_j为第j项指标的权重。

2.4 空间统计分析

探索性空间数据统计分析可以用来描述特定位置现象或事物属性值空间分布格局的相关性特征及规律, 通常采用Moran ′s I指数来进行全局和布局的空间自相关分析, 指数值的范围在-1—1之间, 绝对值越高, 空间相关性越大, 大于0表示呈现空间正相关, 小于0表示呈现空间负相关, 采用G_i*来探究属性值的空间分布格局特征(即冷热点分析)^[27-28]。本研究运用Geoda1.12软件进行空间权重和单变量局部Moran ′s I指数的计算(公式见参考文献^[27]), 并绘制局部空间自相关(Local indicators of Spatial association, LISA) 聚集图, 来分析饶河流域国土空间生态质量评价结果的高值区和低值区的空间分异规律, 聚类结果中“高-高”、“低-低”区域表示评价单元自身与周围的生态质量都较高或较低, 差异较小, 具有较高的空间相关性, “高-低”、“低-高”区域表示自身评价单元与周围生态质量差异较大, 空间相关性较低。运用ArcGIS 10.2软件的空间统计分析工具, 分析四个准则层综合指数的高值和低值区(即冷热点区)格局分布特征。

2.5 国土空间生态质量分区

本研究运用SPSS 24.0软件对评价单元的生态质量评价结果进行K均值聚类分析, 将研究区国土空间生态质量划分为5个等级, 并结合ArcGIS 10.2将结果空间化展示, 划分标准如表 2所示。

2.6 障碍度诊断

定量识别影响国土空间生态质量的主要因素, 可为制定精准的国土空间生态建设管控措施提供科学依据。本研究引入因子贡献度F_j、指标偏离度I_j、障碍度(O_j、U_i)3个指标进行分析诊断, 公式详见参考文献^[29-30]。

3 结果分析 3.1 国土空间生态质量空间分异 3.1.1 准则层评价结果空间分异

研究区国土空间生态质量准则层的评价结果如图 2所示。自然基础综合分值处于0.1017—0.2811之间, 呈现东北西南高、西北东南低的分布态势。其中万年县东部、浮梁县北部、婺源县北部和东北部、鄱阳县南部地区自然条件良好, 基础分值较高；生态结构综合分值处于0.1018—0.1588之间, 高值区主要分布于婺源县、浮梁县、鄱阳县西部、弋阳县北部地区, 总体呈现环绕分布, 这些地区生态用地比例较大, 生态环境优良, 景观类型较为丰富；生态胁迫综合分值处于0.1125—0.2622之间, 乐平市、昌江区、珠山区和浮梁县的分值较高, 东北部的婺源县、南部的弋阳县、万年县胁迫分值相对较低, 主要是受土壤重金属污染以及社会压力的影响, 其中昌江区、珠山区和乐平市主要源于人口密度和社会城镇化率较大, 浮梁县生态压力主要源于草地退化；生态效益综合分值处于0.0085—0.0511之间, 相比于其他三个准则层指标, 该准则层分值偏低, 呈现东北高西南低状分布, 浮梁县和婺源县整体较高, 该区域生态建设和生态效益状况较好, 生态保护意识和自然恢复能力较强, 德兴市、弋阳县、万年县次之, 而昌江区较弱, 主城区需加强粮食安全保障, 增加林草地建设工程, 提高人们的生态保护意识。

3.1.2 生态质量综合指数空间分异

研究区国土空间生态质量综合分值空间分布如图 3所示。从整体来看, 东北部地区的国土空间生态质量高于西南部地区, 生态质量综合指数范围在0.3939—0.6489之间, 均值为0.5425, 最高值栅格位于浮梁县瑶里镇, 最低值栅格位于乐平市清湖乡。其中浮梁县瑶里镇、黄坛乡、鹅湖镇分值较高, 生态质量综合指数均值分别为0.6210、0.6037、0.6001, 婺源县的沱川乡、德兴市的昄大乡的国土空间生态质量分值仅次于浮梁县, 均值分别为0.5954、0.5877。低值区主要集中在弋阳县的清湖乡、乐平市的后港镇和主城区, 均值分别为0.4409、0.4594、0.4602。

3.2 生态质量综合指数空间统计分析 3.2.1 空间自相关分析

基于Geoda 1.12软件的空间分析功能, 分别在Queen邻接、Rook邻接和K-近邻三种不同的空间权重下, 研究区国土空间生态质量四个准则层评价结果的Moran′s I指数如表 3所示。三种空间权重矩阵下, 自然基础、生态结构、生态胁迫、生态效益四个准则层以及生态质量综合分值的Moran′s I指数均处于0.8—1之间, 表明准则层和综合评价结果在空间上均存在较高的正相关性, 具有明显的聚集性特征。准则层Moran′s I指数值由高到低依次是生态效益(0.9400)＞自然基础(0.9096)＞生态胁迫(0.8750)＞生态结构(0.8542), 生态质量综合评价结果的Moran′s I指数为0.9311。

研究区国土空间生态质量综合指数局部空间自相关LISA聚集图如图 4所示, 图中仅出现“高-高”、“低-低”两种聚类结果, 有25.77%的区域为“高-高”聚集的高值区, 总面积为442229.72 hm², 其中近八成分布于浮梁县境内, 还有少量分布于婺源县北部、德兴市东部和万年县南部；“低-低”聚集的低值区占比20.31%, 总面积为348532.62 hm², 主要分布于乐平市、弋阳县以及鄱阳县部分地区。

3.2.2 冷热点分析

准则层和综合评价分值变化冷热点分布结果如表 4、图 5和图 6所示。结果显示准则层和综合评价的冷点和热点区域面积基本持平, 自然基础和生态结构准则层的热点区域面积多于冷点区域面积, 生态胁迫和生态效益准则层与之相反。自然基础准则层的热点区域栅格数为460个, 占比24.14%, 主要集中分布于万年县中南部、婺源县东部和北部地区、浮梁县东部和西部部分地区, 冷点区域栅格数为383个, 占比20.08%, 主要集中分布于乐平市境内；生态结构准则层的热点区域相对于其他准则层占比最大, 占比39.37%, 主要集中分布与婺源县和浮梁县境内；生态胁迫准则层冷热点面积占比最小, 且分布较为分散, 热点地区占比仅为11.55%, 主要分布于昌江区的鲇鱼山镇、万年县的苏桥镇以及浮梁县, 冷点区域分布于弋阳县、万年县的齐埠镇、德兴市和婺源县的部分地区；生态效益准则层冷热点分布相对来说最为集中, 且占比最大, 34.44%为热点区域, 主要集中分布于研究区东北部的浮梁县、婺源县以及德兴市部分地区, 冷点地区占比37.28%, 主要集中分布于研究区西部的鄱阳县、景德镇市和万年县。饶河流域国土空间生态质量综合评价结果的冷热点区域共1341962.14 hm², 其中热点区域占比23.04%, 近六成分布于浮梁县境内, 冷点区域占比55.16%, 主要集中分布于乐平市和弋阳县。

3.3 国土空间生态分区及治理措施

研究区国土空间生态质量分区结果以及障碍度指数测算结果如表 5、表 6、表 7和图 7所示。研究区国土空间生态质量综合指数均值为0.5425, 其中Ⅲ级生态分区面积最大, 占比26.19%, Ⅴ级分区面积最小, 仅占8.71%, 其他等级区占比分别为Ⅰ级22.26%、Ⅱ级22.15%、Ⅳ级20.68%。经测算, 各分区中坡度、有效土层厚度、土壤有机质含量、草地退化率、GDP增长率、景观格局指数指标的障碍度指数较高, 均值分别为13.40、11.31、9.71、8.56、7.64、5.58, 这六个因素是制约饶河流域国土空间生态质量主导因素；自然基础准则层对研究区国土空间生态质量的障碍度都远远高于其他三者, 障碍度指数均值为43.73, 生态结构与生态胁迫对研究区国土空间生态质量的障碍度相近, 生态效益的障碍度最小。

I级生态分区包含424个评价单元, 总面积为381995.87 hm², 占比22.26%, 生态质量综合指数均值为0.5972, 其中66.51%分布于浮梁县境内, 婺源县东北部、德兴市东北部也有少量分布, 占比分别为19.81%、9.67%, 生态质量综合指数最高的栅格位于浮梁县瑶里镇, 分值达到0.6489。该区域拥有植被覆盖率高、规模大、景观斑块连接度好的生态用地, 整体生态质量水平最优, 在涵养水源、水土保持等方面有着很高的生态服务价值, 区域内拥有得天独厚的自然风光, 且已经拥有国内著名的旅游景点, 生态管控措施相对比较完善, 人为活动对于生态环境质量的干扰度较小, 交通便利。可通过招商引资来进一步扩大旅游影响力, 打造知名旅游品牌, 在保障该区域生态环境质量的同时, 保证其生态空间格局, 并通过提供特色生态产品和发展现代生态旅游来提升生态经济效益水平, 严禁毁林开垦、建设污染企业等破坏自然生态质量的工程, 促进一二三产绿色发展。

II级质量生态分区包含422个单元, 总面积为380108.20 hm², 占比22.15%, 生态质量综合指数均值为0.5666, 主要分布于婺源县中部和东北部、德兴市东北部和西南部、鄱阳县南部和北部、万年县南部和东部, 占比分别为37.91%、20.62%、17.54%、10.43%。该区域的森林覆盖面积较广, 生境质量较好, 社会经济压力小, 斑块多以农用地为主, 但该区域生态景观布局紊乱, 综合生态效益水平较低。该区域应注意提升生态空间景观格局, 完善生态廊道建设, 加大林地的养护治理工作, 注重多类型生态用地相互连接的区域保护, 丰富区域的生态结构, 根据实际情况建设生态养护工程, 推广健康农业种植技术, 积极响应乡村振兴战略, 挖掘和发展当地特色农业来推动当地经济发展。

III级生态分区包含499个评价单元, 总面积为449437.19 hm², 所占比重最大, 为26.19%, 生态质量综合指数均值为0.5348, 主要分布于鄱阳县东北和西南、德兴市东部、婺源县西南部、万年县西北部, 占比分别为43.89%、15.83%、13.43%、10.22%；IV级分区包含394个评价单元, 生态质量综合指数均值为0.5019, 主要分布于鄱阳县中部和西北部、乐平市北部和南部、弋阳县北部和南部。III级和IV级分区主要位于城镇集中建设区周围, 自然基础状况欠佳, 景观破碎度较高, 人类活动和工业化影响较为强烈, 生态用地斑块相对较少。该地区在追求社会化发展的同时, 应遵循自然环境发展的客观规律, 因地制宜提出相应的管制措施, 在生态文明相关政策理念的指导下, 积极开展土地整治、生态建设等工程, 确保城镇-农业-生态空间和谐统一发展, 努力提升生态空间效益, 营造人民安居乐业的城镇空间氛围, 打造产出高效的绿色农业生产格局。

V级分区的占比最小, 仅占总面积的8.72%, 生态质量综合指数均值为0.4617, 主要分布于乐平市中部和西部、弋阳县西南部, 其中弋阳县清湖乡分值最低。该区域主要分布于主城区及工业园区周边, 建设用地的增加导致土壤环境和植被覆盖率的下降, 生态用地比例严重不足, 人类活动对自然景观和生态环境的干扰极大。该区域应加快生态建设的步伐, 以节约资源、保护环境为前提调整产业结构, 加强工业三废治理力度, 制定生态管控措施, 加快退化植被的修复工作, 积极统筹并开展山水林田湖草综合治理工程, 增加绿地公园、湿地公园等自然生态景观建设, 确保区域粮食安全, 建设生态宜居的城镇生活环境, 逐步提高区域生态效益水平, 协调好人类经济社会发展的需要与生态环境保护之间的关系。

4 结论与讨论 4.1 结论

基于多种空间分析软件, 着眼于农业、城镇、生态三类空间的现实需求, 选取自然基础属性、生态结构、生态胁迫和生态效益作为四个准则层, 构建饶河流域国土空间生态质量评价指标体系, 对饶河流域国土空间生态质量状况进行综合评估及空间分异研究, 结论如下：

(1) 研究区东北部丘陵地区的国土空间生态质量高于西南部平原地区, 其中浮梁县瑶里镇分值较高, 乐平市和弋阳县部分地区的分值较低。准则层分值和生态质量综合指数分布均具有高度空间自相关性, 生态效益准则层的冷热点分布空间聚集性最大。

(2) 研究区国土空间生态质量I级分区生态质量均值为0.5972, 占比22.26%, 有近七成分布于浮梁县；V级生态分区面积最少, 仅占8.71%, 生态质量均值为0.4617, 主要分布于乐平市和弋阳县境内, 工业化的飞速发展导致开发建设强度相对较大, 附近地区受到不同程度的辐射影响, 植被受到破坏, 土壤环境污染程度相对较高, 此区域的生态问题最为突出。

(3) 结合综合评价和障碍度测算结果, 研究区应依据地形地势、土壤等自然基础条件, 调整社会经济结构, 加大建设生态项目工程力度, 在生态质量状况急需改善的区域, 应以绿色发展为前提推进工业健康转型, 统筹并落实生态修复工程的实施, 合理调整土地用途, 强化耕地保护, 确保农业空间的高效绿色运作, 增加自然景观建造, 优化城镇空间人居环境, 整体提高饶河流域的生态效益水平。

4.2 讨论

健康的土地生态系统是维持经济社会稳步发展的环境基础, 也是国土空间规划项目落地实施首要考虑的环境因素, 生态质量的优劣受到内在本底属性和外部因素的多重影响。本研究在系统地收集和分析研究区自然基础条件和社会经济建设资料的基础上, 立足农业、城镇和生态三类空间的结构、功能、效应和胁迫特征, 选取土地利用/覆被、气候、地形坡度、土壤、生产能力、社会压力等因子, 构建国土空间生态质量综合评价指标体系, 探究研究区国土空间生态质量空间分异特征及主导障碍因素, 针对不同生态分区的现实情况提出不同的质量提升策略, 研究视角较有新意, 研究结果符合实际, 可为研究区下一步开展国土空间规划、生态修复、山水林田湖草综合治理等工作提供基础数据支撑。但仍存在以下不足之处：(1)指标选取稍有欠缺, 除了社会压力外, 频发的自然灾害也对生态质量带来威胁, 但由于数据获取存在困难, 没有将气象灾害、水土保持等自然压力纳入其中, 仅考虑了社会压力, 所使用的社会经济数据多以县域为单元, 在一定程度上削弱了评价结果的区域差异性, 并且在选取指标时没有进行相关性分析, 多是参考相关研究文献而来, 在后续的研究中需注重评价指标的代表性、差异性和综合性。(2)评价模型方面, 本研究使用层次分析法和德尔菲法计算指标权重存在较强的主观性, 后续研究应采取主客观结合的赋权方法, 来降低单一主观性带来的偏差, 可结合多期综合评价结果进一步运用预测模型, 对区域国土空间生态质量进行模拟预测。(3)本研究未涉及研究结果应用方面内容, 后续研究将在完善评价指标体系的基础之上, 进行多期数据的综合评价, 并预测未来国土空间生态质量趋势, 结合区域规划、生态修复等理论方法, 为研究区实施国土空间规划、生态修复分区等项目提供科学可靠的数据支撑和建议。