朋曲河流域位于中国西南边陲西藏自治区日喀则市南部，地理坐标85.703-88.951°E，27.831-29.128°N，涉及西藏自治区日喀则市萨迦县、聂拉木县、定结县、定日县和岗巴县，南与尼泊尔相邻，面积29 400 km²(图 1).流域大部分区域位于喜马拉雅山脉北侧，海拔多在4 000 m以上，焚风效应明显，气温偏低，降雨偏少，蒸发量大，植被多为草地和灌丛；少部分位于喜马拉雅山脉南侧，平均海拔2 400 m，受印度洋湿暖气流影响，降雨充沛，植被多为森林.流域内地带性土壤为亚高山草原土，而河漫滩和洪积扇缘多为草甸土，属沙质土.区内淡水生态系统主要包括沼泽湿地、草甸、河流和湖泊，面积138 691.81 hm²，仅占流域总面积的4.72%，其中，沼泽湿地及草甸面积最大，主要分布在河流两侧及湖泊周围，面积118 380.36 hm²；其次是湖泊淡水生态系统，主要有丁木措和定结措(包括强左措和共左措)等，面积10 822.41 hm²；还有河流淡水生态系统，主要有朋曲河及支流叶如藏布江等，面积9 489.04 hm².

根据TEEB^[7]生态系统服务分类体系，本研究重点针对朋曲河流域淡水生态系统较重要的供给服务(水资源供给、饲料生产和水力发电)、调节服务(废物处理、水文调节、土壤保持、气体调节和气候调节)、栖息地服务(生物多样性维持)和文化服务(娱乐文化)进行价值评估(表 1).其中，供给服务中的水资源供给、饲料供给和水力发电服务价值评估采用市场价值法^[8-10]，调节服务中的水文调节、气体调节服务价值评估分别采用影子工程法和影子价格法^[11-12]，调节服务中的气候调节、废物处理和土壤保持服务，以及栖息地服务和文化服务的价值评估采用基于单位面积价值当量因子的生态系统服务价值评估方法(简称当量因子法)^[13-14].下面简要介绍当量因子法计算过程.

当量因子法最先由Costanza等^[15]提出.后来，谢高地等^[13-14]结合中国实际情况对生态系统服务单位面积价值当量进行了优化，制定了中国陆地生态系统服务单位面积价值当量表(表 2)，并得到广泛应用^[16-17].本研究中的沼泽湿地及草甸对应表中的湿地，河流湖泊对应表中的水体.

为了尽量反映项目区生态系统服务实际价值，这里的生态系统服务价值当量因子，采用西藏莽措湖流域1 hm²农田每年自然粮食平均产量的价值177.32元/(hm²·a)^[18].当量因子再乘以生态系统服务单位面积价值当量，可以获得朋曲河流域各生态系统服务的单位面积价值(表 3)，如公式1所示.然后，把各类生态系统面积与其各生态系统服务单位面积价值的乘积再相加，可以计算出朋曲河流域各类生态系统服务的总价值，如公式2所示.

式中：P_ij为j类生态系统i种生态系统服务单位面积价值，单位为元/(hm²·a)，结果如表 2所示；p_ij为j类生态系统i种生态系统服务单位面积价值当量(表 1)；C为项目区农田每年自然粮食平均产量的价值，这里取177.32元/(hm²·a)；V为区域生态系统服务总价值，单位为元；S_j为j类生态系统面积，单位为hm²；n为生态系统服务类型；m为生态系统类型.

朋曲河流域淡水生态系统水资源丰富，每年向周边居民提供多种生产用水和生活用水，具体包括生活用水、灌溉用水和牲畜饮水等. 2014年朋曲河流域5县有城镇人口22 287人、农村人口13 8851人，人均日用水量分别按照2014年全国平均水平213L和81L估算，则朋曲河流域年生活用水量583.78万m³.生活用水价格取3.0元/m^3[19]估算，则朋曲河流域年生活用水供给价值为1 751.35万元.朋曲河流域耕地面积16 475.97 hm²，灌溉用水量以青稞单位面积需水量3 750 m³/hm^2[20]进行估算，年灌溉耗水量约6 178.49万m³.这里农田灌溉水成本取0.25元/m^3[19]，则朋曲河流域年灌溉水供给价值为1 544.62万元. 2014年朋曲河流域5县共出栏牛43 992头、羊424 676只、猪711只，日需水量分别按50，5，25 L/d^[21]估算，则年牲畜饮水158.44万m³.这里牲畜饮水成本取0.25元/m^3[19]，则朋曲河流域年牲畜饮水供给价值为39.61万元.由上可知，淡水生态系统水资源年供给服务价值为3 335.58万元.

朋曲河淡水生态系统，尤其沼泽湿地与草甸地势较低，草本植物长势好，是良好的四季牧场和冬季牧场.为了合理利用区域草场资源，朋曲河流域农户实施轮牧制度.夏季把羊群赶到高山牧场放牧，冬季则回到宽谷沼泽湿地与草甸放牧.经调查，宽谷沼泽湿地与草甸放牧利用时间约1个月，而牛马等大牲畜则一年四季在沼泽湿地和草甸牧场放牧.朋曲河流域沼泽湿地和草甸单位面积可食干草产量按日喀则平均水平0.50 t/hm^2[22]估算，则沼泽湿地及草甸年生产可食干草58 835.04 t.牧草价格取拉萨市场价0.40元/kg^[23]，则沼泽湿地及草甸年饲料供给服务价值为2 353.40万元.

朋曲河干流落差3 370 m，水能蕴藏量较丰富，但境内干流水能资源尚未开发，目前已开发的主要是农村小水电.这里仅收集到定结县和定日县水力发电数据，其中，定结县正在运行的水电站3座，分别是长所水电站、卡达水电站和曲当水电站，年发电量为1 440万kw/h.定日县正在运行电站有6座，年发电量为360万kw/h.电力售价取居民用电价格0.50元/(kw·h)(实地调查)，则朋曲河流域现有水电站年水力发电服务价值为900万元.

朋曲河沿线集镇居民点两污设施建设严重滞后，尚未建立污水处理设施，生活污水未经处理直接排放，大部分依靠附近及下游淡水生态系统自然净化.此外，放牧时牲畜粪便直接排入沼泽湿地和草甸，通过淡水生态系统自然降解.由于缺乏朋曲河流域淡水生态系统净化能力和污染负荷数据，这里采用当量因子法估算流域淡水生态系统废物处理服务价值.由表 3可知，朋曲河流域沼泽湿地单位面积废物处理价值2 553.41元/(hm²·a)，湖泊河流单位面积废物处理价值2 633.20元/(hm²·a)，则朋曲河流域淡水生态系统年废物处理服务价值为35 575.75万元.

由于年降雨分配不均，朋曲河流域雨季多次发生洪水灾害，冲毁农田；而旱季则降水较少，气候干旱.而沼泽湿地及草甸和湖泊生态系统均具有较强的水源涵养和洪水调蓄能力，雨季能够消纳部分面山径流，降低洪涝灾害风险；干季则可以缓慢释放雨季储蓄水资源，补充周边地下水与河水，缓解干旱.其中，沼泽湿地及草甸水源涵养按8 100 m³/hm^2[24]估算，则沼泽湿地与草甸蓄水能力为95 888.09万m³；湖泊水资源调蓄能力参考青藏高原湖泊水量调节能力评价模型^[25]，则湖泊蓄水能力18 249.45万m³；由上可知，朋曲河流域以上二类淡水生态系统年水资源调蓄能力合计为114 137.54万m³.以西藏地区修建水库投资约6元/m^3[23]计，朋曲河流域淡水生态系统年水文调节服务价值为684 825.26万元.

朋曲河流域土壤属沙质，再加上年降雨集中，风力大，土壤风蚀水蚀现象较突出.区内面山普遍存在深浅不一的冲蚀沟，部分冲蚀沟深达一到两米；还有些区域的沙尘在风力作用下形成连片沙丘，吞没植被和农田，影响大气质量.而沼泽湿地及草甸植被覆盖度高，能够一定程度缓解降低区域土壤侵蚀影响；湖泊则能淤积泥沙.由于缺乏土壤保持相关数据，这里采用当量因子法估算流域淡水生态系统土壤保持服务价值.由表 3可知，朋曲河流域沼泽湿地及湿地单位面积土壤保持价值352.87元/(hm²·a)，湖泊单位面积土壤保持价值72.70元/(hm²·a)，则朋曲河流域淡水生态系统年土壤保持服务价值为4 255.93万元.

朋曲河流域地处青藏高原地区，大气氧浓度低，因此生态系统释放氧气对人类生存具有重要意义.淡水生态系统植物生长良好，初级生产力较高，能不断吸收大气中的CO₂，释放O₂，维持大气中O₂浓度.这里采用光合作用方程法估算氧气释放量.沼泽湿地及草甸和产草量按0.945 t/hm^2[23]估算，则朋曲河流域淡水生态系统年产草量为111 869.44 t.根据光合作用方程，每产生1g干草会释放1.2 g O₂，由此可知朋曲河流域淡水生态系统年释放O₂ 1 34 243.33 t.以拉萨市医用氧气制造价格900元/t^[23]估算，则朋曲河流域淡水生态系统年释放O₂的价值为12 081.90万元.淡水生态系统吸收CO₂的价值同时反映在气候调节中，这里不再重复核算.

朋曲河流域淡水生态系统植被生长过程中，会通过光合作用不断固定CO₂，部分还会形成泥炭，长期固定，进而缓解全球气候变化.同时，淡水生态系统还会通过蒸、散发改变局部湿地、温度和降水状况，使之适宜人类居住.由于缺乏朋曲河流域气候调节相关研究，这里采用当量因子法估算流域淡水生态系统气候调节服务价值.由表 3可知，朋曲河流域沼泽湿地及草甸单位面积气候调节价值2 402.69元/(hm²·a)，湖泊河流单位面积气候调节价值365.28元/(hm²·a)，则朋曲河流域淡水生态系统年气候调节服务价值为29 185.02万元.

朋曲河流域淡水生态系统生境优越，适宜多种动植物生存.据现场初步调查，朋曲河流域淡水生态系统内分布有维管束植物62种、鸟类71种、哺乳类5种和鱼类3种，其中包括黑颈鹤等多种珍稀保护物种.这里采用当量因子法估算流域淡水生态系统生物多样性维持服务价值.由表 3可知，朋曲河流域沼泽湿地及草甸单位面积生物多样性维持价值654.31元/(hm²·a)，湖泊河流单位面积生物多样性维持价值608.21元/(hm²·a)，则朋曲河流域淡水生态系统年生物多样性维持服务价值为8 981.11万元.

朋曲河流域位于青藏高原西南部，拥有丰富的河流、沼泽湿地、冰川、山体、草原、湖泊、沙丘等自然景观资源及神秘的藏族文化资源，可以为人们提供美学体验与休闲等精神享受，旅游价值较高.但由于资金缺乏和配套设施不完善，区内旅游建设开发严重滞后，除定日县珠峰大本营和定结县陈塘镇旅游开发外，其他旅游资源基本未开发.此外，朋曲河属于国际河流，淡水生态系统维持较完整，生态安全意义突出，科学研究价值较高.这里采用当量因子法估算流域淡水生态系统娱乐文化服务潜在价值.由表 3可知，朋曲河流域沼泽湿地及草甸单位面积娱乐文化价值831.63元/(hm²·a)，湖泊河流单位面积娱乐文化价值787.30元/(hm²·a)，则朋曲河流域淡水生态系统年娱乐文化服务价值为11 444.00万元.

综合以上10项生态系统服务价值，朋曲河流域淡水年生态系统服务总价值792 937.95万元(表 4)，其中，调节服务价值最大，为765 923.86万元，占总价值的96.59%，远大于其他3类生态系统服务价值；其次是文化服务和栖息地服务，分别为11 444.00万元和8 981.11万元，占总价值的1.44%和1.13%；供给服务价值最小，为6 588.98万元，占服务总价值的0.83%(图 2).

朋曲河流域近2/3的区域属于珠穆朗玛峰国家级自然保护区，是我国青藏高原生态屏障的重要组成部分，其淡水生态系统的调节服务及栖息地服务十分重要.经测算，不足流域面积5%的淡水生态系统年调节和栖息地服务价值774 904.97万元，相当于5县GDP总和的4倍多.但是，流域内居民往往更重视供给服务的开发利用，忽视其他生态系统服务效益保护.比如，畜牧产业是朋曲河流域的传统产业，是当地居民的重要收入来源.牧业发展利用主要是为流域内淡水生态系统的饲料提供供给服务，该服务价值仅占生态系统服务总价值的0.30%.为了提高收入水平，过去牧民大幅增加牲畜数量，甚至出现超载放牧，导致淡水生态系统出现退化.现场调查发现，部分沼泽湿地及草甸区域已发现有棘豆、狼毒等有毒有害物种分布，甚至局部区域出现棘豆优势群落，这可能与过度放牧有关.因此，朋曲河流域地方政府应加强淡水生态系统保护和管理，努力协调不同生态系统服务间的关系.

供给服务与其他生态系统服务之间主要是竞争关系^[26].过度放牧只会短期内带来经济收入的增长，长期必然造成大面积淡水生态系统退化，影响的不仅是饲料供给服务，还会威胁其他生态系统服务效益，产生一系列不良后果.西藏生态系统十分脆弱，尤其在气候变化背景下，超载放牧等人为活动对淡水生态系统破坏往往具有加速和倍增作用^[27].因此，朋曲河流域应加强淡水生态系统牧业管理，制止挖草皮、湿地开挖、湿地征占、湿地垦殖等一系列破坏淡水生态系统的不良行为，保护淡水生态系统健康安全. 2010年，项目区开始实施草原生态保护补助奖励机制和草畜平衡管理制度，有助于协调朋曲河流域淡水生态系统供给服务与调节服务、栖息地服务和文化服务之间的冲突.但需要注意的是，目前实施的草畜平衡模式是否能充分保护调节服务、栖息地服务和文化服务效益还需要进一步的研究.

文化服务与供给服务以外的其他生态系统服务之间存在较大协同关系^[26].因此，在传统牧业发展受到一定调控时，朋曲河流域可逐步有序开发淡水生态系统旅游业，充分发挥其文化服务效益，打造新的经济增长点，从而促进地方经济社会可持续发展.

此外，朋曲河流域淡水生态系统的调节服务、栖息地服务和文化服务还具有很强的外溢效应(表 4).据研究，科西河流域水旱灾害频发，已成为当地居民生活贫困的一个重要原因^[4].保护朋曲河流域淡水生态系统，维持并提高其水文调节服务效益，将有助于降低下游科西河洪涝灾害及干旱风险.朋曲河年径流量49.2亿m³，扣除流域内水资源供给消耗的6 920.71万m³及其他生态用水，绝大部分水资源无偿流入下游尼泊尔和印度境内的科西河，为下游提供大量水资源.还有，栖息地服务和文化服务的受益对象为全球，是全人类的宝贵财富.因此，朋曲河流域淡水生态系统保护工作还应积极争取国际资金，尤其来自下游地区的生态效益补偿.

西藏朋曲河流域淡水生态系统位于珠穆朗玛峰国家级自然保护区，生态保护意义十分重要，同时又是当地居民生产生活的重要基础和资源.本研究评估结果表明，项目区淡水生态系统具有较高的生态系统服务价值，尤其是调节服务价值.因此，当地政府应当加强流域内淡水生态系统保护，控制对调节服务影响较大的畜牧产业及湿地草甸破坏行为，适当加大开发对调节服务影响较小的文化服务.