人工草地是指利用综合农业技术, 在完全破坏了天然植被的基础上, 通过播种建植新的人工草本群落, 以生产补充饲料, 解决饲料季节不平衡的牧用草地^[1]。为解决草畜矛盾和畜牧业发展的瓶颈制约, 在农区或者半农半牧区的边际土地和低产耕地也常被用于建立人工草地。按照牧草的生长年限, 可划分为一年生人工草地和多年生人工草地。目前, 发展人工草地已成为缓解牧区天然草地放牧压力和提升草地生产力的有效途径, 其高效集约的生产方式, 不仅为牧区畜牧业的稳定发展提供所需的优质牧草, 而且也有助于保护和恢复退化的天然草地^[2-5]。未来畜牧业可持续发展的一项重要措施就是建立集约化经营的人工草地。在畜牧业发达的澳大利亚、新西兰和欧洲等国家, 人工草地占草地总面积的比例已超过50%, 人工草地集约化经营已成为低成本, 高效益的畜牧典范模式。而中国的人工草地仅占草地总面积的约5%, 西藏的人工草地占比就更低了。据2018年西藏自治区农牧厅的最新发布数据, 人工草地面积9.9万hm^2[6], 占可利用草地面积7713万hm²的0.13%^[7], 且主要集中在海拔低于4000 m的藏南农牧结合区域^[8-9]。现阶段我国人工草地建设中依然面临种植规模小、牧草产量低、持续性和稳定性差等问题^[10], 且大部分人工草地建设的实践研究局限于内蒙古和新疆等低海拔牧区^[11-14], 而对于高寒牧区人工草地的布局建设方面探讨较少, 目前开展较多的实验研究主要涉及青藏高原人工草地群落特征、稳定性、持续性和管理技术等^[15-19]。

西藏北部(藏北高原)是西藏自治区重要的牧区, 目前草地畜牧业的发展依然是以天然草地为饲草来源的传统的粗放型畜牧业模式为主。作为西藏今后畜牧业的重点发展区域, 牲畜对饲草料的需求巨大。但由于水热条件的限制, 以高寒草地生态系统为主的天然草地生态系统脆弱, 生产力不稳定, 易受气候变化的影响, 尤其是在降水量季节变化差异大的干旱、半干旱区, 天然草地生产力的不稳定性严重制约着草地载蓄量^[20-23], 加之近些年受气候变化和放牧压力的影响, 草地植被群落结构破坏和生物量减少, 草地呈现普遍退化状态^[24-27], 面临天然草地草畜平衡调控能力严重失调的巨大压力。解决西藏天然草地草畜平衡失调问题是遏制草地退化和构建高原生态安全屏障的关键所在, 要缓解天然草场压力, 除加强畜牧业生产体系自身机制建设外, 还应多途径、多渠道地组织饲草饲料生产^[28]。例如建立人工草地、修建牲畜暖棚并进行牛羊舍饲育肥、退化草地恢复与重建、优化畜群结构和放牧强度等, 都是实现高寒草地畜牧业经营持续发展的有效途径^[29]。通过建立人工草地和发展草地畜牧业, 可以使因长期过牧而普遍退化的天然草地得以退牧还草, 恢复退化天然草地的生态与生产功能^[30-31]。但是, 人工草地建设发展需注重人工草地区域布局、牧草种植区划、人工草地与家畜饲养耦合和人工草地经营管理等战略问题^[32]。

藏北高原由于寒、旱气候条件的约束, 在牧区发展大规模的人工草地是不切实际的, 适宜建设人工草地的土地资源也极其有限, 而且, 目前部分现有人工草地存在选址布局不合理, 种植模式单一, 稳定性差, 缺乏有效管理等问题。因此, 不同于低海拔区域人工草地建设, 藏北高原人工草地布局考虑的条件应更为严格, 不仅要考虑牧草品种生物学特性, 而且还需评价人工草地选址的适宜性及其潜在的环境影响。根据野外调查和青藏高原目前人工草地建设情况分析^{[15-18, 29]}, 藏北高原人工草地布局应该选择在海拔相对较低、土壤和水热条件较好、靠近乡镇点便于管理, 以及严重退化的草地作为潜在的建立区域, 不能盲目破坏天然草地和能够通过改良恢复的退化草地。此外, 根据牧草品种生理特性选择合适的牧草品种和种植模式同样是布局建设中需要重点考虑的问题。

因此, 本研究基于遥感数据、土壤数据、气象数据、地形和土地利用数据, 结合野外人工草地调查数据, 在分析了现有人工草地建设立地条件的基础上, 针对藏北不同区域人工草地布局建设的适宜性进行了探讨, 识别出区域适宜建设人工草地的潜在区域, 并与现有人工草地进行对比分析和验证, 以期为未来人工草地建设布局提供科学依据。

1 研究区概况

藏北高原又称“羌塘高原”(图 1), 南北宽达760 km, 东西长约1200 km, 面积59.5万km², 约占西藏自治区土地面积的一半, 该区域拥有广阔的天然草地资源, 约占西藏天然草地面积的59%, 是西藏主要的牧业发展区。藏北高原大部分地面高度在海拔4600—5100 m之间, 地势总倾向为西北高、东南低, 气候寒冷干燥, 大部分地区年平均气温低于0 ℃, 最暖月均温不及14 ℃, 最冷月均温-10℃以下。区域年平均降水量从东南部的300 mm左右, 逐渐减至西北部的50 mm以下, 且降水主要集中在6—8月, 夏季日照时间长, 太阳辐射强^[33]。藏北高原土壤主要以莎嘎土、冷漠土和寒漠土为主, 有机质含量低。从东往西, 草地植被类型包括以高山嵩草(Kobresia pygmaea)、矮嵩草(Kobresia humilis)和青藏苔草(Carex moorcroftii)为优势种的高寒草甸；以紫花针茅(Stipa purpurea)、丝颖针茅(Stipa capillacea)和座花针茅(Stipa subsessiliflora)为优势种的高寒草原；以及以羽柱针茅(Stipa subsessiliflora)、沙生针茅(Stipa glareosa)和驼绒藜(Ceratoides lateens)为优势种的荒漠草原^[34]。

2 高寒地区人工草地布局原则和方法 2.1 高寒地区人工草地布局原则

考虑到藏北地区特殊的地理环境条件, 本研究选择适宜高寒地区开展人工种草项目建设的区域需满足以下原则(表 1)。

2.2 野外调查

于2017年7月起对藏北各县域现有人工草地进行野外调查, 调查内容主要包括现有人工草地面积、位置、海拔高度、立地条件、种植牧草品种及其长势等。为分析现有人工草地的立地条件, 在研究区内不同区域选取了20块面积规模在60 hm²以上的人工草地(图 1)筛选潜在人工草地适宜性条件。

2.3 数据来源与处理

人工草地布局建设中, 要根据区域气候条件、水热状况、土壤条件、人工草地类型以及社会经济条件等来确定人工草地类型及其面积^[35]。因此, 结合野外调查和高原人工草地布局原则, 本研究选取土壤、地形、气候、植被、土地利用类型五个方面的影响因素(表 2), 各个参数通过ArcGIS 10.2进行处理和提取。

2.4 植被覆盖变化趋势

人工草地建设不能破坏天然草地, 只能以地势较低的退化草地建设。采用2000—2015年的MOD13A3植被指数来反演植被覆盖度, 为反映植被覆盖变化趋势, 利用趋势变化公式, 如下所示, Slope表示影像中每一个栅格单元的植被覆盖变化趋势, > 0为增长, ＜0为退化趋势, n为研究时间段(2000—2015年, 16), NDVI_i表示研究期间第i年的植被覆盖度。本研究选择以Slope＜0植被处于重度退化状态的区域。

在具体操作中, 建立了NDVI与群落盖度之间的相关关系, 选择退化且群落盖度＜20%的重度退化草地。

2.5 空间分析法

本研究根据选取的影响因素的各个参数作为人工草地布局建设适宜性分析的指标, 首先利用ArcGIS 10.2的空间分析模块提取和分析每个参数指标所确定的人工草地布局建设适宜的空间分布, 然后基于各个参数的适宜性条件(表 3), 将各个参数指标的适宜性空间分布图进行叠加, 得到综合的人工草地在藏北的适宜性空间分布图, 最后根据县级单元的边界提取人工草地潜在的适宜性面积, 和现有人工草地分布位置和面积进行对比分析。

3 结果分析 3.1 藏北高原现有人工草地建设立地条件分析

根据实地调查的20块60 hm²规模以上的现有人工草地, 提取了海拔、坡度、土地利用类型、土壤有机质、土壤碎石百分比、水源距离等参数值(图 2)。结果表明藏北地区现有人工草地均位于较高的海拔地区, 海拔最高位于4761 m, 最低位于4249 m, 平均海拔4546 m。考虑灌溉条件, 现有人工草地均选择离水源较近的位置, 或者离县城较近方便灌溉的区域, 最远有40 km, 最近1 km, 平均距离13 km。适宜建立人工草地的土壤属性主要考虑0—30 cm土层的土壤碎石百分比和土壤有机质, 现有人工草地土壤有机质差别不大, 均值1.29 g/100g, 而土壤碎石比最高27.8 g/100g, 最小值1.6 g/100g, 均值16.6 g/100g, 从土壤条件的选择可以看出, 现有人工草地的建立通常选择土壤碎石比较低且利于开垦的区域, 而土壤有机质均不高与区域地质条件有关。此外, 提取土地利用类型和坡度发现现有人工草地均建立在草地覆盖度低的平坦地势区域, 坡度均小于25°, 以利于耕作和灌溉保墒。因此, 筛选人工草地潜在区域时主要考虑选择低覆盖度的退化草地和地势平坦的位置。

按照牧草生育期的≥10℃积温最低的燕麦500℃^[37], 人工牧草的建植的主要区域在尼玛县及其以西的地区(图 3), 而青稞、箭筈豌豆、紫花苜蓿和披碱草的≥10℃积温都高于1500℃, 分别是1500℃、1700℃^、1700℃和1700℃^[38-40], 藏北高原的活动积温能满足这些牧草种的地方微乎其微。因此, 本文以活动积温要求最低的燕麦为例来匡算人工草地建植的潜在最大面积。

据2017年调查, 藏北高原人工草地普遍种植的是一年生燕麦, 而多年生的披碱草和饲料用的青稞种植面积非常少。那曲县、尼玛县和普兰县种植的燕麦鲜草产量普遍达到22500 kg/hm²以上。而那曲县那木切乡披碱草和燕麦混播的产草量约22500 kg/hm², 青稞鲜草产量为18000 kg/hm²；班戈县披碱草鲜草产草可达到27000 kg/hm², 青稞18000 kg/hm², 而普兰县的披碱草鲜草产量只有11250 kg/hm²。可见, 藏北高原燕麦产草量通常高于披碱草和青稞等牧草, 由于披碱草人工草地积温需求较高, 种子不能成熟, 还存在越冬的冻害等问题, 且3年以上的披碱草草地容易退化, 因此, 藏北普遍种植的是燕麦人工草地。尽管一年生燕麦种植需要的积温要求比其他牧草要低许多, 每年种植需要较高的劳动力投入, 但产草量高且比较稳定, 牧民还是普遍愿意种植燕麦。

3.2 藏北高原潜在适宜人工草地建设区域

根据人工草地建立适宜性条件分析, 确立藏北人工草地建立的基本条件(表 4), 将各个参数指标筛选出的适宜区进行叠加, 综合得出藏北潜在适宜建设人工草地的区域面积(图 4)。在本研究野外调查期间, 发现藏北种植最多的牧草品种为燕麦, 以≥10℃积温需求最低的燕麦为例, 适宜藏北区域人工草地建设的土地极其有限, 只有五个县域(尼玛县、改则县、革吉县、日土县、噶尔县)有具备人工草地建设条件的土地, 且绝大部分位于阿里地区, 那曲地区只有尼玛县有少数符合建设条件的土地, 五个县域中以尼玛县的面积最多, 改则县次之, 噶尔县最少。根据面积统计, 区域符合人工草地建设条件的潜在面积为15333 hm², 尼玛县最多6067 hm², 改则县2933 hm², 最少的噶尔县为1333 hm²。此外, 各县域潜在人工草地与现有人工草地面积对比(图 5)表明, 除了噶尔县现有人工草地超过潜在人工草地面积外, 其他4个县域均存在多余尚待开发的面积。另外, 若以其他四种常种牧草(青稞、箭筈豌豆、紫花苜蓿和披碱草)的活动积温需求提取, 只有札达县和噶尔县少数区域符合条件, 但这些区域却不符合人工草地建立的其他需求。因此, 不同牧草品种活动积温需求是影响藏北人工建立的关键限制因素。

3.3 藏北高原人工草地布局建设适宜性验证——以改则县为例

通过野外调查和现有人工草地分布数据收集, 对比分析了改则县现有人工草地和潜在人工草地分布区域的适宜性(图 6)。改则县现有人工草地共1000 hm², 地块分布较为零散, 主要分布在改则镇、洞措乡、物玛乡和麻米乡, 最大的面积为148 hm², 最小的面积1.78 hm², 面积超过60 hm²的人工草地只有7个地块, 其余的均为规模较小的地块。此外, 现有人工草地只有少数几块分布在海拔4500 m较低的位置, 其余大部分海拔在5000 m左右, 且离水源地较远。相反地, 根据适宜性分析筛选出的潜在人工草地分布主要集中在海拔4500 m较低的区域, 离水源地近, 有道路可以通达, 靠近乡镇点便于管理, 且大部分地块面积规模较大, 分布集中。相较于现有的人工草地, 筛选出的潜在人工草地分布位置更适宜作为今后发展建设的区域。

4 讨论

本文利用气象、土壤、地形、土地利用、遥感数据并结合野外调查, 从可利用土地资源的角度考虑, 综合分析藏北高原牧区人工草地建设布局的适宜性, 识别出适宜人工草地建设的潜在分布位置和面积。研究发现, 影响该区人工草地布局最重要的因子是水热条件和海拔, 且在综合考虑灌溉条件、交通的可达性和低盖度退化草地等人工草地布局综合约束条件下, 能满足人工草地建设条件的适宜区域极其有限。

影响人工草地布局的基础为土地适宜性, 其中热量和水分两大因子尤其重要^[41]。低温和生长季节短是限制高寒草地植被生产力的主要限制因子^[42-43], 对牧草的生长物候和生产力具有重要影响。本研究以燕麦的积温需求为例选取温度适宜区, 发现藏北高原大部分区域, 尤其是那曲地区无法满足人工牧草种植的积温需求, 只有阿里地区的五个县能满足燕麦生长的活动积温要求。水分条件也是影响牧草生产力的关键因子, 区域草地的产草量与降水量直接相关^[44]。西部的阿里地区虽满足积温需求但由于降水量低, 气候干旱, 满足灌溉条件的水源点也较少。因此, 从水热条件来看, 藏北高原牧区发展人工草地面临二难境地, 即东部的那曲地区降水较丰富, 但积温较低, 而西部的阿里地区虽然积温达到要求, 但降水显得严重不足。水热条件的不协调, 都会限制人工草地的发展及其生产力。藏北高原地处高寒牧区, 气候呈现寒、旱的特征, 水热条件的制约成为人工草地的首要约束条件。从目前对积温要求较低的燕麦人工草地种植分布区分析来看, 藏北高原牧区种植面积潜在面积只有15333 hm², 难以达到西藏自治区初步在该区每个县建立人工草地1000 hm²的目标, 而且适宜燕麦种植的还在阿里的5个县, 这些县的降水量少, 灌溉条件比较难以满足, 需要通过修建水渠或者引夏季融雪的雪水灌溉, 成本很高。因此, 从人工草地经济可行性和成本效益角度来说, 藏北高原人工草地建设还面临诸多问题, 难以大面积发展人工草地。比较可行的是, 藏北牧区建设以家庭圈窝种草为主要形式的牧户家庭种草模式。同时在农区或者农牧交错带, 在保障粮食安全基础上建设冬饲料牧草生产基地, 通过区域农牧耦合的方式解决藏北牧区饲草问题。

在牧草种类选择方面, 披碱草虽然长期被认为是高寒牧区人工草地建设的首选牧草物种, 但从牧草品种的适应性来看, 该牧草的品种需要的积温较高, 在海拔普遍超过4500 m的藏北高原地区, 面临冻害的风险较高。而且, 披碱草人工植被在自然演替状态下短期内就会发生严重退化^[45-46], 在退化草地补播或者翻耕后种植的垂穗披碱草根系比较浅, 不耐践踏和啃食, 不能形成长期稳定的放牧群落。其他牧草种类也同样面临类似问题。相比之下, 燕麦虽然是一年生品种, 但在6—9月的生长季作为饲料作物种植, 能获得较高的产草量, 在藏北地区受到牧民欢迎, 种植面积也比其他牧草品种广得多。当然, 种植一年生牧草, 频繁地耕作、投入成本、产出效益和可持续性问题都是需要慎重考虑的。否则, 人工草地生产力下降至天然草地水平, 投入和产出效益不成正比这都是没有实践价值的。

适宜不同牧草品种生长所需的积温和水分条件不同, 因此, 未来人工草地建设布局时需考虑牧草品种的抗寒性和抗旱性, 应结合牧草种植区划, 选取适合的种植模式。另外, 本研究选取植被指数反演的覆盖度作为反映草地退化的指标, 植被覆盖单独作为反映草地退化的指标太单一, 草地退化不单是覆盖度的变化, 还包括群落结构的退化, 今后需要结合土壤和植被群落结构变化来综合反映草地退化^[25], 选择严重退化草地作为发展建设人工草地的区域, 而不是可改良和恢复的退化草地。

总得来说, 藏北高原牧区由于地理环境条件的限制, 草地生产力十分有限, 且适宜发展人工草地的潜在区域也极其有限, 从草地生态系统服务功能和国家生态安全建设角度来看, 其供给服务发展潜力有限, 国家也将藏北高原定位于高原生态安全保障。要解决藏北草畜间的矛盾, 关键是草料的供给, 应选择海拔低, 水热条件好的农区作为发展人工草地的重点区域。西藏“一江两河”(雅鲁藏布江及其支流拉萨河、年楚河)流域的农区, 由于较低的海拔和良好的水热条件, 土地资源相对丰富, 与藏北牧区人工草地建设适宜性相比, 有着巨大的潜力。但农区由于粮食保障的需求, 土地资源长期以来以耕地利用为主, 人工草地建设未得到足够重视, 粮草争地的矛盾长期存在, 还未充分发挥农区草地的最大效益^[47]。面对西藏不同区域草地生态系统服务供给和需求存在的空间弹性差异, 只有因地制宜地调控供给和需求差异才能使得生态系统服务效用最大化。因此, 如何合理有效地将以藏北为代表的牧区和“一江两河”流域的农区相结合, 发展区域农牧耦合是今后西藏畜牧业发展的重要战略任务^[48-49]。

5 结论

藏北高原人工草地布局最重要的影响因子是水热条件和海拔, 决定了人工牧草种植的适宜性；其次, 灌溉条件和交通的可达性是人工草地建设可行性需要重点考虑的要素；此外, 人工草地只能建设在低盖度的严重退化草地上, 不能破坏天然草地和可以改良恢复的轻度退化草地。在多因素的约束下, 藏北高原能满足人工草地建设条件的适宜区域极其有限, 以生长季较短、活动积温受限最低的一年生燕麦种植为例, 其潜在的适宜区总面积约15333 hm², 而且主要集中分布在降水条件较差的阿里地区, 考虑引水灌溉的建设经济成本, 要想建立大规模成片的人工草地是非常难的。对于水热因子限制条件来说, 东部的那曲地区降水较充沛, 主要受积温不足的限制；而西部的阿里地区则主要受到水分的限制, 需要加强灌溉设施建设。因此, 要在平均海拔超过4500 m的藏北高寒牧区建立人工草地必须慎重, 人工草地建设应该重点考虑海拔较低的“一江两河”地区。今后人工草地种植规划还需要加强牧草的抗寒性和抗旱性研究, 在人工草地的管理方面要特别关注已建人工草地的可持续性、稳定性、生产力和牧草质量, 防止出现草地退化和沙化等问题。