草地生态系统是中国面积最大的陆地生态系统类型, 占国土面积的29.98%, 也是中国重要的生态安全屏障区, 在土壤保持、防风固沙和水源涵养的方面发挥着重要作用^[1]。草地也是我国北方、西北和青藏高原等地区畜牧业发展的基础, 是当地居民生产生活的重要资源。由于社会历史的原因, 草地生态系统遭到掠夺式开发, 造成草地退化, 草地载畜量下降, 居民承受巨大的经济损失^[2]。因此充分了解各地区草地生态系统的承载状况及演变趋势, 对于我国现阶段及未来的草地生态系统恢复和科学规划畜牧业发展格局具有重要意义。

青藏高原是我国重要牧区之一, 全国264个牧区与半牧区县中有133个位于该地区^[3]。同时, 该地区作为我国“两屏三带”生态系统格局的重要组成部分^[4], 对于区域和国家的生态安全战略具有重要意义。草地生态系统也是青藏高原主要的生态系统类型, 对青藏高原草地承载力展开研究, 可充分了解青藏高原畜牧业的发展状况及潜力, 并对该地区的草畜平衡、草地生态系统保护和区域乃至国家的生态安全提供必要支撑^[2]。

草地承载力是指区域的草地资源在可持续发展的前提下, 所能承载的最大牲畜数量, 是草地牧业发展的重要指标, 对区域草地的开发与保护具有重要意义, 具体表现为草地的理论载畜量等^[5]。测算产草量是计算草地承载力的基础。草地承载力的研究需要通过草地产草量间接计算牲畜承载量。产草量的计算方法包括实测产量法^[6]、模型评估法^[7—8]和遥感评估法。其中, 实测产量法的结果准确, 适用于局部地区的研究使用, 现有研究多将其用于验证其他两种方法^[9—10]。模型评估法包括气候模型和过程模型, 适用于定性评价大尺度区域的产草量地带性分布规律或已有长期试验站区域的研究^[11]；遥感评估法包括NDVI法^[12—13]和NPP法^{[5, 14—15]}, 两种方法均适用于大尺度范围的研究^[16—17], 其中NPP法是利用改进的CASA模型计算产草量, 该模型考虑了光能利用率, 机理性较强, 适用于缺乏大量实测数据的地区^{[5, 18]}。

目前针对青藏高原草地承载力的研究多以局部地区为主^[19—27], 如曹娟等^[25]发现在草畜平衡政策实施之前的青海刚察县处于超载状态, 政策实施后则处于草畜平衡状态；Zhang等^[26]发现, 三江源东部和南部县虽然承载力较强, 但牲畜存栏量过高, 这些地区仍处于超载状态；赵卫等^[27]发现, 西藏地区整体处于超载状态, 且未超载县的草地承载力盈余越来越少。为充分了解青藏高原全域草地承载力的空间变化特征并对其进行预警, 根据青藏高原独特的地理条件, 且缺乏大范围实测数据的情况, 本文选择NPP法展开研究。

本文以多年NPP数据为基础, 以青藏高原全域为研究区域, 选取理论载畜量和草畜平衡指数, 展开草地承载力评价, 主要为了(1)分析青藏高原2000—2015年间草地生产力空间格局及其变化特征, (2)结合牲畜存栏量分析县域草地承载状况, 并对其进行预警, 旨在为青藏高原草地保护的优先级划分, 为畜牧业的发展规划等提供科学依据。

1 研究地区与研究方法 1.1 研究区概况

青藏高原位于亚洲内陆, 是地球上最独特的地质-地理-生态单元, 是中国面积最大, 世界海拔最高的高原, 平均海拔4000 m以上, 被称为“世界屋脊”^[28], 位置介于26°00′—39°47′N, 73°19′—104°47′E之间, 总面积约278.68万km²。草地生态系统总面积158.82万km², 占青藏高原总面积的56.99%, 是该地区主要的生态系统类型(图 1)。青藏高原内的县级行政区有244个, 除新疆边界的8个不完全县外, 共有完整县级行政区236个。

1.2 研究方法

在已有计算方法的基础上^[5], 本研究进一步将计算公式进行提炼, 方法如下：

(1) 理论载畜量：

(1)

式中, CA为理论载畜量(羊单位); G_ij为行政区i第j类草地的全年干草总量(kg); C_ij为行政区i范围内牲畜对第j种草地所产牧草的利用率^[29]; UG为每个羊单位每天需要的干草量(kg/d), 数值为1.8; D为一年的天数, 数值为365。

(2) 全年标准干草产量

(2)

式中, NPP_ij为行政区域i范围内第j种草地年总NPP(kg)；A_ij为行政区域i范围内第j种草地总面积(km²)；S_ij为行政区域i范围内第j种草地的生物量到NPP的转换系数^[30]；SG_ij为行政区域i范围内第j种草地地下部分和地上部分生物量比例系数^[31]。

(3) 产草量模型测算精度分析

选取已发表论文中与研究时间相近的青藏高原地区产草量数据, 与本论文进行对比, 结果见表 1。本研究结果均在相应实测值的误差范围内, 表明研究方法适用于本研究区。

(4) 草畜平衡指数^[35]

(3)

BGLI是草畜平衡指数, A为牲畜存栏量, CA为理论载畜量。

当BGLI≤10%, 表示不超载；当10%＜BGLI≤15%时, 表示临界超载；当BGLI＞15%时, 表示超载。BGLI越高, 表示超载情况越严重。

(5) 草畜平衡预警

在草畜平衡指数分级的基础上, 以各县2000—2010年、2010—2015年两个时间段的平衡指数变化为依据, 分为持续改善型、波动改善型、波动恶化型和持续恶化型, 再将承载状况进行预警等级划分, 未超载状态下的持续改善型和波动改善型分为绿色无警区, 波动恶化型和持续恶化型分为蓝色预警区；临界超载状态下的持续改善型和波动改善型分为黄色预警区, 波动恶化型分为橙色预警区, 持续恶化型分为红色预警区；超载状态下的持续改善型和波动改善型分为橙色预警区, 波动恶化型和持续恶化型分为红色预警区, 具体方法见表 2。

1.3 数据来源

(1) 净初级生产力数据：来源于中国科学院空天信息创新研究^[36], 数据基于CASA模型生产, 为空间栅格数据, 空间分辨率为250 m, 时间分辨率为月。

(2) 生态系统格局数据：来源中国生态系统评估与生态安全数据库(http://www.ecosystem.csdb.cn), 空间分辨率90 m。

(3) 牲畜存栏量数据：来源于青海^[37]、西藏^[38]、新疆^[39]、甘肃^[40]、四川^[41]、云南统计年鉴^[42]和中国县域统计年鉴^[43], 数据以县为单位。

2 结果与分析 2.1 青藏高原草地理论载畜量空间格局及变化趋势

草地生产力的评价指标有：干草产量和理论载畜量等^[5], 本文以理论载畜量作为草地生产力的代表性指标。

青藏高原理论载畜量整体呈现东高西低格局(图 2)。2015年总产草量达到11393.32万吨, 相应的总理论载畜量为9451.53万羊单位, 草地单位面积载畜量为59.51羊单位/km²。高寒草原和高寒草甸是青藏高原草地承载力的主要组成部分, 其理论载畜量占全区的64.99%, 分别为3429.08万羊单位和2713.12万羊单位；温性草原和温性草甸理论载畜量有1467.08万羊单位, 占全区干草总产量的15.52%(图 3), 其单位面积理论载畜量高于高寒草原和高寒草甸(图 4)。

与2000年相比, 2010和2015年的草地生产力均呈增长趋势, 总理论载畜量由8614.89万羊单位上升到9117.54万和9451.53万羊单位, 增长区域以青藏高原东北部和中东部为主, 下降区域以青藏高原中南部地区为主(图 2)。干草总产量由2000年的10403.42万吨到2015年的11393.32万吨。

2.2 青藏高原牲畜存栏量及变化趋势

2015年青藏高原牲畜存栏量合计15846.99万羊单位, 存栏量密度(牲畜存栏量与草地面积的比值)为99.78羊单位/km²。有46个县(19.49%)超过100万羊单位, 主要是传统的牧区与半牧区县, 分布在青藏高原中东部及西部诸县；有53县(22.46%)低于25万羊单位, 分布在青藏高原边缘诸县(图 5)。

与2000年相比, 2010和2015年牲畜存栏量呈先升后降的趋势(图 6), 由14490.52万羊单位增至16385.26万羊单位, 再降至15846.99万羊单位, 增长区域为青藏高原中北部、北部边缘和东南部边缘诸县(图 5), 下降区域为中东部和西南部诸县。

2.3 青藏高原草畜平衡指数及变化趋势

2015年青藏高原整体的草畜平衡指数为67.91%, 处于超载状态。2000和2010年草畜平衡指数也一直处于超载状态(图 7), 草畜平衡指数由67.88%增至79.90%, 再降至67.91%。2010年的超载情况最为严重, 主要原因是2000—2010年的牲畜存栏量增长过快(13.08%), 而理论载畜量仅增长了5.83%, 至2015年草畜平衡指数有所下降, 主要是因为牲畜存栏量下降(-3.28%)以及同期草地承载力上升(3.66%)。

2015年超载的县有191个(80.93%), 超载县集中分布在青藏高原东部和中南部诸县。与2000相比, 2010和2015年超载的县数量呈先稳定后降低的趋势, 由前两个年份的196个下降至191个(图 8)。2000—2010年和2010—2015年两个时间段中, 草畜平衡指数均增大的县有43个(18.22%), 集中在青藏高原中南部(以西藏为主)；草畜平衡指数均减小的县有42个(17.80%), 集中在青藏高原东部(图 7)。超载县在牧区与半牧区中有98个, 与其他县相近(表 3)。

2015年仅有41个县未超载(16.39%), 分布在北部、西部和南部部分县。与2000相比, 2010和2015年超载的县数量呈先减后增的趋势, 从39个减至34个, 再增至41个。未超载的县主要以牧区与半牧区县为主, 其他县的比重较低(表 3)。

2.4 青藏高原草地承载力预警

青藏高原草畜平衡指数预警结果如图所示(图 9)。与表 1对比发现, 实际的预警结果中, 临界超载状态只包括黄色预警区, 橙色和红色预警区均属于超载状态。

因青藏高原县域面积的差异巨大, 红色预警区诸县的面积共43.84万km², 绿色无警区诸县面积共86.90万km², 分别占青藏高原总面积(完整县域)的17.28%和34.25%, 前者约是后者一半, 但红色预警区的县域数量合计72个, 绿色无警区仅有32个, 前者是后者的2.25倍。且红色预警区多分布于青藏高原中南部(以西藏为主)和东部, 绿色无警区多分布于青藏高原西部及北部, 中东部地区的草地承载力远高于西部及北部, 呈现出草地承载力越强, 超载情况越严重的空间分布格局。绿色无警区主要分布在牧区与半牧区县(73.17%), 橙色预警区和红色预警区在两类县中的数量基本相当(表 4)。

红色预警区已经处于超载状态, 且超载状态还在加剧, 在未来畜牧业的管理中亟需重点关注；橙色预警区和黄色预警区的超载状况在改善中, 未来需要进一步保持这一趋势, 使之最终转为绿色无警区；蓝色预警区需要也需要得到相应的关注, 避免恶化为橙色或红色预警区。青藏高原绿色无警区, 是未来畜牧业调整发展的主要实施区域。

3 讨论与结论 3.1 讨论

基于CASA模型的草地净初级生产力研究目前已成熟, 在不破坏草地的情况下, 能够大范围的获得多年的草地生产情况, 在中国的实地验证效果较好^[44—45], 尤其在青藏高原地区, 其拟合优度达0.83^[46]和0.92^[47], 结合表 2的多个地区产草量的验证结果, 表明该模型及其结果适用于青藏高原地区。

青藏高原草地生产力分布格局与其自然本底条件和人口分布格局相一致^{[15, 48—49]}, 并在研究期间内相对稳定。青藏高原东部的草地生产力远高于西部及北部, 该地区与中国的中亚热带、北亚热带和南温带相毗邻^[50], 平均海拔低于4000 m, 平均气温为6 ℃, 年降水为650 mm ^[51]；西部及北部主要由羌塘高原、三江源西部和柴达木盆地等组成, 平均海拔超过4600 m, 平均气温差异较大, 降水量低于300 mm^[52]。集中分布于东部地区的高生产力草地面积为26.32万km², 而集中分布于青藏高原西部及北部地区的低生产力草地面积为53.49万km², 前者不足后者的一半, 使得该区域草地承载力相对较低。青藏高原草地生产力在研究时间内保持增长趋势, 主要原因可能与秦大河等人提出的在全球气候变化背景下, 中国西部地区由“暖干化”趋向于“暖湿化”^[53]有关, 但青藏高原草地生产力的格局仍然相对稳定。随着气温和降雨量的上升, 未来青藏高原地区草地生产力将进一步提升, 相应的草地承载力也会进一步增长。因此可充分利用该气候变化趋势, 积极扩展草地面积, 提高草地质量, 避免草地承载压力进一步加重。

青藏高原草地超载状况要得到有效缓解乃至消除, 必须依靠多种措施的结合。青藏高原的牲畜存栏量远高于本地区草地承载力, 其中牲畜存栏量是决定该地区草畜平衡指数的主要因素。青藏高原草畜平衡指数在研究时间内均处于超载状态, 且十分严峻, 呈现出草地承载力越强, 超载情况越严重的格局。虽然草地承载力增长了9.71%, 但同期牲畜存栏量增长了9.36%, 两者的增长幅度相近, 因此仅通过依靠气候变化以提高草地生产力而提高草地承载力的方法, 无法彻底解决青藏高原地区的超载状况。未来随着生态保护工程在青藏高原地区的进一步实施及居民收入要求的进一步提高, 将对该地区的草畜平衡提出更严峻的挑战^[54]。鉴于2000、2010和2015年未超载的县还有39、34和41个, 可以调整未超载县的牲畜存栏量, 优化畜牧产业布局。但应该优先控制处于红色预警区的牲畜存栏量, 避免其超载状况进一步恶化。

目前的研究是以县域为单位的草地承载力空间分布和预警, 但牲畜并非是均匀分布在所有草地中, 要更加精确的研究栅格尺度的草地承载力状况, 还需要结合人口密度、坡度和净初级生产力等因素将实际牲畜存栏量进行栅格化^[55]。同时本研究将草地的所有植物物种纳入了研究结果, 但并非所有物种均可被牲畜采食, 因此有效区分不可采食区域, 对于进一步提高理论载畜量计算精度, 以便于提出更加精细化的管理措施, 是未来研究的重要方向。

3.2 结论

本文核算了青藏高原地区的2000—2015年的草地生产力, 结合相应年份的县域牲畜存栏量, 获得了青藏高原全域的草畜平衡指数, 并根据其变化趋势对青藏高原草地承载力进行预警, 结论如下：

(1) 青藏高原草地承载力整体呈东高西低的格局, 其中高寒草原和高寒草甸是该地区草地承载力的主要组成部分；

(2) 青藏高原地区草地承载力呈上升趋势, 理论载畜量由8613.35万羊单位增至9437.70万羊单位；

(3) 青藏高原整体处于超载状态, 2000—2010年超载状况加剧, 至2015年超载状况稍有缓解, 草畜平衡指数由67.88%增至79.90%, 再降至67.91%。其中红色预警区的县有72个, 亟需优先调控其牲畜存栏量。