免费文献传递   相关文献

Remotely sensed dynamics monitoring of grassland aboveground biomass and carrying capacityduring 2001-2008 in Gannan pastoral area

2001-2008年甘南牧区草地地上生物量与载畜量遥感动态监测



全 文 :书2001-2008年甘南牧区草地地上
生物量与载畜量遥感动态监测
梁天刚1,崔霞2,冯琦胜1,王莺1,夏文韬1
(1.兰州大学草地农业科技学院 农业部草地农业生态系统学重点实验室,甘肃 兰州730020;2.兰州大学资源与环境学院,甘肃 兰州730000)
摘要:草地地上生物量监测是草地资源空间格局动态研究的重要内容,也是草畜平衡综合分析的基础。利用2006
-2008年甘南牧区草地调查资料和Terra/MODIS每日地表反射率产品 MOD09GA,建立了草地地上生物量遥感
反演模型,模拟分析了甘南州及各县市草地资源在2001-2008年期间的各旬、月和年的生物量及理论载畜量变化
动态。研究结果表明,MODIS增强型植被指数EVI的乘幂函数可以较好地模拟甘南牧区草地地上的生物量鲜
重,拟合模型平均估产精度为76.7%,可很好地模拟牧草生长状况较好时期(5-10月)的草地地上生物量变化动
态。甘南牧区草地生长主要集中在5月上旬-10月下旬期间,草地旬最大地上生物量数字图像可以较客观地反映
草地植被生长发育的总体规律。但是,个别旬生物量受大范围长时间阴雨多云天气状况及放牧家畜数量有较大变
动情况的严重影响。2001-2008年不同草地类型的月最大生物量动态变化曲线均呈单峰抛物线形式,全州平均最
大生物量均出现在7月,但不同年份产量达到最大值的月份有所变化,主要集中在7-8月。甘南州草地地上总生
物量年度之间存在较大的差异。全州8年平均总地上生物量为109.31亿kg。2005年全州草地植被年总生物量最
高,达129.1亿kg,其次为2006,2007和2002年,分别为113.2,110.7和109.0亿kg。由于气候条件和不同县市
草地面积及生长状况等存在较大差异,其理论载畜量也存在显著差别。
关键词:甘南牧区;草地地上生物量;载畜量;遥感监测
中图分类号:S127;S812  文献标识码:A  文章编号:10045759(2009)06001211
 1999年12月和2002年2月发射成功的美国新一代对地观测系统(EarthObservingSystem,EOS)卫星
Terra和Aqua均带有中分辨率成像光谱仪(MODerateresolutionImagingSpectroradiometer,MODIS),具有观
测从太阳短波到热红外,光谱范围从0.415~14.235μm的36个光谱波段,地面分辩率为250,500和1000m,
可以反映陆地、云特征、生物地理、地表温度等特征信息。近年来,EOS/MODIS资料已广泛应用于植被生态、土
地利用和土地覆盖、自然灾害等方面的监测和研究[1]。
利用 MODIS资料开展天然草地资源的遥感动态监测,已成为当前国际草地科学研究中的前沿课题,对快速
确定合理载畜量和加强草场科学管理等均具有重要的意义。植被指数是表征地表植被特征的重要指标,MODIS
植被指数可以更好地反映植被的时空变化特征,已成为当前植被及其变化动态宏观研究的主要遥感资料和热点
问题。MODIS光谱分辨率较高,波幅窄,避免了几个大气吸收带,在计算植被指数时有更严格的去云算法和比较
彻底的大气校正[1,2]。MODIS-NDVI(NormalizedDifferenceVegetationIndex)是已有20年积累的NOAA-
NDVI系列的延续,增强型植被指数EVI(EnhancedVegetationIndex)利用 MODIS辐射仪的优点,订正地表反
射率以提高对高生物量区的敏感性,并通过叶冠背景信号耦合和减少大气影响来提高植被监测精度。这2个植
被指数可以在研究全球植被、提高植被变化的检测和提取叶冠生物物理参数方面相互补充[3]。使用 MODIS资
料已进行的相关研究有草地植被指数的空间变化规律分析[4~7]、草地地上部分生物量及初级生产力反演[8~17]、草
地退化状况、草地植被盖度、草地利用强度、土地覆盖变化监测等[18~23]。但是,在现有的草地地上生物量建模方
面,还存在选用哪种植被指数及地面实测草地样方同 MODIS植被指数产品在时间尺度上难以匹配,导致监测模
12-22
2009年12月
   草 业 学 报   
   ACTAPRATACULTURAESINICA   
第18卷 第6期
Vol.18,No.6
 收稿日期:20081222;改回日期:20090105
基金项目:“863”计划:国家高技术研究发展专项经费(2007AA10Z232)资助。
作者简介:梁天刚(1967),男,甘肃崇信人,教授。
通讯作者。Email:tgliang@lzu.edu.cn
型精度及稳定性较差等问题。
针对以上问题,本研究的目标包括2个方面:1)利用高时间分辨率的Terra/MODIS的每日地表反射率产品
MOD09GA,通过建立甘南州2001-2008年的日、旬、月和年NDVI和EVI数据库,改进草地地上生物量遥感动
态监测模型及其精度;2)模拟分析2001-2008年草地植被的地上生物量和载畜量时空变化特征,为甘南牧区草
地资源监测和草畜动态平衡研究及决策提供科学依据。
1 材料与方法
1.1 研究区概况
甘南牧区位于甘肃省西南部,青藏高原东北缘,地处黄河、长江源头,是甘肃省重要的畜牧业生产基地,也是
全国主要的少数民族集聚的草原牧区之一。地理坐标为33°06′~35°44′N,100°46′~104°44′E,平均海拔3000
m。气候属典型的高原大陆性气候,高寒阴湿,多年平均降水量400~800mm,年均温仅1~3℃,≥10℃积温持
续期仅有2个多月,全年平均日照时数2200~2400h。据1982年草地资源详查的结果[12],甘南州主要草地类
型有高寒灌丛草甸、高寒草甸、温性草甸草原、温性草原、低平地草甸、暖性草丛和沼泽。全州天然草地面积
260.25万hm2,占全州土地总面积的67.64%,其中可利用草地249.50万hm2,占草地总面积的95.86%[24,25]。
自20世纪80年代以来,甘南州生态环境不断恶化,草地退化严重,草地产草量下降了25%~45%,90%的
天然草地存在不同程度的退化,中度以上的退化草地占全州天然草地面积的80%,沼泽面积由1982年的8万
hm2 减少到2万hm2,黄河补充水量减少15%左右。由于长期依赖于牧业生产的牧民占全州总人口的82.65%,
牧区草地超载、低投入和缺乏科学管理等原因,甘南草地资源的退化问题愈演愈烈,面积不断扩大,已成为制约草
地畜牧业及社会经济可持续发展的重要因素[26~30]。
1.2 草地外业调查
外业调查主要集中在2006-2008年草地生长季期间。根据草地类型的空间分布特点和面积大小,对甘南地
区不同类型的草地设计了不同数目的样地。在2006-2007年7-8月期间,主要对分布面积较大的高寒灌丛草
甸和高寒草甸进行了调查。2006年设置样地46个,包括66个样方;2007年设置样地47个,包括133个样方。
2008年进一步扩大了草地调查的范围和类型,设置样地43个,调查样方172个。
样地设置主要选择在草地植被空间分布比较均一,可以代表较大范围的典型地段。在植被群落一致性较好
的样地内设置1或2个0.5m×0.5m的样方,在植被比较复杂、分布不均匀的区域设置3或4个1.0m×1.0m
的样方。样地基本特征调查主要记录样地所隶属行政区、草地类型、地形、坡度、坡向、土壤质地、地表特征、季节
利用方式、利用状况等。利用GPS测定样方的经纬度和海拔,同时在样方内采用常规植被调查法测定植被盖度、
草层平均高度、植物种数、群落组成及优势物种名称、不可食牧草种数、总产草量和可食产草量等指标。
1.3 每日Terra/MODIS地表反射率数据
MODIS数据来自NASAMODIS陆地产品组开发的Terra卫星的每日地面反射率产品 MOD09GA。该产
品包含 MODIS地面分辨率为500m和1km的1~7通道每日地表反射率观测数据和地学定位统计数据,已对
大气、气溶胶及薄卷云的影响进行了校正。本研究使用了空间分辨率为500m且空间位置在全球正弦投影系统
中编号为h26v05的数字图像,数据格式为EOS-HDF,时间序列为2005-2007年1-12月、2001-2004年和
2008年4-10月的每日数据,数字图像值域范围为-100~16000。
1.4 MODIS植被指数合成数据
MODIS归一化差值植被指数NDVI的计算公式为:
NDVI=(NIR-Red)/(NIR+Red) (1)
式中,NIR和Red分别为红外波段和红光波段经过大气校正的地面反射率。
MODIS增强型植被指数EVI是对AVHRR-NDVI的继承和改进,具体计算公式为:
EVI=2.5(NIR-Red)/(NIR+犆1×Red-犆2×Blue+犔) (2)
式中,NIR、Red和Blue分别为经过大气校正的近红外、红光和蓝光通道的反射值;犔=1,为土壤调节参数,参数
犆1和犆2分别为6.0和7.5。
31第18卷第6期 草业学报2009年
每日Terra/MODIS植被指数NDVI和EVI的具体计算方法如下:1)图像投影坐标转换。使用 MRT(MO
DISReprojectionTools)软件,将下载的 MOD09GA产品500m分辨率的1~7通道地表反射率数据的地图投影
SIN转换为 WGS84/Albers地图投影系统,把 HDF格式转化为Tiff格式;2)图像数据格式转换。利用ArcGIS
软件,将1~7通道地表反射率数据的Tiff格式转换为ArcGIS的grid格式;3)每日植被指数图像的合成。利用
ArcGIS的GRID模块及公式(1)和(2),分别计算每日Terra卫星的NDVI和EVI数字图像;4)多日 MODIS最
大植被指数图像的合成。采用最大值合成方法(MVC:MaximumValueComposite),生成研究区2005-2007年
1-12月和2001-2004年及2008年4-10月的NDVI和EVI旬、月、年最大植被指数合成数字图像。
1.5 草地植被地上生物量遥感动态监测模型
在建立草地植被地上生物量遥感动态监测模型时,为了减少气象条件(如降水、阴天多云等)对植被指数的影
响,并且使地面实测样方资料同多日 MODIS合成数据在时间尺度上能够更好地相匹配,本研究对数据进行了如
下处理:1)基础数据库的建立。利用ArcGIS9.1软件,建立了甘南地区2006-2008年地面调查样点空间分布、
行政分区、草地资源类型等数据库。2)每日地面反射率值的提取。利用ArcGIS软件和2006-2008年地面调查
样方点数据库,分别提取各样方点在7-8月的每日地面反射率值;3)样方点10日最大植被指数的计算。利用
MSExcel统计分析软件,计算各样方点每日植被指数NDVI和EVI值。依据各样方点外业具体观测日期及其
之前的9d的植被指数数据,合成10日最大植被指数NDVI和EVI值。4)遥感监测模型的建立。利用SPSS统
计分析软件、各样方点实测的草地鲜草产量数据及相应的10日最大植被指数NDVI和EVI值,建立 NDVI和
EVI与草地地上生物量遥感监测模型,比较分析和选择最优模型。5)模型精度评价。利用误差统计公式[式
(3)],对遥感监测模型的精度进行综合评价。
犘=100×|犔实测值-犔反演值|/犔实测值 (3)
式中,犘为平均估产误差(%),犔实测值和犔反演值分别代表地面实测的草地鲜草产量和遥感监测模型反演的数值。
1.6 草地地上生物量与载畜量变化动态
利用研究区草地生物量遥感动态监测模型及ArcGIS软件空间建模等分析模块,主要研究:1)甘南全州及各
县市不同草地类型在2001-2008年冷季(1-4月、11-12月)和牧草主要生长时期(5-10月)的旬、月和年最大
生物量及生物量密度变化动态;2)利用5-10月期间甘南及各县市草地月最大地上生物量鲜重平均值,结合各类
草地的面积、可利用草地面积系数、可食牧草系数、草地放牧利用系数数据,对比分析2001-2008年期间可食产
草量和理论载畜量的变化特点。
在全年放牧条件下,理论载畜量的计算公式如下:
犆=


犻=1
(犃犻×狔犻×犽犻,1×犽犻,2×犽犻,3)
犐×犇
(4)
式中,犆为某区域内各类型草地总载畜量(SHU,标准羊单位);犃犻为第犻类草地面积(hm2),犻=1,2,3,……,狀,代
表不同的草地类型;狔犻为第犻类草地的牧草单产(kg/hm2);犽犻,1、犽犻,2和犽犻,3分别为第犻类草地的可利用面积系数、
可食牧草系数和草地放牧利用率;犐为1只标准成年绵羊日采食量(kg/D),本研究按5kg鲜草计算;犇为草地放
牧利用的天数,按365d计算。
2 结果与分析
2.1 草地地上生物量遥感监测模型
地上生物量是草地资源动态监测的重要指标之一。根据各样方的牧草生物量和其对应的NDVI、EVI等植
被指数值,通过SPSS11.0统计软件,采用一元线性回归模型、等比级数曲线模型(Growth生长曲线模型)、指数
函数、对数函数、乘幂函数和多项式函数模拟分析植被指数同牧草生物量之间的相关关系,研究结果表明,乘幂函
数可以较好地模拟 MODIS/NDVI、EVI与草地鲜重之间的相关关系。EVI与草地鲜重的相关性大于NDVI与
草地鲜重的相关性,相关系数达0.797(表1、图1)。
41 ACTAPRATACULTURAESINICA(2009) Vol.18,No.6
表1 甘南牧区草地植被指数与地上部分生物量鲜重回归模型
犜犪犫犾犲1 犚犲犵狉犲狊狊犻狅狀犿狅犱犲犾狊犫犲狋狑犲犲狀狏犲犵犲狋犪狋犻狅狀犻狀犱犲狓犪狀犱犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊
犻狀犌犪狀狀犪狀狆犪狊狋狅狉犪犾犪狉犲犪(犘<0.001,犖=199)
回归模型 Model 植被指数 VI 模拟方程Equation 犚2 犉
1、线性模型Linearmodel EVI 狔=15860狓-3361 0.584 276.37
狔=犪狓+犫 NDVI 狔=17470狓-7202 0.308 87.55
2、指数模型Exponentialmodel EVI 狔=953.7e2.94狓 0.631 337.21
狔=犪e犫狓 NDVI 狔=421.6e3.38狓 0.361 111.27
3、生长模型 Growthmodel EVI 狔=e(6.86+2.94狓) 0.631 337.21
狔=e(犪+犫狓) NDVI 狔=e(6.04+3.38狓) 0.361 111.27
4、对数模型Logarithmmodel EVI 狔=10880+8835ln狓 0.571 261.72
狔=犪+犫ln狓 NDVI 狔=9637+12660ln狓 0.303 85.65
5、乘幂模型Powermodel EVI 狔=13583狓1.665 0.636 344.64
狔=犪狓犫 NDVI 狔=10996狓2.469 0.362 111.90
6、多项式模型Polynomialmodel EVI 狔=-3943+16330狓-400.6狓2 0.584 137.49
狔=犪+犫狓+c狓2 NDVI 狔=37.7-2390狓+13440狓2 0.309 43.89
图1 甘南牧区草地地上生物量与犖犇犞犐和犈犞犐回归模型比较
犉犻犵.1 犆狅犿狆犪狉犻狊狅狀狅犳狋犺犲狉犲犵狉犲狊狊犻狅狀犿狅犱犲犾狊犫犲狋狑犲犲狀犵狉犪狊狊犾犪狀犱犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊,犖犇犞犐犪狀犱犈犞犐犻狀犌犪狀狀犪狀狆犪狊狋狅狉犪犾犪狉犲犪
从回归分析的结果可以看出(表1),基于 MODIS-EVI的草地植被地上生物量鲜重最优拟合模型为:
狔=13583狓1.665  (犚2=0.636,犘<0.001,犖=199) (5)
式中,狔为草地鲜重(kg/hm2),狓为 MODIS的EVI值。
统计分析的结果表明,式(5)平均估产精度为76.7%,平均高估误差(负向误差)大于低估误差(正向误差)约
10%。当EVI大于0.4时,估产精度可达70.2%~82.5%,但EVI小于0.4时估产精度较低(表2)。由此可见,
反演模型可很好地模拟牧草生长状况较好时期(5-10月)的地上生物量变化动态。对该时期生物量的准确估
算,是进行全年草地理论载畜量、草畜平衡动态研究的关键。
2.2 甘南州草地地上生物量旬变化动态
草地生物量是衡量草地生产力大小的标准之一,也是衡量植物生长发育的一个重要指标。对2001-2008年
甘南州旬草地最大地上总生物量的分析表明,近8年草地旬生物量具有如下变化规律:从4月下旬草地植被开始
生长,4月下旬-8月上旬生物量逐渐增大,7月上旬-8上旬达到产量高峰值;8月中旬-11月下旬生物量逐渐
减小;12月上旬-次年1-4月中旬,草地植被基本停止生长(图2)。可见,甘南牧区草地生长主要集中在5月上
旬-10月下旬,草地旬最大地上生物量数字图像可以客观地反映草地植被生长发育的总体规律。但是,在草地
51第18卷第6期 草业学报2009年
表2 草地地上生物量鲜重最优拟合模型精度评价
犜犪犫犾犲2 犃犮犮狌狉犪犮狔犲狏犪犾狌犪狋犻狅狀犳狅狉狋犺犲狅狆狋犻犿狌犿狊犻犿狌犾犪狋犲犱犿狅犱犲犾狅犳犵狉犪狊狊犾犪狀犱犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊(狔)
EVI值
EVIvalue
样本数
No.sites
EVI均值
Mean
EVI标准差
Std.error
狔均值
Mean狔
狔标准差
Std.error
平均精度
Meanaccuracy
平均误差
Meanerror
正误差(低估误差)
Positiveerror
负误差(高估误差)
Negativeerror
<0.4 14 0.367 0.031 2965.5 1099.0 58.9 41.1 28.7 -68.9
0.4~0.5 25 0.443 0.030 3503.5 1030.7 70.2 29.8 20.8 -38.2
0.5~0.6 63 0.545 0.032 4983.2 1488.0 74.6 25.4 16.0 -35.6
0.6~0.7 57 0.652 0.028 7317.4 2007.6 82.5 17.5 17.7 -17.3
0.7~0.8 31 0.739 0.022 8342.7 2008.6 81.1 18.9 16.5 -21.0
>0.8 5 0.831 0.026 9470.0 2235.8 80.4 19.6 19.8 -19.6
0~0.1 199 0.603 0.126 6200.6 2608.1 76.7 23.3 18.4 -28.6
图2 甘南州2001-2008年草地植被平均旬生物量变化动态
犉犻犵.2 犃狏犲狉犪犵犲10犱犱狔狀犪犿犻犮狊狅犳犵狉犪狊狊犾犪狀犱狏犲犵犲狋犪狋犻狅狀犫犻狅犿犪狊狊犱狌狉犻狀犵2001-2008犻狀犌犪狀狀犪狀犘狉犲犳犲犮狋狌狉犲
F:上旬Firstthirdofthemonth;M:中旬 Middlethirdofthemonth;L:下旬Lastthirdofthemonth
植被主要生长期间的个别旬生物量(如2002年9月中旬、2007年10月中旬)出现较低的异常现象(图3)。在冷
季期间(1月上旬-4月中旬及11月上旬-12月下旬),生物量旬变化也出现中旬比下旬小的异常情况(如2005
年1-2月和12月、2006年1-3月和12月、2007年1-3月)。研究结果表明,出现这种情况的主要原因是部分
旬最大植被指数图像受大范围长时间阴天多云天气状况的严重干扰,其次是因为在部分旬期间放牧家畜数量的
较大变动(如集中出栏等)引起采食量较大幅度的减小。因此,利用草地植被主要生长期(5月上旬-10月下旬)
各月内的旬最大生物量值(或月最大生物量值),可以避免个别旬生物量的异常数值,计算该区全年草地资源生物
量的空间分布动态,为草畜动态平衡和草地资源的综合评价提供科学依据。
2.3 甘南州草地生物量月度变化动态
虽然草地旬最大地上生物量数字图像的时间分辨率高(上旬和中旬为10d、下旬为8~11d),但旬生物量反
演结果受天气状况的较大干扰,因此研究甘南牧区月度生物量变化动态具有重要意义。对2001-2008年草地的
月最大生物量及其密度的统计分析表明,不同草地类型8年1-12月的生物量均值的月度动态变化均呈单峰抛
物线形式,8年平均最大生物量均出现在7月,其平均最大生物量密度达6370.96kg/hm2,但不同年份产量达到
最大值的月份有所变化,主要集中在7-8月。不同草地类型的生物量在冷季(1-4月及11-12月)和生长季节
的暖季时期(5-10月)具有显著差异。在草地植被生长季期间,平均生物量密度最高的为沼泽类,达4885.26
kg/hm2,其次分别为高寒灌丛草甸类、高寒草甸类、温性草甸草原类、低平地草甸类、暖性草丛类和温性草原类。
61 ACTAPRATACULTURAESINICA(2009) Vol.18,No.6
图3 甘南州2001-2008年草地植被旬生物量变化动态
犉犻犵.3 犜犺犲10犱犱狔狀犪犿犻犮狊狅犳犵狉犪狊狊犾犪狀犱狏犲犵犲狋犪狋犻狅狀犫犻狅犿犪狊狊犳狉狅犿2001狋狅2008犻狀犌犪狀狀犪狀犘狉犲犳犲犮狋狌狉犲
F:上旬Firstthirdofthemonth;M:中旬 Middlethirdofthemonth;L:下旬Lastthirdofthemonth
71第18卷第6期 草业学报2009年
在草地植被几乎停止生长的冷季时期,沼泽类的平均地上生物量密度仍然最高,达894.57kg/hm2,其次分别为
高寒灌丛草甸类、高寒草甸类、暖性草丛类、温性草甸草原类、低平地草甸类和温性草原类(表3)。
2.4 甘南州草地地上生物量年度变化动态
对比分析甘南州草地地上生物量2001-2008年度的动态变化特征,可以发现,2005年各县及全州的草地植
被年总生物量最高,全州达129.13亿kg,其次为2006,2007和2002年,分别为113.17,110.66和109.00亿kg。
全州8年平均总地上生物量为109.31亿kg(表4)。分析甘南州8个气象台站2001-2007年的降水量等气象数
据、地上生物量和畜牧业统计资料表明,不同年度间的降水量存在较大差异是引起草地地上生物量年度变异的主
要原因。
2.5 甘南州草地理论载畜量变化动态
表5和6是2001-2008年甘南州及各县市草地主要生长期间的月最大理论载畜量和年最大理论载畜量统
计分析结果。在全年放牧条件下,不同县市因草地面积、产量等存在较大差异,其载畜量在不同月份的差别也较
大(表5)。玛曲、碌曲和夏河3县在5-10月的平均载畜量分别占全州相同月份总载畜量的37.3%,16.7%和
19.3%。由于气候条件和草地生长状况的较大差异,同一地区在不同年度之间的理论载畜量也存在明显的差异。
例如,在2001-2008年,各县市的载畜量在2005年均达到最大值,全州达395.32万SHU,超过8年平均载畜量
约18.9%(表6)。此外,甘南牧区具有比较明确的冬春场和夏秋场,因此除计算出全年不分季节的载畜量外,还
表3 甘南州2001-2008年草地月平均最大地上生物量密度统计结果
犜犪犫犾犲3 犚犲狊狌犾狋狅犳狋犺犲犿狅狀狋犺犾狔犿犪狓犻犿狌犿犵狉犪狊狊犾犪狀犱犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犱犲狀狊犻狋狔犱狌狉犻狀犵
2001-2008犻狀犌犪狀狀犪狀犘狉犲犳犲犮狋狌狉犲 kg/hm2
项目
Item
低平地草甸
Lowland
meadow
暖性草丛
Warmtemperate
herbosa
沼泽
Marsh
温性草原
Temperate
steppe
温性草甸草原
Temperatemeadow
steppe
高寒灌丛草甸
Alpinebrush
meadow
高寒草甸
Alpine
meadow
平均值
Average
value
1月Jan. 553.51 490.17 750.23 472.69 526.47 669.91 627.78 662.99
2月Feb. 537.40 539.21 754.11 491.99 558.41 706.02 652.91 697.67
3月 Mar. 505.56 631.80 757.71 495.45 540.82 711.09 651.01 702.19
4月 Apr. 624.49 1128.15 1025.67 566.85 719.93 843.69 793.94 837.42
5月 May 1408.63 2531.24 2642.29 1407.15 2383.64 1946.96 1863.10 1940.72
6月Jun. 3106.22 3601.47 5658.17 3089.31 4797.93 4462.38 4468.47 4455.25
7月Jul. 4920.01 3896.35 7294.27 4219.68 5410.40 6398.72 6400.20 6370.96
8月 Aug. 4391.90 3313.02 6603.86 3367.80 4252.39 5957.76 5808.60 5905.97
9月Sep. 3417.86 2686.82 4708.31 2474.36 2877.08 4262.89 4049.79 4215.04
10月 Oct. 1942.52 1643.62 2404.68 1384.95 1630.19 2128.65 1975.21 2102.78
11月 Nov. 785.75 938.97 1115.13 628.98 795.18 941.78 889.36 934.18
12月 Dec. 610.15 580.25 819.29 501.49 587.06 709.98 677.57 704.51
5-10月均值
MeanofMaytoOct.
3197.86 2945.42 4885.26 2657.21 3558.60 4192.89 4094.23 4165.12
1-4月均值
MeanofJan.toApr.
555.24 697.33 821.93 506.74 586.41 732.68 681.41 725.06
11-12月均值
MeanofNov.toDec.
697.95 759.61 967.21 565.23 691.12 825.88 783.46 819.34
冷季均值
Meanofcoldseason
626.59 728.47 894.57 535.99 638.76 779.28 732.44 772.20
 1-3月和11-12月的数据为2005-2007年的均值,其余数据均为2001-2008年8年的平均值。ThedatainJanuary,February,March,No
vemberandDecemberaretheaveragevaluesduring2005-2007,andtherestdataaretheaveragevaluesduring2001-2008.
81 ACTAPRATACULTURAESINICA(2009) Vol.18,No.6
应分别计算冷季(1-4月和11-12月)草场和暖季(5-10月)草场的载畜量,这将是今后工作的重点研究内容之
一。
表4 2001-2008年甘南州年地上生物量鲜重统计结果
犜犪犫犾犲4 犢犲犪狉犾狔犳狉犲狊犺狑犲犻犵犺狋狅犳犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犱狌狉犻狀犵2001-2008犻狀犌犪狀狀犪狀犘狉犲犳犲犮狋狌狉犲 ×108kg
项目Item 2001 2002 2003 2004 2005 2006 2007 2008 8年平均 Mean
玛曲县 Maqu 38.64 41.21 37.92 39.54 47.71 42.53 41.53 37.78 40.86
碌曲县Luqu 17.33 18.73 16.85 17.67 20.88 19.35 18.66 17.38 18.36
夏河县Xiahe 18.58 20.45 19.74 20.54 25.41 22.02 22.12 20.86 21.21
合作市 Hezuo 6.35 6.79 6.43 6.77 8.26 7.02 6.93 6.77 6.91
临潭县Lintan 1.84 2.12 1.98 1.94 2.38 2.09 2.04 2.12 2.07
卓尼县Zhuoni 11.96 12.71 11.86 11.89 15.56 13.00 12.38 12.42 12.72
舟曲县Zhouqu 2.13 1.97 2.14 2.01 2.60 2.18 2.14 2.24 2.18
迭部县 Diebu 4.87 5.02 4.71 4.56 6.34 4.97 4.87 4.68 5.00
全州 Total 101.70 109.00 101.64 104.92 129.13 113.17 110.66 104.25 109.31
 草地主要生长期(5-10月)各月最大地上生物量鲜重平均值。Thedataaretheaveragevaluesofthemonthlymaximumabovegroundbiomassdur
ingthemajorgrowingperiodofgrassland(fromMaytoOctober).
表5 甘南州2001-2008年草地生长期间(5-10月)的月平均载畜量变化动态
犜犪犫犾犲5 犇狔狀犪犿犻犮狊狅犳犿狅狀狋犺犾狔犮犪狉狉狔犻狀犵犮犪狆犪犮犻狋狔犻狀犵狉犪狊狊犾犪狀犱犵狉狅狑犻狀犵狆犲狉犻狅犱(犳狉狅犿犕犪狔狋狅犗犮狋狅犫犲狉)
犱狌狉犻狀犵2001-2008犻狀犌犪狀狀犪狀犘狉犲犳犲犮狋狌狉犲 ×104SHU
月份
Month
玛曲县
Maqu
碌曲县
Luqu
夏河县
Xiahe
合作市
Hezuo
临潭县
Lintan
卓尼县
Zhuoni
舟曲县
Zhouqu
迭部县
Diebu
全州
Total
5 54.18 24.06 28.83 9.59 3.62 19.18 5.37 8.87 153.69
6 126.27 58.61 68.79 22.92 7.73 44.30 8.00 17.99 354.61
7 189.23 84.68 99.81 32.53 9.81 60.35 9.28 22.92 508.62
8 180.04 80.30 91.65 29.95 8.46 53.97 8.03 19.80 472.19
9 129.88 57.14 64.56 20.80 5.56 37.94 6.14 14.86 336.88
10 64.32 27.83 31.37 10.27 2.91 18.93 3.59 8.43 167.66
5-10月均值Average 123.99 55.44 64.17 21.01 6.35 39.11 6.74 15.48 332.27
表6 甘南州2001-2008年草地生长季(5-10月)年载畜量统计结果
犜犪犫犾犲6 犇狔狀犪犿犻犮狊狅犳狔犲犪狉犾狔犮犪狉狉狔犻狀犵犮犪狆犪犮犻狋狔犻狀犵狉犪狊狊犾犪狀犱犵狉狅狑犻狀犵狆犲狉犻狅犱(犳狉狅犿犕犪狔狋狅犗犮狋狅犫犲狉)
犱狌狉犻狀犵2001-2008犻狀犌犪狀狀犪狀犘狉犲犳犲犮狋狌狉犲 ×104SHU
项目Item 2001 2002 2003 2004 2005 2006 2007 2008 8年平均 Mean
玛曲县 Maqu 109.15 126.23 116.16 121.09 146.04 130.28 127.20 115.72 123.99
碌曲县Luqu 50.95 56.77 51.07 53.54 63.28 58.65 56.55 52.68 55.44
夏河县 Xiahe 53.30 62.23 60.07 62.51 77.37 67.00 67.32 63.55 64.17
合作市 Hezuo 18.22 20.77 19.68 20.71 25.27 21.47 21.22 20.72 21.01
临潭县Lintan 5.50 6.55 6.12 6.00 7.34 6.46 6.30 6.53 6.35
卓尼县Zhuoni 34.03 39.46 36.83 36.92 48.30 40.36 38.43 38.55 39.11
舟曲县Zhouqu 6.58 6.10 6.63 6.23 8.04 6.76 6.61 6.94 6.74
迭部县 Diebu 14.77 15.57 14.63 14.15 19.67 15.43 15.10 14.53 15.48
全州 Total 292.50 333.69 311.18 321.15 395.32 346.42 338.74 319.21 332.27
91第18卷第6期 草业学报2009年
3 结论
1)MODIS增强型植被指数EVI的乘幂函数可以较好地模拟甘南牧区草地地上生物量鲜重与EVI之间的
相关关系,拟合模型为狔=13583狓1.665(kg/hm2)。该模型平均估产精度为76.7%,可很好地模拟牧草生长状况
较好时期(5-10月)的草地地上生物量变化动态。
2)甘南牧区草地生长主要集中在5月上旬-10月下旬期间,草地旬最大地上生物量数字图像可以较客观地
反映草地植被生长发育的总体规律。但是,草地植被旬生物量受大范围长时间阴天多云天气状况及放牧家畜数
量有较大变动时期的严重影响。
3)不同草地类型的生物量在冷季和暖季时期具有显著差异。2001-2008年不同草地类型的平均月最大生
物量动态变化曲线均呈单峰抛物线形式,不同年份产量达到最大值的月份有所变化,主要集中在7-8月。
4)在全年放牧条件下,甘南州不同县市因气候条件和草地类型、面积、生长状况等存在较大差异,其草地地
上总生物量及载畜量在不同年度之间也存在明显的差别。研究多年草地植被生物量和载畜量变化动态,是草畜
动态平衡研究和草地资源综合评价的基础。
参考文献:
[1] 王正兴,刘闯,HueteA.植被指数研究进展:从AVHRR-NDVI到 MODIS-EVI[J].生态学报,2003,23(5):979986.
[2] 王正兴,刘闯,陈文波,等.MODIS增强型植被指数EVI与NDVI初步比较[J].武汉大学学报,2006,31(5):407427.
[3] HueteA,DidanK,MiuraT,犲狋犪犾.OverviewoftheradiometricandbiophysicalperformanceoftheMODISvegetationindices[J].
RemoteSensingofEnvironment,2002,83:195213.
[4] 卢丽萍,赵成义.基于 MODIS数据不同荒漠植被指数的空间变化研究———以新疆三工河流域为例[J].干旱区地理,2005,
28(3):381385.
[5] 邢旗,刘爱军,刘永志.应用 MODIS-NDVI对草原植被变化检测研究[J].草地学报,2005,13(增刊):1519.
[6] 赵冰茹,刘闯,王晶杰,等.锡林郭勒草地 MODIS植被指数时空变化研究[J].中国草地,2004,26(1):18.
[7] 毕晓丽,王辉,葛剑平.植被归一化指数(NDVI)及气候因子相关起伏型时间序列变化分析[J].应用生态学报,2005,
16(2):284288.
[8] 王正兴,刘闯,赵冰茹,等.利用 MODIS增强型植被指数反演草地地上生物量[J].兰州大学学报,2005,41(2):1015.
[9] 徐斌,杨秀春,陶伟国,等.中国草原产草量遥感监测[J].生态学报,2007,27(2):405413.
[10] YangYH,FangJY,PanYD,犲狋犪犾.AbovegroundbiomassinTibetangrasslands[J].JournalofAridEnvironments,2009,
73:9195.
[11] ZhaoMS,HeinschFA,NemaniRR,犲狋犪犾.ImprovementsoftheMODISterrestrialgrossandnetprimaryproductionglobal
dataset[J].RemoteSensingofEnvironment,2005,95:164176.
[12] TurnerDP,RittsWD,CohenWB,犲狋犪犾.EvaluationofMODISNPPandGPPproductsacrossmultiplebiomes[J].Remote
SensingofEnvironment,2006,102:282292.
[13] LiZG,YuGR,XiaoX M,犲狋犪犾.ModelinggrossprimaryproductionofalpineecosystemsintheTibetanPlateauusing
MODISimagesandclimatedata[J].RemoteSensingofEnvironment,2007,107:510519.
[14] YangF,Ichi K,WhiteM A,犲狋犪犾.Developingacontinentalscalemeasureofgrossprimaryproductionbycombining
MODISandAmeriFluxdatathroughSupportVectorMachineapproach[J].RemoteSensingofEnvironment,2007,110:
109122.
[15] SimsDA,RahmanAF,CordovaVD,犲狋犪犾.AnewmodelofgrossprimaryproductivityforNorthAmericanecosystems
basedsolelyontheenhancedvegetationindexandlandsurfacetemperaturefrom MODIS[J].RemoteSensingofEnviron
ment,2008,112:16331646.
[16] ChasmerL,BarrA,HopkinsonC,犲狋犪犾.ScalingandassessmentofGPPfromMODISusingacombinationofairbornelidar
andeddycovariancemeasurementsoverjackpineforests[J].RemoteSensingofEnvironment,2009,113:8293.
02 ACTAPRATACULTURAESINICA(2009) Vol.18,No.6
[17] PrietoBA,NorthPR,BarnsleyMJ,犲狋犪犾.SatelitedrivenmodelingofNetPrimaryProductivity(NPP):Theoreticalanal
ysis[J].RemoteSensingofEnvironment,2009,113:137147.
[18] 姜立鹏,覃志豪,谢雯.基于单时相 MODIS数据的草地退化遥感监测研究[J].中国草地学报,2007,29(7):3943.
[19] 韩爱惠.用 MODIS植被指数的时间序列分析提取荒漠化敏感区域的方法[J].林业资源管理,2004,(1):5760.
[20] 刘勇洪,牛铮.基于 MODIS遥感数据的宏观土地覆盖特征分类方法与精度分析研究[J].遥感技术与应用,2004,19(4):
217224.
[21] 刘爱军,韩建国.利用遥感技术监测锡林郭勒天然草原利用强度方法初探[J].中国草地学报,2007,29(2):7074.
[22] 宋小宁,赵英时.应用 MODIS卫星数据提取植被-温度-水分综合指数的研究[J].地理与地理信息科学,2004,20(2):
1317.
[23] CleughHA,LeuningR,MuQZ,犲狋犪犾.RegionalevaporationestimatesfromfluxtowerandMODISsatelitedata[J].
RemoteSensingofEnvironment,2007,106(3):285304.
[24] 甘肃省草原总站.甘肃草地资源[M].兰州:甘肃科学技术出版社,1999.
[25] 中华人民共和国农业部畜牧兽医司全国畜牧兽医总站.中国草地资源[M].北京:中国科学技术出版社,1996.
[26] 陈建华.环境恶化:甘南州如何拉生态建设一把[J].执政决策,2005,(2):1819.
[27] 洛桑·灵智多杰.青藏高原甘南生态经济示范区研究[M].兰州:甘肃科学技术出版社,2005.
[28] 王兮之,杜国祯,梁天刚,等.基于RS和GIS的甘南草地生产力估测模型构建及其降水量空间分布模式的确立[J].草业
学报,2001,10(2):95102.
[29] 赵雪雁.高寒牧区草地退化的人文因素研究———以甘南牧区玛曲县为例[J].草业学报,2007,16(6):113120.
[30] 仁青吉,罗燕江,王海洋,等.青藏高原典型高寒草甸退化草地的恢复———施肥刈割对草地质量的影响[J].草业学报,
2004,13(2):4349.
12第18卷第6期 草业学报2009年
犚犲犿狅狋犲犾狔狊犲狀狊犲犱犱狔狀犪犿犻犮狊犿狅狀犻狋狅狉犻狀犵狅犳犵狉犪狊狊犾犪狀犱犪犫狅狏犲犵狉狅狌狀犱犫犻狅犿犪狊狊犪狀犱犮犪狉狉狔犻狀犵犮犪狆犪犮犻狋狔
犱狌狉犻狀犵2001-2008犻狀犌犪狀狀犪狀狆犪狊狋狅狉犪犾犪狉犲犪
LIANGTiangang1,CUIXia2,FENGQisheng1,WANGYing1,XIAWentao1
(1.KeyLaboratoryofGrasslandAgroecologySystem,MinistryofAgriculture;ColegeofPastoral
AgricultureScienceandTechnology,LanzhouUniversity,Lanzhou730020,China;2.Colege
ofEarthandEnvironmentalScience,LanzhouUniversity,Lanzhou730000,China)
犃犫狊狋狉犪犮狋:Monitoringgrasslandabovegroundbiomass(GAB)isnotonlyanimportantcontentfortheresearch
ofspatialpatternsofgrasslandresources,butisalsoabasisfortheintegrationanalysisofthebalancebetween
grasslandforagesupplyandlivestockdemand.Grasslandinvestigationdatafrom2006-2008andTerra/
MODIS(ModerateresolutionImagingSpectroradiometer)dailysurfacereflectanceproductofMOD09GAin
Gannanpastoralarea,wereusedtoestablishbothamonitoringmodelofGAB,andthe10day,monthlyand
yearlydynamicsofbiomassandtheoreticallivestockcarryingcapacityduring2001-2008inGannanPrefecture.
Andeachcountywassimulatedandanalyzed.ResultssuggestedthatthepowerfunctionofMODIS-EVI(En
hancedVegetationIndex)canbestsimulatethegreenyieldofGABinGannanpastoralarea.Themeanpreci
sionofthemonitoringmodelwas76.7%,whichmightsimulatetheGABdynamicsintheperiodfromMayto
Octoberwhenthegrasslandwasatagoodgrowingstage.Itmainlyconcentratedintheperiodfromthefirst10
daysofMaytothethird10daysofOctoberduringthegrasslandgrowthtime.Duringthisperioddigitalimages
ofthemaximumGABina10dayperiodmayobjectivelyreflectthegeneralruleofgrasslandgrowth,butafew
ofthedigitalimagesofthemaximumGABwereseverelyaffectedbyrainyorcloudyweatherconditionwhich
lastedalongtimeandcoveredalargescale,withaconsiderablechangeinthequantityofgrazinglivestock.The
dynamiccurveofthemonthlymaximumGABforeachgrasslandtypeduring2001-2008wascharacterizedbya
unimodalparabolaform.TheaveragemaximumGABappearedinJuly,butthemonththattheGABreached
themaximumvaluevariedfromyeartoyear,mainlyintheperiodfromJulytoAugust.Thereweregreat
differencesintheoveralbiomassesindifferentyearsintheGannanpastoralarea.Theaverageoveralbiomass
duringthoseyearsintheoveralPrefecturewas109.31×108kg.Thehighestvalueappearedin2005(129.1×
108kg),folowedby2006and2007,whichreached113.2×108,110.7×108,and109.0×108kg,respectively.
Duetothedifferencesofweathercondition,grasslandareaandgrowthsituation,therewasasignificantdiffer
enceforthetheoreticallivestockcarryingcapacityineachcounty.
犓犲狔狑狅狉犱狊:Gannanpastoralarea;grasslandabovegroundbiomass;grasslandcarryingcapacity;remotely
sensedmonitoring
22 ACTAPRATACULTURAESINICA(2009) Vol.18,No.6