论文部分内容阅读
随着全球气候变暖成为当今世界各国政府和科学界关心的热点问题,对森林碳计量的研究也备受关注。森林碳计量是指在一定的时限内和给定的地域内,对不同森林类型碳贮存量与碳流通量进行估算,其中森林碳储量及碳汇能力是碳计量研究中的关键点。中国是世界上最大的发展中国家,同时也是温室气体排放大国,气候变化对我国生态环境和社会经济造成的负面影响也日益突出。因此,积极主动地参与林业碳计量并进行相关研究对我国的生态、经济、社会的发展具有十分重要的意义。杉木(Cunninghamia lanceolata)是我国重要的用材树种,栽培历史悠久,不仅为我国经济发展提供了大量的商品用材林,而且在调节碳平衡、减缓大气中二氧化碳(CO2)等温室气体浓度上升方面也发挥着重要的作用。迄今为止,我国已经积累了大量有关杉木的基础研究数据,但是由于数据来源和研究方法的不同,对杉木碳计量的研究结果还存在较大差异。众所周知,森林碳计量的测定非常困难,而且耗时费力,确定一种行之有效而又准确的调查预估方法是十分重要的。因此利用已有的广泛可靠的杉木碳储量(生物量)研究数据,尝试建立精度高,适用广的碳计量模型或方法能大大减少测定立木碳含量和林分碳储量的外业工作。只要通过测定一定数量样木、样地碳储量(生物量)的数据,建立模型后就可以在同类的林分中结合相应的林分调查数据来估计单株碳含量或是整体林分碳储量。不仅如此,在大尺度的研究上还能充分利用这些模型结合已经有的森林资源清查体系,通过尺度换算外推对区域内杉木林总碳储量进行估测,这也是研究我国杉木整体碳汇功能及其动态变化的有效途径。鉴于此,本研究通过收集大量不同来源的杉木生物量(碳储量)相关研究数据,包括课题组实测数据、已发表的前人研究成果及公开的全国森林资源清查资料,以福建省杉木为例,建立不同尺度下杉木碳计量模型,并对各尺度模型之间的耦合换算进行探讨。同时,利用多种动态预估的方法对福建省杉木碳储量和碳汇能力进行模拟预测。最后基于学科交叉的角度,结合实际生产问题,对杉木碳计量模型的应用进行实证研究。主要研究结果如下:1.通过对个体尺度下杉木单木各器官及全树碳含量的建模研究,发现以胸径(D)和树高(H)两个维量组合(DH.D2H和D&H)作为自变量建模的拟合结果优于仅以胸径(D)建模的结果,且最佳的自变量组合形式是D&H(分离的二元变量)。本研究新提出的二元通用模型f(x,y)=(α·xβ·yε+γ·eη·x·y)θ(式中e为自然对数的底,α,β,γ,η,θ为待求参数)对叶、枝、地上部分及全树碳含量的拟合结果最好,而干和根最优拟合结果所使用的基础模型形式分别为幂函数和多项式函数。在进行区域数据外推拟合过程中,通用模型也表现出了较高的预测精度和较灵活的适用性。同时,通过理论分析和实例检验,利用比例平差法和非线性似乎不相关回归法都能较好的解决全树模型与各器官模型之间的碳含量预测结果不相容的问题。2.在林分尺度下,以福建省金森林业股份有限公司所经营的杉木人工林400个样地数据为基础,建立了两类不同理论基础的林分碳储量模型:(1)将林分蓄积量方程:ln M=h1+d1 ln N+g1 ln (Db1Hc1)与碳储量方程:C=a2(Db2-Hc2)Nd2联立,并根据二元材积方程约束b1=b2且c1=c2,便能得到与材积兼容的林分碳储量模型(式中C为林分碳储量、M为林分平均蓄积、N为林分密度、D为林分平均胸径、H为林分平均树高,b1、c1、d1、g1、h1、a2、b2、c2、d2为待求参数),所得到的林分碳储量模型对建模样本的拟合精度R2达到0.9165,验证样本的预测精度R2为0.8997;(2)以林龄、平均胸径、平均树高、林分密度、地位指数作为输入向量,林分碳储量为输出向量,采用BP神经网络和支持向量机2种机器学习方法进行建模,结果表明所建立的机器学习模型对建模样本和验证样本的拟合预测结果的精度R2都高于0.94,且在同样的优化算法(遗传算法)下,用支持向量机建立的模型比利用BP神经网络建立的模型精度更高,结果更为稳定。3.收集不同生长区域的杉木林分碳储量(生物量)研究数据,利用地理加权回归模型(GWR)建立区域尺度下的杉木林分碳储量模型。结果表明,GWR模型的拟合结果优于利用全局线性回归模型(OLR)建立的结果。此外,通过对马尾松(Pinus massoniana Lamb.)、落叶松(Larix gmelinii)、油松(Pinus tabuliformis Carriere)等主要针叶树种的建模,辅助验证的结果也同样显示GWR模型的精度优于OLR模型。由此可见,将数据的地理位置纳入回归参数中后所得到的拟合结果更接近于区域内的现实情况,因此在区域尺度下对林分碳储量进行研究时,可以尝试运用GWR模型进行建模。4.利用课题组调查的54株杉木单木含碳率实测结果,首次建立了福建省杉木全树平均含碳率模型,结果显示Peal-Read模型的拟合精度最高,二次曲线次之,且R2都在0.95以上。结合福建省杉木林生长的平均状态,计算得到各龄组的平均含碳率分别为:幼林龄(0.5388)、中龄林(0.5095)、近熟林(0.4865)、成熟林(0.4840)、过熟林(0.4867)。同时,利用不同的区域碳储量估算方法对福建省杉木林2009~2013年间的总碳储量进行估算,得到结果范围是38706.15~45413.37万t。5.以第八次(2009~-2013)全国森林资源清查数据基础,结合同期全国林业统计年鉴及福建省统计年鉴数据,利用C02FIX模型对福建省未来60年(2016~2075)杉木林碳汇潜力的模拟结果显示,福建省杉木林的累计净固碳量在2075年时应该处于1.6×108t至2.1×108t之间,占现有(2009~~2013)杉木林总碳储量的35%~50%,表明了杉木林拥有巨大的碳汇潜力,在福建省的森林碳汇作用中占有重要的地位。6.以福建省金森林业股份有限公司所经营的杉木人工林为实例,利用差分自回归移动平均(ARIMA)模型和多维时间序列(CAR)模型对林分碳储量进行时间序列分析,结果表明2种时间序列模型在短期内对林分碳储量的预测准确度都较高。同时,利用年均净现值法和林地期望价法探讨了碳汇木材复合经营目标下杉木林经济收益及最优轮伐期的确定,并比较分析了林分最优轮伐期和最大复合经济收益对立地质量、碳价格、利率这三个因素变化的响应,结果发现:(1)在目前木材价格远远高于碳价格的情况下,增加碳汇目标不会改变原来的采伐策略,但复合经济收益的增加是十分明显的;(2)复合经营目标下最优轮伐期随着地位指数上涨而提前,最大复合经济收益也随之显著提高;(3)碳价格在较长变化范围内(0~350元/t)对最优轮伐期的影响并不显著,只有当碳价格提高一定程度(350~600元/t),最优轮伐期才有提前的趋势,且随着碳价格的提高,复合经济收益明显增加;(4)随着利率增加,碳汇木材复合经营目标下杉木林最优轮伐期不断提前,但最大复合经济收益明显降低,因此高利率不利于杉木林的经营。7.利用收集的福建省杉木林分生物量(碳储量)研究数据,建立了福建省杉木净生产力与群落生长量和年凋落量之间的关系,结果表明在2009~2013年间,福建省杉木林乔木层净生产力为16.2022x 106t/a,平均净生产力达到11.9406t/(a·hm2)。结合建立的福建省杉木含碳率模型,计算得到福建省杉木林乔木层的固碳量为8.2299x 106t/a,其中幼龄林、中龄林、近熟林、成熟林和过熟林固碳量的比例分别为33.65%、23.38%、22.10%、18.56%、2.41%。全省杉木林乔木层的平均固碳能力为6.07t/(a·hm2),中龄林平均固碳能力最强,过熟林最弱。同时,计算得到2009~2013年间福建省杉木林乔木层固碳能力对同期化石能源碳排放的碳抵消率为8.37%,其中幼龄林的碳抵消贡献最大。