-
岷江是长江主要支流,生态地位突出,而岷江干旱河谷区域受“焚风效应”严重影响,导致辐射强、气温高、风速大、蒸发量大于降水量,土壤瘠薄、干燥,植被以灌、草丛为主,并形成稳定的灌木群落[1, 2]。灌木在岷江干旱河谷区水土保持和植被恢复中具有重要的生态地位,而生物量是灌木群落结构和功能的主要测度指标之一,研究灌木生物量对于研究灌木生长发育规律以及在生态系统中的作用和地位等都具有重要意义。现已有研究证明了该区几种主要灌丛地上生物量的分布规律[3, 4]。岷江干旱河谷区有丰富的蔷薇资源,而峨眉蔷薇(Rosa omeiensis)是优良的抗旱树种和绿化景观树种,其中分布广,适应性强,抗病性强,其根和果实有止血,止痢等功效,多用于治疗吐血,衄血,崩漏,带下病,赤白痢疾。本文在大量野外调查的基础上,研究了岷江干旱河谷区不同海拔梯度峨眉蔷薇单株地上生物量及估测模型,目的是分析峨眉蔷薇的生物量与环境因子的关系,找到该区域峨眉蔷薇最适合的单株生物量估测模型,从而为岷江干旱河谷区植被保护与恢复工作提供科学理论依据。
-
本研究在2016年9月下旬植物停止生长时进行调查,分别于岷江干旱河谷区(理县甘堡乡)两边的阴坡和阳坡由山脚至山顶沿海拔各设3个梯度,海拔间隔约为100~150 m,在阴、阳坡各海拔梯度人畜干扰少的地段选择有代表性的植被类型,设置1—2个20 ×20 m 2的标准样地,在样地内对峨眉蔷薇进行每木检尺,测定基径(D)、树(H)、冠幅(C)等指标。按不同径级,结合峨眉蔷薇分布密度,每个样地选择5—10株样木,采取直接收获法,刈割全部灌木地上部分称其总鲜重,再按干、枝、皮、叶、花(果)等不同部位,分别称鲜重并取样,将样品带回实验室在85 ℃的通风干燥箱内烘干至绝对干重,称重并计算出各部位干重及单株地上部分总干重(W干样)。同时采用GPS全球定位仪调查海拔、坡度、坡向等环境因子(见表1)。
表 1 岷江干旱河谷区样地特征表
Table 1. Characteristics of sample plots in the dry valley of Minjiang River
样地号 海拔/m 坡向 坡度 坡位 主要伴生植被 样本数/株 盖度/% 阳01 1980 SE118° 35° 中下坡 鞍叶羊蹄甲 Bauhinia brachycarpa,刺旋花 Convolvulus tragacanthoides,白刺花 Sophora davidii,小蓝雪花 Ceratostigma minus,铁杆蒿 Artemisia gmelinii 13 42 阳02 2160 SE129° 33° 中坡 鞍叶羊蹄甲,白刺花,刺旋花,铁杆蒿,黄花亚菊 Ajania nubigena,光果莸 Caryopteris tangutica 13 39 阳03 2320 SE123° 28° 上坡 鞍叶羊蹄甲,多花蔷薇 Rosa multiflora,鲜黄小檗 Berberis diaphana ,白刺花,铁杆蒿 16 48 阴01 1815 N 0° 30° 下坡 刺旋花,白刺花,甘肃矮探春 Jasminum floridum,金花蚤草 Pulicaria chrysantha,黄花亚菊,光果莸,华北驼绒藜 Krascheninnikovia arborescens,小蓝雪花 12 45 阴02 1950 NW 335° 36° 中坡 虎榛子 Ostryopsis davidiana,橿子栎 Quercus baronii,鲜黄小檗,散生栒子 Cotoneaster divaricatus 15 42 阴03 2060 NE 9° 37° 上坡 虎榛子,橿子栎,多花蔷薇,球花荚迷,散生栒子,华北驼绒藜 Viburnum betulifolium 13 52 -
估测灌木地上生物量,选择自变量十分重要。既可采用由D、H、C等单一变量,也可以采用D、D2、C、C2、C2H、D2H、HC、HD、DC等多因子相组合的复合变量。本文采用EXCEL对峨眉蔷薇的基径、树高、冠幅进行统计,采用SPSS13.0的Bivariate Correlation分析对基径、树高、冠幅与峨眉蔷薇生物量的相关性进行分析,最终选择预测模型自变量,同时建立峨眉蔷薇生物量预测回归方程,选出拟合性好,相关系数高的模型作为峨眉蔷薇生物量预测模型。可选择许多种线性、对数、指数、幂函数以及多项式等回归模型来进行灌木生物量预测模型的建立[7, 8],常用的回归方程有一元方程y=a+bx、幂函数方程y=axb、指数函数y=aebx等。
-
表2为阴阳坡不同海拔梯度样地峨眉蔷薇地上生物量分配比例表,从表中可以看出,不论阴阳坡,不同径级的峨眉蔷薇茎和枝在地上生物量(即干重)中占比最大,阳坡为0.41,阴坡为0.39,其中茎占比略大于枝,阳坡为0.35,阴坡为0.34,叶占比最小,阴阳坡均为0.11,而皮占比略大于叶,阳坡为0.13,阴坡为0.16;而从地上总生物量对鲜重的占比来看(生物量占鲜重比越大,说明植株含水率越低,木质化程度高),阳坡中,呈现出中下坡(0.67)>上坡(0.66)>中坡(0.61)的趋势,阴坡中,上坡(0.71)>下坡(0.69)>中坡(0.61)的趋势,说明不论阴阳坡,中坡位的峨眉蔷薇植株干鲜重比例都是最低的。
表 2 阴阳坡不同海拔梯度峨眉蔷薇地上生物量分配比例表
Table 2. Aboveground biomass allocation ratio of Rosa omeiensis at different altitude gradients on shady and sunny slopes
样地号 基径
区间/mm生物量
区间/g茎干重/
鲜重占总干
重比皮干重/
鲜重占总干
重比枝干重/
鲜重占总干
重比叶干重/
鲜重占总干
重比总干重/
鲜重阳01 0.61~1.45 19.22~234.51 0.67 0.32 0.76 0.15 0.68 0.39 0.56 0.13 0.67 阳02 0.87~1.49 39.72~383.14 0.63 0.53 0.66 0.10 0.63 0.24 0.52 0.12 0.61 阳03 0.51~2.9 14.90~418.55 0.68 0.36 0.71 0.13 0.64 0.42 0.55 0.08 0.66 平均 0.66 0.41 0.71 0.13 0.65 0.35 0.55 0.11 0.65 阴01 0.87~2.15 32.94~446.05 0.68 0.40 0.81 0.14 0.71 0.34 0.58 0.12 0.69 阴02 0.66~1.18 16.57~106.8 0.66 0.38 0.68 0.18 0.57 0.30 0.54 0.14 0.61 阴03 0.65~1.53 17.94~160.92 0.74 0.39 0.75 0.15 0.70 0.38 0.58 0.07 0.71 平均 0.69 0.39 0.75 0.16 0.66 0.34 0.57 0.11 0.67 -
从表3可以看出,峨眉蔷薇地上部分生物量W与D、H、C、D2H、HC、HD、DC等多因子或多因子组合的复合变量有着极显著的相关关系。在所有因子中,以D与W相关性最高(0.854),而C与W相关性最低(0.691);在单因子中,也以D与W相关性最高(0.854),H与W相关性次之(0.733),而C与W相关性相对较低(0.691);在复合因子中,以DC与W相关性最高(0.848),而HC与W相关性相对较低(0.70)。在单因子之间,D与C之间也有着较为明切的关系(0.762),而H与C和D之间相关性均不太显著。
表 3 峨眉蔷薇自变量因子和地上生物量相关关系分析表
Table 3. Correlation analysis results between independent variable factors and aboveground biomass of Rosa omeiensis
项目 W D H C HC HD DC D2H W 1 0.854** 0.733* 0.691** 0.700** 0.804** 0.848** 0.808** D 0.854** 1 0.428* 0.782** 0.762** 0.909** 0.973** 0.946** H 0.733* 0.428* 1 0.375* 0.775** 0.735** 0.388* 0.487** C 0.691** 0.782** 0.375* 1 0.854** 0.721** 0.869** 0.711** HC 0.700** 0.762** 0.775** 0.854** 1 0.897** 0.804** 0.775** HD 0.804** 0.909** 0.735** 0.721** 0.897** 1 0.886** 0.941** DC 0.848** 0.973** 0.388* 0.869** 0.804** 0.886** 1 0.946** D2H 0.808** 0.946** 0.487** 0.711** 0.775** 0.941** 0.946** 1 **代表在P<0.01时极为显著。*. 代表在P<0.05时极为显著。 通过选择与生物量密切相关的测树因子和复合因子为自变量,进行曲线拟合,筛选出估测峨眉蔷薇单株地上生物量的统计模型(见表4)。从筛选出的模型来看,预测生物量最优模型多以三次曲线和幂函数为佳,其中以单因子D和HD、DC、D2H复合因子为自变量对峨眉蔷薇地上生物量预测模型最佳,其相关系数R2均为0.95以上,又以单因子C为自变量对峨眉蔷薇地上生物量预测模型较差,其相关系数R2为0.8919。
表 4 峨眉蔷薇地上生物量模型
Table 4. Aboveground Biomass Model of Rosa omeiensis
对应关系 因变量 回归方程 相关系数 标准差 F P 样本数 D-W D y = 646.37x3 − 1 449.9x2 + 1 211.8x − 296.84 R2 = 0.9546 27.65 70.02 <0.001 30 y = 116.1169x3.211 R2 = 0.935 30.17 172.73 <0.001 30 H-W H y = 286.7x3 − 594.99x2 + 581.45x − 147.55 R2 = 0.9386 32.15 50.95 <0.001 30 C-W C y = 2 451.4x3 − 4 564.1x2 + 2 879.7x − 519.18 R2 = 0.8919 42.65 27.49 <0.001 30 HD-W DH y = 126.88482x1.4191 R2 = 0.969 20.85 375.04 <0.001 30 y = 14.773x3 − 11.985x2 + 131.2x − 8.7881 R2 = 0.9323 21.56 117.34 <0.001 30 DC-W DC y = 246.14x3 − 513.59x2 + 463.98x − 42.172 R2 = 0.9673 23.45 98.65 <0.001 30 HC-W HC y = 115.01x3 − 238.68x2 + 321.87x − 23.58 R2 = 0.9196 18.52 160.19 <0.001 30 D2H-W D2H y = 9.2257x3 − 43.052x2 + 155.91x + 2.9289 R2 = 0.9511 20.88 125.17 <0.001 30 y = 124.3118x0.9831 R2 = 0.9536 21.52 351.19 <0.001 30
Study on Aboveground Biomass and Model of Rosa omeiensi in Dry Valley of the Minjiang River
More Information-
摘要: 岷江干旱河谷区是植被恢复困难地带,研究峨眉蔷薇生物量及模型,可对岷江干旱河谷区植被保护与恢复工作提供科学理论依据。研究结果表明:①不同径级的峨眉蔷薇各器官生物量分配大小,在阴阳坡均表现为干生物量>枝生物量>皮生物量>叶生物量,占比最大的主干部分对地上生物量总量贡献较大;②峨眉蔷薇对海拔梯度的响应,在阴阳坡均表现为中坡位的峨眉蔷薇植株地上生物量占鲜重比是最低的,说明在该区域的中坡位,受土壤含水率低的影响,峨眉蔷薇植株生长受到影响,木质化程度低,在相同鲜重的情况下,地上部分生物量明显低于其他坡位;③冠幅(C)、基径(D)和树高(H)都是和灌木生物量有密切相关关系的自变量因子;而峨眉蔷薇生物量估测模型筛选结果显示,幂函数模型和三多项式模型都有令人满意的R2值,最优模型多以三次多项式为佳,因为三次多项式模型R2值更高,但由于灌木形态各异,所以应根据实际情况进行自变量因子的选择。Abstract: Vegetation restoration is very difficult in the dry valley of the Minjiang River. Research on the biomass and model of Rosa omeiensis can provide scientific theoretical basis for vegetation protection and restoration in the dry valley of the Minjiang River. The results showed that: (1) Regardless of the slopes, the biomass distribution of different organs of Rosa omeiensis with different diameters was in the order: dry biomass > branch biomass > bark biomass > leaf biomass. The main stem with the largest proportion contributed greatly to the total aboveground biomass. (2) The ratio of aboveground biomass to fresh weight of Rosa omeiensis plants in the middle slope position was the lowest, which indicated that the growth of Rosa omeiensis plants was affected by the low soil moisture content in the middle slope position in this area, and the degree of lignification was low. Under the same fresh weight, the aboveground biomass was significantly lower than other slope positions. (3) The crown width (C), basal diameter (D) and tree height (H) were all independent variables closely related to shrub biomass. The screening results of Rosa omeiensis biomass estimation model showed that both power function model and triple polynomial model had satisfactory correlation coefficient values. The optimal model was mostly the cubic polynomial, because the R2 value of cubic polynomial model was higher. Considering the different shrub morphology, the independent variable factors should be selected according to the actual conditions.
-
Key words:
- Biomass;
- Model;
- Correlation relationship;
- Independent variable factor
-
表 1 岷江干旱河谷区样地特征表
Tab. 1 Characteristics of sample plots in the dry valley of Minjiang River
样地号 海拔/m 坡向 坡度 坡位 主要伴生植被 样本数/株 盖度/% 阳01 1980 SE118° 35° 中下坡 鞍叶羊蹄甲 Bauhinia brachycarpa,刺旋花 Convolvulus tragacanthoides,白刺花 Sophora davidii,小蓝雪花 Ceratostigma minus,铁杆蒿 Artemisia gmelinii 13 42 阳02 2160 SE129° 33° 中坡 鞍叶羊蹄甲,白刺花,刺旋花,铁杆蒿,黄花亚菊 Ajania nubigena,光果莸 Caryopteris tangutica 13 39 阳03 2320 SE123° 28° 上坡 鞍叶羊蹄甲,多花蔷薇 Rosa multiflora,鲜黄小檗 Berberis diaphana ,白刺花,铁杆蒿 16 48 阴01 1815 N 0° 30° 下坡 刺旋花,白刺花,甘肃矮探春 Jasminum floridum,金花蚤草 Pulicaria chrysantha,黄花亚菊,光果莸,华北驼绒藜 Krascheninnikovia arborescens,小蓝雪花 12 45 阴02 1950 NW 335° 36° 中坡 虎榛子 Ostryopsis davidiana,橿子栎 Quercus baronii,鲜黄小檗,散生栒子 Cotoneaster divaricatus 15 42 阴03 2060 NE 9° 37° 上坡 虎榛子,橿子栎,多花蔷薇,球花荚迷,散生栒子,华北驼绒藜 Viburnum betulifolium 13 52 表 2 阴阳坡不同海拔梯度峨眉蔷薇地上生物量分配比例表
Tab. 2 Aboveground biomass allocation ratio of Rosa omeiensis at different altitude gradients on shady and sunny slopes
样地号 基径
区间/mm生物量
区间/g茎干重/
鲜重占总干
重比皮干重/
鲜重占总干
重比枝干重/
鲜重占总干
重比叶干重/
鲜重占总干
重比总干重/
鲜重阳01 0.61~1.45 19.22~234.51 0.67 0.32 0.76 0.15 0.68 0.39 0.56 0.13 0.67 阳02 0.87~1.49 39.72~383.14 0.63 0.53 0.66 0.10 0.63 0.24 0.52 0.12 0.61 阳03 0.51~2.9 14.90~418.55 0.68 0.36 0.71 0.13 0.64 0.42 0.55 0.08 0.66 平均 0.66 0.41 0.71 0.13 0.65 0.35 0.55 0.11 0.65 阴01 0.87~2.15 32.94~446.05 0.68 0.40 0.81 0.14 0.71 0.34 0.58 0.12 0.69 阴02 0.66~1.18 16.57~106.8 0.66 0.38 0.68 0.18 0.57 0.30 0.54 0.14 0.61 阴03 0.65~1.53 17.94~160.92 0.74 0.39 0.75 0.15 0.70 0.38 0.58 0.07 0.71 平均 0.69 0.39 0.75 0.16 0.66 0.34 0.57 0.11 0.67 表 3 峨眉蔷薇自变量因子和地上生物量相关关系分析表
Tab. 3 Correlation analysis results between independent variable factors and aboveground biomass of Rosa omeiensis
项目 W D H C HC HD DC D2H W 1 0.854** 0.733* 0.691** 0.700** 0.804** 0.848** 0.808** D 0.854** 1 0.428* 0.782** 0.762** 0.909** 0.973** 0.946** H 0.733* 0.428* 1 0.375* 0.775** 0.735** 0.388* 0.487** C 0.691** 0.782** 0.375* 1 0.854** 0.721** 0.869** 0.711** HC 0.700** 0.762** 0.775** 0.854** 1 0.897** 0.804** 0.775** HD 0.804** 0.909** 0.735** 0.721** 0.897** 1 0.886** 0.941** DC 0.848** 0.973** 0.388* 0.869** 0.804** 0.886** 1 0.946** D2H 0.808** 0.946** 0.487** 0.711** 0.775** 0.941** 0.946** 1 **代表在P<0.01时极为显著。*. 代表在P<0.05时极为显著。 表 4 峨眉蔷薇地上生物量模型
Tab. 4 Aboveground Biomass Model of Rosa omeiensis
对应关系 因变量 回归方程 相关系数 标准差 F P 样本数 D-W D y = 646.37x3 − 1 449.9x2 + 1 211.8x − 296.84 R2 = 0.9546 27.65 70.02 <0.001 30 y = 116.1169x3.211 R2 = 0.935 30.17 172.73 <0.001 30 H-W H y = 286.7x3 − 594.99x2 + 581.45x − 147.55 R2 = 0.9386 32.15 50.95 <0.001 30 C-W C y = 2 451.4x3 − 4 564.1x2 + 2 879.7x − 519.18 R2 = 0.8919 42.65 27.49 <0.001 30 HD-W DH y = 126.88482x1.4191 R2 = 0.969 20.85 375.04 <0.001 30 y = 14.773x3 − 11.985x2 + 131.2x − 8.7881 R2 = 0.9323 21.56 117.34 <0.001 30 DC-W DC y = 246.14x3 − 513.59x2 + 463.98x − 42.172 R2 = 0.9673 23.45 98.65 <0.001 30 HC-W HC y = 115.01x3 − 238.68x2 + 321.87x − 23.58 R2 = 0.9196 18.52 160.19 <0.001 30 D2H-W D2H y = 9.2257x3 − 43.052x2 + 155.91x + 2.9289 R2 = 0.9511 20.88 125.17 <0.001 30 y = 124.3118x0.9831 R2 = 0.9536 21.52 351.19 <0.001 30 -
[1] 刘兴良,慕长龙,向成华,等. 四川西部干旱河谷自然特征及植被恢复与重建途径[J]. 四川林业科技,2001,22(2):10−17. [2] 叶延琼,陈国阶,樊宏. 岷江上游脆弱生态环境刍论[J]. 长江流域资源与环境,2002,11(4):383−387. [3] 刘国华,张洁瑜,张育新,等. 岷江干旱河谷三种主要灌丛地上生物量的分布规律[J]. 山地学报,2003,21(1):24−32. [4] 刘国华,马克明,傅伯杰,等. 岷江干旱河谷主要灌丛类型地上生物量研究[J]. 生态学报,2003,23(9):1757−1764. [5] 朱炳海. 中国之气候[M]. 北京: 科学出版社, 1962. [6] 刘文彬. 岷江上游半干旱河谷灌木种群的主要类型[J]. 山地研究,1994,2(1):27−31. [7] 朱志诚,贾东林. 陕北黄土高原铁秆蒿群落生物量初步研究[J]. 生态学报,1993,13(3):243−251. [8] Mariessii, A., Betule, E., Toshihicoko, H. Growth patterns of tree height and stem diameter in populations of Abies veitchi[J]. Ecology, 1991, 79: 1085−1095. doi: 10.2307/2261100 [9] 黎燕琼,郑绍伟,龚固堂,等. 不同年龄柏木混交林下主要灌木黄荆生物量及分配格局[J]. 生态学报,2010,30(11):2809−2818. [10] 郭威星,蓝登明,王玉婕,等. 荒漠草原6种灌丛地上生物量分析[J]. 内蒙古林业调查设计,2019,42(6):87−89, 93. [11] 严代碧,岳永杰,郑绍伟,等. 岷江上游干旱河谷区土壤水分含量及其动态[J]. 南京林业大学学报(自然科学版),2006,30(4):64−68. [12] 曾慧卿,刘琪璟,冯宗炜,等. 红壤丘陵区林下灌木生物量估算模型的建立及其应用[J]. 应用生态学报,2007,18(10):2185−2190. [13] 张峰,上官铁梁,李素珍. 关于灌木生物量建模方法的改进[J]. 生态学报,1993,12(6):67−69.