以青藏高原为核心的第三极^[1]以及受其影响的东亚、南亚、中亚、西亚、中东欧等泛第三极地区是地球上生态环境最脆弱和人类活动最强烈的地区, 面积广达2000多万km^2[2], 其中喀喇昆仑山-喜马拉雅山脉地区则是泛第三极的重要组成部分。这个地区位于青藏高原的边缘地带, 受大断裂控制, 由一系列高大山系和山脉组成^[3], 生态环境极其独特, 集聚了全球所有14座海拔8000 m以上的高峰, 是地球上最壮观的高山景观, 平均海拔都在4500 m以上, 众多世界大河或者国际河流皆出于此。这个区域涵盖了全球34个生物多样性热点地区中的3个, 属于多种地理成分的动植物交汇过渡地带^[4]。

复杂的生态环境、过渡地带特征、政治边界的分割, 使得喀喇昆仑山-喜马拉雅山脉地区作为一个地理区域在全球尺度上难以区划。Schultz将其北部划为干旱中纬度的荒漠和半荒漠, 南部则划为夏季湿润热带或中年湿润热带, 未能对内部结构细分^[5]；Bailey只能将其粗略划分为山地草原、山地荒漠和山地稀树草原^[6]。Olson等提出的以多套生物群区为基础的, 自下而上的全球陆地生态区划^[7-8], 较好地将生物地理分区与生态功能单元相衔接, 着重强调地理单元中生态系统特征分异, 形成具有相对同质性的以植被为基础的生态系统组合-生态区(ecoregion)的概念, 被世界野生生物基金会和大自然协会采用, 为本地区生态系统保护发挥了关键的空间指引作用^[9]。

从20世纪30年代开始, 特别是80年代以来, 中国以及周边国家在此建立了多个国家公园和自然保护区, 有效地促进了喀喇昆仑山-喜马拉雅山脉地区生态环境的保护。随着本地区跨国合作和大尺度景观保护的需求增加, 连续性保护理念在跨国项目中体现为跨界保护理念^[10], 跨界自然保护区作为跨界资源管理的一种典型方式正成为研究热点^[11], 加强国际间或国家内不同行政区之间的合作来保证资源的可持续利用可能是必然的选择。喀喇昆仑山-喜马拉雅山脉地区包含东喜马拉雅山、干城章嘉等跨境景观^[12], 当前主要研究包括整合地理空间数据进行生物多样性和物种信息保护、应用跨界景观方法保护文化遗产、总结跨界保护合作机制及管理经验等^[13-20]。尽管在部分热点区域已进行跨界研究, 但目前缺少从喀喇昆仑山-喜马拉雅山脉地区的角度对涉及的7个国家进行全面评估的研究, 因此厘清本区域生态区划特征和主要类型, 阐明保护地现状及其空间分布规律, 掌握保护地与生态区之间的相互关系, 有利于未来更好地确定喀喇昆仑山-喜马拉雅山脉地区保护地跨境合作的重点区域和建设路径, 促进泛第三极(Pan-TPE)地区的可持续发展。

1 区域简介

喀喇昆仑山-喜马拉雅山脉地区位于24°0′31.76″—39°5′38.97″N, 71°16′3.23″—99°56′41.77″E之间, 涵盖了中国、阿富汗、巴基斯坦、印度、尼泊尔、不丹和缅甸等7个国家的边境山地, 全长约2500 km, 宽200—300 km^[21]。喀喇昆仑山脉是世界山岳冰川最发达的高大山脉, 位于中国、塔吉克、阿富汗、巴基斯坦和印度等国的边境上, 包括东延部分(西藏高原的羌臣摩山和潘顿山)在内的喀喇昆仑山, 宽度约为240 km, 长度为800 km, 平均海拔超过5500 m, 是世界上高山和高纬度之外最长的冰川最集中的地区。喜马拉雅山脉西起巴属克什米尔境内的南迦帕尔巴特峰, 山麓有印度河上游大峡谷围绕, 东至南迦巴瓦峰, 有雅鲁藏布江下游围绕作急转弯南流而形成的举世闻名的大拐弯峡谷, 主体由四列大致平行的山脉组成, 呈向南凸出的弧形, 从南至北被命名为外喜马拉雅山脉、小喜马拉雅山脉、大喜马拉雅山脉和西藏喜马拉雅山脉。本地带北接青藏高原, 南连恒河、印度河平原, 喜马拉雅山连绵成群的高峰挡住了从印度洋上吹来的湿润气流, 喜马拉雅山脉南坡雨量充沛, 植被茂盛, 而北坡的雨量较少, 植被稀疏, 形成鲜明的对比。

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

本文所使用的生态区数据来源于Olson等基于生物地理区划的全球陆地生态区划方案^[7], 数据来自于世界自然基金会网站(https://www.worldwildlife.org/biome-categories/terrestrial-ecoregions)。该数据对全球陆地、淡水和海洋三大主要生境类型进行生物多样性区域分析得出陆地包含26种栖息地类型(Habitat types), 用于描述全球具有相似的环境条件、栖息地结构和生物复杂性模式, 并包含相似的群落和物种适应性的不同区域。为了代表世界各大洲和海洋的独特动植物, 每种生境类型都进一步细分为7个生物地理区：埃塞俄比亚界、澳新界、东洋界、新北界、新热带界、大洋界、古北界。最后在每个生物地理领域内, 将物种丰富度、分类学唯一性、生态或进化现象、稀有性等作为参数, 选择确定代表主要生境类型的867个生态区。

国家公园和自然保护区位置、面积、设立时间、级别等数据来自于研究区涉及的7个国家政府网站的公开数据, 包括中国生态环境部(http://www.mee.gov.cn/)、阿富汗野生动物保护组织(https://afghanistan.wcs.org/)、巴基斯坦野生动物网(http://www.wildlifeofpakistan.com/)、联合国教科文组织世界遗产中心(https://whc.unesco.org/en/list/)、印度政府环境、森林和气候变化部自治机构(https://wii.gov.in/)、尼泊尔国家森林与环境部(http://www.dnpwc.gov.np/)、不丹森林和公园服务部(http://dofps.gov.bt/)、世界保护区数据库(https://www.protectedplanet.net/)等。

2.2 研究方法 2.2.1 空间分布均衡性分析

(1) 地理集中指数

用来衡量保护地在各个国家之间的集中程度。计算公式为^[22]：

(1)

式中, E为地理集中指数, P_i为i国家内保护地的个数, S为保护地的总数, n为研究区内国家的个数, 即n=7。E值越大, 保护地分布越集中。若研究区中的保护地完全平均分布于各个国家内, 则其集中指数为：

(2)

(2) 不平衡指数

不平衡指数采用来衡量保护地各规模级分布的齐全程度, 采用罗伦兹曲线中计算集中指数的公式, 计算式为^[23]：

(3)

式中, B为不平衡指数；C_i为各研究单元保护地的数量占总样本数量的比重, 把各个研究单元的比重从大到小排列后求累计百分比, C_i为排在第i位的累计百分比；n为研究单元的个数。如果保护地平均分布在不同等级中, 则B=0；如果保护地分布极不平衡, 集中在一个规模级, 则B=1。

2.2.2 空间核密度估计法

核密度估计法的基本思想在于地理事件在空间点密度大的区域发生的概率大, 在空间点密度低的区域发生的概率低^[24]。核密度估计法的原理是以P点为圆心, 阈值r为半径, 统计以P、r为圆的范围内保护地分布的数量, 并除以圆的面积。一般以P点为圆心的保护地的核密度为:

(4)

式中, n表示距离尺度范围内所包含的内保护地分布的数量；K(·)表示核密度函数。h表示距离阈值, 即核密度估计法的尺度, d(x, x_i)表示两点之间的欧式距离。

本文的核密度估计借助ArcGIS 10.7的密度分析工具, 对保护地空间分布进行分析。

3 结果与分析 3.1 保护地与生态区关系紧密, 跨界特征明显

喀喇昆仑山-喜马拉雅山脉地区是全球两大重要生物地理区, 古北界和东洋界的过渡地带。根据Olson等对陆地生态区的分类^[7], 在古北界中, 主要为温带针叶林和山地草原和灌丛二大群落, 合计有7个生态区, 面积为572058.50 km², 占总面积的63.36%。在东洋界中, 则有热带和亚热带湿润阔叶林、热带和亚热带针叶林、温带阔叶和混交林、温带针叶林等四大群落, 合计7个生态区, 面积为330785.26 km², 占总面积的36.64%。在合计面积902843.76 km²的范围内分布有55处国家公园和自然保护区(表 1, 表 2, 图 1)。高海拔造成植被垂直带分带明显, 自下而上从热带、亚热带、温带以至寒带等各种景观交替出现；在水平方向上, 喜马拉雅山脉的南北向山谷使得印度洋暖湿气流得以通达, 南北动植物也借此迁移运动, 形成交汇过渡地带。高海拔的山地使得本地区植被类型极其丰富, 为各类动植物生存栖息提供了多样化的生境, 成为全球重要的生物多样性热点地区。

喀喇昆仑山-喜马拉雅山脉地区是全球保护地的集中区域。对中国、阿富汗、巴基斯坦、印度、尼泊尔、不丹和缅甸等7国在此地区的保护地进行统计分析, 可以看到这个地区已经形成一个举世闻名的保护地集聚带(表 3, 图 1)。初步统计共有55处国家公园和自然保护区, 合计面积为159063.30 km², 占总面积的17.62%。这55个保护地全部覆盖14处生态区, 其中有22处涵盖2—5个生态区, 占总数的40%, 某些保护地更是跨越了古北界和东洋界两大生物地理区单元, 充分显示本区域生态环境复杂多样性, 过渡性强的特征。从各保护地中生态区的数量统计看(表 4), 同时跨5个生态区有尼泊尔安纳普尔纳保护区, 其正处于东、西喜马拉雅生态区交汇处；跨4个生态区的有印度干城章嘉国家公园、不丹旺楚克百年国家公园和吉格·辛格·旺楚克国家公园；跨3个生态区的保护地数量最多, 达到10个, 其中尼泊尔数量最多达到6个；跨2个生态区的保护地为8个, 中国占其中3个。

喀喇昆仑山-喜马拉雅山脉地区55处保护地中, 有25处与其他国家毗邻(表 5), 占总数的45.45%, 其中中国塔什库尔干野生动物自然保护区和阿富汗瓦罕国家公园2个保护地同时与3个国家相邻；印度干城章嘉国家公园、尼泊尔阿皮南帕自然保护区、尼泊尔干城章嘉自然保护区、不丹吉格梅·凯萨尔严格自然保护区等4个保护地同时与二个国家相邻；19个保护地与一个国家相接, 其中中国的珠穆朗玛峰自然保护区同时与尼泊尔的萨加玛塔国家公园、马卡鲁贝润国家公园、朗塘国家公园、高里三喀自然保护区和马纳斯鲁自然保护区等5个保护地接壤, 形成由6个保护地组成的规模巨大的世界第一高峰保护地集群。

3.2 保护地在不同生态区内密度差异较大

喀喇昆仑山-喜马拉雅山脉地区自然环境的复杂多样, 造就了多个生态区跨越不同国家, 不同国家在这些生态区中典型地段布局保护地的格局(表 6)。在14个生态区中, 共有3个布局了10处及以上保护地, 只有3个功能区仅有1处保护地。从每万平方千米生态区内保护地数量来看, 东喜马拉雅亚高山针叶林(IM0501)生态区密度最高, 包括印度、尼泊尔、不丹3个国家的14个保护地, 每万平方千米达到5.10个；喜马拉雅亚热带阔叶林(IM0115)生态区和东喜马拉雅阔叶林(IM0401)生态区密度次之, 均达到1.57个/万hm²；东喜马拉雅高山灌丛和草地(PA1003)生态区则涵盖5个国家的16处保护地, 数量最多, 密度达到1.32个/万hm²。每万平方千米内保护地不到0.5个的生态区分别为：东北喜马拉雅亚高山针叶林(PA0514)生态区覆盖中国雅鲁藏布大峡谷自然保护区和察隅慈巴沟自然保护区, 密度0.43个/万hm²；雅鲁藏布干旱草原(PA1022)生态区覆盖中国珠穆朗玛峰自然保护区, 密度0.17个/万hm²；怒江澜沧江峡谷高山针叶林和混交林(PA0516)生态区覆盖中国高黎贡山自然保护区, 密度0.12个/万hm²。这些保护地密度过低的生态区应该是今后着重布局保护地的区域。

3.3 保护地空间结构呈现集聚性分布

将研究区中的保护地按照建立时间共分为3个批次, 分别是第一批次(1936—1984年)、第二批次(1985—1999年)和第三批次(2000—2014年), 分别计算喀喇昆仑山-喜马拉雅山脉地区的国家公园、自然保护区和保护地的地理集中指数以及不均衡指数, 结果如表 7所示, 表中的保护地的数量为国家公园与自然保护区数量的总和。

国家公园与自然保护区的地理集中指数(E)都大于完全平均分布值(37.796), 说明研究区各类保护地分布不均衡, 且呈集中分布的趋势, 但自然保护区的集中指数都在65以上, 分布不均衡度大于国家公园。而第一批次的自然保护区地理集中指数更是达到了峰值100, 因为1936—1984年期间只有中国建立了2个自然保护区。

从不均衡指数(B)看, 保护地的不平衡指数大多数在0.36—0.50波动, 说明保护地在喀喇昆仑山-喜马拉雅山脉地区空间分布不均衡。第一批次的自然保护区不均衡指数为1, 是因为在第一批次中仅中国建立了2个自然保护区, 与其他5个国家比不均衡；第二批次的自然保护区B值为0, 因为中国和尼泊尔各建立了3个自然保护区, 平均分布在两个国家。在喀喇昆仑山-喜马拉雅山脉地区, 国家公园主要分布在印度、尼泊尔两个国家, 自然保护区主要分布在中国和尼泊尔。总体看来, 研究区中的不均衡指数为0.45, 还是呈集中分布的趋势。

图 2中显示了喀喇昆仑山-喜马拉雅山脉地区保护地的空间分布情况。将保护地空间分布的核密度分析共分为9个级别, 从绿色到红色集中程度逐渐增加, 其中绿色区域为密度最小地区, 红色区域为密度最高地区。从图中可以明显看出喀喇昆仑山-喜马拉雅山脉地区的保护地主要集中分布在三个区域(核密度在6.0512×10^-5—9.0768×10^-5之间), 分别是不丹境内, 共包含6个保护地(5个国家公园和1个自然保护区)；尼泊尔-中国边境环珠峰地区, 包含11个保护地(7个国家公园和4个自然保护区)；印度-中国边境普兰、札达县, 包含12个保护地(11个国家公园和1个自然保护区)。而缅甸与中国藏南及云南部分所分布的保护地核密度最低。巴基斯坦、阿富汗与中国新疆塔什库尔干县所在地区中的保护地的核密度中等。总的来说, 喀喇昆仑山-喜马拉雅山脉地区保护地呈集中分布的趋势, 具有很明显的不均衡性。

3.4 各国采取灵活策略, 设置多样化的保护地

1936年, 印度成立喀喇昆仑山-喜马拉雅山脉地区第一个国家公园-吉姆库柏国家公园(, 从20世纪60年代开始, 各国在边境地区纷纷开始建立各类保护地, 至今已有55处(表 1), 其中国家公园有41处, 占总数的74.55%, 其余14处为自然保护区。在41处国家公园中, 最多的是印度, 计有15处, 尼泊尔有11处, 巴基斯坦有8处, 不丹有5处, 阿富汗和缅甸各有1处, 中国尚没有国家公园。截止2020年3月, 该地区共有5处世界遗产, 涉及6个国家公园。除印度干城章嘉国家公园是自然和文化双重遗产外, 印度大喜马拉雅山脉国家公园、花谷国家公园、楠达戴维国家公园和尼泊尔萨加玛塔国家公园、奇特旺国家公园均为世界自然遗产。

中国是本地区最早成立自然保护区的国家, 重点保护青藏高原特有的高寒自然生态系统和生物多样性。1983年在新疆喀喇昆仑山和帕米尔高原地区建立了塔什库尔干野生动物国家级自然保护区, 1988年在西藏建立了跨越定日县、定结县、聂拉木县和吉隆县等四县的珠穆朗玛峰国家级自然保护区。此后又建立了主要保护山地亚热带森林生态系统及扭角羚(Budorcas taxicolor)、孟加拉虎(Panthera tigris tigris)等濒危动物的察隅慈巴沟自然保护区, 主要保护喜马拉雅山脉东端南坡山地垂直带带谱及野生动植物的雅鲁藏布大峡谷自然保护区^[25]。在现有的8处自然保护区中, 国家级自然保护区有6处。虽然数量不多, 但总面积达到70227.93 km², 是喀喇昆仑山-喜马拉雅山脉地区保护地面积最大的国家。

印度的15处保护地成立于1936—2010年, 形成了以“国家公园为核心区, 野生动植物保护区、保护预留地、社区保护地、森林保护地为缓冲区或廊道”的保护地体系^[26]。截止2019年7月, 印度已建有104个国家公园、551个野生动物保护地、88个保护区、127个社区保护地, 共计870个自然保护地^[27], 本论文梳理了印度在喀喇昆仑山-喜马拉雅山脉地区的国家公园。印度国家公园总面积为14329.24 km ², 位列喀喇昆仑山-喜马拉雅山脉地区第四, 占9.01%, 其中还有4个被列入联合国教科文组织世界遗产名录。

自1973年成立第一个国家公园——奇特旺国家公园以来, 尼泊尔建立了20个由“国家公园—野生动物保护地—自然保护区—狩猎保护区”组成的保护地网络。喀喇昆仑山-喜马拉雅山脉地区16个保护地由11处国家公园和5处自然保护区组成, 面积为26910.02 km², 占16.92%。其中萨加玛塔国家公园、奇特旺国家公园为世界自然遗产。尼泊尔保护地设立的特点就是集中在尼泊尔和中国边境的喜马拉雅山区, 相互毗邻, 单个面积不大, 便于灵活管理。既相互衔接, 又独立管理, 与珠峰北侧中国境内覆盖四个县的珠穆朗玛峰自然保护区相比, 类型多样, 既有以严格保护为主的自然保护区, 也有可供开展生态旅游的国家公园, 同时由于面积都不大, 对周边城镇影响较小, 自然保护具有更好的操作性。

4 结论

通过上述分析, 可以看到喀喇昆仑山-喜马拉雅山脉地区各国经过多年的建设, 已经初步建成类型多样, 设置灵活的自然保护地体系, 中国-尼泊尔交界处的珠穆朗玛峰地区, 以及中国-巴基斯坦交界处的喀喇昆仑山地区是本地区自然保护地分布较为集中的区域, 具备较好的保护地跨境合作基础。2010年颁布的《中国生物多样性保护战略与行动计划(2011—2030年)》中, 提出中国将在新疆红其拉甫山口、西藏珠峰等地区研究建立跨国界保护区。2011年, 为了应对喀喇昆仑-帕米尔跨界景观保护面临的挑战, 保护马可波罗盘羊(Ovis ammon polii)这一帕米尔高原的特有种和旗舰种, 中国和巴基斯坦学者共同提出《促进中巴跨界合作的喀喇昆仑山-帕米尔景观规划纲要草案》, 以保护生活在两国毗邻保护地中野生动物^[28]。

2019年中国政府发布《关于建立以国家公园为主体的自然保护地体系的指导意见》, 明确将从2020年设立国家公园, 并将自然保护地按生态价值和保护强度高低依次分为国家公园、自然保护区和自然公园, 形成完整的自然保护地体系。中国也希望借鉴国际经验, 注重与国际自然保护体系对接, 积极参与全球生态治理, 共谋全球生态文明建设^[29]。同年, 西藏宣布启动地球第三极国家公园群建设, 青藏高原将规划布局一批以保护具有国家代表性的自然生态系统的国家公园群, 对于加强青藏高原生态安全屏障建设, 提高当地居民福祉水平、促进区域绿色发展、增强可持续性能力, 促进青藏高原进入高质量的自然和人文生态系统良性循环过程具有重要意义^[30]。因此未来中国将在边境地区加强国家公园建设, 也有利于进一步推进国家公园跨境合作。通过跨境合作, 进一步加强喀喇昆仑山-喜马拉雅山脉地区各国的生物多样性和生态环境保护, 更好地应对全球气候变化, 更好地落实联合国2030年可持续发展目标, 促进当地社会经济发展。