1. 问题背景
习近平总书记在2005年提出“绿水青山就是金山银山”,这句话强调了森林对于人们的重要性。近年来,随着温室气体(如CO2、CH4等)的大量排放,全球气候变暖,引发了冰川融化、土地荒漠化等严重的生态危机,生态环境问题开始备受关注。图1展示了人均碳排放量。森林是减缓气候变化的一个重要组成,森林可以将二氧化碳封存在植物和树木生产的产品中,如图2所示,这个过程被称为碳封存。
然而,随着森林的成长,释放的二氧化碳将与封存的二氧化碳量相等,甚至更多。此时森林将会变为碳源,不利于环境保护。通过砍伐森林,可以防止这种情况的产生,因此,如何制定合理的砍伐策略具有重要的意义 [1]。
2. 森林及其产品固碳量衡量模型
本文将模型分为三个部分:生物量模块、土壤模块和砍伐模块。生物量模块可以计算树木的固碳量,土壤模块用于计算凋落物的固碳量,砍伐模块用于计算木制品的固碳量。
2.1. 生物量模块
树木的固碳量是整片森林固碳量的主要部分,因此衡量树木的固碳量是最重要的。由于不同树木的生物特征差别较大,因此我们选择最有代表性的落叶林和针叶林作为树木的种类 [2]。
生物量模块使用森林面积和碳萃取量作为输入数据。将一片森林分为落叶林区域和针叶林区域,分别计算固碳量,以增加结果的合理性。模块计算两个指标:
净初级生产力(NPP)
净初级生产力表示树木所固定的有机碳中扣除本身呼吸消耗的部分 [3]。它可以通过下式进行估计 [4]
(1)
其中r表示森林的内增长率,B表示现存的固碳量,K为环境承载力。
固碳量变化量(B)
固碳量的变化情况与一段时间内的凋零量有关。这个变化情况可以使用微分方程表示 [4]
(2)
其中
代表树木不同部分的比例,
代表树木不同部分的自然损失率,具体数值见表1。
Table 1. (a) Proportion of different organs; (b) Natural losses from biomass to litter
表1. (a) 树木不同部分所占比例;(b) 不同部分的自然损失率
根据上述方程可以计算出
,它代表一年的固碳量增长量,通过迭代这个过程,可以得到一段时间内的固碳量变化情况。
2.2. 土壤模块
一片森林中含有大量的凋落物,包括落叶和枝干等,这些凋落物中封存了大量的碳,但会随着时间分解而释放,因此模拟这一过程对估算森林固碳量也十分重要。
本文将凋落物分为四种:树干类、主根和枝干类、细根类、树叶类 [5]。我们也将砍伐导致的废料算入凋落物中。将分解物分为两类:快速腐殖质和慢速腐殖质。整个过程如图3所示。
Figure 3. Processes of carbon sequestration in litter
图3. 凋落物的固碳量变化过程
凋落物转变为快速腐殖质
每类凋落物都有着各自的分解率,见表2。细根类和树叶类的分解率与气温和湿度有关,可以用如下式衡量 [6]:
(3)
其中T代表年平均气温,M代表五月到九月的平均湿度值与PET的差。
Table 2. Decay rate of different litter compartments
表2. 不同凋落物的分解率
腐殖质的变化过程
腐殖质分为两个部分:快速腐殖质和慢速腐殖质。凋落物分解后转化为快速腐殖质,快速腐殖质分解后变为慢速腐殖质。凋落物分解后一部分转化为二氧化碳,一部分转化为快速腐殖质,见表3。
Table 3. Fraction of litter decay converted into fast humus
表3. 凋落物分解后转化为快速腐殖质的比例
快速腐殖质的分解率取决于年平均气温,它通常用下式计算 [4]
(4)
慢速腐殖质的分解率大约为快速腐殖质分解率的0.007倍,并且大约有0.0033的快速腐殖质转变为慢速腐殖质。
通过模拟上述过程,我们可以得到一段时间内凋落物的固碳量变化情况。
2.3. 砍伐模块
现实中,砍伐的树木一部分用于燃烧,一部分用于制造木制品。在实际中,由于焚烧、加工过程的损耗等,木制品的固碳量更加接近于理论下限,而不是上限,因此,我们设砍伐的树木有四分之一做成了木制品。
我们将砍伐树木抽象为从森林中抽取固碳量,一次砍伐相当于从森林中取走了一定的固碳量。这些取走的碳,一部分作为木制品,一部分作为残留物开始腐烂的过程(残留物比例见表4),一部分作为燃料,一次性释放所有碳。
Table 4. Losses of different organs
表4. 不同部分的砍伐残留物比例
木制品在寿命结束后释放所有碳(通常是降解或燃烧),然而,木制品的寿命难以衡量。木制品的使用寿命大约为5~60年,且集中分布在10~20年间。因此,我们假设使用寿命服从正态分布
,如图4所示。
通过上述过程,我们可以模拟木制品在一段时间内的固碳量变化情况。
2.4. 模型测试与分析
通过将上文三个模块的结果相加,我们就可以得到一段时间内森林固碳量的变化情况。为了更好的评估砍伐的影响,我们将一个地区的所有树木视为一篇森林。
首先,我们不考虑木制品的固碳量变化。在一定时间后,执行砍伐操作。然后,我们观察森林固碳量的变化情况,如图5所示。我们可以看到砍伐后,森林固碳量立刻减少,但是随着时间的推移,森林逐渐恢复到与未砍伐时的固碳量一致。这与现实相符合,说明了模型的合理性。
Figure 5. The change in the amount of carbon sequestration with no felling vs felling
图5. 不砍伐和砍伐时固碳量的变化情况
接着,我们考虑木制品。我们让砍伐后森林固碳量减去不砍伐时森林的固碳量,如果结果大于0,则说明砍伐后森林固碳量增加了,反之,如果结果小于0,则说明砍伐后森林固碳量减少。如图6所示,我们将结果分为三个阶段:
1) 由于砍伐导致的浪费和燃烧等造成的固碳量减少。
2) 由于森林再生固碳量逐渐增加。最终,由于木制品的固碳量原因,总的固碳量可能会大于自然生长的森林固碳量。
3) 由于木制品使用寿命结束,释放二氧化碳,导致森林固碳量逐渐变少,但随着时间推移,最终趋于与自然生长的森林一致。
这三个阶段与现实相符,也从一方面说明了模型的合理性。
Figure 6. The change in the amount of carbon sequestration with products
图6. 考虑木制品时固碳量的变化情况
然而,当我们不断增大砍伐量和缩短砍伐周期时,固碳量变化情况如图7所示,我们可以看出砍伐量过大和砍伐周期过短时,都会产生负面影响。因此,通过模型的结果,我们可以制定最有益的砍伐策略。
Figure 7. Impact of felling strategies on carbon sequestration
图7. 不同砍伐策略对固碳量的影响
3. 模型在法国地区的应用
我们收集了法国各个地区的相关数据,为了对未来的固碳量进行模拟,我们对法国历年气温进行了拟合,年平均气温可以通过下式估计:
通过模拟我们得到了法国各个地区未来100年内固碳量的变化,如图8所示。
Figure 8. Changes of carbon sequestration in different regions of France
图8. 法国各个地区固碳量变化情况
从图8可以看出,并不是所有地区都适合砍伐,当该地区林地面积较小时,砍伐可能会导致森林固碳量降低,带来负面影响。当森林成长率较大时,砍伐可以带来更大的收益。
我们以Isère-Drôme-Ardèche地区为例,搜索得到最优砍伐率为落叶林1.47%、针叶林6.21%,砍伐周期为5年。
因此,通过模型结果,我们可以制定更加有益的砍伐策略,不仅提高了收益,还能有效提高固碳量,从而减缓温室效应。
4. 总结
本文建立的模型通过模拟森林自然变化情况、凋落物自然分解过程以及砍伐后制成木产品的使用过程来估算一个地区内森林固碳量的变化情况。根据固碳量的变化,可以制定更加合理的砍伐策略,以提高收益,并且减缓温室效应带来的影响,对森林管理和环境保护具有一定的实际意义。