免费文献传递   相关文献

Heterogeneity and Error Estimation of Grassland Biomass Based on Cyclic Sampling Approach

大样地循环采样的草地生物量空间异质性及误差分析



全 文 :第20卷 第2期
 Vol.20  No.2
草 地 学 报
ACTA AGRESTIA SINICA
   2012年  3月
  Mar.  2012
大样地循环采样的草地生物量空间异质性及误差分析
邴龙飞1,2,邵全琴1*,王军邦1,刘纪远1
(1.中国科学院地理科学与资源研究所,北京 100101;2.中国科学院研究生院,北京 100049)
摘要:利用2008年和2010年青海省草地生物量野外采样结果,用变差函数探讨地下和地上生物量的空间变异情
况,二者存在空间相关性的阈值分别为3720.52和2258.7 5m。结果表明:玛多采样点前后2年地上生物量的变
异程度都高于海北;2010年各样地地上生物量变化相对较小。珍秦/巴塘平均地下生物量高于海北地下生物量,玛
多平均地下生物量最小;海北和玛多前后2年平均地下生物量年际间变异程度低于珍秦/巴塘样地。对海北、玛
多、巴塘和珍秦的地下、地上生物量的误差进行分析,地下生物量和地上生物量误差变异不同。研究不同采样距离
对误差的影响发现,随着采样距离的加大,不论地上还是地下生物量,其误差总体呈增加的趋势,但在2500m采样
间距下,地下生物量、地上生物量以及二者的比值都出现较为明显的低值。
关键词:循环采样;生物量;误差;样方调查;空间异质性
中图分类号:Q-9;Q948;N33    文献标识码:A     文章编号:1007-0435(2012)02-0257-11 
Heterogeneity and Error Estimation of Grassland Biomass
Based on Cyclic Sampling Approach
BING Long-fei 1,2,SHAO Quan-qin1*,WANG Jun-bang1,LIU Ji-yuan1
(1.Institute of Geographic Sciences and Natural Resources Research,CAS,Beijing 100101,China;
2.Graduate University of Chinese Academy of Sciences,Beijing 100049,China)
Abstract:From cyclic sampling results of 2008and 2010in Qinghai,variogram parameters were used to de-
scribe the heterogeneity of above-ground biomass and below-ground biomass.The spatial correlation
threshold was 2258.75and 3720.52m,respectively.The variation of above ground biomass was greater in
Maduo than in Haibei for two years.The variation of above-ground biomass for al the sampling points was
smaler in 2010.The below-ground biomass was larger in Zhenqin and Batang than in Haibei for 2008and
2010.The smalest below-ground biomass was in Maduo.The interannual variation of mean below-ground
biomass was smaler in Haibei and Maduo than in Zhenqin and Batang.Estimation errors of both above-
ground biomass and below-ground biomass were quite different.The estimation errors increased for both
above-ground biomass and below-ground biomass as the sampling distance increased.The lowest value of
above-ground biomass,below-ground biomass and their ratios were al at the sampling distance of 2500
meters.
Key words:Cyclic sampling;Biomass;Error;Quadrat survey;Spatial heterogeneity
  草地是我国最主要的陆地植被类型之一,约占
我国国土面积的41.7%,对于保障我国的国土生态
安全和可持续发展具有至关重要的意义。在我国近
4亿hm2 的天然草原中,有90%的可利用草原已有
不同程度退化,草原生态环境持续恶化,草地过牧超
载是重要原因之一。为此,自2000年以来,国家和
有关地方政府相继实施了“退牧还草”、“天然草原保
护”等工程,出台了《草畜平衡管理办法》、《草原保护
补助奖励机制》等政策和措施。草畜平衡是指为保
持草原生态系统良性循环,在一定时间内,草原使用
者或承包经营者通过草原和其他途径获取的可利用
饲草饲料总量与其饲养的牲畜所需的饲草饲料量保
持动态平衡。要实现草畜平衡,草地产草量的准确
估算非常重要。
收稿日期:2011-10-13;修回日期:2011-12-20
基金项目:国家科技支撑项目(2009BAC61B01)资助
作者简介:邴龙飞(1979-),男,山东青岛人,博士,研究方向为区域水热平衡,E-mail:binglf@lreis.ac.cn;*通信作者 Author for corre-
spondenc,E-mail:shaoqq@lreis.ac.cn
草 地 学 报 第20卷
  草地产草量的估算方法主要有生物量调查法、
遥感估产法、模型模拟方法[1~3]。生物量调查多在
植被生长茂盛的时候,通过刈割和根钻取样,然后烘
干、称重,分别获得地上、地下部分干重[4,5],然后估
算出区域草地生物量。遥感方法通常是利用
AVHRR或 MODIS的 NDVI数据与同期、相应地
点的地面样方生物量建立相关关系,然后利用ND-
VI来估算草地生物量[6~10]。模型模拟方法是利用
遥感数据驱动,利用生态模型模拟初级生产力,然后
再利用地上地下生物量比和根系周转率等估算草地
生物量[11,12]。其中,生物量调查法是草地产草量估
算的传统方法,其结果通常被认为可以代表地面真
实情况,用来为遥感方法和模型方法建立经验关系
和精度验证。
但生物量调查法工作量大,野外作业辛苦,通常
一个样点只设3~6个1m×1m或50cm×50cm、
甚至是25cm×25cm草地样方[1,4,13~15]。农业部草
原监理中心发布的《全国草原监测技术手册》提出草
地样地和样方设置的原则,若样地内只有草本、半灌
木及矮小灌木植物,样方的面积一般为1m2,若样地
植被分布呈斑块状或者较为稀疏,应将样方扩大到
2~4m2,一个样地内不少于3个样方。面积大、地
形复杂、生态变异大,应多设样方。然而,遥感数据
和模型数据的分辨率多是500m×500m,或1km
×1km,因此,传统的草地生物量调查方法获得的
数据,用来为遥感和模型建立经验关系和结果验证
存在着不确定性。
为此,本研究采用大样地循环采样方法,分别于
2008年和2010年在青海海北、玉树、玛多、称多等
地开展了5个大样地的调查,并利用变差函数讨论
生物量的空间变异和误差情况。对比传统方法和大
样地循环采样方法之间的差异,探索不同采样步长
对生物量估算误差的影响,为消减草地生物量估算
的不确定提供参考方法。
1 材料与方法
大样地循环采样方法是由 Burrows等人在
2002年提出,循环采样方案在采样频率较低的情况
下,仍能够在不同的步长(距离)上保持采样的密度,
此方法最早是用于时间序列数据变异特征的研
究[16],后来将其扩展到了二维的空间上[17~20],是研
究叶面积指数、植被生产力有效的采样方法[21]。
本研究分别于2008年和2010年在青海省门源
回族自治县海北通量观测塔附近、玛多县、玉树巴
塘、称多县珍秦设置4个5km×5km或3km×3
km的样地(其中海北重复2次,共5个),进行循环
采样;并于2010年分别在甘德县、格尔木市、共和
县、玛沁县、囊谦县、曲麻莱县、同德县、玉树县、杂多
县、治多县,采用传统方法进行样方采样(表1)。
大样地的采样方法(以海北为例):以通量塔为
中心,沿东南西北四个方向使用4/13循环的方法进
行采样,采样的基本步长为250m。4/13循环采样
就是循环的长度为13个基本步长,只在0,1,3,9,
13基本步长的整数倍处设置采样样方。在每个区
内从南到北和从东到西共设置18条采样线,最南和
最西的采样线记为1,最北和最东的采样线记为9。
地上生物量的采样样方设置在线与线的交叉点处
(图1),这样,一个采样区的总样点数为81个。在每
个采样点设置5个50cm×50cm样方,采集地上生
物量,地下生物量则利用根钻重复3次,分3层采集。
图1 青海海北大样地循环采样点空间分布图
Fig.1 Distribution of cyclic sampling
points in Haibei,Qinghai
  玛多采样区群落类型为高寒草原;巴塘/珍秦采
样区为高寒草甸/草原化草甸类型;海北样点主要植
被类型为高寒草甸、高寒灌丛和沼泽化草甸,土壤为
高山草甸土、高山灌丛草甸土和沼泽土。具体的物
种类型主要有嵩草(Kobresia Wild)、火绒草(Leon-
topodium leontopodioides)、针茅(Stipa)、早熟禾
(Poa)、羊茅(Festuca ovina)等,各样地详细的物种
分布情况如表2所示。
852
第2期 邴龙飞等:大样地循环采样的草地生物量空间异质性及误差分析
表1 采样点基本情况
Table 1 Basic information of sampling points
日期
Date
地点   
Site/station   
生态系统      
Ecosystem      
观测内容  
Observing content  
样方数/采样数
Quadrat number/Sampling number
2008/7/27-8/7
海北站
Haibei
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  81×5=405
地下生物量 Under-ground biomass  81×3×2层=486
LAI  81×3=243
土壤取样Soil sampling  81×3×2层=486
2008/8/16-20
玉树巴塘
Batang,Yushu
高寒草甸/草原化草甸
Alpine meadow/Steppe meadow
地上生物量 Above-ground biomass  81×5=405
地下生物量 Under-ground biomass  81×3×2层=486
土壤取样Soil sampling  81×3×3层=729
2008/8/21-24
玛多县
Maduo county
高寒草原
Alpine steppe
地上生物量 Above-ground biomass  81×5=405
地下生物量 Under-ground biomass  81×3×2层=486
土壤取样Soil sampling  81×3×3层=729
2010/8/12、14
称多县
Chenduo county
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  63×5=315
地下生物量 Under-ground biomass  63×3×2层=378
2010/8/8
甘德县
Gande county
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  1×5=5
地下生物量 Under-ground biomass  1×3×2层=6
2010/8/26、29
格尔木市
Golmud city
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  7×5=35
地下生物量 Under-ground biomass  7×3×2层=42
2010/8/30
共和县
Gonghe county
高寒草甸/草原
Alpine meadow/Alpine steppe
地上生物量 Above-ground biomass  2×5=10
地下生物量 Under-ground biomass  2×3×2层=12
2010/8/9、11
玛多县
Maduo county
高寒草原
Alpine steppe
地上生物量 Above-ground biomass  21×5=105
地下生物量 Under-ground biomass  21×3×2层=126
2010/8/6-7
玛沁县
Maqin county
高寒灌丛/高寒草甸
Alpine shrubs/Alpine meadow
地上生物量 Above-ground biomass  3×5=15
地下生物量 Under-ground biomass  3×3×2层=18
2010/7/23-24
海北站
Haibei
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  77×5=385
地下生物量 Under-ground biomass  67×3×2层=402
2010/8/22
囊谦县
Nangqian county
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  4×5=20
地下生物量 Under-ground biomass  4×3×2层=24
2010/8/25
曲麻莱县
Qumarleb county
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  3×5=15
地下生物量 Under-ground biomass  3×3×2层=18
2010/8/6
同德县
Tongde county
草原
Alpine steppe
地上生物量 Above-ground biomass  4×5=20
地下生物量 Under-ground biomass  4×3×2层=24
2010/8/23-24
玉树县
Yushu county
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  2×5=10
地下生物量 Under-ground biomass  2×3×2层=12
2010/8/23
杂多县
Zaduo county
高寒草甸
Alpine meadow
地上生物量 Above-ground biomass  1×5=5
地下生物量 Under-ground biomass  1×3×2层=6
2010/8/24
治多县
Zaduo county
高寒草原/高寒草甸
Alpine steppe/Alpine meadow
地上生物量 Above-ground biomass  3×5=15
地下生物量 Under-ground biomass  3×3×2层=18
  首先对地上生物量和地下生物量的空间变异进
行研究,得到变量自相关性影响的范围,在此基础上
运用统计方法研究不同采样间距下误差的分布。
对生物量空间变异的研究采用变异函数描述。
变异函数又称半变差函数、半变异矩,以区域化变量
理论为基础,分析自然现象的空间变异和空间相关
特征,被认为是描述生态学特征空间异质性的有效
方法[22~24]。通过变异函数可以得到3个重要的参
数:块金值(Nugget)、基台值(Sil)和变程(Range)。
随着计算间距h的增加,变异函数γ(h)逐渐趋于稳
定。常数(C0+C)称为基台值,表示系统属性或区
域化变量最大变异,其值越大变量的空间异质性程
度越高。而块金值则表示最小抽样尺度以下变量的
变异性及测量误差,表示随机部分的空间异质性。
上面2个参数都无法比较不同变量之间的随机性方
面的差异,在比较不同变量的空间变异性时,多采用
952
草 地 学 报 第20卷
表2 各样地物种类型
Table 2 Species of each sampling site in Qinghai
海北     
Haibei     
拉丁名   
Latin name   
玛多   
Maduo   
拉丁名    
Latin name    
巴塘/珍秦
Batang/Zhenqin
拉丁名   
Latin name   
矮嵩草 Kobresia humlis 小嵩草 Kobresia pygmaea 小嵩草 Kobresia pygmaea
小蒿草 Kobresia pygmaea 紫花针茅 Stipapurpurea 紫花针茅 Stipapurpurea
线叶嵩草 Kobresia duthiei 冰草 Agropyron cristatum 火绒草 Leontopodium leontopodioides
藏嵩草 Kobresia tibetica 点地梅 Androsace umbellata 雪白委陵菜 Potentilla nivea
高山唐松草 Thalictrum alpinum 细叶亚菊 Ajania tenuifolia 多裂委陵菜 Potentilla multifida
金露梅灌丛 Potentilla fruticosa shrub 火绒草 Leontopodium leontopodioides 二裂委陵菜 Potentilla bifurca
异针茅 Stipa aliena 二裂委陵菜 Potentilla bifurca 黄帚橐吾 Ligularia virgaurea
早熟禾 Poa annua 黄芪 Astragalus membranceus 高山唐松草 Thalictrum alpinum
美丽风毛菊 Saussurea pulchra 早熟禾 Poa annua 细叶亚菊 Ajania tenuifolia
羊茅 Festuca ovina 矮嵩草 Kobresia humlis
短枝发草 Deschampsia littoralis 杜鹃 Rhododendron simsii
钝苞雪莲(瑞苓草) Saussurea nigrescens 藏嵩草 Kobresia tibetica
垂穗披碱草 Elymus nutans
块金与基台值的比值(C0/(C0+C)),表示系统变量
空间自相关变异所占的比例,隐含系统变量的空间
相关性的程度,如果该比值高,说明样本间的变异更
多的是由随机因素引起的。基台值对应的采样距离
称为变程,超过变程后,各采样点统计属性的空间自
相关性消失[23]。
在得到变量空间变异特征的基础上,统计各采
样点地上生物量和地下生物量的平均值、极大值和
极小值,得到生物量总体的变化情况。
野外工作条件往往限制采样的频率和时间,大
样地循环采样方法虽然考虑了样点之间的自相关
性,但是采样的工作量大,对样地破坏性也较大。本
研究对2008年和2010年2年的采样点按500,
1000,1500,2000,2500和3000m间距处理,落于对
应大小网格内的点,取网格内各样点生物量的均值
作为此栅格所代表的生物量(图2),以大样地循环
采样的结果为真值,依次讨论扩大空间采样距离后
相对误差的变化规律,计算各采样点地上和地下生
物量与上述真值的差,以此作为地上/地下生物量的
估算误差,以误差下限2%,5%,10%,15%,20%,
25%和30%为横坐标,以落于此范围内的样点个数
为纵坐标绘制出误差变异曲线,目的是研究适合三
江源地区的,介于传统采样和大样地循环采样之间
的折衷的采样方法。
图2 大样地循环采样样点数处理方式示意图
Fig.2 Processing methods of raw data by cyclic sampling in Qinghai
2 结果与分析
2.1 生物量变化的空间异质性
对所有采样点得到的地上生物量和地下生物量
分别计算不同间距下的变异函数,发现地下生物量
变程高于地上生物量,地下生物量变程约为3720
m,而地上生物量约为2258m,说明采样点间距在
3720m以上时,地下生物量不再具有空间自相关
性,而地上生物量的这个距离为2258m。
由块金值和基台值的比值来看,地下生物量对
应的这个比值也大于地上生物量,说明地下生物量
的空间变异受随机因素的影响更强,在变程内二者
都具有中等强度的空间自相关性。这也隐含了在变
程范围内,似乎可以探索其他采样间距,使得在此间
距下,既维持了一定的采样精度,又能在某种程度上
减少采样的工作量。
2.2 各样地的生物量变化
统计发现,2008年和2010年海北平均地上生
物量分别为284.48和225.04g·m-2左右。2008
062
第2期 邴龙飞等:大样地循环采样的草地生物量空间异质性及误差分析
年海北样地最大地上生物量为363.78g·m-2,
2010年最大为394.83g·m-2。2008年海北样地
最小地上生物量为224.12g·m-2,2010年最小为
69.20g·m-2。
玛多县2008年采用大样地循环采样的方法,而
2010年则主要是沿路采样,其地上生物量实测值极
值分布范围较大,2008年最大值和最小值分别为
136.14和13.85g·m-2;2010年最大值和最小值
分别为323.74和53.78g·m-2。与海北相比,不
论平均地上生物量还是最大地上生物量,玛多采样
点前后2年地上生物量的变异程度都高于海北采样
点(表3和图3-a,b)。
图3 2008年和2010年各样地平均地上、地下生物量变化
Fig.3 Comparison of mean above-and below-ground biomass variation between 2008and 2010in Qinghai
  2010年玉树县内仅有2个采样点,且距2008
年采样点较远,因此,选择距离2008年巴塘(玉树)
采样点最近的珍秦采样点进行分析,后者在2010年
采用了大样地循环采样法,由表3可知,2008年巴
塘平均地上生物量为119.57g·m-2,2010年珍秦
平均地上生物量为171.73g·m-2。
不同样地估算的平均地下生物量也有明显差
异,但是2010年各样地平均地下生物量变化程度比
2008年小。从前后2年的计算结果来看,珍秦/巴
塘平均地下生物量高于海北对应年份地下生物量,
而玛多对应年份平均地下生物量最小。从年际间变
异来看,2010年各样地地下生物量变化相对较小;
海北和玛多前后2年平均地下生物量年际间变异程
度低于珍秦/巴塘样地(表3和图3-b)。
表3 2008年和2010年生物量变化情况
Table 3 Variation of above-ground and below-ground biomass of 2008and 2010in Qinghai,respectively
项目Content           年份Year 海北 Haibei 玛多 Maduo 珍秦/巴塘Zhenqin/Batang*
平均地上生物量 2008  284.48  49.24  119.57
Mean above-ground biomass/g·m-2  2010  225.04  133.33  171.73
最大地上生物量 2008  363.78  136.14  213.88
Maximum above-ground biomass/g·m-2  2010  394.83  323.74  456.14
最小地上生物量 2008  224.12  13.85  50.84
Minimum above-ground biomass/g·m-2  2010  69.20  53.78  49.50
平均地下生物量 2008  2077.94  957.61  3382.27
Mean under-ground biomass/g·m-2  2010  1575.81  1048.42  1759.34
最大地下生物量 2008  4769.87  3135.23  5590.46
Maximum under-ground biomass/g·m-2  2010  6979.50  2670.66  2856.67
最小地下生物量 2008  1019.69  172.14  1945.17
Minimum under-ground biomass/g·m-2  2010  561.74  219.54  859.31
  注:*:珍秦为称多县样地,巴塘为玉树县样地。下同
Note:*:Zhenqin:sampling site of Chengduo county.Batang:sampling site of Yushu county.The same as below
2.3 不同样地相对误差分布
2.3.1 海北 2008年海北各采样点估算的地上生
物量相对误差变异极大,相对误差较大的点集中在
中心点以外,2008年地上生物量最大相对误差为
26.93%,最小相对误差为1.83%(图4-a)。2010
年地上生物量相对误差最大值和最小值分别为
77.06%和0.46%(表4),相对误差低于10%的采
样点既有中心点又有周围部分;南部采样点相对误
162
草 地 学 报 第20卷
差较大,多数高于10%(图4-b)。由相对误差分布
情况可以看出,多数采样点地上生物量相对误差在
30%以下。在误差区间[0,10]内,随着值的增加,相
对误差落于此区间内的采样点数迅速增加,但是增
加的速率逐渐趋缓(图5-a,b)。
图4 2008年(a)和2010年(b)海北地上生物量相对误差
Fig.4 Relative errors of above-ground biomass in Haibei of 2008(a)and 2010(b)
表4 各样地地下和地上实测生物量相对误差
Table 4 Relative errors of above-ground and below-ground measured biomass
年份
Year
项目
Content
海北地下
Under-ground
biomass in Haibei
海北地上
Above-ground
biomass in Haibei
玛多地下
Under-ground
biomass in Maduo
玛多地上
Above-ground
biomass in Maduo
巴塘地下
Under-ground
biomass in Batang
巴塘地上
Above-ground
biomass in Batang
2008
平均值 Mean -7.09 -1.50 -56.46  0.00  13.78  17.16
标准差SD  27.66  12.31  123.10  50.83  22.81  37.54
最大值(绝对值)
Maximum(absolute)
103.78  26.93  69.65  74.53  73.88  135.16
最小值(绝对值)
Minimum(absolute)
0.17  1.83  0.13  0.09  0.09  0.83
2010
平均值 Mean -0.44  0.92  0.00  0.00 -9.55  21.61
标准差SD  83.47  30.38  66.50  46.18  30.61  72.50
最大值(绝对值)
Maximum(absolute)
340.98  77.06  154.73  142.82  165.68  223.01
最小值(绝对值)
Minimum(absolute)
0.10  0.46  3.14  2.89  0.95  3.33
图5 2008年(a)和2010年(b)海北不同误差区间采样点个数
Fig.5Sampling numbers of different error range in Haibei of 2008(a)and 2010(b)
  地下生物量相对误差分布情况和地上部分有很
大差别,2008年地下生物量相对误差最大值和最小
值分别为103.78%和0.17%(图6-a)。2010年多
数采样点平均地下生物量都低于真值,少数的几个
高于真值的采样点离群程度都很大(图6-b),地下生
物量相对误差的最大值和最小值分别为340.98%和
0.10%(表4),导致采样点间地下生物量出现较大
变异。从相对误差的分布情况来看,2008年和2010
262
第2期 邴龙飞等:大样地循环采样的草地生物量空间异质性及误差分析
年,落于误差区间内的采样点个数随区间下限的增
大而迅速增加,但是增加的速率不同,2008年,当相
对误差超过30%以后,落于此误差范围内的采样点
个数趋缓(图7-a);而2010年采样结果则正好相
反,误差限内采样点数随下限的增加而迅速增加(图
7-b)。
2.3.2 玛多 图9-a,b所示为前后2年玛多地上生
物量相对误差的分布情况。2008年部分采样点相对
误差偏正,这些点主要分布在中心部分,而中心点以
外的采样点相对误差偏负,2008年玛多地上生物量
相对误差的最大值和最小值分别为74.53%和
0.09%。2010年各采样点相对误差分布较为均匀,主
要为沿路采样,涉及的空间范围较大,且样点个数较
少(图8-a),采样点间距超过空间自相关的最大距离。
2008年玛多各采样点地下生物量相对误差变
异程度高于同年地上生物量相对误差情况(图9-a、
图10-a),其中地下生物量相对误差偏负的点,其绝
对值也相对较大。2010年采样结果显示,多数采样
点地下生物量较平均地下生物量偏低,而数个偏正
的点,其误差都较大(图10-b)。
图8 玛多(a)和巴塘(b)采样点空间分布情况
Fig.8 Distribution of sampling points in Maduo(a)and Batang(b)of Qinghai
362
草 地 学 报 第20卷
2.3.3 巴塘和珍秦 从2008年地上生物量的计算
结果看出(图11-a),中心点相对误差偏大,离中心点
越远,相对误差越小,至外围,大部分点的相对误差
为负,2008年地上生物量相对误差最大值和最小值
分别为135.16%和0.83%(表4);由相对误差的分
布结果(图11-c)可以看出,巴塘采样点在25%的误
差下限内,随着值的增加,落于此限内的采样点个数
迅速增加,至25%趋于稳定。
2010年巴塘采样点较少(图8-b),因此选择邻
近的称多县珍秦采样点分析,珍秦各样点地上生物
量估算误差相对均匀(图11-b)。随着误差区间下
限的增加,落于此限内的采样点个数也随之迅速增
加,不存在采样点个数比较稳定的区间下限(图
11-d)。
图11 巴塘2008年(a,c)和珍秦2010年(b,d)各采样点地上生物量相对误差
Fig.11 Relative errors of above-ground biomass in Batang of 2008(a,c)and in Zhenqin of 2010(b,d)
462
第2期 邴龙飞等:大样地循环采样的草地生物量空间异质性及误差分析
  2008年,巴塘地下生物量估算的相对误差分布
比较均匀(图12-a),从各误差区间样点分布情况可
以看出,样点个数随误差区间的增加几乎呈线性关
系增加(图12-b)。由2010年采样结果可以看出,
珍秦各点地下生物量相对误差分布亦比较均匀(图
13-a),样点个数随误差区间的增大亦几乎呈线性关
系增加(图13-b)。
2.4 采样距离对相对误差的影响
在500,1500和2500m处,地上生物量相对误
差出现低值,皆低于1000m采样间距地上生物量
相对误差大小(图14-a和图15-a)。而地下生物量
相对误差低值主要分布在500和2500m采样间距
下(图14-b和图15-b)。
图14 海北不同间距地上(a)/地下(b)生物量相对误差
Fig.14 Relative errors of above- (a)and below-ground(b)biomass under different sampling distance in Haibei of Qinghai
  如果以大样地循环采样的结果作为真值,与不
同采样间距下相对误差的分布进行对比,可以看出
随着采样距离的增加,不论地上生物量还是地下生
物量,其相对误差总体呈增加趋势,但是在2500m
采样间距下,地下生物量、地上生物量以及二者的比
值,其误差都出现较为明显的低值,相对误差多数低
于30%(图14和图15)。
562
草 地 学 报 第20卷
图15 珍秦不同间距地上(a)/地下(b)生物量相对误差
Fig.15 Relative errors of above-ground(a)and below-ground(b)biomass
under different sampling distance in Zhenqin of Qinghai
3 讨论与结论
影响生物量的主要因素有环境因子和人为因
子,前者包括水分、温度和光照条件等,而人为因子
则包括放牧、刈割等利用方式以及灌溉、施肥等抚育
措施,不同的样点布局方式对生物量计算结果影响
极大。
由青海省的采样结果可知,不同样点地上生物
量差异较大,2008年和2010年玛多各采样点平均
地上生物量和最大地上生物量估算结果的变异程度
都高于海北样点,巴塘和珍秦前后2年地上生物量
估算的相对误差则无太大差异。不同样地估算的平
均地下生物量亦有明显差异,但是2010年各采样点
地下生物量变化程度比2008年小。珍秦和巴塘平
均地下生物量高于海北对应年份地下生物量,对应
年份平均地下生物量最小的是玛多。从年际间变异
来看,2010年各样地的地下生物量变化相对较小;
海北和玛多前后2年平均地下生物量年际间变异程
度低于珍秦和巴塘样地。
2008年海北不同采样点地上生物量估算相对
误差差异极大,误差较大的点集中在中心点以外;
2010年误差分布相对均匀。2008年海北地下生物
量误差分布亦相对均匀,但2010年则变化较大。
2008年玛多地上生物量,中心点相对误差偏正,中
心点以外的采样点相对误差偏负;地下生物量估算
的相对误差高于同时段地上生物量估算的相对误
差。2008年巴塘地上生物量估算,中心点相对误差
偏大,离中心点越远,相对误差越小,外围大部分点
的相对误差为负;在25%的误差下限内,随着值的
增加,落于此限内的采样点个数迅速增加,至25%
~30%趋于稳定。珍秦各样点地上生物量估算误差
较为均匀,随着误差区间下限的增加,落于此限内的
采样点个数也随之增加,不存在采样点个数比较稳
定的误差区间。巴塘和珍秦地下生物量估算的相对
误差分布则相对均匀。
受环境异质性和采样方法的限制,即使在同一
研究地域,其研究结果仍然会有较大差异。而且在
大量实验的时候,由于成本和可达性的限制,不得不
人为舍弃一些较典型的采样点,即使随机采样也无
法完全避免人为的干扰因素。在研究区空间范围较
大的时候,不得不加大采样间距,以减少采样次数,
这样做不可避免地会增加误差,大样地循环采样兼
顾了变量的空间异质性和空间自相关性,但是在环
境严苛的地区仍显复杂。在分析青海省地上和地下
生物量空间异质性特征的基础上,本研究对2年的
采样结果按500,1000,1500,2000,2500和3000
m间距处理,对地上和地下生物量估算误差的分析
发现:随着采样距离的加大,其相对误差呈增加的趋
势,但是在2500m采样间距下,地下生物量、地上
生物量以及二者的比值都出现较为明显的低值,其
相对误差多数低于30%。
在进行生物量调查分析的时候,大样地循环采
样能够兼顾采样密度和采样精度,应该作为野外研
究的首选方案。如果在不同时段内多次取样,在一
定的误差允许范围内,可以加大采样间距,以减少工
作量,这样做的前提是要至少有前期一年,此年内样
点的空间分布范围较大,该方法假定通过采样方法
估算的生物量能够代表该年度研究区域内生物量的
实际情况。而针对三江源的研究发现,2500m间距
采样可以明显降低采样频数,同时其相对误差又能
很好地控制在30%以内。
参考文献
[1] 安卯柱,高娃,朝鲁.内蒙古第四次草地资源调查草地生产力
测定及计算方法简介[J].内蒙古草业,2002,14(4):20-21
[2] 徐斌,杨秀春,陶伟国,等.中国草原产草量遥感监测[J].生态
学报,2007,27(2):405-413
[3] 杨秀春,徐斌,朱晓华,等.北方农牧交错带草原产草量遥感监
测模型[J].地理研究,2007,26(2):213-221
662
第2期 邴龙飞等:大样地循环采样的草地生物量空间异质性及误差分析
[4] 樊江文,杜占池,钟华平.红三叶、鸭茅生物量和叶面积时空结
构特征[J].草地学报,2004,12(3):204-208,213
[5] 宋璐璐,樊江文,钟华平,等.重庆红池坝地区草地群落生物量
及物种丰富度沿海拔梯度格局变化[J].草地学报,2010,18
(2):160-166
[6] Skinner R H,Wylie B K,Gilmanov T G.Using normalized
difference vegetation index to estimate carbon fluxes from
smal rotationaly grazed pastures[J].Agronomy Journal,
2011,103(4):972-979
[7] Chen J,Gu S,Shen M,et al.Estimating aboveground bio-
mass of grassland having a high canopy cover:an exploratory a-
nalysis of in situ hyperspectral data[J].International Journal
of Remote Sensing,2009,30(24):6497-6517
[8] 陈世荣,王世新,周艺.基于遥感的中国草地生产力初步计算
[J].农业工程学报,2008,24(1):208-212
[9] 罗玲,王宗明,任春颖,等.基于 MODIS数据的松嫩草原产草
量遥感估算模型与空间反演[J].农业工程学报,2010,26(5):
182-187
[10]毛飞,张艳红,侯英雨,等.藏北那曲地区草地退化动态评价
[J].应用生态学报,2008,19(2):278-284
[11]Fan J,Zhong H,Harris W,et al.Carbon storage in the
grasslands of China based on field measurements of above-and
below-ground biomass[J].Climatic Change,2008,86(3/4):
375-396
[12]樊江文,邵全琴,刘纪远,等.1988-2005年三江源草地产草
量变化动态分析[J].草地学报,2010,18(1):5-10
[13]陈生云,赵林,秦大河,等.青藏高原多年冻土区高寒草地生物
量与环境因子关系的初步分析[J].冰川冻土,2010,32(2):
405-413
[14]王启基,周兴民.高寒小嵩草草原化草甸植物群落结构特征及
其生物量[J].植物生态学报,1995,19(3):225-235
[15]王长庭,曹广民,王启兰,等.青藏高原高寒草甸植物群落物种
组成和生物量沿环境梯度的变化[J].中国科学C辑:生命科
学,2007,37(5):585-592
[16]Clinger W,Van Ness J W.On unequaly spaced time points in
time series[J].The Annals of Statistics,1976,4(4):736-745
[17]Clayton M K,Hudelson B D.Confidence intervals for autocor-
relations based on cyclic samples.[J].Journal of the American
Statistical Association,1995,90(430):753-757
[18]Miler T F,Mladenoff D J,Clayton M K.Old-growth north-
ern hardwood forests:spatial autocorrelation and patterns of
understory vegetation[J].Ecological Monographs,2002,72
(4):487-503
[19]Liang L,Schwartz M D,Fei S.Validating satelite phenology
through intensive ground observation and landscape scaling in a
mixed seasonal forest[J].Remote Sensing of Environment,
2011,115(1):143-157
[20]李正泉.陆地生态系统生产力的多尺度分析与尺度转换方法
研究[D].北京:中国科学院研究生院,2006:66-68
[21]Burrows S N,Gower S T,Clayton M K,et al.Application of
geostatistics to characterize leaf area index(LAI)from flux
tower to landscape scales using a cyclic sampling design[J].E-
cosystems,2002,5(7):667-679
[22]Rossi R E,Mula D J,Journel A G,et al.Geostatistical tools
for modeling and interpreting ecological spatial dependence[J].
Ecological Monographs,1992,62(2):277-314
[23]李哈滨,王政权.空间异质性定量研究理论与方法[J].应用生
态学报,1998,9(6):651-657
[24]Robertson G P.Geostatistics in ecology:Interpolating with
known variance[J].Ecology,1987,68(3):744-748
(责任编辑 刘云霞)
762