以动态监测样地(340 m×300 m)的网格(5 m×5 m)为基本单元,选择5个群落结构属性变量,采用分形几何学和地统计学方法对祁连山大野口流域青海云杉林群落结构的空间异质性进行研究.结果表明: 青海云杉林5个结构属性的变异程度大小依次为:密度>平均冠幅>显著度>盖度>平均树高,变异系数为43.7%~79.6%.Moran I系数表明,各个结构属性均有不同程度的空间自相关性,自相关大小顺序为:密度>平均树高>盖度>平均冠幅>显著度,变化范围为-0.047~0.382.指数理论变异函数模型能很好地拟合不同结构属性的空间变异,变程为24.6~68.1 m,各属性变量的结构比除盖度表现为中等空间相关,其他指标均为强烈空间相关,各属性指标分维数接近2,空间依赖性较小.植被密度和盖度空间分布呈带状结构和斑块结构叠加的特点,其他指标呈较强的斑块状空间结构,密度和盖度对平均冠幅、显著度和平均树高有一定的空间依赖性.青海云杉林群落结构空间异质性适宜的采样间隔和采样面积分别为10 m和0.5 hm2.
We selected the grid of 5 m × 5 m in a dynamic monitoring plot (340 m × 300 m) as the sampling unites and chose 5 structural characteristics (density, average crown breadth, coverage, conspicuousness and average height) to study the spatial heterogeneity of community structure of Picea crassifolia forest in Dayekou Basin of Qilian Mountains by the fractal geometry and geostatistics methods. The results showed that the order of spatial variation in these characteristics was: density > average crown breadth > conspicuousness > coverage > average height, with the variation coefficient ranging from 43.7% to 79.6%. Moran’s I index indicated that the structural variables had different degrees of spatial autocorrelation, and the order of autocorrelation was density > average height> coverage > average crown breadth > conspicuousness, with the range of -0.047-0.382. The exponential semivariation model well fitted the spatial variability in different structural features, and the range was 24.6-68.1 m. The variables displayed moderate spatial autocorrelation except for coverage, while the other variables had strong spatial autocorrelation, and the fractal dimension of the variables was close to 2, indicating a low spatial dependence among variables. The variables presented a superposing characteristic of zonal and patchy structures except for density and coverage, while the other variables presented strong patchiness property. Density and coverage had a certain spatial dependence on average crown breadth, conspicuousness and average height. Density and coverage for the spatial heterogeneity of community structural of P. crassifolia forests were 10 m and 0.5 hm2, respectively.
全 文 :祁连山青海云杉林群落结构的空间异质性∗
赵维俊1,2,3 刘贤德1,2,3∗∗ 敬文茂1,2,3 徐丽恒4 牛 赟2,3 齐 鹏2,5 赵永宏1
( 1甘肃农业大学林学院, 兰州 730070; 2甘肃省祁连山水源涵养林研究院甘肃省森林生态与冻土水文水资源重点实验室, 甘
肃张掖 734000; 3甘肃张掖生态科学研究院甘肃省祁连山生态科技创新服务平台, 甘肃张掖 734000; 4庆阳市林业局, 甘肃庆
阳 745000; 5 甘肃农业大学资源环境学院, 兰州 730070)
摘 要 以动态监测样地(340 m×300 m)的网格(5 m×5 m)为基本单元,选择 5 个群落结构
属性变量,采用分形几何学和地统计学方法对祁连山大野口流域青海云杉林群落结构的空间
异质性进行研究.结果表明: 青海云杉林 5个结构属性的变异程度大小依次为:密度>平均冠
幅>显著度>盖度>平均树高,变异系数为 43.7%~79.6%.Moran I系数表明,各个结构属性均有
不同程度的空间自相关性,自相关大小顺序为:密度>平均树高>盖度>平均冠幅>显著度,变
化范围为-0.047~0.382.指数理论变异函数模型能很好地拟合不同结构属性的空间变异,变
程为24.6~68.1 m,各属性变量的结构比除盖度表现为中等空间相关,其他指标均为强烈空间
相关,各属性指标分维数接近 2,空间依赖性较小.植被密度和盖度空间分布呈带状结构和斑
块结构叠加的特点,其他指标呈较强的斑块状空间结构,密度和盖度对平均冠幅、显著度和平
均树高有一定的空间依赖性.青海云杉林群落结构空间异质性适宜的采样间隔和采样面积分
别为 10 m和 0.5 hm2 .
关键词 青海云杉林; 群落结构; 空间异质性; 地统计学; 祁连山
∗国家自然科学基金项目(41461004)、甘肃省基础研究创新群体项目(145RJIG337)和甘肃省科技创新服务平台项目(144JTCG254)资助.
∗∗通讯作者. E⁃mail: liuxiande666@ 163.com
2014⁃11⁃11收稿,2015⁃05⁃05接受.
文章编号 1001-9332(2015)09-2591-09 中图分类号 Q948.1 文献标识码 A
Spatial heterogeneity of community structure of Picea crassifolia forest in Qilian Mountains,
China. ZHAO Wei⁃jun1,2,3, LIU Xian⁃de1,2,3, JING Wen⁃mao1,2,3, XU Li⁃heng4, NIU Yun2,3, QI
Peng2,5, ZHAO Yong⁃hong1 ( 1College of Forestry, Gansu Agricultural University, Lanzhou 730070,
China; 2Gansu Province Key Laboratory of Hydrology and Water Resources of Forest Ecology and
Frozen Soil, Gansu Academy of Water Resources Conservation Forests in Qilian Mountains, Zhangye
734000, Gansu, China; 3Gansu Province Science and Technology Innovation Service Platform of
Ecology in Qilian Mountains, Zhangye Academy of Ecology Science, Zhangye 734000, Gansu, Chi⁃
na; 4Forestry Bureau of Qingyang City, Qingyang 745000, Gansu, China; 5College of Resources
and Environmental Sciences, Gansu Agricultural University, Lanzhou 730070, China) . ⁃Chin. J. Ap⁃
pl. Ecol., 2015, 26(9): 2591-2599.
Abstract: We selected the grid of 5 m × 5 m in a dynamic monitoring plot (340 m × 300 m) as
the sampling unites and chose 5 structural characteristics ( density, average crown breadth, cove⁃
rage, conspicuousness and average height) to study the spatial heterogeneity of community structure
of Picea crassifolia forest in Dayekou Basin of Qilian Mountains by the fractal geometry and geosta⁃
tistics methods. The results showed that the order of spatial variation in these characteristics was:
density > average crown breadth > conspicuousness > coverage > average height, with the variation
coefficient ranging from 43.7% to 79.6%. Moran’s I index indicated that the structural variables had
different degrees of spatial autocorrelation, and the order of autocorrelation was density > average
height> coverage > average crown breadth > conspicuousness, with the range of -0.047-0.382. The
exponential semivariation model well fitted the spatial variability in different structural features, and
the range was 24.6-68.1 m. The variables displayed moderate spatial autocorrelation except for co⁃
verage, while the other variables had strong spatial autocorrelation, and the fractal dimension of the
variables was close to 2, indicating a low spatial dependence among variables. The variables presen⁃
应 用 生 态 学 报 2015年 9月 第 26卷 第 9期
Chinese Journal of Applied Ecology, Sep. 2015, 26(9): 2591-2599
ted a superposing characteristic of zonal and patchy structures except for density and coverage,
while the other variables presented strong patchiness property. Density and coverage had a certain
spatial dependence on average crown breadth, conspicuousness and average height. Density and
coverage for the spatial heterogeneity of community structural of P. crassifolia forests were 10 m and
0.5 hm2, respectively.
Key words: Picea crassifolia forest; community structure; spatial heterogeneity; geo⁃statistics;
Qilian Mountains.
森林群落是一个具有镶嵌性质的异质体[1] .群
落通过与环境因素相互作用,导致群落结构分布的
空间异质性,特别是在特殊的环境条件下,群落内的
关键种或优势种的空间分布反映了群落为维持自我
稳定而做出对环境干扰的生态格局和过程,对维持
群落稳定起着关键作用[2] .然而,由于群落结构的分
布格局对尺度依赖性,使得其在空间上出现新的变
化和特点[3] .如何将所有尺度上的群落结构和其可
能发生的生态学特征联系起来,对其结构及维持机
制有更加深入和完整的认识,分形几何学和地统计
学为解决这一问题提供了切实可行的工具.因为它
注重了变量因子的空间过程,考虑了其空间分布特
征和空间自相关[4],在国内林业科学研究领域得到
了广泛应用.例如,王政权等[5]对阔叶红松(Pinus
koraiensis)林的研究得出土壤水分具有明显的空间
异质性的结论,为土壤空间异质性与根系结构和生
长提供了参考;刘少冲等[6]对小兴安岭阔叶红松混
交林一定面积林隙土壤养分的空间异质性进行了研
究,为林木更新奠定了基础;彭晚霞等[7]对喀斯特
常绿落叶阔叶混交林植被的空间异质性研究,得出
适宜的采样间隔为 10 m 左右,采样范围大于 10
hm2,对该区域森林长期动态监测样地的建立具有
指导意义.
祁连山区作为青藏高原向北东方向扩展生长的
前缘地区之一,其北、西两面受腾格里、巴丹吉林和
塔克拉玛干沙漠的影响和高差悬殊及距海洋远近的
不同,受多种因素的影响,总体气候属大陆性高寒半
干旱、半湿润气候[8] .分布在其阴坡、半阴坡的关键
种或优势种青海云杉林受山地特殊环境因素、气候、
地形、土壤、干扰和生物学特性的综合影响,森林生
态系统环境较为脆弱,生态系统的稳定性和抗干扰
性差[9-10] .由于人类活动的长期干扰,该区生态系统
退化严重,服务功能降低,逆向演化显著,进程加
快[11],给区域生态环境建设和经济发展带来了负面
影响.本文利用国家级祁连山森林生态站建立的青
海云杉林 10.2 hm2长期动态监测样地开展试验,从
青海云杉(Picea crassifolia)林群落 5 个空间结构属
性即密度、平均冠幅、盖度、显著度和平均树高等变
量出发,采用经典统计学、分形几何学及地统计学的
方法[12-14],定量分析青海云杉林群落结构属性的空
间异质性、分形特征和空间分布格局,了解青海云杉
林群落结构属性的空间分异规律并探讨其维持机
制,为青海云杉林生态系统的恢复与生态重建提供
科学依据.
1 研究地区与研究方法
1 1 研究区概况
研究区位于祁连山国家级自然保护区北麓中段
西水林区的大野口河流域,流域面积 68.06 km2,海
拔 2650~4600 m,平均海拔为 3330 m,低山地带坡
度为 20° ~30°,高山地带坡度在 40°左右.属大陆性
高寒半湿润山地气候,年均气温 5.4 ℃,最低月平均
气温-12.5 ℃,最高月平均气温 19.6 ℃;年降水量为
300~500 mm,多集中在 6—9 月,年蒸发量为 1488
mm[15] .海拔从低到高,植被类型依次为:山地荒漠
植被、山地草原植被、山地森林草原植被、亚高山草
甸植被、高山冰雪植被,土壤类型依次为:山地灰钙
土、山地栗钙土、山地灰褐土、亚高山灌丛草甸土、高
山寒漠土.在各类土壤中山地森林灰褐土和亚高山
灌丛草甸土是生长森林的土壤,山地灰褐土分布在
海拔 2400~3300 m 地带,是乔木林的主要分布带;
亚高山灌丛草甸土分布在海拔 3300 ~ 4000 m 亚高
山地带,是湿性灌木林的主要分布带[16] .建群种青
海云杉呈斑块状或条状分布在海拔 2400 ~ 3300 m
阴坡和半阴坡地带,与阳坡草地交错分布;祁连圆柏
(Sabina przewalskii)呈小块状分布于阳坡、半阳坡;
灌木优势种有:金露梅(Potentilla fruticosa)、箭叶锦
鸡儿(Caragana jubata)、吉拉柳( Salix gilashanica)
等;草本植物主要有:珠芽蓼(Polygonum vivipamm)、
黑穗苔(Carex atrata)和针茅(Stipa)等[17] .
1 2 样地概况
动态监测样地位于祁连山西水林区的大野口河
流域[18-19],2010—2011年设立了一块面积为 340 m×
300 m的青海云杉林动态监测样地(图1) ,样地位
2952 应 用 生 态 学 报 26卷
图 1 祁连山大野口流域植被类型分布及青海云杉林样地
位置
Fig. 1 Distribution of vegetation types of Dayekou basin in
Qilian Mountains and location of Picea crassifolia forest plots.
A:针叶林 Coniferous forest; B:水体Water; C:稀疏植被 Sparse vege⁃
tation; D: 草地 Grassland; E: 高山草甸 Alpine meadow; F: 灌丛
Scrub; G: 裸土 Bare soil; H: 裸岩 Bare rock; I: 样地 Plot.
于阴坡青海云杉林相对平坦的林地中(38°31′48″ N,
100°15′00″ E),平均海拔 2835 m,坡度 13°,坡向
WN.林分起源为天然次生林,土壤为山地森林灰褐
土,年均降水量 368 mm,年均气温 0.5 ℃,平均相对
湿度 60%(1994—2004 年),水热条件较好,空气湿
度较大[20],样地内乔木树种组成单一,只有青海云
杉,DBH≥1 cm 的木本植物活个体数有 19927 株,
群落成熟稳定,林木更新良好,郁闭度为 0.67,林下
光照不足.灌木层生长发育不良,种类较少,主要由
蔷薇科金露梅和银露梅(Potentilla glabra)等落叶灌
木组成.草本层主要有披针苔草(Carex lanceolata)、
珠芽蓼、马先蒿(Pedicularis reaupinanta)等典型高山
草甸植物,物种较为丰富.苔藓层发育良好,盖度达
到 80%以上,厚度约 10 cm,主要由山羽藓(Abietinel⁃
la abietina)、欧灰藓(Hypnum cupressiforme)、提灯藓
(Mnium cuspidatum)等组成[21] .
1 3 样地的设置与调查
利用全站仪将整个样地划分成 255 个 20 m ×
20 m的样方,以样地东南角为原点,东西向为纵轴
(y),南北向为横轴(x),每隔 20 m 节点埋设石桩,
石桩上标有编号,编号前两位为横轴坐标,后两位为
纵轴坐标,再将每个 20 m×20 m样方分隔为 16 个 5
m×5 m的样格,然后按“N”字形调查,悬挂铝制标牌
进行个体标记,以待再次复查[21] .
调查对象为样方内所有胸径(DBH)≥1 cm 的
活立木木本植物,DBH≥5 cm 的个体用胸径尺测量
胸径,1 cm≤DBH<5 cm 的个体用数显游标卡尺测
量.调查内容包括:样方编号、样格编号、三维坐标以
及样地内物种组成、生长状况等,活立木木本植物的
平均胸径、平均树高、平均冠幅、平均枝下高、个体相
对坐标(x,y),记录 DBH<1.0 cm的活立木木本植物
株数和高度,及灌木层、草本层和苔藓层的物种组成
和生长状况等.
1 4 区域化变量设置及数据分析
将每个 20 m × 20 m 的样方用插值法细分成
1020(255×4)个 10 m×10 m 的小样方,统计青海云
杉的密度(株·m-2)、平均冠幅(m2)、盖度(%)、显
著度(%,即所有林木的胸高断面积之和占样方面积
的比例)和平均树高(m)等 5 个群落结构属性指
标[7] .经典统计学分析采用 SPSS 17.0 软件,半变异
函数数据分析、模型优化模拟、自相关Moran I指数和
分维数 D 的计算均在 GS+软件中完成,采用 ArcGIS
9.2软件绘制 Kriging等值线图,Moran I指数和半方差
函数及分维数的分析方法参见文献[14,22-24].
2 结果与分析
2 1 群落结构属性的描述性统计分析
从表 1 可以看出,青海云杉林 5 个植被结构属
性的变异程度依次为:密度( 79 .6%) >平均冠幅
表 1 群落结构属性的描述性统计分析
Table 1 Descriptive statistics analysis of community structure properties (n=1020)
项目
Item
最小值
Minimum
value
最大值
Maximum
value
平均值
Mean
value
标准差
SD
变异系数
CV
(%)
偏度
Skewness
峰度
Peak
value
K⁃S值
K⁃S
value
密度 Density ( ind·m-2) 0.000 1.510 0.195 0.155 79.6 2.513 10.706 4.417
密度 Density#(ind·m-2) 0.000 3.521 1.523 0.544 35.7 0.026 0.569 1.196∗
平均冠幅 Average crown breadth (m2) 0.000 53.848 8.400 5.742 68.3 2.261 9.047 4.160
平均冠幅 Average crown breadth#(m2) 0.000 4.697 2.177 0.738 33.9 0.003 0.222 0.125∗
盖度 Coverage (%) 0.000 3.778 1.194 0.539 45.1 0.592 1.472 1.499
盖度 Coverage#(%) 0.000 3.550 1.498 0.525 35.1 0.066 0.858 0.893∗
显著度 Conspicuousness (%) 0.000 0.890 0.290 0.132 45.5 0.002 0.006 1.237∗
平均高度 Average height (m) 0.000 23.400 8.448 3.689 43.7 0.714 1.034 2.060
平均高度 Average height#(m) 0.000 7.601 4.155 1.158 27.9 0.015 0.211 1.383∗
# 数据经过 Box⁃Cox转换 Data after Box⁃Cox conversion. ∗ 正态分布 Normal distribution.
39529期 赵维俊等: 祁连山青海云杉林群落结构的空间异质性
(68 3%)>显著度(45.5%) >盖度(45.1%) >平均高
度(43.7%),变异系数为 43.7% ~ 79.6%,达到强变
异程度[25],波动较为剧烈,表明青海云杉林更新良
好,结构和稳定性较低,可塑性和成长空间较大.这
与金铭等[26]的研究结论一致.总体来看,不同结构
属性均存在一定的空间异质性,其原因可能主要是
森林植被对许多基本生态学过程和物理、化学过程
在时间和空间连续系统上响应和反馈的长期作用的
结果[27] .因此,需要进一步用地统计学分析其空间
异质性.通过对偏度、峰度的观察和单样本 K⁃S 法进
行非参数检验,在 5%的检验水平下,只有显著度数
据的偏度和峰度均偏离数值 0较近,接近正态分布,
采用地统计学原理插值时,为避免比例效应,其他指
标数据均经过 Cox⁃Box 转换,转换后的数据经 K⁃S
检验服从正态分布或近似正态分布(P<0.05).
2 2 群落结构属性的空间自相关分析
从图 2 可以看出,青海云杉林 5 个植被结构指
标均有一定程度的空间依赖性,空间自相关的大小
依次为:密度 ( 0. 382) >平均树高 ( 0. 253) >盖度
(0 239)>平均冠幅(0.174) >显著度(0.138),Moran
I系数为-0.047~0.382,随着滞后距的增加,Moran I
系数呈快速下降的趋势,到距离 50 m 处,出现明显
的拐点,随着滞后距的增加,Moran I 系数减小或者
增加,但是变化幅度不大,基本在横轴附近呈波动性
变化,表明青海云杉林植被空间分布具有较好的整
体性,植被结构破碎度很小,种群维持稳定的能力较
图 2 群落结构属性的自相关分析和半变异函数
Fig.2 Spatial correlation analysis and variograms of community structure properties.
a) 密度 Density ( ind·m-2); b) 平均冠幅 Average crown breadth (m2); c) 盖度 Coverage (%); d) 显著度 Conspicuousness (%); e) 平均树高
Average height (m). 下同 The same below.
4952 应 用 生 态 学 报 26卷
表 2 群落结构属性半变异函数理论模型的相关参数
Table 2 Parameters of the best variogram models for community structure properties (n=1020)
项目
Item
模型类型
Model type
块金值
Nugget
(C0)
基台值
Sill
(C0+C)
结构比
C / (C0+C)
变程
Range
(m)
R2 分维数
Fractal
D
密度 Density∗(ind·m-2) 指数 0.038 0.305 0.887 53.1 0.842 1.923
平均冠幅 Average crown breadth∗(m2) 指数 0.059 0.549 0.893 25.8 0.661 1.970
盖度 Coverage∗(%) 指数 0.136 0.294 0.537 68.1 0.681 1.947
显著度 Conspicuousness (%) 指数 0.002 0.018 0.880 24.6 0.664 1.972
平均树高 Average height∗(m) 指数 0.137 1.351 0.899 27.6 0.614 1.970
∗ 数据经过 Cox⁃Box转换后达正态分布或近似正态分布 Data obey normal distribution or normal distribution after Cox⁃Box conversion.
强.同时,Moran I系数在滞后距 50~200 m处在横轴
上下微小浮动,表明自相关范围均在采样范围内.
2 3 群落结构属性的空间变异特征
青海云杉林 5个植被结构属性指标实验半方差
函数及理论拟合模型拟合结果显示,各结构属性指
标的拟合模型均为指数理论变异函数模型(图 2,表
2),而且决定系数 R2为 0.61~0.84,模型拟合效果较
好,表明理论模型能很好地反映植被的空间结构特
征.由表 2可以看出,植被结构各属性变量的结构比
为 0.537~0.899,比较高,说明本研究采样的最小步
长即采样密度相对合理.以 Cambardella 等[28]对空
间结构比率的研究作为判断,盖度的空间结构比率
为 25.0%~75.0%,属中等空间相关,表明在该研究
尺度上随机变异占的比率较高,可能存在更小尺度
上的控制过程,也可能是最大有效距离不够,其变程
最大,为 68.1 m.密度、平均冠幅、显著度和平均树高
的空间结构比率均大于 75.0%,为强烈的空间自相
关,由空间自相关部分引起的空间异质性程度较大,
其变异主要由林地小地形、土壤理化性质、地被物苔
藓枯落物层等结构性变异所致,变程的变幅较大,为
24.6~53.1 m.其中,显著度的变程较小,为 24.6 m;
平均冠幅和平均树高的变幅也较小,变程大小基本
相同(25.8 和 27.6 m).这 3 个指标的空间连续性范
围较小,其分布和形成过程对青海云杉林生态系统
其他结构属性起着决定性作用,尤其是显著度.密度
的变程为 53.1 m,所有林木结构属性变程均在采样
范围之内,采样范围包括了样地地形、土壤、林下地
被物等生态因子的综合控制尺度.另外,半变异函数
模型的变程和 Moran I 系数自相关范围结果基本一
致.分形维数可以对不同变量间空间自相关的强度
进行比较,其值越接近 2,表明小尺度上的变异越显
著[29-30] .青海云杉林木结构属性空间分布的分形分
析结果显示(表 2),林木结构属性分形维数 D 值接
近 2,分别为 1.923、1.970、1.947、1.972 和 1.970,表
明在此尺度范围内植被结构空间依赖性较小,均一
性程度较好.
2 4 群落结构属性的空间分布格局
从图3可以看出,5个结构属性空间分布均表
图 3 群落结构属性的空间分布
Fig.3 Spatial distribution of vegetation structure properties.
59529期 赵维俊等: 祁连山青海云杉林群落结构的空间异质性
现为围绕多个中心呈聚集分布状态,但其复杂程度
以及斑块等级不同,通过分形维数可以确定斑块的
复杂程度,密度和盖度的分形维数相对较小,格局变
异的空间依赖性较强,具有明显的带状异向性空间
结构的轮廓,又凸显出斑块状的特点,总体呈现出带
状结构和斑块结构叠加的特点,而平均冠幅、显著度
和平均树高的分形维数相对较大,且分形维数值接
近,其格局是由于林木对水分、养分等资源在小尺度
上的竞争以及中小尺度上的空间互补等因素引起
的,在图中表现为明显的斑块状空间结构,斑块呈多
而复杂的特点,盖度和密度对平均高度、显著度和平
均冠幅有一定的空间依赖性.
3 讨 论
3 1 青海云杉林群落结构属性的空间变异和分布
目前对植被空间结构属性的研究多采用植物的
物种、冠幅、盖度、高度、生物量等单个指标或多个指
标作为空间结构属性[25,31-32],也有用 Simpson 生态
优势度和 Herlbert多样性指数作为植被结构的空间
变量[27],还有用物种多样性或去势对应分析(DCA)
排序值等综合参数作为区域化变量描述植被空间格
局特征[24,33] .本研究对象为青海云杉纯林,主要选择
了能够表征青海云杉林群落结构的 5 个重要参数:
密度、平均冠幅、盖度、显著度和平均树高作为描述
植被分布的空间变量,分析群落结构属性的空间变
异与分布.描述性统计分析表明,青海云杉林 5 个结
构属性变量整体变异较大,变异系数为 43. 7% ~
79 6%,除了与自身的生物学和生态学特性有关外,
也与群落环境条件密切相关.张立杰等[34]研究表
明,青海云杉种群分布格局与海拔、含水量、坡向和
坡度等因子有关.本研究还发现,青海云杉林幼苗和
小树多分布在坡度较大的小生境,而中树和大树多
分布在坡度相对平缓的坡面.
青海云杉林 5个结构属性指标均有不同程度的
空间相关性且规律性较强,均可用指数函数模型进
行拟合,能很好地反映群落结构的空间结构特征,属
中等或强烈的自相关,表明 10 m的采样间隔相对合
理,仅盖度指标表现为中等程度的空间自相关,这是
由结构性和随机性因素共同决定的,表明在当前观
测尺度上,可能与干扰、试验误差有关.而其他指标
表现为强烈的空间自相关,说明随机因素对这些指
标空间分布贡献较小,其变异与青海云杉林生长环
境有很大的关系,即结构性变异引起的.半变异函数
模型和 Moran I 系数自相关分析基本一致,密度的
变程最大,为 68.1 m,其他指标的变程为 24.6 ~ 53.1
m.另外,分形维数接近 2,指示着大部分空间变异发
生在较小的尺度上.
3 2 青海云杉林群落结构属性与维持机制
植被的分布格局与维持机制一直是生态学的研
究热点[35] .本研究中,青海云杉林不同结构属性空
间分布格局图分析表明,青海云杉种群的平均冠幅、
显著度和平均树高的空间分布和形成过程决定了其
密度和盖度的空间分布和形成过程,进而影响了群
落结构的空间格局.冠幅、显著度和树高的空间异质
性是青海云杉林对提高光合效能以及适应环境而做
的一种生态对策,包括林木生长对几何空间的差异
性利用、基质的异质性及与大气界面的微气候相互
作用等.青海云杉通过聚集分布格局对这种异质性
做出反应[15,36],一方面青海云杉种子大部分分散在
母树周围,不同树高下的冠层相互重叠,遮阴条件使
得耐阴、喜湿的幼苗、幼树聚集在一起提高了抵抗
力,随着林木幼体的生长,对光的需求增加,通过竞
争在林隙中生长并逐渐进入群落上层.另一方面,植
被异质性与局部地形、气候、土壤基质异质性等为云
杉更新生长创造了条件.青海云杉具有生长慢、寿命
长、耐阴性强以及种子扩散能力弱等特点,能够维持
自我稳定.这与大多研究认为顶极群落物种的自我
维持机制一致[37-38] .青海云杉稳定从根本上讲依赖
于个体的生活史特征[39],从大范围来看,青海云杉
分布在高海拔地区,适宜的热量和水分条件相结合,
加上地形、地貌、小生境等特点,为青海云杉生长创
造了相对开阔的群落环境[40-41] .
有研究表明,干扰是群落结构形成的主要原
因[42-43],而林隙干扰是维持森林群落结构的必要成
分[44-46],林隙通过干扰对林内光照、温湿度、土壤和
微地形产生影响,增强了干扰生境的异质性.林隙干
扰状况对青海云杉幼苗的持续更新和建立创造了条
件,进而影响林木结构的空间异质性[47] .青海云杉
群落的生活史、生长环境及林隙干扰,加上青海云杉
的个体死亡及降雪、间歇性大风等自然扰动,可能是
青海云杉林在祁连山能够形成稳定群落的重要机
制,在一定程度上也支持了生态位理论,即生态学上
相同物种不能共存,如要共存,则必须有生态位的分
化[48] .
3 3 青海云杉林群落结构属性空间分布格局对取
样设计的启示
物种的空间分布格局不仅依赖于潜在的生态机
理,还受不同研究方法和研究尺度的影响,其中研究
6952 应 用 生 态 学 报 26卷
尺度对物种空间分布格局有着重要的影响[49],林地
环境因子诸如地形、光照条件、土壤养分、干扰等在
不同尺度上具有明显的异质性[50] .因此,选择合理
的采样间隔和采样范围可以更清楚地判断物种的空
间分布格局,还有利于更准确地探讨植被格局形成
的机制或生态学过程.本研究中,通过自相关分析、
半方差模型的变程和 Kriging 等值线图综合分析得
出,青海云杉植被结构空间分布格局研究的合适采
样间隔和采样范围为 10 m 和 0.5 hm2 .本研究也存
在一定局限性:一是研究区青海云杉林分布在海拔
2400~3300 m,而本研究的模拟取样仅仅基于一个
海拔点(2835 m 左右)进行的,随着海拔梯度变化,
环境因子特别是地形、温度、降水等因子会发生显著
变化.二是本研究对象林木植被结构密度、平均树
高、盖度、显著度和平均冠幅都是群落结构静态参
数,与群落结构的静态参数相比,群落动态参数往往
表现出更大的变异性[25,51] .因此,在不同海拔开展青
海云杉林木结构长期动态监测相关研究,可以提高
研究结果的普适性,例如,可以沿海拔梯度建立宽为
70 m、长为若干米的样带,对青海云杉林进行长期动
态监测.
参考文献
[1] Ou Y⁃D (区余端), Su Z⁃Y (苏志尧). Spatial hetero⁃
geneity dynamics of canopy structure in a montane ever⁃
green broadleaved forest following a natural disturbance
in North Guangdong. Plant Science Journal (植物科学
学报), 2012, 30(3): 223-229 (in Chinese)
[2] Sun R⁃Y (孙儒泳), Li Q⁃F (李庆芬), Liu C⁃J (刘翠
娟), et al. Fundamentals of Ecology. Beijing: High
Education Press, 2002: 152-173 (in Chinese)
[3] Du H⁃Q (杜华强), Tang M⁃P (汤孟平), Zhou G⁃M
(周国模). Linking species diversity with spatial scale
dependence to spatial patterns using multifractal analy⁃
sis. Acta Ecologica Sinica (生态学报), 2007, 27
(12): 5038-5049 (in Chinese)
[4] Wang H (王 蕙), Zhao W⁃Z (赵文智), Chang X⁃X
(常学向). Spatial variability of soil moisture and vege⁃
tation in desert⁃oasis ecotone in the middle reaches of
Heihe River Basin. Acta Ecologica Sinica(生态学报),
2007, 27(5): 1731-1739 (in Chinese)
[5] Wang Z⁃Q (王政权), Wang Q⁃C (王庆成). The spa⁃
tial heterogeneity of soil physical properties in forests.
Acta Ecologica Sinica (生态学报), 2000, 20 ( 6):
945-950 (in Chinese)
[6] Liu S⁃C (刘少冲), Duan W⁃B (段文标). Spatial he⁃
terogeneity of soil nutrition in the mixed broad⁃leaved
Korean Pine forest gap. Journal of Soil and Water Con⁃
servation (水土保持学报), 2011, 25(3): 142- 146
(in Chinese)
[7] Peng W⁃X (彭晚霞), Song T⁃Q (宋同清), Zeng F⁃P
(曾馥平), et al. Spatial heterogeneity of vegetation in
Karst mixed forest of evergreen and deciduous broadleaf.
Acta Botanica Borcali⁃Occidentalia Sincia (西北植物学
报), 2011, 31(4): 815-822 (in Chinese)
[8] Wang J⁃Y (王金叶), Chang X⁃X (常学向), Ge S⁃L
(葛双兰), et al. Vertical distribution of the vegetation
and water and heat conditions of Qilian Mountain north⁃
ern slope. Journal of Northwest Forestry University (西北
林学院学报), 2001, 16(suppl.): 1-3 (in Chinese)
[9] Peng S⁃Z (彭守璋), Zhao C⁃Y (赵传燕), Zheng X⁃L
(郑祥霖), et al. Spatial distribution characteristics of
the biomass and carbon storage of Qinghai spruce (Picea
crassifolia) forests in Qilian Mountains. Chinese Journal
of Applied Ecology (应用生态学报), 2011, 22(7):
1689-1694 (in Chinese)
[10] Yang G⁃J (杨国靖), Xiao D⁃N (肖笃宁). Spatial pat⁃
tern analysis of forest landscape in low coteau of Middle
Qilian Mountains. Chinese Journal of Applied Ecology
(应用生态学报), 2004, 15(2): 269-272 ( in Chi⁃
nese)
[11] Liu P⁃X (刘普辛), Sun J (孙 婧), Chen F⁃H (陈
发虎), et al. Changes of landscape pattern in forest
ecosystem in Sunan mountainous region at upper Heihe
Stem River. Journal of Desert Research (中国沙漠),
2008, 28(1): 89-93 (in Chinese)
[12] Zhang W⁃T (张文彤). SPSS 11 Statistical Analysis Tu⁃
torial. Beijing: Hope Electronic Press, 2002 ( in Chi⁃
nese)
[13] Zu Y⁃G (祖元刚), Ma K⁃M (马克明), Zhang X⁃J
(张喜军). A fractal method for analysing spatial hetero⁃
geneity of vegetation. Acta Ecologica Sinica (生态学
报), 1997, 17(3): 333-337 (in Chinese)
[14] Wang Z⁃Q(王政权). Geostatistics and Its Application
in Ecology. Beijing: Science Press, 1999: 150-156 (in
Chinese)
[15] Li X⁃X (李效雄), Liu X⁃D (刘贤德), Zhao W⁃J (赵
维俊). Population structure and spatial distribution pat⁃
tern of Picea crassifolia in Dayekou Basin of Qilian
Mountains. Arid Land Geography (干旱区地理 ),
2012, 35(6): 960-967 (in Chinese)
[16] Wang J⁃Y (王金叶), Wang Y⁃L (王艺林), Jin B⁃W
(金博文), et al. Studies on regulating function of forest
hydrology and microclimate in arid and semi⁃arid area of
west China. Scientia Silvae Sinicae (林业科学), 2011,
37(5): 120-125 (in Chinese)
[17] Niu Y (牛 赟), Liu X⁃D (刘贤德), Wang L (王
79529期 赵维俊等: 祁连山青海云杉林群落结构的空间异质性
立), et al. Feature analysis on stand structure of Picea
crassifolia and its water and temperature of soil in
Dayekou Basin of Qilian Mountains. Ecology and Envi⁃
ronmental Sciences (生态环境学报), 2014, 23(3):
385-391 (in Chinese)
[18] State Forestry Bureau (国家林业局). Observational
Methods of Standard System for Long⁃Term Located Ob⁃
servation of Forest Ecosystems: Forestry Industry Stan⁃
dard of the People’ s Republic of China LY / T 1952 -
2011. Nanjing: Phoenix Press, 2011
[19] Condit R, Ashton PS, Baker P, et al. Spatial patterns
in the distribution of tropical tree species. Science,
2000, 288: 1414-1418
[20] Jiang C (蒋 冲), Wang F (王 飞), Mu X⁃M (穆
兴民), et al. Characteristics of temperature and humidi⁃
ty in forest ecosystem of Heihe River basin. Bulletin of
Soil and Water Conservation (水土保持通报), 2012,
32(3): 96-100 (in Chinese)
[21] Li X⁃X (李效雄), Liu X⁃D (刘贤德), Zhao W⁃J (赵
维俊). Community structure of a dynamical plot of Picea
crassifolia forest in Qilian Mountains,China. Journal of
Desert Research (中国沙漠), 2013, 33(1): 94-100
(in Chinese)
[22] Goovaets P. Geostatistical tools for characterizing the
spatial variability of microbiological and physico⁃chemi⁃
cal soil properties. Biology and Fertility of Soils, 1998,
27: 314-315
[23] Zhu Y (朱 岩), Zhang Y (张 岩), Zhu Q⁃K (朱
清科), et al. Spatial distribution patterns of vegetations
on natural grassland in Loess Plateau interfluve. Chinese
Journal of Ecology (生态学杂志), 2011, 30 ( 6):
1063-1069 (in Chinese)
[24] Shang Z⁃H (尚占环), Yao A⁃X (姚爱兴), Long R⁃J
(龙瑞军). Spatial heterogeneity of the north slope’s
plant communities in the mountain desert grassland of
the arid region. Acta Ecologica Sinica (生态学报),
2005, 25(2): 312-318 (in Chinese)
[25] Yan L (颜 亮), Zhou G⁃S (周广胜), Zhang F (张
峰), et al. Spatial heterogeneity of vegetation coverage
and its temporal dynamics in desert steppe, Inner Mon⁃
golia. Acta Ecologica Sinica (生态学报), 2012, 32
(13): 4017-4024 (in Chinese)
[26] Jin M (金 铭), Li Y (李 毅), Liu X⁃D (刘贤
德), et al. Spatial structure characteristic of Picea cras⁃
sifolia in Qilian Mountains. Arid Land Geography (干旱
区地理), 2012, 35(4): 587-593 (in Chinese)
[27] Qian Y⁃B (钱亦兵), Wu Z⁃N (吴兆宁), Yang H⁃F
(杨海峰), et al. Vegetation spatial heterogeneity across
longitudinal dunes in the southern Gurbantunggut Desert.
Journal of Desert Research (中国沙漠), 2011, 31(2):
420-427 (in Chinese)
[28] Cambardella CA, Moorman TB, Novak JM, et al. Field⁃
scale variability of soil properties in central Iowa soils.
Soil Science Society of America Journal, 1994, 58:
1501-1511
[29] Shao F⁃L (邵方丽), Yu X⁃X (余新晓), Yang Z⁃J
(杨志坚), et al. Spatial heterogeneity of topsoil mois⁃
ture and litter lay in natural Populus davidiana-Betula
platyphylla secondary forests. Journal of Soil and Water
Conservation (水土保持学报), 2012, 26(3): 199-
204 (in Chinese)
[30] Burrough PA. Multiscale sources of spatial variation in
soil. Ⅰ. The application of fractal concepts to nested
levels of soil variation. Journal of Soil Science, 1983,
34: 577-597
[31] He Z⁃B (何志斌), Zhao W⁃Z (赵文智). Characteris⁃
tics of spatial heterogeneity of the desert riparian forest
in the Hei River Basin. Journal of Glaciology and Geco⁃
cryology (冰川冻土), 2003, 25(5): 591 - 596 ( in
Chinese)
[32] Augustine DJ. Spatial heterogeneity in the herbaceous
layer of a semi⁃arid savanna ecosystem. Plant Ecology,
2003, 167: 319-332
[33] Yang X⁃L (杨小林), Zhao K⁃T (赵垦田), Ma H⁃P
(马和平), et al. Ecological studies on vegetation quan⁃
tity in the semi⁃arid valley region of Lasa. Scientia Silvae
Sinicae (林业科学), 2010, 46(10): 17-24 ( in Chi⁃
nese)
[34] Zhang L⁃J (张立杰), Zhao W⁃Z (赵文智), He Z⁃B
(何志斌). Characteristics in Picea crassifolia forest
fractal dimension and its influencing factors. Acta Eco⁃
logica Sinica (生态学报), 2008, 28(4): 1383-1389
(in Chinese)
[35] Agustin R, Adrian E. Small⁃scale spatial soil⁃plant rela⁃
tionship in semi⁃arid gypsum environments. Plant and
Soil, 2000, 220: 139-150
[36] Zhao L⁃Q (赵丽琼), Huang H⁃G (黄华国), Liang D⁃
S (梁大双), et al. Spatial distribution pattern of Picea
crassifolia population in Dayekou, Gansu Province. Jour⁃
nal of Beijing Forestry University (北京林业大学学
报), 2010, 32(4): 59-64 (in Chinese)
[37] Huston M, Smith T. Plant succession: Life history and
competition. American Naturalist, 1987, 130: 168-198
[38] Guo Y⁃X (郭垚鑫). Characteristics and Stability of
Batula albosinensis Communities in Qinling Mountains.
Master Thesis. Yangling: Northwest A&F University,
2013 (in Chinese)
[39] Laurence DM, Joshi A, Borash DJ. Does population sta⁃
bility evolve? Ecology, 2000, 8: 1273-1285
[40] Liu X⁃C (刘兴聪). Picea crassifolia. Lanzhou: Lanzhou
8952 应 用 生 态 学 报 26卷
University Press, 1992 (in Chinese)
[41] Peng S⁃Z (彭守璋), Zhao C⁃Y (赵传燕), Xu Z⁃L
(许仲林 ), et al. Potential distribution of Qinghai
spruce and assessment of its growth status in the upper
reaches of the Heihe River in the Qilian Mountains of
China. Chinese Journal of Plant Ecology (植物生态学
报), 2011, 35(6): 605-614 (in Chinese)
[42] Denslow JS, Spies T. Canopy gaps in forest ecosystems:
An introduction. Canadian Journal of Forest Research,
1990, 205: 642-648
[43] Lemmans I. Simulation and future protection of succes⁃
sion in a Swedish broad⁃leaved forest. Forest Ecology
and Management, 1992, 48: 305-319
[44] Zang R⁃G (臧润国), Liu J⁃Y (刘静艳), Dong D⁃F
(董大方). Gap Dynamic and Forest Biodiversity. Bei⁃
jing: China Forestry Press, 1999: 73-200 (in Chinese)
[45] Yan S⁃J (闫淑君), Hong W (洪 伟), Wu C⁃Z (吴
承祯), et al. Gaps and their natural disturbance charac⁃
teristics in mid⁃subtropical evergreen broad⁃leaved forest
in Wanmulin. Chinese Journal of Applied Ecology (应用
生态学报), 2004, 15(7): 1126-1130 (in Chinese)
[46] Luo D⁃Q (罗大庆), Guo Q⁃S (郭泉水), Xue H⁃Y
(薛会英), et al. Characteristics and disturbance status
of gaps in subalpine fir forest in Southeast Tibet. Chinese
Journal of Applied Ecology (应用生态学报), 2002, 13
(7): 777-780 (in Chinese)
[47] Song X⁃Z (宋新章), Xiao W⁃F (肖文发). Research
advances of microsites and regeneration within canopy
gap. Scientia Silvae Sinicae (林业科学), 2006, 42
(5): 114-119 (in Chinese)
[48] Zhang D⁃Y (张大勇). The Theory of Ecology. Beijing:
Higher Education Press, 2000 (in Chinese)
[49] Song H⁃J (宋厚娟), Ye J (叶 吉), Lin F (蔺
菲), et al. Effects of sampling scales on spatial patterns
of woody species in forests. Chinese Science Bulletin (科
学通报), 2014, 59(24): 2388-2395 (in Chinese)
[50] Wu H, Li ZL. Scale issues in remote sensing: A review
on analysis, processing and modeling. Sensors, 2009,
9: 1768-1793
[51] Wagner F, Rutishauser E, Blanc L, et al. Effects of plot
size and census interval on descriptors of forest structure.
Biotropica, 2010, 42: 664-671
作者简介 赵维俊,男,1981年生,博士研究生. 主要从事森
林生态和湿地生态研究. E⁃mail: zhaoweijun1019@ 126.com
责任编辑 孙 菊
99529期 赵维俊等: 祁连山青海云杉林群落结构的空间异质性