我国地域辽阔, 地貌复杂, 湖泊分布广泛、类型多样且演变历史复杂^[1-2]。高原湖泊位于高原地带, 具有对外无连通性、水体交换率低及自净能力差等特点^[3-4]。由于地理位置特殊, 高原湖泊具有极高的生态价值, 在维护高原地区生态系统及我国大江大河生态安全稳定中扮演着重要角色。与平原湖泊相比, 高原湖泊结构更加复杂, 生态功能更为重要, 具有一定的特殊性^[5-7]。我国高原湖泊主要分为青藏高原湖区、云贵高原湖区和蒙新高原湖区。在全球气候变暖的大背景下, 除青藏高原以外, 其他高原湖泊普遍产生了湖泊水位下降、湖面萎缩、湖水污染等问题, 并且这些问题所产生的影响已经从局部扩展为流域性乃至区域性^[8-10]。高原湖泊流域是内蒙古地区重要的水源区域, 在防汛抗旱、气候调节、生物多样性维持等方面发挥着重要作用^[3]。近几十年, 粗放的经济发展模式及人类生产生活的影响给内蒙古高原湖泊带来了巨大的破坏。据内蒙古高原湖泊演变历史的研究表明, 近几十年, 内蒙古高原湖泊面积减少了1223 km², 下降约22%, 加之湖水污染及周围生态系统的破坏产生了严重的生态环境问题^[11-12]。一方面, 流域内居民生存及产业发展越来越依赖流域生态系统；另一方面, 受到破坏的生态系统所能提供的价值却在减少, 已经严重威胁到流域可持续发展及资源可持续利用^[13]。目前对于高原湖泊的研究主要针对青藏高原和云贵高原地区^[14-20], 对内蒙古高原湖泊流域生态环境变化及生态保护对策仍需进一步深入研究, 进而减缓内蒙古高原湖泊流域生态环境退化趋势。

对目前国内外相关文献的梳理可以发现, 对湖泊及其流域进行综合治理与保护是未来研究的主要发展方向, 然而对自然条件较为复杂、生态环境更为脆弱及生态保护难度较大的内蒙古地区高原湖泊的研究依然不多^[21-25]。作为内蒙古四大湖泊之一, 达里湖是北方地区比较典型的高原湖泊, 其流域地处干旱/半干旱地区, 气候干燥, 降雨稀少。近年来, 在气候变化及人类活动双重影响下, 湖面萎缩、湖泊污染等现象严重, 流域草原退化、土地沙化等问题持续发生, 严重威胁着内蒙古地区的生态平衡^[26]。目前对达里湖流域的研究多集中在达里湖水域变化、水体污染等方面, 鲜有以流域为研究对象进行生态功能分区及保护措施的研究, 亟需针对流域特点及问题开展科学有序的保护措施研究^[27-30]。围绕“如何针对地理位置重要、生态环境脆弱的高原湖泊流域进行生态分区并实施有效地保护与管理”这一科学问题, 本文选取内蒙古高原湖泊达里湖流域为研究对象, 在明晰其生态环境状况及主要问题的基础上, 开展生态系统敏感性和生态系统功能重要性评价, 找出达里湖流域生态脆弱区域及生态保护重点区域, 在此基础上提出具有针对性的分区保护对策, 研究结果可以指导流域的综合规划、治理和保护, 并为流域的可持续发展提供科学依据, 在实现区域人与自然和谐发展具有重要的科学意义。

1 研究区概况及主要环境问题

达里湖是我国北方典型的干旱—半干旱地区草原性湖泊, 其流域位于内蒙古自治区赤峰市克什克腾旗西部, 地理坐标为116°22′—117°33′E, 42°33′—43°45′N(图 1)。流域面积约为4675 km²。流域属于温带半干旱大陆性季风气候, 多年平均气温约为2—4℃, 年降水量大约为300 mm, 年蒸发量是降水量的5倍多, 约为1600 mm。流域内景观类型多样, 动植物资源较为丰富。常住人口约为1.33万人, 以牧民为主, 且多数为少数民族。流域以草原农牧业和旅游业为主要支柱产业, 是经济发展欠发达地区^{[26, 31-32]}。流域内于1987年建立了达里诺尔自然保护区, 并于1997年经国务院批准为国家级自然保护区, 重点保护珍稀鸟类及其生境, 在流域内发挥了重要的生态系统功能^[33]。

长期以来, 在气候变化及人类活动的双重影响下, 流域的生态系统结构和功能都遭受到严重影响, 主要包含以下几个方面：

(1) 达里湖水质水量问题。由于入湖河流多次出现断流, 使得达里湖湖面面积萎缩明显, 从1980年的220 km²下降到2018年的197 km², 减少了近12%；储水量也从15亿m³下降到13亿m³。达里湖属于封闭湖泊, 由于湖水持续浓缩, 导致湖中各污染物浓度持续增加, 水质为劣V类。

(2) 草原退化问题。流域内草原生态系统受损明显, 退化草场面积达到1773 km², 占总面积的38%；盐渍化草原面积达到200 km², 占总面积的4%。草地的持续退化使得植物多样性及优质建群植物衰退明显, 多年生牧草由19种减少到8种, 一年生牧草由1种上升到7种, 草群高度由28 cm减少到5 cm。

(3) 水土流失问题。流域所在地为北方典型的风蚀水蚀交错区, 由于长期的侵蚀, 流域内水土流失面积高达2704 km², 占流域一半的面积, 局部地区甚至成为流动沙丘群。

(4) 土地沙化问题。流域位于干旱/半干旱生态系统脆弱区域, 降水量持续减少、气温持续升高, 再加上生态资源使用不合理, 使得流域内土地沙化面积多达1646 km², 严重影响了流域内生态环境的健康持续发展。

2 数据来源及研究方法 2.1 数据来源

通过前期与地方政府紧密合作, 开展了达里湖流域野外调研及数据资料收集工作, 获取了流域较为全面的基本情况信息。本文所选用的数据主要来源于以下途径：

(1) 土地利用/覆盖数据来源于中国科学院资源环境科学数据中心(http://www.resdc.cn)。

(2) 流域内降雨等气象数据来源于中国气象数据网(http://data.cma.cn/)及地方政府。

(3) 流域内土壤数据来源于国家地球系统科学数据共享服务平台土壤科学数据(http://soil.geodata.cn/)及地方政府。

(4) 流域内道路数据来源于全国地理信息资源目录服务系统(http://www.webmap.cn/)及地方政府。

(5) 流域内DEM(Digital Elevation Model, 数字高程模型)数据下载于地理空间数据云(https://www.gscloud.cn/search)上的ASTER GDEM V2数据源。

2.2 研究方法 2.2.1 生态系统敏感性评价方法

根据研究区的生态环境特征及状况, 并借鉴同类型湖泊流域生态系统敏感性分析已有案例^[34-37], 本文选取流域内容易导致生态危机的水土流失敏感性、土地沙化敏感性和土壤侵蚀敏感性三项指标对流域进行整体分析。参考相关文献, 3个指标的计算过程如下^{[35, 38-41]}：

(1) 水土流失敏感性：基于2015年遥感影像数据, 使用ArcGIS中的spatial analyst中的叠加分析, 对流域内沟谷密度、地形坡度、土壤可蚀性和土地利用类型进行加权叠加得到, 并将水土流失敏感性分为极敏感、高度敏感、中度敏感、轻度敏感和不敏感5个等级。

(2) 土地沙化敏感性。

(1)

式中, D_i为i评价区土壤沙化敏感性指数, I_i为干燥度指数, W_i为起沙风天数, K_i为土壤质地, C_i为植被覆盖度, 同样将结果分为5个敏感等级。

(3) 土壤侵蚀敏感性：选取植被覆盖度、地形起伏度、土壤质地和降雨侵蚀力R值使用栅格计算器叠加得到。

(2)

式中, SS_j为j空间单位土壤侵蚀敏感性指数, C_i为因子敏感性等级值, n为j空间单元所包含的敏感性因子指标个数, 并将结果分为5个敏感等级。

结合以上3个指标的评价结果, 通过使用ArcGIS平台的空间分析叠加模块, 对3个指标的分布图进行空间叠加, 最终得到达里湖流域生态系统敏感性分布图, 并将敏感性分为3个等级, 即一般敏感、敏感和极敏感。

2.2.2 生态系统功能重要性评价

依据研究区主要生态环境状况, 并充分参考同类型流域的相关研究文献^[42-46], 本文选取达里湖流域能提供的重要生态系统水源涵养、土壤保持和生境质量重要性3个指标对流域进行整体评价。根据相关参考文献, 3个指标的计算过程如下^[47-56]：

(1) 水源涵养重要性：在计算出产水量的基础上, 依据Bai等人提出的公式, 即各地类的水源涵养量等于各地类产水量减去各地类的地表径流量。本文选取流域内年均降水量、潜在蒸散量、植被可利用含水量、土壤深度、根系深度等数据, 通过公式3—5的计算, 得到达里湖流域水源涵养重要性, 并将结果分为极重要、高度重要、中等重要、较重要及一般重要5个等级。

(3)

(4)

(5)

式中, WR_ij是LUCC_j上像元i的年持水量, Runoff_ij是LUCC_j上像元i的年地表径流量, P_ij表示LUCC_j上像元i的年均降水量, C_j表示第j类土地利用类型的地表径流系数。

(2) 土壤保持重要性：选取降雨侵蚀力、土壤可蚀性、地表覆盖管理、土壤保持措施等因子, 通过InVEST模型中的RUSLE(Revised Universal Soil Loss Equation)模型对研究区土壤保持量进行估算得到, 并将结果分为5个等级。公式如下：

(6)

(7)

(8)

式中, A是流域内任意栅格x的实际土壤侵蚀量(t), A_p和A_c分别为流域内任意栅格x的潜在土壤侵蚀量和土壤保持量。R_x、K_x、C_x、P_x分别是栅格的降雨侵蚀力、土壤可蚀性、地表覆盖管理因子和土壤保持管理因子, L_x、S_x为坡长和坡度因子。

(3) 生境质量重要性: 选取生境退化度、生境质量指数、生态胁迫源、生境类型对胁迫源的敏感性和生态胁迫源因子图层等数据, 运用InVEST模型中的生境质量(Habitat quality)模块和ArcGIS 10.3软件处理进行分析得到, 结果同样分为5个等级。

结合以上3个指标的评价结果, 通过使用ArcGIS平台的空间分析叠加模块, 对3个指标的分布图进行空间叠加, 最终得到达里湖流域生态系统功能重要性分布图, 将流域按一般重要、重要和极重要等级划分。

2.2.3 分区方法

达里湖流域的分区遵循可持续发展原则、整体性原则、主导性原则与可操作性原则, 将生态系统敏感性及生态系统功能重要性的最终评价结果按照表 1的规则进行叠加分区, 将整个流域分为极重要敏感区、一般重要敏感区和低重要敏感区。

3 研究结果 3.1 生态系统敏感性评价

达里湖流域的3个生态系统敏感性指标评价结果如图 2所示, 水土流失、土地沙化和土壤侵蚀敏感性分布有一定差异。水土流失极敏感地及高敏感地区主要分布在流域的南部部分地区, 其余大部分地区为水土流失敏感性较低地区；土地沙化极敏感及高度敏感地区主要分布在流域的西北部地区及东北部零星地区；土壤侵蚀极敏感及高度敏感地区主要分布在流域的东部。

从流域生态系统敏感性分布图(图 2)来看, 生态极敏感地区主要分布在流域的西北部及东北和南部的零星地区, 面积大约为1577 km², 约占流域总面积的34%, 该区域土地覆盖类型主要为林地、耕地及部分草地。敏感地区主要分布在流域的中部, 零星分布在流域东北部及南部, 面积大约为2363 km², 约占流域总面积的51%。一般敏感区主要分布在流域的中部地区, 面积大约为732 km², 约占流域总面积的15%。

3.2 生态系统功能重要性评价

水源涵养、土壤保持和生物多样性维持是达里湖流域所能提供的重要生态系统功能, 对维持流域内生态环境及经济社会发展起到重要作用。通过图 3可以看出, 水源涵养中等重要区域主要分布在流域的北部, 较重要地区主要分布在流域的南部；土壤保持重要地区主要分布在流域的东北林地地区, 其余大部分地区为一般重要区域；生境质量极重要地区分布在流域内的湖区及东北地区, 中等重要地区分布在流域的北部, 较重要地区主要分布在流域的南部。

将水源涵养、土壤保持和生境质量3个结果叠加后, 可以得到达里湖流域生态系统功能重要性分布图(图 3)。如图 3所示, 达里湖流域生态系统功能极重要区域主要分布在达里湖区及东北部的林地区域, 零星分布在部分草地, 面积大约为964 km², 约占整个流域面积的20%。生态系统功能重要区域主要分布在流域北部的耕地及部分草地, 面积约为2912 km², 约占整个流域面积的62%。生态系统功能一般重要区域主要分布在流域南部的建设用地、未利用地与部分草地, 面积约为824.65 km², 约占流域总面积的18%。

3.3 流域分区结果

将生态系统敏感性分布与生态系统功能重要性分布结果进行空间叠加, 并按照极重要敏感区、一般重要敏感区及低重要敏感区划分, 并依据每个区域的生态环境问题与保护特点, 同时考虑到流域内达里诺尔国家级自然保护区的重要性, 将流域分为禁止开发区、重点保护区、质量提升区与潜在威胁区, 得到达里湖流域的最终分区图, 如图 4所示。

由于流域内达里诺尔国家级自然保护区具有的重要生态意义, 将其单独划分为禁止开发区进行保护与管理。禁止开发区位于流域中间部位, 面积约为1216 km², 为整个流域面积的26%。该区域湖水水质水量问题较为严重, 盗猎野生动物现象时有发生。重点保护区涵盖流域内生态极重要敏感区, 面积约为762 km², 为流域面积的16%。重点保护区涵盖了流域东北部的林地、北部部分草地及南部零星草地。该区域内林地资源开发无序, 且工作人员保护意识较低, 生态保护工作效率不高。质量提升区包括一般重要敏感区, 面积约为1599 km², 为流域面积的34%。该区域主要涵盖流域东北部草原地带、北部耕地及南部部分草地地区。质量提升区内草地退化现象严重, 耕地结构亟需调整, 存在农业废水现象。潜在威胁区包括低重要敏感区, 面积约为1099 km², 为流域面积的24%。该区域主要涵盖南部的建设用地与未利用地, 以及周边零星草地。潜在威胁区是流域内人类活动最活跃地区, 生态环境受到较大影响, 且有逐渐影响其他区域的可能。

3.4 分区保护对策

达里湖流域具有地理位置重要、生态文化浓厚、少数民族聚居等特征, 是一个以旅游业和草原农牧业为支柱产业的特殊经济区域。流域所在地是生态脆弱区, 也是经济欠发达地区。在对达里湖流域进行生态分区分析后, 以保护生态系统和提升生态功能为重点, 通过对达里湖流域内4个区域的特点及问题分析, 提出针对每个区域特色的保护对策, 为流域的可持续发展提供科学依据。

禁止开发区涵盖了达里湖自然保护区, 应建立生态、环境、资源一体化的管理模式。达里湖保护区内划分了5个核心区、3个缓冲区和1个实验区, 保护对策应在此基础上开展。核心区是保护区内的精华所在, 要严格禁止一切人类活动, 并采取联合执法严厉打击破坏自然资源、偷猎等现象发生。缓冲区位于核心区周围地带, 在核心区保护中起到缓冲作用。要充分利用智能化、现代化技术, 对缓冲区内的自然资源及人类活动进行实时监测, 实现对缓冲区内珍稀物种和人类活动的实时监管。其他区域为实验区, 在该区域内应该依据实际情况, 有计划、有目的的进行相关基础设施建设以及开展科研、教学、参观、考察等活动^[57]。

重点保护区内的林地、草地具有重要的生态价值, 并且是流域内敏感地区, 对外界干扰的抵御能力较低, 需要加强生态功能重要及敏感地区的生态保护建设, 包括加强天然林的保护以及生态林和经济林的建设, 以加速被毁林地修复以及林地生态环境的改善。要科学、合理地利用区域内的自然资源, 确保生态环境为当地居民提供足够的生存资源, 以实现人与自然和谐发展。

质量提升区以重要敏感区的部分草地、耕地为主, 该区域内可以在生态环境承载力允许的情况下, 适当的进行资源开发利用, 同时对有害的生产生活方式要绝对避免。对于该地区的草地资源, 要严格实施草蓄平衡政策, 指导牧民进行蓄群结构调整, 要尽可能的防止其退化、沙化及盐碱化, 并根据草地的分布特点及生态环境特征, 进行因地制宜的分类规划保护。在该区域的耕地布局中, 要大力发展节水型农业, 积极发展节水技术、更改农作物结构, 提高水的利用率和水的生产效率, 以保护流域内的水资源, 提升流域生态环境质量^[58]。

潜在威胁区以低敏感重要的建设用地与未利用地为主, 零星包含周边部分草地。与其他地区相比, 潜在威胁区的人类活动相对密集, 生态环境较容易遭受破坏, 且影响范围有向周边区域扩散的可能。在此区域内, 应在生态环境承载力允许的条件下, 加强控制和合理规划资源的开发利用, 适度进行基础设施建设, 并科学、合理的进行农业、工业发展, 以减少对生态环境造成的破坏。要鼓励推进矿工企业进行规模化生产, 对于效益低下、环境破坏强的小企业进行合并, 加强对高污染企业的监督, 避免发生环境破坏行为, 同时鼓励环境修复及保护的行为^[59]。坚持生态优先、绿色发展理念, 推进生态旅游产业的健康发展, 打造体现当地少数民族特色文化的旅游业, 在不破坏生态环境的前提下, 增加当地居民的收入^[60]。

4 结论

高原湖泊是我国极为宝贵的自然资源, 是影响高原地区生态环境状况的关键因素之一, 具有较高的生态价值。达里湖是内蒙古典型的干旱-半干旱草原性湖泊, 其流域是一个以草原农牧业和旅游业为支柱产业的特殊经济区域, 具有地广人稀、生态文化浓厚、地理位置重要、少数民族聚居等特征。近十几年来, 达里湖流域内生态功能退化严重, 使流域内经济、社会和生态系统矛盾进一步恶化, 保护与治理内蒙古高原湖泊及其流域已成为社会各界关注的重大问题。为了达到扭转内蒙古高原重要湖泊当前严重而紧迫的萎缩现状、缓解生态系统功能降低等问题、保持流域生态系统健康和完整的目标, 本文针对达里湖流域开展了保护对策研究, 基于生态系统敏感性及生态系统功能重要性评价结果, 将达里湖流域分为禁止开发区、重点保护区、质量提升区与潜在威胁区, 并根据每个区域的生态特性, 制定合理的保护对策。研究有利于丰富流域生态保护的基础理论与技术方法, 可为研究区及其他同类型湖泊流域的生态保护建设提出具有可操作性的建议及对策, 同时为管理者制定生态保护策略提供科学依据, 对恢复生态脆弱地区的环境及资源可持续利用具有指导作用。