全 文 :书青海省草原鼠害区域草地生物量及盖度动态监测研究
陈梦蝶1,黄晓东1,侯秀敏2,冯琦胜1,于惠1,郭正刚1,梁天刚1
(1.兰州大学草地农业科技学院 草地农业生态系统国家重点实验室,甘肃 兰州730020;2.青海省草原总站,青海 西宁810008)
摘要:本研究利用Terra/MODIS植被指数产品 MOD13A3以及2005-2006年青海省草地调查资料,建立了草地
地上生物量和植被盖度遥感反演模型,模拟分析了青海省鼠害发生区域、危害区域和正常区域草地资源在2001-
2010年期间的各月及年际生物量及植被盖度变化动态。结果表明,1)EVI指数函数模型为模拟生物量鲜重最优模
型,而对数函数模型为植被盖度模拟最优模型,拟合精度分别为70.18%和77.43%;2)研究区不同草地类型鼠害
发生区域、危害区域的地上生物量和植被盖度随着鼠害危害状况的加重而呈下降趋势,在盛草季表现最明显,这种
趋势近10年基本保持不变。
关键词:鼠害;MODIS;植被指数;青海省
中图分类号:S812.6 文献标识码:A 文章编号:10045759(2013)04024710
犇犗犐:10.11686/cyxb20130430
草地是陆地生态系统的重要组成部分,在生态环境中起着举足轻重的作用[1,2]。它不仅是发展国民经济的
物质基础,还是维护陆地生态环境的天然屏障[3]。近年来,生物灾害连年持续发生,造成草地严重沙化、退化,生
产力逐年下降,草地生态系统恶性循环,直接威胁广大农牧民群众的生存环境,严重制约着畜牧业经济的稳步发
展[4]。因此,对草原生产力的监测,尤其是对鼠害区域生产力的监测对于研究陆地生态系统的碳循环、草地畜牧
业生产以及荒漠化研究与治理具有重要的意义[5]。
草地生物量和植被盖度是衡量天然草地生产力高低的重要参数,也是指示生态系统变化的重要指标[6,7]。
20世纪80年代初期,国内外开始利用“3S”技术进行草地植被研究,并且取得了一些成就[811]。随着美国Terra
和Aqua卫星的发射,中分辨率成像光谱仪 MODIS(moderateresolutionimagingspectroradiometer)数据在草
地资源监测方面得到了广泛的应用。梁天刚等[12]利用2006-2008年甘南牧区草地调查资料和 MODIS/Terra
每日地表反射率产品 MOD09GA,建立了草地地上生物量遥感反演模型,分析了甘南州及各县草地资源在
2001-2008年期间的各旬、月和年的生物量及理论载畜量变化动态。王正兴等[13]利用 MODIS增强型植被指数
反演草地地上生物量,对比了4种草地类型和2种植被指数,建立了VIANPP(abovenetprimaryproductivity)
的线性模型和幂模型。除多等[14]根据西藏高原中部地面观测的植被覆盖度和同期接收的EOS/MODIS数据,
分别建立了250m分辨率归一化植被指数(NDVI)、土壤调节植被指数(SAVI)与地面观测的植被覆盖度之间的
相关关系,并验证了Carlson和Ripley植被覆盖度算法的精度。王国成等[15]利用遥感和气象数据,采用CASA
(carnegieamesstanfordapproach)模型和统计分析方法,研究了1981-2001年内蒙古自治区草地7-8月植被
净初级生产力(netprimaryproductivity,NPP)时空变化特征。
纵观国内外相关研究,基于“3S”技术的草地资源生产力监测的研究已经取得了不少研究成果,然而对鼠害
发生区域的生产力监测以及对鼠害与草地生产力之间的相互关系等方面的研究却鲜有报道。本研究利用 MO
DIS数据,分别建立青海省草地生物量和植被盖度估测模型,深入探讨近10年(2001-2010年)青海省鼠害区域
的草地生产力的变化动态,并且与该省非鼠害发生区域的植被生产力状况进行对比分析,最终明确鼠害的发生对
鼠害区域植被生产力各项指标的影响,以期为利用“3S”技术监测草原鼠害提供理论依据。
第22卷 第4期
Vol.22,No.4
草 业 学 报
ACTAPRATACULTURAESINICA
247-256
2013年8月
收稿日期:20110918;改回日期:20120111
基金项目:国家自然科学基金项目(31172258)和教育部科技创新工程重大项目培育资金项目(708089)资助。
作者简介:陈梦蝶(1989),女,浙江台州人,在读硕士。Email:chenmd11@lzu.cn
通讯作者。Email:huangxd@lzu.edu.cn
1 材料与方法
1.1 研究区概况
青海省位于青藏高原东北部,地理位置介于 N31°39′~39°19′,E89°35′~103°04′。全境天然草地面积
3.64494×107hm2,其中可利用草地面积3.16103×107hm2,是我国发展草地畜牧业的主要基地之一[16]。但
是,青海也是鼠害发生较为严重的省份之一。全省有草地害鼠约30余种,常见的有高原鼠兔(犗犮犺狅狋狅狀犪犮狌狉狕狅狀犻
犪犲)、高原鼢鼠(犕狔狅狊狆犪犾犪狓犫犪犻犾犲狔犻)、高原田鼠(犘犻狋狔犿狔狊犻狉犲狀犲)、长尾仓鼠(犆狉犻犮犲狋狌犾狌狊犾狅狀犵犻犮犪狌犱狌狋狌狊)及五趾跳鼠
(犃犾犾犪犮狋犪犵犪狊犻犫犻狉犻犮犪)等。其中以高原鼠兔分布最广,数量最多,危害也最为严重。1998年调查资料显示,全省发
生鼠害面积971.67×104hm2,其中严重危害面积736.19×104hm2;由于害鼠的挖掘和啃食活动,致使草地植被
破坏,造成大量水土流失,形成大片次生裸地和寸草不生的“黑土滩”。据调查,全省由鼠害造成的“黑土滩”型退
化草地面积达333.33×104hm2;同时,害鼠大量啃食优良牧草,与家畜争食,加剧了草畜矛盾,降低了草地的载
畜能力,给畜牧业生产造成了极大的危害[17,18]。因此,利用“3S”技术对青海省鼠害区域的监测,不但能够及时准
确地了解草地生物量的时空分布状况,而且能够更科学地量化草地生产和系统内动物活动的关系,从而才能科学
合理地利用、管理草地,保证畜牧业生产持续稳定发展和生态环境的有效改善。
1.2 卫星遥感资料
本研究所使用的遥感数据为NASA(NationalAeronauticsandSpaceAdministration)提供的EOS/MODIS
MOD13A3数据。该数据空间分辨率为1km,为月最大值合成图像,数据格式为EOSHDF,空间位置在全球正
弦曲线投影SIN(sinusoidalprojection)系统中编号为h26v05和h25v05的两景图像,时间序列为2001-2010
年,共计240幅EVI(enhancedvegetationindex)数字图像。近年来其应用得到极大的推广,除了对草地生产力
的监测,还应用于自然灾害监测[19]、土地覆盖与土地利用[20]、生态环境动态监测[21]等方面。
1.3 地面采样数据
野外调查时间分别为2005年和2006年的7-9月,采样点分布于青海省各县(图1)。总计调查样地563个、
包括地面样方1167个。其中2005年设置样地172个,包括583个样方;2006年设置样地391个,包括584个样
方。样地设置尽量选择草地植被空间分布比较均一的区域。每个样地内设置1~3个1m×1m的样方,在植被
群落一致性较好的样地内设置1~2个,当样地植被分布不均匀时,设置2~3个,个别极复杂样地设3~5个。记
录每个样点的隶属行政区、草地类型、坡度、坡向,利用状况评价等,同时利用GPS测定样方的经纬度和海拔,并
在样方内采集其盖度、高度、植物种数、不可食牧草种数、总产草量和可食牧草产量等数据。
1.4 鼠害区域划分方法
根据青海省1998年草地害鼠普查结果,综合青海草地害鼠种群数量、挖掘能力以及采食特点,本研究将青海
省草地划分为如下3个区域:1)正常区域:草地害鼠种群密度未达到防治指标。2)发生区域:草地害鼠种群密度
超过防治指标。3)危害区域:草地害鼠种群密度达到防治指标2倍(含2倍)以上。
1.5 植被指数提取及处理
目前所使用的植被指数大约有20多种[22],如比值植被指数RVI(ratiovegetationindex)、归一化差值植被指
数NDVI(normalizeddifferencevegetationindex)[23]、差值植被指数DVI(differencevegetationindex)等,本研究
选取增强型植被指数EVI来进行研究。增强型植被指数EVI(enhancedvegetationindex)是对AVHRRNDVI
的继承和改进,相对于AVHRRNDVI,MODISEVI对基础数据进行了全面的大气校正,采用“大气抵抗植被指
数”对残留的气溶胶做了进一步校正。因此具有以下一些特点:第一,具有较高的空间分辨率,可详细地反映地表
植被特征。第二,红光和近红外探测波段的范围设置更窄,不仅提高了对稀疏植被探测的能力,而且减少了水汽
的影响,同时,引入了蓝光波段对大气气溶胶进行了修正[24]。
犈犞犐= 2.5
(ρNIR-ρred)
(ρNIR+犆1-犆2+犔)
(1)
式中,ρNIR和ρred为经过大气校正的近红外和红光通道的反射值;犔=1,为土壤调节参数;参数犆1和犆2分别为6.0
和7.5。
842 ACTAPRATACULTURAESINICA(2013) Vol.22,No.4
图1 研究区鼠害区域及地面样方分布图
犉犻犵.1 犜犺犲犱犻狊狋狉犻犫狌狋犻狅狀狅犳狉狅犱犲狀狋狉犲犵犻狅狀犪狀犱犻狀狊犻狋狌狅犫狊犲狉狏犪狋犻狅狀狊犻狀狋犺犲狊狋狌犱狔犪狉犲犪
利用 MODISTools工具软件将下载的 MODIS原始影像进行格式转换,由 HDF转换为TIFF格式,同时完
成空间拼接和投影转换;在ArcMap软件中使用ConversionTools把影像转换成Grid格式;用双线性插值方法
提取各样方点对应的 MODISEVI值,并作归一化处理。
2 结果与分析
2.1 草地生物量、盖度模型构建及精度分析
利用一元线性回归和曲线回归方法,分析了研究区2006年草地生物量和植被盖度同EVI值之间的相关关
系并进行了统计分析(表1,表2)。结果表明,1)犉检验和犜 检验的统计分析结果表明,MODISEVI值和生物
量、植被盖度之间呈显著相关。2)MODISEVI与生物量的4种模型中,指数函数模型最优,决定系数犚2 为
0.5024,明显优于其他模型;3)MODISEVI与植被盖度的4种模型精度都较高,其中对数函数模型最优,决定系
数犚2 为0.5669。因此,本研究将利用指数函数模型和对数函数模型来分别反演研究区草地地上生物量和植被
盖度。
利用2005年的实测数据,采用交叉验证法,对模拟数据进行了验证(图2,图3)。结果表明,模拟生物量和植
被盖度与实测值吻合度较高。同时,通过计算该组数据的平均误差和绝对标准差进行精度分析。
相对误差=
实测值-估测值
实测值 ×100% (2)
平均误差=∑|相对误差|/犖 (3)
式中,实测值为样地实际生物量和植被盖度测量值,估测值为依据模型计算生成的模拟值,犖 为样本数。结果表
明,生物量模型的绝对标准差为1375.75kg/hm2,平均误差系数为29.82%,模型精度为70.18%;盖度模型的绝
对标准差为11.25%,平均误差系数为22.57%,模型精度为77.43%。2个模型的精度都比较高,能较好地模拟
研究区草地地上生物量和植被盖度。
942第22卷第4期 草业学报2013年
表1 青海省草地植被指数与生物量和植被盖度回归模型
犜犪犫犾犲1 犜犺犲狉犲犵狉犲狊狊犻狅狀犿狅犱犲犾犳狅狉狋犺犲狏犲犵犲狋犪狋犻狅狀犻狀犱犲狓犪狀犱犫犻狅犿犪狊狊/狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉犪犵犲
统计分析指标
Statisticalanalysisindicators
回归模型
Regressionmodel
模型方程
Modelequation
决定系数犚2
Decisioncoefficient
犉
草地生物量 一元线性回归模型Alinearregressionmodel 狔=8996.2狓-340.51 0.3693 340.764
Grasslandbiomass 指数模型Exponentialmodel 狔=657.13e3.3817狓 0.5024 587.524
对数模型Logarithmsmodel 狔=2278ln(狓)+5567 0.2910 239.846
乘幂模型Powerfunctionmodel 狔=6475狓0.918 0.4560 489.538
植被盖度 一元线性回归模型Alinearregressionmodel 狔=108.55狓+32.401 0.5109 599.701
Vegetationcover 指数模型Exponentialmodel 狔=29.91e2.128狓 0.4350 443.135
对数模型Logarithmsmodel 狔=32.534ln(狓)+109.15 0.5669 751.222
乘幂模型Powerfunctionmodel 狔=142.4狓0.689 0.5640 742.971
注:表示相关极显著(犘<0.001)。下同。
Note:indicatessignificantcorrelation(犘<0.001).Thesamebelow.
表2 基于犈犞犐的生物量和植被盖度最优模型的估计参数
犜犪犫犾犲2 犈狊狋犻犿犪狋犲犱狆犪狉犪犿犲狋犲狉狊狅犳狋犺犲狅狆狋犻狅狀犪犾犿狅犱犲犾犳狅狉犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犪狀犱狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉犪犵犲犫犪狊犲犱狅狀犈犞犐
指标Index 模型 Model
参数估计与检验Estimationandtestingoftheparameters
参数Parameters 估计值Estimatedvalue 犜
生物量Abovegroundbiomass 指数Exponent b 3.3817 24.239
a 657.13 17.299
植被盖度Vegetationcover 对数Logarithm b 109.15 21.051
a 32.534 -11.624
图2 草地地上生物量实测值和估测值对比
犉犻犵.2 犆狅犿狆犪狉犻狊狅狀犫犲狋狑犲犲狀犳犻犲犾犱狊狌狉狏犲狔犻狀犵犪狀犱犲狊狋犻犿犪狋犲
狏犪犾狌犲狅犳犵狉犪狊狊犾犪狀犱犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊
图3 植被盖度实测值和估测值对比
犉犻犵.3 犆狅犿狆犪狉犻狊狅狀犫犲狋狑犲犲狀犳犻犲犾犱狊狌狉狏犲狔犻狀犵犪狀犱
犲狊狋犻犿犪狋犲狏犪犾狌犲狅犳狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉
2.2 青海省草地生产力变化动态分析
2.2.1 青海省各区域草地生产力的月季动态变化分析 青海省草地包含的草地类型共有13种,其中受鼠类危
害的草地类型主要包括5种,分别为高寒草甸类、低地草甸类、山地草甸类、温性草原类和高寒草原类。不同草地
类型草地4月份开始返青,7、8月份因日照时数最长,雨水量充沛,进入盛草期。2010年7月份高寒草甸月平均
生物量达5103kg/hm2,植被盖度达92%。9月份之后牧草进入枯黄期,生物量逐渐下降,到11月中旬达到最小
值。从11月中旬到次年的4月份是草地的枯萎期,这些时段研究区草地生物量均维持在1000kg/hm2,植被盖
052 ACTAPRATACULTURAESINICA(2013) Vol.22,No.4
度也都在50%之下。除了低地草甸,其他4种草地类型在牧草枯萎期间,鼠害发生区域、危害区域的地上生物量
和正常区域的地上生物量差异不大,正常区域和发生区域的植被盖度略高于危害区域。5月中旬到10月中旬,
草地生长旺盛,鼠类也开始猖獗,原生植被受到破坏,鼠害区域草地地上生物量和植被盖度的增加量相对正常区
域明显减少,并且随着鼠害程度的加剧,植被受影响程度随之加大。以受鼠害影响面积最大的高寒草甸类草原为
例,2010年7月,青海省正常区域的高寒草甸草地月平均生物量为5103kg/hm2,鼠害发生区域高寒草甸草地月
平均生物量为4862.5kg/hm2,而鼠害危害区域高寒草甸草地月平均生物量仅为3363.9kg/hm2,明显低于非
鼠害区域草地生物量,并且危害区域的草地生物量相对于正常区域和发生区域的生物量,差异较大(表3)。说明
草原啮齿类动物在盛草期最为肆虐,此时对草原生产力的影响也是最大的。此时,正常区域的植被盖度为92%,
发生区域为91%,危害区域为84%(表4),说明随着鼠害危害程度上升,啮齿动物的采食行为以及鼠兔打洞和鼢
鼠堆土丘的行为加重了对植被的破坏和次生裸地产生,植被覆盖度下降。但是相对于生物量来说,下降程度减
少,这可能是盖度下降到一定程度,一些阔叶毒杂草侵入,植被盖度不再下降甚至有所回升[25]。
表3 青海省各个草地类型地上生物量的月季动态变化
犜犪犫犾犲3 犕狅狀狋犺犾狔犱狔狀犪犿犻犮狊狅犳犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊狅犳狏犪狉犻狅狌狊犵狉犪狊狊犾犪狀犱狋狔狆犲狊犻狀犙犻狀犵犺犪犻 kg/hm2
草地类型
Grassland
type
区域
Zone
生物量Biomass
1月
Jan.
2月
Feb.
3月
Mar.
4月
Apr.
5月
May
6月
Jun.
7月
Jul.
8月
Aug.
9月
Sep.
10月
Oct.
11月
Nov.
12月
Dec.
平均
Average
温性草原类
Temperate
steppe
正常Rodentfree 984 989 1013 1045 1402 2863 4091 3606 2397 1332 1009 976 1809
发生Rodentoccurred 834 826 812 805 811 838 882 895 885 860 848 863 847
危害Rodentdamaged 840 838 864 921 1156 1415 1506 1289 1215 1068 903 845 1072
高寒草原类
Alpine
steppe
正常Rodentfree 957 961 974 1000 1236 2126 3376 3295 2218 1260 982 966 1613
发生Rodentoccurred 878 884 889 898 993 1682 1916 1814 1405 1037 904 905 1184
危害Rodentdamaged 856 860 856 868 951 1290 1601 1564 1245 976 888 884 1070
低地草甸类
Lowland
meadow
正常Rodentfree 905 900 918 955 1123 1760 2322 2345 1910 1182 998 978 1358
发生Rodentoccurred 915 882 901 927 1126 1843 2359 2232 1816 1162 952 939 1338
危害Rodentdamaged 850 868 865 873 902 1297 1528 1465 1167 933 883 886 1043
山地草甸类
Mountain
meadow
正常Rodentfree 951 975 1014 1061 1391 2554 4035 3778 2595 1436 1000 972 1813
发生Rodentoccurred 949 965 982 1010 1461 2147 3148 2808 1996 1369 1004 958 1566
危害Rodentdamaged 905 918 929 950 1250 1832 2682 2392 1664 1207 943 915 1382
高寒草甸类
Alpine
meadow
正常Rodentfree 1015 1030 1060 1091 1398 2867 5103 4765 2963 1514 1057 1025 2074
发生rodentoccurred 1004 1027 1067 1101 1383 2969 4862 4606 2935 1570 1080 1046 2054
危害rodentdamaged 955 980 1006 1029 1293 2534 3364 3332 2266 1342 1005 984 1674
2.2.2 青海省各区域草地生产力状况及年际动态变化比较分析 青海省不同草地类型一般在每年的7月或8
月份生物量达到最高,为草地的生长旺季,对草地生产力具有典型的代表性,因此本研究选取2001-2010年8月
份的平均生物量和植被盖度来反映整个青海省草地生产力概况以及2001-2010年草地生物量和植被盖度变化
动态。2001-2010年期间,不同草地类型8月份平均生物量的差异极其明显,高寒草甸类和山地草甸类的生物
量较高,正常区域月平均地上生物量>4500kg/hm2;高寒草原类和低地草甸类生物量相对较低,正常区域月平
均产量<3000kg/hm2(表3);植被盖度的差异情况大致和平均生物量相同,只是各个草地类型植被盖度的差距
较小(表4)。
从全境范围的草地生产力来看,2001-2008年各个草地类型草地生物量和植被盖度曲线都十分平缓,但是
2008-2010年高寒草甸类、山地草甸类、高寒草原类和温性草原类的平均生物量和植被盖度的曲线均出现明显
152第22卷第4期 草业学报2013年
上升,低地草甸类略微有所不同,在2008-2009年低地草甸生物量和植被盖度降低,2009-2010年回升。说明
了研究区在2008-2010年间的草地生产力整体有所回升。这可能跟近几年研究区气候状态的改变相关,也可能
是近年来人们对生态环境的保护举措起到了一定的效果,具体原因需要进一步探究。
而对于鼠害发生区域、鼠害危害区域和正常区域的各种草地类型的生物量和植被盖度分别进行分析,结果表
明,1)2001-2010年,各个区域草地生物量和植被盖度动态变化趋势和全境的动态变化趋势一致。2)由于高原
鼠兔或高原鼢鼠等啮齿类动物的影响,2001-2010年各个草地类型的地上生物量和植被盖度都有不同程度的降
低,并且都遵循正常区域>发生区域>危害区域这个规律,说明青海省的鼠害发生严重危害了草地生产,这种情
况延续近10年一直未见改善,需要强烈重视和采取措施来解决这个问题,从而从根本上防治草地次生演替和退
化。3)其中温性草原类和高寒草原类草地,鼠害发生区域和危害区域的生物量和植被盖度相对于正常区域,都有
很大程度的降低,说明对于青海省来说,鼠害对温性草原类和高寒草原类影响较大,是今后鼠害防治和治理工作
的重点区域(图4~8)。
表4 青海省各个草地类型植被盖度的月季动态变化
犜犪犫犾犲4 犕狅狀狋犺犾狔犱狔狀犪犿犻犮狊狅犳狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉狅犳狏犪狉犻狅狌狊犵狉犪狊狊犾犪狀犱狋狔狆犲狊犻狀犙犻狀犵犺犪犻 %
草地类型
Grassland
type
区域
Zone
植被盖度 Vegetationcover
1月
Jan.
2月
Feb.
3月
Mar.
4月
Apr.
5月
May
6月
Jun.
7月
Jul.
8月
Aug.
9月
Sep.
10月
Oct.
11月
Nov.
12月
Dec.
平均
Average
温性草原类
Temperate
steppe
正常Rodentfree 39 40 42 44 57 79 88 85 76 57 41 39 57
发生Rodentoccurred 22 21 19 17 19 23 28 30 29 26 25 27 24
危害Rodentdamaged 22 22 27 33 49 58 58 52 51 44 31 22 39
高寒草原类
Alpine
steppe
正常Rodentfree 34 36 37 39 49 68 78 78 69 51 38 37 51
发生Rodentoccurred 28 29 29 30 35 58 63 62 54 39 30 29 40
危害Rodentdamaged 24 25 25 26 32 49 56 56 48 35 28 26 36
低地草甸类
Lowland
meadow
正常Rodentfree 30 28 32 35 43 58 65 66 61 46 37 36 45
发生Rodentoccurred 33 30 29 33 43 60 69 69 64 46 35 34 45
危害Rodentdamaged 25 27 27 28 31 53 60 59 49 32 29 30 37
山地草甸类
Mountain
meadow
正常Rodentfree 35 36 41 44 57 75 85 83 75 59 39 37 56
发生Rodentoccurred 36 37 39 41 58 73 83 80 71 58 41 37 55
危害Rodentdamaged 29 32 33 35 48 64 74 72 61 50 34 30 47
高寒草甸类
Alpine
meadow
正常Rodentfree 42 44 46 47 59 80 92 91 83 63 45 43 61
正常Rodentfree 40 42 46 48 58 79 91 90 82 64 46 44 61
发生Rodentoccurred 35 39 41 43 53 77 84 84 75 57 41 39 56
3 结论与讨论
目前,“3S”技术在草地生产力以及灾害研究方面的应用得到了极大的推广[17],广泛应用的 MODIS植被指
数产品也多样化[13]。前人在草原鼠害监测方面的研究做了一系列的工作,鼠害动态监测的常规方法和内容一直
在不断探讨和改进[26,27]。本研究在前人利用“3S”技术对草地生产力和鼠害监测研究的基础上,利用 MODIS植
被指数产品,结合地面调查资料,利用一元线型回归和曲线回归方法,建立了青海省草地生物量和植被盖度遥感
监测模型,并对研究区2001-2010年草地地上生物量和植被盖度进行了模拟。同时利用GIS手段将研究区域
分成鼠害发生区域、鼠害危害区域和正常区域分别进行研究和对比分析,研究结果表明:1)青海省共有草地类型
13种,其中受鼠类危害较严重的主要包括5种,分别为高寒草甸类、低地草甸类、山地草甸类、温性草原类和高寒
252 ACTAPRATACULTURAESINICA(2013) Vol.22,No.4
图4 青海省温性草原各区域生物量、植被盖度的年季动态变化
犉犻犵.4 犃狀狀狌犪犾犱狔狀犪犿犻犮狊狅犳犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犪狀犱狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉狅犳狋犲犿狆犲狉犪狋犲狊狋犲狆狆犲犻狀犙犻狀犵犺犪犻
图5 青海省高寒草原各区域生物量、植被盖度的年季动态变化
犉犻犵.5 犃狀狀狌犪犾犱狔狀犪犿犻犮狊狅犳犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犪狀犱狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉狅犳犪犾狆犻狀犲狊狋犲狆狆犲犻狀犙犻狀犵犺犪犻
图6 青海省低地草甸各区域生物量、植被盖度的年季动态变化
犉犻犵.6 犃狀狀狌犪犾犱狔狀犪犿犻犮狊狅犳犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犪狀犱狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉狅犳犾狅狑犾犪狀犱犿犲犪犱狅狑犻狀犙犻狀犵犺犪犻
草原类。受高原鼠兔和高原鼢鼠等啮齿类动物的影响,2001-2010年各个草地类型鼠害区域地上生物量和植被
盖度都有不同程度的降低,遵循正常区域>发生区域>危害区域这个规律,说明随着危害等级的上升,鼠害对草
地的破坏更加严重,与前人的研究结果一致[28]。进一步说明青海省的鼠害严重危害草地生产力,这种情况延续
近10年一直未见改善,需要强烈重视和采取措施来解决这个问题,从而从根本上防治草地次生演替和退化。2)
王建柱[28]研究的桑科草原不同鼠类群落主要栖息于高寒草甸草场、亚高山草甸草场、草原草场以及草原草甸草
场等,武晓东[29]研究的内蒙古荒漠半荒漠区草地不同鼠类群落主要集中于整个荒漠半荒漠草原,本研究结果表
352第22卷第4期 草业学报2013年
图7 青海省山地草甸各区域生物量、植被盖度的年季动态变化
犉犻犵.7 犃狀狀狌犪犾犱狔狀犪犿犻犮狊狅犳犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犪狀犱狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉狅犳犿狅狌狀狋犪犻狀犿犲犪犱狅狑犻狀犙犻狀犵犺犪犻
图8 青海省高寒草甸各区域生物量、植被盖度的年季动态变化
犉犻犵.8 犃狀狀狌犪犾犱狔狀犪犿犻犮狊狅犳犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犪狀犱狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉狅犳犪犾狆犻狀犲犿犲犪犱狅狑犻狀犙犻狀犵犺犪犻
明高寒草甸类、低地草甸类、山地草甸类、温性草原类和高寒草原类5种草地受鼠类危害较重,草地大类与王建
柱[28]在桑科草原的研究结果一致。其中温性草原类和高寒草原类草地鼠害发生区域和危害区域的生物量和植
被盖度相对于正常区域,都有很大程度的降低,说明对于青海省来说,可能因其区域性的鼠害种群特点和分布规
律,鼠害对温性草原类和高寒草原这两类草地的影响较大,是今后鼠害防治的重点区域。3)本研究利用 MODIS
植被指数产品,建立了草地生物量和盖度遥感监测模型,分析了青海省草原鼠害对草地生产力的影响,研究结果
证明了鼠害对草地生产力的破坏作用,以期为合理制定鼠害防治提供科学依据,对青海省草地畜牧业的可持续发
展具有重要的研究意义。由于本研究只采用 MODISEVI数据对整个青海省的数据进行统一建模,并未对每种
草地类型进行分类型建模,没有充分考虑不同草地类型相对独立的自然特点和规律,因此建立的模型具有一定的
不确定性。另一方面,植被生长过程受到很多因素(环境条件、人为因素等)的影响,在今后的研究中,可探索通过
建立害鼠种群密度与植被盖度、生物量以及 MODISEVI的关系,进一步分析鼠害对草地生产力的影响,以期得
到更加科学、合理的结果,从而为今后的鼠害监测和防治提供指导。
参考文献:
[1] 王新欣,朱进忠,范燕敏,等.基于 MODISNDVI的天山北坡中段草地动态估产模型研究[J].草业科学,2009,26(7):24
27.
[2] 刘兴元,龙瑞军,尚占环.草地生态系统服务功能及其价值评估方法研究[J].草业学报,2011,20(1):167174.
[3] 李霞.北疆地区草地生物量遥感监测[D].兰州:兰州大学,2008.
[4] 徐秀霞.青海省草地生物灾害现状及应对措施[J].草业与畜牧,2006,129(8):3638.
452 ACTAPRATACULTURAESINICA(2013) Vol.22,No.4
[5] 辛晓平,张保辉,李刚,等.1982—2003年中国草地生物量时空格局变化研究[J].自然资源学报,2009,24(9):15821592.
[6] 俞联平,李发弟,李新媛,等.基于3S技术的甘肃省甘州区荒漠草地生产力评价[J].草原与草坪,2008,130(5):3639.
[7] 张云霞,李晓兵,陈云浩,等.草地植被盖度的多尺度遥感与实地测量方法综述[J].地球科学进展,2003,18(1):8593.
[8] TaylorBF.DeterminationofseasonalandinterannualvariationinNewZealandpasturegrowthfromNOAA7Data[J].Re
moteSensingofEnvironment,1985,18:177192.
[9] ParueloJM,EssteiHE,LauenrothWK,犲狋犪犾.ANPPestimatesfromNDVIforthecentralgrasslandregionoftheUnited
States[J].Ecology,1997,78(3):953958.
[10] 刘海启,金敏毓,龚维鹏.美国农业遥感技术应用状况概述[J].中国农业资源与区划,1999,20(2):5660.
[11] 赵冰茹,刘闯,刘爱军,等.利用 MODISNDVI进行草地估产研究———以内蒙古锡林郭勒草地为例[J].草业科学,2004,
21(8):1215.
[12] 梁天刚,崔霞,冯琪胜,等.2001-2008年甘南牧区草地地上生物量与载畜量遥感动态监测[J].草业学报,2009,18(6):
1222.
[13] 王正兴,刘闯,AlfredoH.植被指数研究进展:从AVHRR-NDVI到MODIS-EVI[J].生态学报,2003,23(5):979986.
[14] 除多,次仁多吉,王彩云,等.利用 MODIS数据估算西藏高原地表植被覆盖度[J].遥感技术与应用,2010,25(5):707
713.
[15] 王国成,张稳,黄耀.1989-2001年内蒙古草地净初级生产力时空变化特征[J].草业科学,2011,28(11):20162025.
[16] 鲍文.青藏高原草地资源发展面临的问题及战略选择[J].农业现代化研究,2009,30(1):2023.
[17] 张生合,任程,陈国民,等.青海省草地鼠害防治及今后设想[J].青海草业,2001,10(2):2224.
[18] 韩立辉,尚占环,任国华,等.青藏高原“黑土滩”退化草地植物和土壤对秃斑面积变化的响应[J].草业学报,2011,20(1):
16.
[19] 高懋芳,覃志豪,刘三超.MODIS数据在林火监测中的应用研究[J].国土资源遥感,2005,(2):6064.
[20] 柳海鹰,高吉喜,李政海.土地覆盖及土地利用遥感研究进展[J].国土资源遥感,2001,(4):712.
[21] 李聪,肖继东,宫恒瑞,等.MODIS数据在乌鲁木齐地区植被景观动态监测中的应用[J].新疆气象,2005,28(3):2729.
[22] 田庆久,阅祥军.植被指数研究进展[J].地球科学进展,1998,13(4):327333.
[23] 刘睿,孙九林,张金区,等.中国北方草地覆被的 HJ星NDVI校正研究[J].草业学报,2011,20(1):189198.
[24] 张春桂.用 MODIS植被指数研究福州城区空间扩展变化[J].气象,2006,32(10):2026.
[25] 聂学敏,石红霄,李志强.鼠类种群密度变化及其对高寒草地植物群落的影响[J].湖南农业科学,2010,(2):7981.
[26] 杨东生,赵新春,熊琳,等.新疆草原鼠害的动态监测[J].新疆畜牧业,2004,(6):5860.
[27] 王玮,冯琦胜,于惠,等.“3S”技术在草地鼠虫害监测与预测中的应用[J].草业科学,2010,27(3):3139.
[28] 王建柱.基于3S技术系统鼠害监测的基础研究[D].兰州:甘肃农业大学,2003.
[29] 武晓东.内蒙古半荒漠与荒漠区啮齿动物群落研究、区域性危害区划及GIS分析[D].呼和浩特:内蒙古农业大学,2004.
552第22卷第4期 草业学报2013年
犇狔狀犪犿犻犮犿狅狀犻狋狅狉犻狀犵狅犳犫犻狅犿犪狊狊犪狀犱狏犲犵犲狋犪狋犻狅狀犮狅狏犲狉犪犵犲犻狀狉狅犱犲狀狋犱犪犿犪犵犲犱犵狉犪狊狊犾犪狀犱
狉犲犵犻狅狀狊狅犳犙犻狀犵犺犪犻犘狉狅狏犻狀犮犲,犆犺犻狀犪
CHENMengdie1,HUANGXiaodong1,HOUXiumin2,FENGQisheng1,
YUHui1,GUOZhenggang1,LIANGTiangang1
(1.ColegeofPastoralAgricultureScienceandTechnology,LanzhouUniversity,StateKey
LaboratoryofGrasslandAgroecosystem,Lanzhou730020,China;2.Qinghai
ProvinceGrasslandsStation,Xining810008,China)
犃犫狊狋狉犪犮狋:CombiningTerra/MODISvegetationindexproductsMOD13A3andfieldsurveyingdatafrom2005-
2006,remotesensinginversionmodelsofgrasslandabovegroundbiomass(GAB)andvegetationcoverage
(VC)inQinghaiProvincewereestablished.Theywereusedtoseparatelysimulateandanalyzethemonthly
andannualdynamicchangesofbiomassandvegetationcoverageduring2001-2010inrodentinfestedregions,
harmfulregionsandrodentfreeregionsinthestudyarea.Rodentshadasignificantinfluenceongrasslandpro
ductivityinthestudyarea.1)AnexponentialmodelofMODISEVI(enhancedvegetationindex)andaloga
rithmicmodelweretheoptimalmodelstoestimatethegreenyieldofGABandVC,respectivelyandhadapre
cisionof70.18%and77.43%,respectively.2)TheGABandVCfordifferenttypesofgrasslandsinthestudy
areadecreasedwiththeaggravationofrodentpests,especialyinJulyandAugustofeachyear,andthistend
encyremainedstableduringthewholestudyperiod.
犓犲狔狑狅狉犱狊:rodentpest;MODIS;vegetationindex;QinghaiProvince
652 ACTAPRATACULTURAESINICA(2013) Vol.22,No.4