湿地是陆地与水体系统相互作用形成的特殊的自然综合体^[1-2], 具有生产、生态和水文等多项功能^[3], 尤其是对于保护生物多样性具有重要作用。在全球变化的大背景下, 随着人口的增加和经济的发展, 湿地系统面临的退化和萎缩的风险与日俱增^[4]。大面积分布于世界第三极-青藏高原的高寒湿地, 是在高寒气候背景条件下形成的包含湖泊、河流、沼泽和水库等特殊湿地类型^[5], 被称为“高原之肾”。青藏高原高寒湿地不仅是重要的生物多样性保护区域, 也是中国乃至亚洲生态安全屏障区。青藏高原高寒湿地由于其地理区位的特殊性和独特的气候土壤特点, 使其在全球变化研究中占有特殊的重要地位^[6-10]。

拉萨河流域是藏民族的主要聚居区, 人口密度相对较大, 尤其是拉萨市辖区范围内, 农业、建筑业和旅游业相对发达, 水资源消耗和污染问题不容忽视。麦地卡湿地是拉萨河的源头, 也是黑颈鹤、赤麻鸭等珍稀水禽的栖息地, 在维护拉萨河流域气候稳定、调蓄水量、净化水质和保护水源地等方面具有极其重要的作用^[11]。随着气候变暖和人类活动的加剧, 麦地卡湿地出现了湿地退化、草场沙化等生态问题, 严重威胁当地生态安全和生物多样性。但是目前关于麦地卡湿地的研究较少, 已发表文献仅研究了麦地卡湿地保护区内种子植物资源多样性、环境保护现状以及湿地生态系统经济价值等内容, 对湿地保护区景观格局和生态系统服务价值尚未有深入研究。因此, 本文希望利用多期遥感数据获得麦地卡湿地保护区土地覆被变化情况, 利用修正的生态服务价值评估模型计算湿地保护区生态系统服务价值, 分析气候因素对麦地卡湿地的影响, 以期对湿地保护提出有效措施, 从而加强湿地保护区建设, 维护高原生态系统安全。

1 研究区概况与数据源 1.1 研究区概况

麦地卡湿地保护区位于西藏那曲地区嘉黎县北部措拉乡境内(92°45′55″—93°19′25″E, 30°51′04″—31°09′44″N), 距拉萨东北方向230 km, 那曲县东南方向88 km(图 1), 保护区面积89237 hm², 平均海拔4900 m。气候属高原亚寒带半湿润气候, 年均气温-1.7—0.7℃, 年均降水量700 mm, 干湿季节分明^[11]。地貌属高原湖盆谷地平原, 多湖盆、谷地、山地, 区内分布着240多个面积不等的湖泊, 其中面积最大的湖泊彭措为拉萨河源头。麦地卡湿地土壤类型多样, 主要包括高山寒漠土、高山亚高山草甸土、沼泽土以及新积土和粗骨土, 并孕育了丰富的动植物资源, 共记载维管束植物304种, 动物332种。麦地卡湿地是藏北地区典型的高原湖泊沼泽草甸湿地, 其湿地类型丰富, 包含永久性河流、洪泛平原湿地、永久性淡水湖、草本沼泽、沼泽化草甸等5类, 不仅是西藏地区的国际重要湿地、国家级自然保护区, 也是黑颈鹤(Grus nigricollis)、赤麻鸭(Tadorna ferruginea)、斑头雁(Anser indicus)等珍稀水禽冬季迁徙的重要栖息与繁殖地^[12]。

1.2 数据源与预处理

以1988、1990、1995、2000、2004、2015年生长季(6—10月)Landsat30 m遥感影像为基础数据源(数据来源于美国地质调查局USGS官网)。结合土地利用分类现状(GB/T 21010-2017)和全国湿地资源调查技术规程, 将研究区土地利用类型划分为湿地、草地、未利用地三大类, 并进一步分为河流湿地、湖泊湿地、沼泽湿地、天然牧草地、裸岩石砾地、冰川积雪等6种二级分类。2016年夏季对研究区进行了实地调查, 获取了研究区6种土地利用类型解译标志和57个GPS采样点。原始数据在ENVI软件做大气、几何校正后, 结合研究区地形、植被分布状况, 利用面向对象方法进行解译。解译结果使用混淆矩阵方法检验, 结果显示分类精度均大于80%。

选取研究区附近11个站点气象资料进行ANUSPLIN插值(气象数据来源于中国气象科学数据共享服务网), 根据研究区边界提取气象数据平均值。气象资料选取1985年以来的年平均气温、降水、地表温度和相对湿度。

2 研究方法 2.1 气候要素分析方法

利用回归统计方法分析麦地卡湿地保护区4种气候要素随时间的变化状况, 将时间t作为自变量, 气象因素x为因变量, 利用二阶方程进行拟合。同时, 采用非参数Spearman相关分析, 计算气象因素与土地覆被之间的相关性。

2.2 生态服务价值系数修正与价值计算

(1) 生态服务价值系数修正

标准生态系统生态服务价值当量因子经济价值量指1 hm²全国平均产量的农田每年自然粮食产量的经济价值, 结合生态系统服务当量因子, 可以表征和量化不同类型生态系统对生态服务功能的潜在贡献能力^[13-14]。根据谢高地的研究, 1个标准单位生态系统生态服务价值当量因子价值量相当于单位面积农田提供的食物生产服务经济价值的1/7, 参考王源等得到的纳木错流域当量因子生态系统服务价值628.6元^[15], 以及谢高地得到的西藏自治区农田生态系统生物量修正因子^[16], 计算得到麦地卡湿地保护区单位面积农田生态系统当量因子的价值量为471.45元。

采用千年生态系统评估(MA)方法, 将生态系统服务分为四大类11小类(表 1)。在求得研究区单位农田生态系统生产服务价值后, 结合“单位面积生态系统服务价值当量表”生态系统二级分类^[17], 最终得到研究区6类生态系统服务价值系数。

(2) 生态系统服务价值计算

牧业是研究区所属嘉黎县的主要产业, 研究区内除分散少数村庄外主要是草地景观, 缺少农田、森林生态系统。考虑研究区实际状况, 河流湿地、湖泊湿地当量值参考文献[17]表 1中的水系, 沼泽湿地参考湿地, 裸岩石砾地参考荒漠。生态系统服务价值计算公式：

(1)

(2)

式中, ESV和ESV_f为生态系统服务价值(元)和第f项功能的服务价值(元)；A_i为研究区第i类土地利用/覆被类型面积(hm²)；VC_i和VC_fi为第i类土地利用/覆被类型单位面积生态系统服务价值系数(元/hm²)和第i类土地利用/覆被类型所对应生态系统第f项功能的服务价值系数(元/hm²)。

3 结果与分析 3.1 土地利用类型空间分布特征

麦地卡湿地保护区土地利用类型空间分布如图 2。从一级分类看, 分布面积草地>湿地>未利用地, 从二级分类看, 分布面积天然牧草地>沼泽湿地>裸岩石砾地>湖泊湿地>河流湿地>冰川积雪。结合实地调查资料可知, 天然草地大面积分布在研究区西部和中部, 并在中部地形较为平坦区域与沼泽湿地相互交错。沼泽湿地和河流湿地主要分布在保护区中部、北部以及东北部海拔低于5100 m的沟谷和低洼地段。湖泊湿地主要分布在西南部高山脚下海拔5200 m左右的地带。未利用地主要分布在保护区南缘-东南缘高山地区, 冰川积雪主要分布在海拔高于5500 m的地区, 冰雪带至山脚为植被稀疏的裸岩石砾地。总体来看, 麦地卡湿地保护区土地利用分布与海拔和地形有关。

3.2 土地覆被动态变化

1988—2015年麦地卡湿地保护区土地利用类型变化如表 2。结合图 2、表 2可知, 麦地卡湿地保护区土地利用面积呈现明显的草地持续扩张而湿地连年减少的特点。研究期间, 草地占总面积的比例分别为53.4%、51.8%、54.3%、53.5%、61.1%、59.8%, 27年增长了5727 hm²。草地年均增长率在1988—2000年间为4.8 hm²/a, 而2000年之后激增为377.9 hm²/a, 主要是由湿地干旱演替成草地和裸岩石砾地向草地扩张造成。湿地面积占总面积的比例为32.4%、35.3%、31.1%、29.6%、26.9%、25.8%, 27年减少了5825 hm², 年均减少率为215.8 hm²/a。湿地面积减少主要以沼泽湿地为主, 减少了20%。通过叠加6期土地利用类型图可知, 沼泽湿地减少主要集中在保护区中心区和东北部彭措湖周边。27年间未利用地面积变化不大, 仅增加98.5 hm²。

3.3 格局演变的驱动因素

分析1985—2015年气象资料(图 3)可知：近30年来保护区气温和地表温度分别增加了1.2℃和2.5℃, 增温幅度为0.4℃/10a和0.8℃/10a。在1985—2000年期间, 年平均降水量由706 mm增加至773 mm, 之后到2015年减少至675 mm左右, 是30年间的历史最低水平。相对湿度变化规律与降水相似。

麦地卡湿地保护区沼泽湿地和天然牧草地面积发生明显变化, 将面积与4种气象要素做Spearman非参数相关分析, 探索麦地卡湿地保护区面积变化与气候要素之间相关性。年平均气温、降水、地表温度和相对湿度多重共线性检验结果为VIF=4.568、2.020、3.346、1.559, 表明气象要素之间不存在严重共线性, Spearman相关分析可以准确反应面积与气象要素的关系。

相关分析结果可知(表 3)：4种气象要素与沼泽湿地面积变化为负相关, 与天然牧草地面积变化为正相关。气温、地表温度和沼泽湿地、天然牧草地相关性大, 分别为0.543、0.468和0.543、0.371, 表明研究区沼泽湿地面积缩小和天然草地面积扩张受气温和地表温度升高的影响较大。结合1985年来气象要素变化可知, 降水减少和温度升高(尤其是地表温度的升高)加剧了土壤和气候的暖干化, 土壤湿度降低, 促使研究区藏嵩草为优势种群落的沼泽湿地演化为藏嵩草和高山嵩草为优势种群落的天然牧草地。同时温度的升高, 导致冰川和积雪融化, 在一定程度上增加了保护区内部分区域的水分供给, 促进了高寒草地向裸岩石砾地的扩张, 进一步增加了天然草地的面积。

3.4 湿地系统主要服务功能评估

如表 4所示, 麦地卡湿地保护区最主要的生态服务功能为调节功能, 占总服务功能的70%。二级生态服务功能中主要表现为水文调节、气候调节、生物多样性, 三项功能合计约占系统服务功能的70%；维持养分循环、食物生产、原材料生产功能较弱, 平均生态系统服务价值为464×10⁴元、1239×10⁴元和1361×10⁴元, 三者合计仅约系统服务功能的8%；其余水资源供给(均值5372×10⁴元)、气体调节(4851×10⁴元)、净化环境(7560×10⁴元)、土壤保持(5908×10⁴元)和美学景观(6647×10⁴元)合计约占系统服务功能的22%。

麦地卡湿地不同覆被类型生态系统服务功能价值变化如图 4所示, 在1988—2015年期间, 麦地卡湿地保护区沼泽湿地生态服务价值最大, 随着沼泽湿地面积的减少, 其生态服务价值由57979×10⁴元减少到44906×10⁴元。天然牧草地生态服务功能逐年增加, 27年增加了3086×10⁴元。湖泊湿地和河流湿地生态服务价值变化不大, 基本稳定在22585×10⁴元和8444×10⁴元。冰川积雪与裸露石砾地生态服务功能价值量很小, 基本维持在1166×10⁴元。

4 讨论

湿地的空间分布与其生态功能相互影响、相互制约, 任何一方改变都会引起另一方的变化^[18]。根据麦地卡湿地土地利用面积和生态系统服务价值变化可知, 沼泽湿地的生态服务价值是麦地卡湿地生态功能的重要体现, 沼泽湿地面积的大幅度减少, 是麦地卡湿地保护区生态服务价值逐年减少的主要原因。这一结果与李晋昌等对若尔盖高原生态系统服务价值的研究结果相同, 其发现若尔盖高原土地利用变化尤其是湿地的退化, 是该区生态系统服务价值减少的主要原因^[19]。

已有研究证明气候尤其是气温的变化是影响青藏高原湿地土地覆被的主要自然因素^[20-22]。在全球性气温升高、降水量重新分配、冰川消融的背景下, 如何减缓湿地面积萎缩, 充分发挥湿地生态调节功能, 有效制定湿地保护、规划管理方案迫切重要。基于以上研究, 针对麦地卡湿地保护区土地利用格局和生态系统服务价值现状, 以缓解湿地萎缩、提高湿地生态服务功能为目的, 提出以下几点建议：

(1) 加强麦地卡湿地保护区核心区域的科学管理, 迁出核心区和缓冲区居民, 极力减少人类影响。健全麦地卡湿地国家级保护区管理条例, 制定严格的保护措施。

(2) 加大麦地卡湿地保护区资金投入, 注重科研基础理论研究, 厘清麦地卡湿地退化的驱动机制和退化过程, 为有效保护、管理高寒湿地奠定理论基础。

(3) 为有效应对气候变化对湿地生态系统的负面影响, 建议对麦地卡湿地保护区采用生态恢复工程措施, 如梳理河道、筑坝蓄水、人工促进湿地植被恢复, 增加湿地土壤保水性, 促进湿地面积恢复。

(4) 土地利用覆盖类型间的斑块、廊道或基质的空间结构也会影响生态系统服务价值^[23-24], 今后可利用景观指数或模型探索麦地卡湿地景观结构、空间配置关系与生态系统服务价值的内在联系, 为湿地保护、修复和管理提供理论依据。

5 结论

本文以麦地卡湿地保护区为研究对象, 利用遥感解译和修正后的ESV模型分析了其土地利用格局和生态系统服务价值变化, 得出以下结论：

(1) 麦地卡湿地土地利用空间分布受海拔和地形影响, 草地和湿地分布在中部和北部的大部分地区, 而未利用地分布在南缘-东南缘高山地区。6类土地覆被类型面积为天然牧草地>沼泽湿地>裸岩石砾地>湖泊湿地>河流湿地>冰川积雪。1988—2015年间, 麦地卡湿地保护区草地面积增加了5727 hm², 湿地面积减少5825 hm², 并以沼泽湿地面积减少最明显。

(2) 近30年来, 麦地卡湿地保护区气温和地温增加显著, 降水无明显变化。相关分析结果显示, 气温、地表温度的升高对草地面积起着正向作用, 对湿地面积起着负向作用, 累积降水、相对湿度对二者影响较小。

(3) 调节服务是麦地卡湿地保护区的主要生态服务功能, 水文调节、气候调节、生物多样性占总生态服务功能的70%。湿地对保护区生态服务贡献最大, 沼泽湿地面积的缩小是麦地卡湿地保护区生态服务价值减小的主要原因。