生态空间是人类赖以生存、生产和生活的重要物质基础。随着社会经济发展对自然资源需求不断增长, 城市过快扩张造成生态用地占用, 生态空间与建设需求的矛盾愈发突显, 生态结构破坏、功能退化、环境恶劣等一系列生态安全问题频发^[1]。自然生态空间保护作为世纪可持续发展的新目标, 是国内外相关学术领域关注的热点和趋势之一。国外许多国家通过空间规划的编制来实现对生态空间的合理分区, 如德国于20世纪初率先将国土空间划分为不同区域并进行空间规划^[2], 欧洲、美洲和亚洲的各个国家紧随其后, 如美国在20世纪初采用土地使用分区制度, 其后针对不同时期的土地问题进行相应的措施制定^[3-7], 到21世纪颁布了《21世纪海洋蓝图》、《美国海洋行动计划》等法令, 对自然生态空间的保护重心向海洋资源倾斜, 促进了海洋资源的保护和可持续开发利用^[8-9]。在欧洲, 关于生态空间分区的研究逐渐趋向定量化, 如荷兰学者Janne Soininen^[10]通过社会空间结构的分析对生态空间进行定量化研究；法国学者Govaert Lynn^[11]以生态和进化对人口和社区的重要性作为生态进化分区的指标。我国对自然生态空间的分区研究始于对耕地的保护, 为提高耕地质量和促进土地资源可持续利用, 我国将土地分为农用地、建设用地和未利用地三种类型^[12-13], 之后在对森林资源^[14]、水资源^[15-16]以及海洋资源^[17]的保护中进行更加细致的分类管理。2017年《自然生态空间用途管制办法(试行)》(国土资发[2017]33号)的颁布, 为自然生态空间用途管制工作提供了科学的依据。自然生态空间合理的管制是我国生态文明建设的重要内容^[18], 目前国内对其研究的内容主要集中在生态空间的识别与结构优化^[19-22]、生态空间的失衡与重构^[23-24]、生态空间的演变与规划^[25-26]等方面, 已有的研究为自然生态空间的管理提供有力的理论基础和技术支持, 能为后续的研究提供较好的借鉴, 但大多处于宏观层面的考量, 较少考虑空间内部生态功能重要性和生态环境敏感性的空间特点。在合理判断自然生态空间生态系统服务功能和生态环境敏感的空间分布差异, 划分空间管制区域, 是实现区域生态社会可持续发展, 维护生态系统环境稳定发展的有效手段。

近年来, 随着快速城市化和工业化的发展, 江西省出现建设用地过快扩张、植被资源破坏、水资源质量下降和生态系统退化等生态问题, 生态环境受到威胁, 严重制约区域经济社会的可持续发展。江西省是我国第一批自然生态空间管制试点区域, 合理划分江西省自然生态空间格局, 是迫切解决当前经济发展与生态保护矛盾的有力保障。基于地理学空间分区方法^[27], 以江西省为研究对象, 基于江西省自然生态空间有效管理的现实需求, 在实地考量研究区生态维护现状情况的基础上^[28], 以生态系统服务功能重要性评价和生态敏感性评价为手段, 进行自然生态空间管制分区研究, 并提出管制建议, 试图解决生态保护与空间开发的合理分配, 以期为自然生态空间管制方案和相应管理措施的制定提供理论依据和技术支持。

1 研究区概况

江西省位于中国东南偏中部长江中下游南岸, 24°29′—30°04′N、113°34′—118°28′E, 与浙江、福建、广东、湖南、湖北、安徽毗邻, 区域总面积为16.69万km², 辖管11个区市。江西省东、西、南面环山, 中北部地势平坦, 形成一个狭长向北开口的盆地。研究区属亚热带季风性气候, 水热资源丰富, 多年平均降雨量1638 mm, 河流以鄱阳湖流域水系为主, 森林覆盖率高, 生物资源类型丰富。江西省作为全国生态大省, 于2017年开展自然生态空间用途管制试点工作。

2 数据与方法 2.1 数据来源与处理

本研究生态系统服务功能重要性和生态敏感性评价所使用的评价指标数据主要包括植被净初级生产力(Net Primary Production, NPP)数据、归一化植被指数(Normalized Difference Vegetation Index, NDVI)数据、DEM数据、气象数据、土地利用类型数据、地质灾害数据、土壤类型数据等。其中, DEM数据来源于地理空间数据云(http://www.gscloud.cn/), 空间分辨率为30 m；2015年NDVI数据来源于资源环境数据云平台(http://www.resdc.cn/)；气象数据为1981—2010年中国地面累年值数据集, 主要包括累年降雨量、累年气温数据, 来源于中国气象科学数据共享服务网；江西省2015年土地利用类型数据、2015年NPP数据与土壤类型数据来源于江西省自然资源厅。对空间数据进行统一的坐标投影, 并重采样至250 m×250 m分辨率的栅格。

数据具体处理方法如下：(1)通过对NPP元数据进行提取、统一坐标、栅格计算等处理, 得到NPP_mean值空间分布图。(2)通过对江西省87个气象站点的平均气温与平均降水数据进行空间克里格插值处理, 并将结果数据归一化得到F_pre的空间分布图, 通过公式计算得到R值。(3)通过ArcGIS 10.2空间分析工具对DEM数据进行坡度、地形起伏度、海拔运算并进行归一化处理, 得到F_slo、LS的空间分布图。(4)对土壤质地数据进行提取、分类, 进行归一化处理得到F_sic值, 并在查询相关文献^[29]和资料的基础上, 对土壤进行可蚀性赋值, 得到K值。(5)对NDVI数据进行处理并归一化得到植被覆盖因子C的空间分布图。(6)在谢高地等^[30]制定的生物多样性服务当量的基础上, 结合研究区实际情况, 得到江西省耕地、园地、林地、草地、水域、湿地等地类的生物多样性服务当量。通过已有研究结果可知, NDVI与生态系统服务存在着较好的正相关性^[31], 因此根据NDVI数据对生物多样性服务当量予以修正, 即计算不同景观类型的生物多样性服务当量和全年NDVI均值的乘积与不同景观类型NDVI的比值, 作为该用地类型的生境重要性指数S_H, 计算得到S_H的取值范围为0—4.51。(7)通过运用Conefor Sensinode 2.2、Conefor Inputs for ArcGIS 10.x以及ArcGIS 10.2软件计算得到研究区湖泊河流景观连通性指数S_R。

2.2 研究方法 2.2.1 研究技术路线

本研究的具体思路如图 1所示。(1)通过对研究区的实地考量, 收集相应的资料, 并运用ArcGIS等软件对数据进行处理, 建立空间属性数据库。为更加精确的表达研究区的空间差异性, 以250 m×250 m的栅格为基本评价单元, 进行生态系统服务功能重要性评价和生态环境敏感性评价, 根据评价结果进行级别分区。(2)根据土地利用类型数据进行自然生态空间识别, 按照《自然生态空间用途管制办法(试行)》(国土资发[2017]33号)中对自然生态空间的有关定义, 提取具有自然属性的土地利用类型, 得到研究区自然生态空间。(3)通过对自然生态空间中生态系统服务功能重要性和生态敏感性空间分布规律的识别分析, 构建二维关联判断矩阵, 得到自然生态空间管制分区结果。

2.2.2 评价方法与指标选取

根据《江西省主体功能区规划》、《江西省生态功能区划》、《自然生态空间用途管制办法(试行)》(国土资发[2017]33号)、《生态保护红线划定技术指南》, 确定江西省主导生态系统服务功能为水源涵养功能、水土保持功能和生物多样性保护功能, 主要的生态环境敏感性类型为水土流失敏感和地质灾害敏感, 生态重要性程度越高, 敏感性越强, 对空间开发的约束能力就越强, 越有利于生态空间的保护与管制。

评价指标因子的选取需要考虑指标的可获取性、稳定性和客观性, 并能充分突显出研究区生态系统服务功能和敏感特点的空间差异性。生态系统服务功能重要性评价通常采用模型评价法^[32], 但由于模型参数数量较多、要求精度较高且不适合大中尺度区域研究, 因此结合研究区域的生态环境实际情况和前人研究结果^[33], 确定水源涵养、水土保持功能适宜采用NPP定量指标评价法^[34], 生物多样性保护功能在参考已有研究的基础上^[35], 采用生物多样性保护综合指数法, 即以植被覆盖指数、生境重要性指数、河流湖泊景观连通性指数为指标进行综合评价。生态系统服务能力指数计算法具有数据参数较少, 且较易获取, 能够更加客观、科学的进行评价的优点。生态敏感性评价主要从水土流失敏感和地质灾害敏感方面, 选用降雨侵蚀因子、土壤可蚀性K值、地形起伏因子、植被覆盖因子、坡度因子以及地质灾害易发性等指标, 通过水土流失方程模型和多因素综合评价法进行相应评价。

由于无法比较3种主导服务功能和两种敏感指数对总体服务功能和生态环境综合敏感度的贡献比重, 因此不涉及具体的权重偏侧, 为保持研究的客观性、中立性, 本研究3种服务功能指数和两种敏感指数的权重一样。指标因子空间分布图(图 2)与评价方法(表 1)如下所示。

2.2.3 二维关联矩阵判断分析

通过对研究区进行生态系统服务功能重要性评价和生态环境敏感性评价, 进行聚类归并和趋同性分析, 将生态系统服务功能按其功能指数划分为极重要、重要、中等重要和一般重要4种类型, 生态环境敏感性则划分为极敏感、敏感、中度敏感和一般敏感4种类型。在对自然生态空间进行生态系统服务功能重要性和生态敏感性识别的基础上, 以此等级类型为基础, 按其高低顺序, 构建二维关联判断矩阵表(表 2)。在分析的过程中, 为了更好地进行研究区自然生态空间分区管制, 维护区域自然生态系统健康, 分区结果遵守“短板效应”原理, 将极重要、重要区域和极敏感、敏感区域确定为高重要区和高敏感区。

3 结果与分析 3.1 生态系统服务功能重要性 3.1.1 水源涵养功能重要性

水源涵养功能是指陆地生态系统如森林、草地、湿地等对水资源调节和供给功能, 主要表现为对地表径流的缓和、补充地下水、减缓河流季节波动等方面, 水源涵养功能是生态系统水文服务的核心^[36]。水源涵养功能重要性指数越高, 对生态系统服务功能的要求和价值就越高, 对建设开发的约束强度就越大。水源涵养功能重要性结果基本呈东部高、西部低的空间分布趋势(图 3), 其中, 高值区主要分布在景德镇市北部、上饶市西北部、鹰潭市东南部和抚州市东部等地区, 面积较小, 约为5341.33 km², 约占研究区总面积的3.20%, 这是由于该区域降雨量高、植被覆盖程度高, 水热均衡条件较好, 对水资源的调节和控制起着重要作用。低值区面积最大, 约为80152.00 km², 约占研究区总面积的48.02%, 主要分布在南昌市西部、九江市北部、吉安市南部和赣州市的大部分地区, 这是由于该区域经济建设活动较为频繁或由于地势较为险峻, 不利于涵养水源。

3.1.2 水土保持功能重要性

水土保持功能是指生态系统过程中, 通过其自身结构和过程的相互作用得以减少水流侵蚀导致的土壤侵蚀的功能^[37]。从评价结果可以看出, 高值区域分布较为零散, 面积约为16377.74 km², 约占研究区域总面积的9.81%, 主要分布在研究区的东北部、中部和东南部, 基本与高坡度走势一致, 坡度越大, 水土流失敏感性越明显。低值区面积约为51372.10 km², 约占研究区域总面积的30.78%, 主要分布在上饶市南部、吉安市东部和赣州市北部等地势较为平缓的区域(图 3)。

3.1.3 生物多样性保护功能重要性

生物多样性保护功能反映生态系统在维护物种、生态系统多样性等方面的能力, 生物多样性保护功能重要性评价是指研究区对生态系统和物种保护重要性的评价^[38]。从总体空间上看, 高值区分布面积最广, 约为86361.58 km², 约占研究区总面积的51.74%, 这是由于高值区域内多森林、湖泊分布, 植被覆盖指数的大小可以有效反映出生物量大小和植被种植特点。一般来说, 植被覆盖度指数越高, 生物活动越频繁, 生物多样性程度越高, 而湖泊作为自然生态空间中的重要景观, 能为生物提供良好的生存环境, 有利于动植物的迁徙、繁衍以及生物物质和能量的转换, 在生物多样性保护中具有不可替代的作用。低值区域面积最小, 约为12608.08 km², 约占研究区总面积的7.55%, 主要为经济建设活动频繁的城区、镇区或靠近人口聚集程度较高的地区(图 3)。

通过将水源涵养功能重要性指数、水土保持功能重要性指数和生物多样性保护功能重要性指数进行空间叠加, 得到研究区生态系统服务功能综合指数, 然后采用聚类分析的方法按其指数的高低划分为极重要、重要、中等重要和一般重要4个等级。总体上, 鄱阳湖核心湖体、景德镇市北部、上饶市北部、鹰潭市南部、抚州市、赣州市东南部等地区生态服务功能指数最高, 因此, 极重要区域分布面积最广。此外, 吉安市、萍乡市、新余市、宜春市和九江市的部分地区也有分布(图 3)。

3.2 生态环境敏感性 3.2.1 水土流失敏感性

反映研究区水土流失的敏感程度和空间分布特征, 采用水土流失方程模型进行综合评价, 结果表明：敏感性高值区主要分布在景德镇市、上饶市、鹰潭市、抚州市, 吉安市、宜春市和九江市有少量分布(图 4), 高值区面积约为16300.19 km², 约占研究区总面积的9.77%。敏感性低值区主要分布在研究区中西部地势较为平坦区域, 以南昌市中部、新余市东部、吉安市东部和赣州市西北部等区域为主, 面积约为55711.6 km², 约占研究区总面积的33.38%。

3.2.2 地质灾害敏感性

研究区地形多为丘陵区, 土壤以红壤为主, 土层较薄, 加之处于亚热带季风气候区, 易发生洪水、滑坡、泥石流等地质灾害。从地质灾害敏感性空间分布结果可以看出, 敏感性高值区分布较广(图 4), 主要分布在易发生洪涝灾害的鄱阳湖核心湖体区和其他坡度较大, 易发生地质灾害风险的地区, 高值区域面积约为71849.97 km², 约占研究区总面积的43.05%。敏感性低值区主要集中在研究区中北部鄱阳湖平原一带和吉安市、赣州市部分区域, 这些区域多为平原, 地势平缓且地质结构稳定, 多为无地质灾害风险区。

通过对水土流失敏感性和地质灾害敏感性的空间叠加分析, 得到生态环境综合敏感性评价结果, 按其生态敏感综合指数的高低划分为极敏感、敏感、中度敏感和一般敏感4个等级。从生态敏感性空间分布结果可以看出, 研究区除中部鄱阳湖平原和南部部分平原的敏感性较低, 极敏感地区基本集中在研究区四周地势起伏较大的区域, 其余地区则呈过渡状分布(图 4)。

3.3 自然生态空间管制分区

根据《自然生态空间用途管制办法(试行)》(国土资发[2017]33号)中对自然生态空间的相关定义, 对土地利用数据中具有自然属性的地类进行提取, 主要为森林、草地、湿地、河流、湖泊、滩涂等自然生态空间类型, 构成研究区自然生态空间(图 5), 面积约为117924.67 km², 约占研究区总面积的70.66%。通过ArcGIS 10.2软件中的Identify Task功能, 得到自然生态空间的生态服务功能重要性和生态敏感性分布结果(图 5), 并构建二维关联矩阵进行分析, 将自然生态空间划分为高重要高敏感、中度重要敏感和低重要低敏感3种类型区。

3.3.1 高重要高敏感区

指生态服务功能重要程度和生态敏感程度都极高的区域, 各市均有所分布, 主要集中在鄱阳湖核心湖体区以及研究区东部地区(图 5), 区域总面积43570.38 km², 占自然生态空间总面积的36.95%。本区生态系统服务功能状况极佳, 水源涵养、水土保持和生物多样性保护功能显著, 涵盖大量国家级自然保护区、一级水源涵养区和重要湿地, 且生态环境极其敏感, 是维护生态安全的核心区域, 因此, 该区内应严禁一切人为因素对自然生态的干扰, 禁止一切与生态保护和提升生态系统质量无关的人为活动^[39], 除保留必要的生态设施之外, 严格禁止与生态保护不符的建设用地占用, 严格控制农业产业活动, 实现“零污染、零排放”, 对区域内现有居民可实施生态移民工程, 建立相应的生态补偿机制, 减少人类活动对该区域的影响。该区未来发展的重点应着眼于加强林地保护、水域保护、湿地保护、水源涵养、生物多样性保护和实施严格的区域准入措施。

3.3.2 中度重要敏感区

指生态服务功能重要程度和生态敏感程度总体处于中等水平的区域。主要集中在研究区中南部, 以赣州市和吉安市为主(图 5), 此外, 基本围绕高重要高敏感区分布, 区域总面积为59316.3 km², 占自然生态空间总面积的50.30%。该区总体生态状况较好, 对高重要高敏感区生态系统起着不可替代的保障功能, 因此, 区域应以生态保护为主, 严格控制建设用地的占用规模和开发强度, 除必要的基础设施、民生工程和生态保护项目建设之外, 应禁止其他建设开发活动, 对区域内现有的居民点, 应明确生产、生活活动范围, 禁止居民点的开发再建。该区未来应以优化森林生态系统结构, 提高林地生态系统功能, 开展流域综合治理、加强生态保护, 依托区域现有生态资源优势, 合理开展以维护、改善生态系统服务功能为主要目的的生态经营活动。

3.3.3 低重要低敏感区

指生态系统服务功能重要程度和生态敏感程度均较低的区域。该区主要分布在自然生态空间靠近建设空间的边缘地带(图 5), 区域总面积为15037.99 km², 占自然生态空间总面积的12.75%。该区与人类活动区域联系较为密切, 受到经济建设活动影响较大, 原则上禁止建设工业项目和与生态保护相悖的建设项目, 允许在生态系统结构和功能不受破坏的前提下, 适度开展一定规模农业生产建设活动, 减少污染排放, 实现生态农业产业集约发展, 增强区域生态功能。

4 讨论与结论 4.1 讨论

目前, 关于自然生态空间涵义的认知主要是依据自然资源部下发的《自然生态空间用途管制办法(试行)》(国土资发[2017]33号), 认为自然生态空间是具有自然属性, 能够提供生态产品或以生态功能为主导功能的空间^[40]。实施自然生态空间管制的目的在于合理保护自然生态资源, 改善生态环境, 实现资源可持续利用, 缓解建设空间过快扩张过程中与生态保护之间的矛盾, 促进城市的可持续发展。因此, 自然生态空间管制分区可以作为生态保护规划、土地利用规划和城市发展规划的有力支撑, 确保研究区在经济转型时期生态环境健康发展。

总体上, 江西省自然生态空间整体生态服务功能重要程度与生态敏感程度较高, 三类管制区中, 森林占比最高(表 3), 因此, 均不适宜大规模的开发建设, 根据不同分区类型, 因地制宜的制定相应的管理措施, 同时加强生态修复和生态保护, 维护自然生态系统结构完整性和功能稳定性。将江西省生态环境厅公布的16个国家级自然保护区与管制分区结果进行空间叠加, 可以发现, 自然保护区基本落入高重要高敏感区域范围(图 6), 由此可见, 自然生态空间管制分区结果较为可靠。同时, 《江西省土地利用总体规划(2006—2020)》要求重点加强国家自然保护区的保护, 加强对生态敏感脆弱区域的保护, 自然生态空间管制分区与江西省土地利用总体规划的相关布局相辅相成, 也进一步验证了自然生态空间分区管制的合理性。研究可为江西省实施生态规划提供相应的理论基础, 有利于建设开发和生态保护之间的协调发展, 也是省级尺度自然生态空间管制分区的有益探索, 但从分区方法看, 还存在许多不足, 自然生态空间是一个复杂的综合体, 合理科学的分区与影响因素选择息息相关, 由于受知识体系和指标可获取性的限制, 评价指标体系有待完善, 特别是在生态系统服务功能评价中还应考虑生态服务产品供给功能和文化服务功能的指标选取, 使生态系统服务功能更加全面, 此外, 有关管制分区的方法还须进一步深入探索与完善。

4.2 结论

自然生态空间管制分区着眼于自然生态环境的改善和自然资源的合理利用, 通过对区域进行级别分区, 为区域生态系统服务和生态环境提供相应的保障。因此, 对自然生态空间进行相应的生态服务功能重要性和生态敏感性评价, 是实现管制分区的重要手段, 但级别分区也并非仅仅是自然生态空间的一种管制方式, 而是在分区的基础上, 关注不同类型区域生态环境问题的发生作用机制, 以维护生态系统结构和功能稳定为原则, 通过人为干预等手段来寻求区域生态环境可持续健康发展的空间途径。

本研究以江西省为研究对象, 运用GIS空间技术手段, 以栅格为评价单元, 进行生态系统服务功能重要性和生态敏感性识别下的自然生态空间管制分区研究。研究结果表明：江西省自然生态空间总面积为117924.67 km², 约占全省总面积的70.66%, 从空间上可划分为高重要高敏感区、中度重要敏感区和低重要低敏感区3种类型区；其中, 高重要高敏感区以生态保护和生态修复为主, 实施严格的区域准入措施；中度重要敏感区可依托区域生态资源优势, 合理开展以维护、改善生态系统服务功能为主要目的的生态经营活动；低重要低敏感区允许在不破坏生态系统结构和功能的前提下, 适度开展一定规模的生产建设活动, 减少污染排放, 增强区域生态功能。基于生态系统服务功能和生态敏感性评价的自然生态空间管制分区相较于传统方法具有更强的客观性和科学性, 并且, 江西省16个国家自然保护区基本落入高重要高敏感区域范围内, 也充分说明基于生态系统服务功能重要性与生态敏感性自然生态空间管制分区方法的适用性。自然生态空间管制分区是实现自然生态空间用途管制的新探索, 然而, 理论能否运用在实践工作中, 能否产生预期效应, 还须经过实践的验证。分区管制的实现不仅需要有关部门进行组织、管理, 还需要建立监督反馈机制, 如何使分区管制得到有效的实施并产生预期的效果, 是研究的重点。