湿地是水陆相互作用形成的独特生态系统, 是一切与水相关的生态系统服务的重要提供者, 在营养循环、气候变化、物种多样性、旅游文化等方面均发挥着重要作用^[1-2]。在全球气候变化和快速城市化的双重影响下^[3], 湿地资源面临着前所未有的压力, 不合理的使用和破坏导致湿地的数量和质量下降^[4-5], 迫切需要各国政府采取相应措施保护和修复湿地资源。湿地生态系统服务评价研究通过明确湿地生态系统与人类福祉之间的联系, 为湿地的规划和可持续发展提供了一个科学可靠的框架。

1997年, Costanza对全球17种生态系统服务进行的经济价值评估引起了巨大的关注, 此后有关生态系统服务的研究日渐丰富, 有关其定义、分类、评估、规划和管理等研究不断增加, 集中在森林^[6-7]、草地^[8]、流域^[9]和湿地^[10-11]等生态系统类型。生态系统服务价值评估很重要, 但是仅仅提供了价值上的评估结果, 对生态系统的格局和过程并不清楚；关于空间格局的研究多以土地利用数据为基础数据, 在模型运算基础上分析得到研究区的服务分布及变化情况；而对小区域的生态系统服务空间分布情况研究较少, 主要是因为小区域土地利用数据要求精度高, 且当前生态系统服务的运算模型多适用于大、中型区域。而且国内对城市湿地这类半自然半人工的生态系统服务研究较少^[12-13], 在评估基础上对空间格局展开的相关研究更是匮乏；此外, 明确生态系统服务的作用边界一直是复杂难以解决的问题, 同一区域内不同湿地的生态系统服务之间的空间溢出和叠加, 相关研究也亟待开展。

断裂点模型主要应用于经济地理研究中, 主要关于城市腹地范围、经济区划分确定和城镇吸引范围等^[14-15]；2017年乔旭宁等^[16]将断裂点模型和场强模型用于研究确定两区域之间生态系统服务价值的相互作用关系, 并确定生态系统服务的作用边界, 为断裂点的应用研究扩展了全新的方向。

本文以位于西宁市的湟水国家湿地公园为例, 在生态系统服务价值评价的基础上, 应用断裂点理论方法, 开展湿地价值和空间分布格局的研究, 解决生态系统服务的作用边界和空间溢出叠加确认问题, 对支持高原城市湿地生态系统的可持续利用具有重要的理论和现实意义。

1 研究区概况和数据来源 1.1 研究区概况

湟水国家湿地公园地处青海省西宁市城区段, 其分布范围为36°33′41″-36°44′42″N, 101°37′06″-101°54′50″E(图 1)。属高原半干旱气候, 气压低, 日照长, 太阳辐射强, 昼夜温差大；年平均气温为7.6℃, 降水的季节分配和年度分配不均匀, 年平均降水量为380 mm^[17]；公园内湿地可划分为2类2型, 为河流湿地和人工湿地2类, 永久性河流和洪泛平原湿地2个湿地型。

湟水国家湿地公园依湟水河及其一级支流北川河而建, 贯穿西宁城市中心, 总面积508.7hm²。以西宁市人民公园T字形水系的河道为中心, 北至北川河康家桥, 南至北川河与湟水河交汇处, 西至湟水河解放渠进水闸, 东至湟水河小峡口闸亭。其中包括海湖湿地公园、火烧沟公园、北川湿地公园、宁湖湿地公园和湟水河主干道公园五个主要公园。本研究以湟水国家湿地公园的五个主要公园为研究区, 具体监测点分布如图 1所示。

1.2 数据来源

数据主要采集于2020年, 包括野外调查与实验分析数据、社会调查, 以及历史数据。在湟水国家湿地公园共布设30个样点, 根据湿地公园与湟水河的进出水断面位置, 在各公园进水口之前和出水口之后的湟水河主干道、北川河河道总共设置了10个样点；5个公园的人工湿地内, 沿进水口、出水口方向均匀布设样点, 总计20个(图 1)。

每月进行一次野外监测, 包括：样点的基本概况、水文参数的监测、各样点收集三个平行水样(1L)；植物调查, 在每个样点附近布设3个50cm×50cm样方, 采用网格法和目测法确定植物盖度, 盖度、植物丰富度、生物量都以3个样方的平均为准；鸟类样方调查记录以各样方的50m为半径的范围内见到和听到的鸟类种类和数量, 每个样方停留不低于10min, 指标包括：物种数、个体数。实验室进行水质的TOC、TN、TP含量测定；其中TOC采用TOC分析仪法, TN和TP分别采用碱性过硫酸钾消解紫外分光光度法和钼酸铵分光光度法。于2020年5月-7月的8点-20点前往各公园进行调查问卷发放, 共发放1300份问卷, 剔除无效问卷和不完整问卷40份, 有效回收率为96.92%；访问对象选择公园内的游客, 进行面对面访谈；访问结果选择分层抽样方法以使样本更具有代表性, 并剔除支付意愿超出收入5%及以上的数据。土地利用面积情况参考《青海西宁湟水国家湿地公园总体规划(2013-2020年)》。部分评估数据参考《中国生物多样性国情研究报告》^[18]、《青海水文手册》^[19], 和中华人民共和国林业行业标准《森林生态系统服务功能评估规范》(LY/T1721-2008)^[20]中各产品的推荐使用价格, 以及TN、TP和TOC的生活污水处理成本^[21], 由于以上各评估对象的市场价格出现年份不同, 通过2020年物价指数修正, 参考以下公式^[22]：

(1)

式中, P₂₀₂₀为评估对象修正后2020年价格(元)；P_i为评估对象原价格(元)；L₂₀₂₀为2020年物价指数；L_i为评估对象所在时间的物价指数。

2 研究方法 2.1 生态系统服务评价方法

(1) 湟水国家湿地公园评估指标体系

根据湟水国家湿地公园各生态系统服务类型特点, 剔除不重要的指标, 选取最适合的评估指标和估算方法, 结合实测数据与不同的土地利用类型, 进行生态系统服务价值量估算。支持服务是供给、调节和文化服务的基础与源泉, 为整个湿地生态系统提供物质循环和能量流动, 在计算过程中容易与其他三大服务出现重复计算^[23], 本研究仅对湿地为鸟类提供栖息地的服务进行估算, 与其他服务类型无明显交集。此外, 文化服务中的休闲娱乐价值指游客为了自己或他人能够在湿地公园休闲放松、锻炼健身而支付的费用；存在价值指没有任何功利目的, 只为了让湿地资源永续存在, 游客自愿支付费用；遗产价值是游客为了将湿地资源作为遗产留给子孙后代, 而自愿支付的费用。最终以千年生态系统评估(MA)的生态系统服务分类为基础, 结合湟水国家公园的生态系统特征, 将各项服务分为2个层级, 第一级分类分为13个服务项目, 第二级具体评估指标对第一级进行细化区分, 保证各层次指标间相互独立；最终筛选出研究区的生态系统服务具体评价指标共15项。具体评价体系见表 1。

(2) 城市湿地生态系统服务评估方法

根据湟水国家湿地公园的生态特点, 依据《湿地生态系统服务评估规范》^[24]理论基础, 以市场价值法为基础, 筛选最优方法对各项服务进行价值评估, 选择最适合、最易操作的计算方法, 结合野外调查数据和实验分析数据, 使研究结果更贴近实际情况。具体采用市场价值法、碳税法、工业制氧价格法、费用支出法、影子工程法、成果参数法和条件价值评估法(CVM)等功能价值法对各项服务进行价值评估, 综合各类方法准确评估城市湿地生态系统服务。具体计算过程见表 2。

2.2 断裂点理论及方法

断裂点理论包括断裂点公式, 用于确定起点城市到断裂点的距离；场强公式, 即起点城市在断裂点的辐射范围；辐射半径由断裂点处场强大小决定起点城市的辐射范围。通过断裂点模型定量计算不同公园的生态系统服务的作用边界^[15-16]。具体过程如下：

(1) 断裂点公式^[30]：

(2)

式中, D_i为公园i到断裂点的距离；D_ij为公园j中心点与公园i中心点间的距离；V_i、V_j分别为公园i和公园j生态系统服务价值量。

(2) 场强公式^[30]：

(3)

式中, F_ik为公园i在断裂点k处的辐射能力大小(转移能力里的强度)；V_i为公园i的生态系统服务价值量；D_ik为公园i中心点到断裂点k的距离。

(3) 辐射半径公式^[30]：

(4)

式中, D_r为公园r生态系统服务作用半径；V_i为公园i的生态系统服务价值量；F为边界场强, 可以根据各公园的综合关系确定F取值。

公式(2)中两个公园的中心点间的距离, 采用交通距离；结合各个公园的生态系统服务总价值量, 计算得到5个湿地公园彼此之间的断裂点到起点公园的距离, 找到两公园的断裂点位置；利用公式(3)计算公园在断裂点处生态系统服务的辐射强度F, 即公园生态系统服务在转移过程中的强度。

各个公园的辐射半径, 即为该公园生态系统服务的辐射影响范围, 利用公式(4)对辐射半径的计算主要在场强的选择上。同一起点公园到其他不同公园, 在平衡点处发挥的生态系统服务辐射强度不同, 选择能发挥最大辐射强度的场强, 进行生态系统服务强辐射半径的计算, 在该半径内公园的生态系统服务能够得到高效率发挥。通过湿地公园与其他四个公园的场强平均值进行弱辐射半径的计算, 采用场强的平均值可以代表该湿地公园对其他公园的平均影响强度。计算得到各公园的强辐射半径和弱辐射半径, 利用ArcGIS9.2的缓冲区分析工具得到各个公园的辐射范围, 强辐射区作为公园生态系统发挥效用的核心区, 弱辐射区为公园的生态系统过渡区。

任意两公园之间的断裂点位置基本位于两公园中间位置, 采用起点公园到断裂点距离值, 在两公园总距离的占比情况, 占比总会大于或小于50%, 以此说明不同公园受起点公园的影响情况, 即由于起点公园不同, 在任意两公园中间的断裂点位置不同, 当占比大于50%时, 说明该方向的公园存在生态系统服务溢出, 与另一公园的生态系统服务叠加。

3 结果与分析 3.1 生态系统服务的价值构成

结合实测数据、调查问卷数据和物价修正后物品单价等, 由表 2的评估方法计算得到湟水国家湿地公园的各项生态系统服务价值量(图 2)。湟水国家湿地公园生态系统提供生态系统服务总价值为6.18亿元, 核心服务为增湿服务、休闲娱乐和遗产价值。5大公园的生态系统服务价值量对比中, 湟水主河道的价值最高, 约占总体价值的37.92%, 达到了23441.37万元/a。其次为北川湿地和海湖湿地, 分别为14008.26万元和8964.78万元；火烧沟公园和宁湖湿地差别较小, 分别为7802.04万元和7600.74万元。

本研究中评估的15项生态系统服务类型价值量由高到低依次为：增湿>休闲娱乐>遗产价值>存在价值>地表水调蓄>河流输沙量>文教科研>温度调节>生物多样性保护>生态环境用水>净化空气价值>释氧>固碳>污染物降解>植物产品价值。调节服务中的增湿服务、文化服务中休闲娱乐、存在价值和遗产价值是湟水国家湿地公园的核心功能, 价值量达56437.47万元, 对总价值的贡献率91.30%。而年价值量低于100万元的服务类型包括净化空气价值、释氧、固碳、植物产品价值和污染物降解, 占总价值的百分比分别为0.144%、0.081%, 0.037%, 0.017%和0.036%；其中污染物降解的价值量最低, 但当水质污染超过自净能力时, 不仅会使水质继续恶化, 还会因水体富营养化影响休闲娱乐功能的发挥, 因此, 加强对水质的管理和净化尤为重要。

3.2 湿地公园的空间分布格局

由断裂点理论计算得到5个湿地公园彼此之间的断裂点到起点公园的距离和起点公园在断裂点处发挥的场强, 结果见表 3。根据表 3计算结果选择各公园的最大场强和场强平均值计算辐射半径；海湖湿地的强辐射半径采用海湖湿地-火烧沟公园的场强, 火烧沟公园、北川湿地、宁湖湿地和湟水主河道的强辐射半径分别采用火烧沟公园-湟水主河道、北川湿地-湟水主河道、宁湖湿地-湟水主河道, 以及湟水主河道-火烧沟公园的场强, 以上两公园的场强均大于其他公园组合。弱辐射半径通过湿地公园与其他四个公园的场强均值进行计算, 海湖湿地、火烧沟公园、湟水主河道、北川湿地和宁湖湿地的平均场强分别为1054.97、2089.95、497.52、2030.80和118.42。根据以上选择和计算得到的场强值, 进一步计算各公园强辐射半径和弱辐射半径, 结果见表 4, 在此基础上缓冲得到各个公园的辐射范围(图 3)。

增湿、休闲娱乐、遗产价值是湟水国家湿地公园的三大核心服务, 根据断裂点模型, 结合各公园的核心服务价值量, 分别计算各个公园核心服务的辐射半径, 利用ArcGIS缓冲区分析得到三大服务各个公园的辐射范围。将其与公园总价值分布图叠加, 得到增湿、休闲娱乐、遗产价值的空间分布图, 并与总价值分布图对比, 相关结果见图 3)。

由图 3明显看出, 宁湖湿地的生态系统服务总价值分布范围最大, 其次为北川湿地和湟水主河道, 最后是海湖湿地和火烧沟公园, 这与各个公园的总价值量情况不同。宁湖湿地的分布区最大, 该湿地的生态系统服务的贡献范围最广, 主要是因为该湿地位于城东区, 西宁市东部地区湿地资源匮乏, 居民需求大, 到公园的可达距离增大；北川湿地的生态系统服务范围基本覆盖城北区, 作为新修复完成的湿地, 服务价值增加的潜力巨大, 未来北川湿地的生态系统服务的分布区有望进一步扩大。湟水主河道选择以北川河汇入处的湟水河为中心点, 服务范围以河道两侧的城北区、城中区为主。海湖湿地和火烧沟公园的核心区和过渡区分布面积较小, 且重合区域大, 海湖湿地和火烧沟公园的生态系统服务面积小的原因可能与两个公园距离较近, 游客竞争力大, 且城西区公园数量较多, 海湖湿地和火烧沟公园的生态系统服务的发挥空间较小。

湟水国家湿地公园的三大核心服务增湿、休闲娱乐、遗产价值的空间分布图与总价值分布图的对比情况, 见图 3。宁湖湿地和火烧沟公园的休闲娱乐、遗产价值的空间分布面积大于总价值分布范围, 而增湿服务的分布面积远小于总价值分布面积；海湖湿地的增湿服务、休闲娱乐的空间分布与总价值分布范围大体一致, 而遗产价值分布面积远小于总价值分布面积；北川湿地、湟水主河道的增湿服务分布面积与总价值分布范围相差较小。此外, 湟水主河道的休闲娱乐、遗产价值的分布区小于总价值服务的分布区；北川湿地的休闲娱乐空间分布面积大于总价值分布范围, 遗产价值分布区与总服务的分布区大致相同。

3.3 生态系统服务重叠和溢出

在湟水国家湿地公园生态系统服务总价值分布图中看到, 海湖湿地、火烧沟公园、湟水主河道和北川湿地分布相对集中, 各公园的空间分布重合区域较多。将图 4各公园的最大分布范围与西宁市区的城市功能区相叠加(图 4), 对湿地公园的空间分布重合情况进行研究分析。三个公园生态系统服务的共同分布区重叠后形成1个优选区, 两个公园生态系统服务的共同分布区重叠后形成5个次优区。

优选区为火烧沟公园、湟水主河道和北川湿地3个公园的服务分布区形成一个交叉区域, 面积为139.33hm²；该区域以住宅区为主。5个次优区中北川湿地和湟水主河道的交叉区面积最大, 为1145.23hm², 主要分布在城北区, 以住宅区为主；火烧沟公园和湟水主河道的重叠区较大, 为473.29hm², 该区域拥有丰富的文化资源和住宅基地；海湖湿地和火烧沟公园的重叠区面积为142.40hm², 该区域为当前西宁市的主要商业分布区, 商业价值高；火烧沟公园和北川湿地的交叉区面积最小, 为44.40hm²；海湖湿地和北川湿地的重叠区域区别于其他4个区域, 该区域远离城市中心且受高速影响, 居民聚集程度低, 总面积74.39hm²。综上, 各大湿地的生态系统服务分布区不同, 其重叠区也不同, 但各重叠区域均为重要的西宁市区住宅、文化和商业集中分布的区域, 人口聚集程度高, 商业价值大, 人类可利用价值得到最大发挥；而且在原各自湿地服务空间中, 相比其他区域占绝对优势, 对人类福祉贡献较大。

生态系统服务也存在溢出的情况。由表 3的占比情况可知, 湟水国家湿地公园中任意两公园之间的断裂点位置基本位于两公园中间位置；同时受各公园价值量不同影响, 占比总会大于或小于50%。如海湖湿地到湟水主河道的断裂点在距海湖湿地2.63km处, 占总距离的38.21%, 而从湟水主河道到海湖湿地的断裂点在距湟水主河道4.26km处, 占总距离的61.79%, 说明在海湖湿地和湟水主河道的生态系统服务发挥功能时, 海湖湿地受湟水主河道的影响较大。

海湖湿地、火烧沟公园和北川湿地的占比存在大于50%, 影响其他公园的情况, 也有小于50%, 受其他公园影响的时候, 类似于生态系统服务的流动“中转站”；而湟水主河道与其他4个公园的占比均大于50%, 则说明湟水主河道的生态系统服务均影响其他公园, 属于生态系统服务的“供给方”, 是为其他公园提供生态系统服务的生态中心, 这与实际情况一致, 海湖湿地、火烧沟公园、北川湿地和宁湖湿地均是依托湟水主河道修建, 并由湟水主河道将其串联。宁湖湿地与以上情况相反, 与其他4个公园的占比均小于50%, 表明宁湖湿地在5个湿地公园中属于“需求方”, 受其他公园影响；但由于距离较远, 受影响强度较小。

4 结论

本文筛选构建了一套适合湟水国家湿地公园的生态系统服务评价体系, 基于野外调查数据和实验分析数据, 采用功能价值法等计算方法计算生态系统服务价值, 并应用断裂点模型对生态系统服务的空间分布和流转情况进行分析, 主要结论如下：

(1) 对2020年湟水国家湿地公园的15项生态系统服务进行评估, 得到总价值为6.18亿元, 四大服务价值量按照价值高低排序依次为文化服务、调节服务、支持服务和供给服务；各分项价值量对比, 由高到低依次为增湿>休闲娱乐>遗产价值>存在价值>地表水调蓄>河流输沙量>文教科研>温度调节>生物多样性保护>生态环境用水>净化空气价值>释氧>固碳>污染物降解>植物产品价值。

(2) 5大湿地区中, 湟水主河道的价值最高, 其次为人工湿地, 依次为北川湿地、海湖湿地、火烧沟公园和宁湖湿地。由于高原城市湿地的特殊性, 湟水国家湿地公园的核心服务类型体现在存在价值、休闲娱乐和小气候调节上。

(3) 湿地价值辐射范围由大到小依次为：宁湖湿地>北川湿地>湟水主河道>海湖湿地>火烧沟公园；对湿地服务辐射范围叠加后, 形成1个由三个公园的服务分布区重叠的优选区和5个两两交叉的次选区。

(4) 湟水国家湿地公园的生态系统服务存在空间溢出和叠加现象, 湟水主河道的生态系统服务流向周边湿地公园溢出, 其他四个湿地区之间生态系统服务存在相互溢出、叠加的现象。

值得注意的是, 受地形、交通、人口分布等因素的影响, 生态系统服务的空间并非均衡的。受数据等条件限制, 本研究对生态系统服务的辐射空间做了均衡处理, 这也是本研究的不足之处, 需要在后续研究中进一步改善。