省域视角下中国森林碳汇空间外溢效应与影响因素

付 伟,李 龙,罗明灿,陈建成,王福利

1 西南林业大学经济管理学院,昆明 650233

2 北京林业大学经济管理学院,北京 100083

全球经济的迅速发展带来了自然生态环境的破坏,CO2在大气中的浓度不断增高,给人类的生存环境带来了切身的危害,低碳减排成为全球共识。我国CO2排放量不断增长,已成为世界第一大碳排放国,低碳减排将成为我国未来发展的重要战略方向。在第七十五届联合国大会上习近平总书记第一次明确提出了“双碳”目标,中央财经委员会第九次会议提出“要把碳达峰、碳中和纳入生态文明建设整体布局”[1]。中国共产党第二十次全国代表大会上进一步提出要推动绿色发展,加快实施重要生态系统保护和修复重大工程。“双碳”目标的提出昭示着我国低碳减排的决心,要在减少“碳源”和增加“碳汇”两方面持续推进温室气体减排[2],与碳减排政策一致,大力推动森林碳汇的发展将是缓解全球气候变化的重要方式[3—4]。森林吸收大气中的CO2并储存就是森林的碳汇功能[5—6]。森林是陆地生态系统中最大的碳库,在“双碳”目标的实现过程中存在着巨大的生态效益[7]。森林的碳汇功能对全球的碳循环有着深刻的影响[8—9]。2020年联合国粮农组织公布的全球森林资源评估报告显示,全球森林碳储量约占全球植被碳储量的77%,森林碳储能力的提高在“双碳”目标的实现中具有重要作用[10]。森林碳汇的发展将成为应对气候变化及“双碳”目标实现的重要途径之一[11—12]。

国际上最早于1960年左右开始进行森林碳汇的研究[13]。马学威等[14]将森林碳汇的研究分为萌芽期和快速发展期两个阶段,2006年之前为萌芽期,2007年至今为快速发展期[15]。早期的研究代表Paul[16]使用森林生长模型将森林集约经营对碳储量的影响进行了探析,学者们在之后的研究中,将自然科学和经济学的融合纳入了动态分析框架和模型,森林碳汇的研究不断走向完善[17—18]。国内对森林碳汇空间溢出效应的研究存在着不同的观点,在森林碳汇的空间溢出效应是否显著及存在正向还是负向的空间溢出效应之间产生分歧[19]。薛龙飞等[20]通过研究我国31个省市的森林碳汇情况,证明了我国森林碳汇之间有着显著的负向空间溢出效应,但大量国内学者还是在省域层面上对森林碳汇进行研究[21—22]。森林资源禀赋类似的地区在森林碳汇量上存在较大差异[23],学者们从经济发展水平[24]、森林管理水平[25—27]以及森林灾害[28]等方面对其影响因素进行研究,这些因素对森林碳汇量都有着重要的影响。在碳汇的空间溢出效应研究中,孙建卫等[29]、孙丽文等[30]和徐敬俊等[31]分别对林业碳汇的经济和区域效应、省域碳锁定的空间外溢效应以及渔业碳汇的空间外溢效应方面进行了研究。

综上所述,国内外对森林碳汇的研究较为广泛,为学者对森林碳汇的研究提供了丰富的材料,但大量学者的研究视角仍关注于单一省份,对空间分布状况的研究较少且局限于一般性分析,对国家层面森林碳汇的空间溢出效应和影响因素未能进一步探究。由于森林资源的分布不是独立的,地理环境的相似性、政策引导的溢出性和生产创新的互补性等会产生较强的区域关联效应,地区间的协调发展问题值得深入探究。本文通过森林蓄积量扩展法测算我国31个省(市、自治区)在1993—2018年六次森林资源清查期间的森林碳汇量,对省域空间森林碳汇总量相关性特征进行对比,利用空间计量模型进行森林碳汇的外溢效应和影响因素分析。为我国各地区差异化森林碳汇政策的制定与实施提供支撑,利用我国林业政策的总体空间规划来综合统筹各区域森林政策,促进林业绿色高质量发展。

1 数据介绍与研究设计

1.1 数据介绍

1.1.1数据来源

选择1993、1998、2003、2008、2013和2018年我国六次森林资源清查的统计资料,数据源自于《全国森林资源清查报告》和《中国林业统计年鉴》,对相关数据取对数来消除异方差的影响。

1.1.2变量选取

本文在相关文献的基础上,进行了以下变量的选取:(1)林业产业发展水平,森林碳汇量与林业产业发展水平之间有着正相关的关系,选取林业产业总产值与森林面积的比值,即森林单位面积产值代表林业产业发展水平,记为F-level。(2)劳动力要素投入,在劳动力素质保持恒定的基础上,劳动者人数在很大程度上影响林业产量,从而对森林碳汇量有着重要影响。文中以林业系统年末从业人员人数代表劳动力要素投入水平,记为F-labor。(3)林业管理水平,造林行为是土地利用变化的方式,可以大幅提高植被碳汇[32]。由于造林行为对森林碳汇的影响存在10年作用的生长期[33]。以被解释变量早10年的造林面积代表林业管理水平,记为F-area。(4)森林采伐面积,因为采伐量对森林碳汇有着重要的关系,选取原木的采伐量作为影响森林碳汇的重要指标,记为F-logs。(5)森林灾害程度,不同地区的森林生态情况不同,森林火灾、病害、虫害对森林碳汇产生一定的影响,以森林火灾、病虫害面积与森林面积的比值代表森林灾害程度,记为F-disaster。(6)森林蓄积水平,森林蓄积量对森林碳汇有着主导影响,以单位面积森林蓄积量代表各地区森林蓄积水平,记为F-stock。

1.2 研究设计

根据空间计量经济学的观点,不同区域的事物和现象之间在空间上相互影响,存在一定的空间相关性。地理空间中各经济事物之间存在着不同相关程度和交互效应,是事物本身所固有的空间经济属性,空间计量经济学用空间外溢效应来定义和分析这种相关性。

本文在对我国31个省(市、自治区)的森林碳汇量进行测算的基础上,运用空间自相关中的全局和局部Moran′sI指数来考虑森林碳汇的时空格局,以空间计量模型分析森林碳汇空间关联性。对三种空间计量模型进行比较,选择最优模型对森林碳汇的空间外溢效应以及影响因素进行研究。

2 研究方法

2.1 森林碳汇量测算

本文以森林蓄积量扩展法核算森林碳汇量,公式中CS、CV、CB和CF依次为不同森林结构的碳汇量,分别为森林土壤碳汇量、林下植被碳汇量、森林生物量碳汇量和森林总体碳汇量[23]。Sij、Cij、Vij分别代表第i类地区第j类森林类型的森林面积、森林碳密度和森林单位面积蓄积量,具体公式如下:

CF=CS+CB+CV

(1)

CF=∑(Sij×Cij)+α∑(Sij×Cij)+β∑(Sij×Cij)

(2)

Cij=Vij×σ×ρ×γ

(3)

式中,α、β、σ、ρ和γ分别为林下植物碳转换系数、森林碳转换系数、微生物含量扩大系数、容积系数和含碳量。各项折算系数均根据联合国政府间气候变化专门委员会的默认参数值取值。其中,α取值0.195,β取值1.244。σ取值1.90,ρ取值0.50 t/m3,γ为取值0.50[34]。

2.2 空间关联性分析

探索性空间数据分析(Exploratory Spatial Data Analysis,ESDA)主要使用地理可视化技术来揭示空间数据的特征。文献中经常使用它来确定空间数据分布模式、聚集热点和空间异质性[35]。ESDA的优势在于能够提取其他方法无法识别的复杂空间现象,并为发现新的研究问题奠定基础[36]。它反映了特定空间中观测数据的相互依赖性。空间自相关方法分为全局空间自相关(Global Moran′sI)和局部空间自相关(Local Moran′sI)。全局Moran′sI指数可以发现空间的集聚或异常值并发现空间分布的差异和关联性。取值为(-1,+1),全局Moran′sI指数>0时,数值越大则空间正关联性越明显;全局Moran′sI指数=0时,则空间分布关联性较弱;全局Moran′sI指数<0时,数值越小则空间差异性越显著。全局空间自相关分析仅用一个值来反映研究区域间的空间差异的平均程度,却不能详细地说明区域间各个对象间具体的空间关联模式,没有考虑到空间异质性,因此无法反映地理单元内的局部空间相关性。所以,有必要使用局部自相关来确定具体集聚情况[37],通过进行局部自相关分析来描述邻近空间的属性值的相关程度。相应公式如下:

(4)

(5)

2.3 空间计量模型

在进行空间计量模型的研究中,先忽略各空间单元之间的相关作用,进行普通最小二乘回归(Ordinary Least Squares regression,OLS),预设回归模型如下:

yit=βixit+μit

(6)

式中,t、i、yit和xit分别为时间、观测样本、n×1维的被解释变量以及n×k维的解释变量,βi为k×1维的xit的影响系数,μit为n×1维的相互独立的随机误差项向量,且μit—(0,σ2)分布。

在最小二乘模型的基础上,检验相邻空间单元要素之间的相关性。若显著相关,则利用拉格朗日乘子检验(Lagrange Multiplier,LM)对三种空间计量模型做抉择:空间滞后模型(Spatial Lagged Model,SLM)、空间误差模型(Spatial Error Model,SEM)以及空间杜宾模型(Spatial Durbin Model,SDM)。

所有的空间计量模型给出以下的形式:

(7)

(8)

空间滞后模型(SLM)、空间误差模型(SEM)和空间杜宾模型(SDM)分别由以上公式的调整而来。

空间滞后模型(SLM)(λ=γ=0):

(9)

空间误差模型(SEM)(ρ=γ=0):

yit=φ+xitβ+μi+ηi+φit

(10)

(11)

空间杜宾模型(SDM)(λ=0):

(12)

3 结果分析

3.1 我国各省(市、自治区)森林碳汇量分析

我国森林碳汇总量丰富,年均碳汇量为1509780.85万t(表1)。根据1993—2018年各省市的碳汇量分布及变化情况,整体上看,我国森林碳汇分为三个梯队:第一梯队为西南省份和东北林区,由于地理因素,这些地区拥有大量的森林面积,为我国森林碳汇总量的主力军。第二梯队为南方林区,在我国森林碳汇中居于中等水平。第三梯队为东部沿海地区,碳汇量相较于前两个梯队较少。2018年,我国森林碳汇总量为1976234.52万t。西藏碳汇量为264438.45万t,居我国首位,占我国碳汇总量的13.38%。云南、四川、黑龙江三省的碳汇总量紧随其后,之间的碳汇量差距不大,四省的森林碳汇总量占我国碳汇总量的46.68%。天津、上海和宁夏由于其森林面积和地理因素的影响,碳汇量较低。2018年,上海碳汇量最低,但是从碳汇量的变化情况来看,上海碳汇增速较快,我国多个省份的碳汇增长幅度都较大,整体处于上升趋势。

表1 我国各省(市、自治区)森林碳汇量汇总/104tTable 1 Summary of forest carbon sequestration in all provinces (cities and autonomous regions) of China

3.2 我国各省(市、自治区)森林碳汇量空间关联性分析

3.2.1我国森林碳汇的全局空间关联性分析

本文利用Geoda 18.0和Stata 15.1对被解释变量-森林碳汇量进行全局Moran检验,选取Rook邻接矩阵将海南处理为与广东、广西邻接,作为空间权重矩阵[39]。结果如表2,1993—2018年,Moran′sI指数均为正值且通过1%的显著性水平检验。显示出我国森林碳汇具有显著的全局空间关联性。

表2 我国森林碳汇量全局Moran′s I指数检验Table 2 Global Moran′ I index test of forest carbon sequestration in China

图1 全局Moran′s I指数趋势图Fig.1 Global Moran′s I index trend chart

如图,Moran′sI指数的变化呈现出不同趋势。1993—2003年期间,森林碳汇Moran′sI指数为倒“V”形的先增加后降低的趋势,在2003年之后,又以较为稳定的趋势上升。Moran′sI指数的变动说明我国各省(市、自治区)之间的森林碳汇在空间上的集聚程度并不恒定。1988年之前,由于我国各地区加大对森林的基础保护并投入大量的资金支持,对各地区Moran′sI指数的提升起到了促进作用。1998—2003年之间,为了保护生态环境,维护生态安全,一系列林业重点工程在各地区实施。其中包括:天然林保护工程、“三北”防护林以及京津风沙源治理工程等,但是由于各地区保护工程实施的规模和进程存在差异,导致了各地区之间空间集聚效应的降低。2003年之后,各地区林权改革和退耕还林工程的实施,有效的保护了森林资源,带动了各地区的森林碳汇发展,Moran′sI指数开始稳步提升。

3.2.2我国森林碳汇的局部空间关联性分析

对森林碳汇的全局关联性分析并没有说明森林碳汇量在各地区之间的局部空间自相关情况。因此,应用局部Moran′sI指数散点图来进行局部的关联性分析。局部Moran′sI指数散点图四个象限表示四种空间集聚模式,其中,第一象限(H-H)表示高森林碳汇量省(市、自治区)被高森林碳汇量的其他省(市、自治区)所包围;第二象限(L-H),表示低森林碳汇量省(市、自治区)被高森林碳汇量的其他省(市、自治区)所包围;第三象限(L-L),表示低森林碳汇量省(市、自治区)被低森林碳汇量的其他省(市、自治区)所包围;第四象限(H-L)表示高森林碳汇量省(市、自治区)被低森林碳汇量的其他省(市、自治区)所包围。第一、三象限表示森林碳汇之间存在着空间正相关作用,第二四象限表示森林碳汇存在着空间负相关作用。

图2为我国森林碳汇局部Moran′sI指数散点图。图中各地区基本落入第一、二、三象限,对比1993和2018年的变化情况。1993年高高相关省份为:西藏和云南。低高相关省份为青海。2018年高高相关省份为西藏和云南。低高相关省份为青海和贵州。高高相关的地区主要集中在我国西南地区,优良的森林资源禀赋,森林管理和林业重点保护工程的实施,推动西南地区林业碳汇集聚发展。

图2 我国森林碳汇Moran′s I指数散点图Fig.2 Moran′s I index scatter diagram of forest carbon sink in China

3.3 我国各省(市、自治区)森林碳汇空间外溢效应与影响因素分析

3.3.1空间计量模型的选择

全局和局部空间相关性分析表明我国森林碳汇的空间外溢效应显著,为此,加入空间因素的影响分析外溢效应。

第一步:LM检验。首先利用Stata 15.1软件对所选取变量的面板数据进行OLS检验,并根据LM检验以及稳健性LM检验结果选取空间计量模型。如表3,林业产业发展水平、森林受灾程度对森林碳汇量存在着显著的负向作用,劳动力要素投入、林业管理水平、森林采伐面积和森林蓄积量水平对森林碳汇量存在显著的正向作用。

表3 OLS回归检验结果Table 3 OLS regression test results

如表4,LM检验和稳健性检验都通过了1%的显著性检验,4个检验均拒绝了原假设,说明本文所选变量兼具空间滞后和空间误差自相关效应。在三种空间计量模型中,SDM模型兼具两种效应,初步判断选择SDM模型来进行空间外溢效应的分析。

表4 LM和Robust LM检验结果Table 4 LM and Robust LM inspection results

第二步:三种空间计量模型的比较。通过豪斯曼检验对固定效应和随机效应进行选择,检验结果中,SDM的豪斯曼检验值为14.73,Prob>chi2 =0.0225,通过了5%的显著性检验。故选用固定效应的空间杜宾模型来进行外溢效应的分析。在空间杜宾模型中引用了解释变量的空间滞后项,并使用极大似然估计法可以有效地消除内生性问题。由表5,在三种空间计量模式的比较中,SDM模型的σ2=0.009,优于SEM和SAR模型,拟合优度R2=0.753也优于SEM模型。综合考虑,SDM模型是最理想的,最终选择固定效应的SDM模型进行分析。

表5 模型回归结果Table 5 Model regression results

第三步:似然比检验(Likelihood Ratio,LR)。对SDM模型进行LR检验,检验空间杜宾模型是否会退化成SLM和SEM模型。如表5,LR检验分别为15.50和25.77,均通过了1%的显著性检验。

综上,本文最终以固定效应的空间杜宾模型分析森林碳汇的空间外溢效应与影响因素。

3.3.2空间外溢效应分析

SDM模型的主要参数ρ的系数为0.238,存在着明显的正向效应并经过了10%的显著性试验,数据表明各省(市、自治区)的森林碳汇量之间具有空间溢出效应,并且相邻区域的森林碳汇水平每变化1%,本区域的森林碳会正向变化0.238%。对邻近地区森林碳汇的溢出效应,充分体现了地域上的空间关联性。林业产业发展水平、劳动力要素投入水平、林业管理水平、森林采伐面积、森林灾害程度和森林蓄积水平的空间效应分别为:-0.0410、-0.0274、0.0576、-0.0380、0.0528、0.522。其中林业产业发展水平、劳动力要素投入和森林采伐面积对临近省份为负向影响。林业管理水平、森林灾害程度和森林蓄积水平对临近省份为正向影响。

3.3.3空间外溢效应影响因素分析

空间杜宾模型很好地解释了各地区之间的空间经济相关性,然而参数估计结果并不能将直接作用效果和空间外溢效应直观表现出来,本文将所选变量对森林碳汇的影响效应分解为直接、间接和总效应来更加清晰地展示影响效果。具体结果如表6所示。

表6 空间效应分解Table 6 Spatial effect decomposition

一是各变量的直接效应。林业产业发展水平、劳动力要素投入、林业管理水平、森林采伐面积、森林灾害程度和森林蓄积水平的直接效应分别为-0.0905、0.0339、0.0207、0.00996、-0.0435和0.855。其中,除了劳动力要素投入和森林采伐面积没有通过10%的显著性检验,其余变量均通过1%的显著性检验。本地区的林业管理水平和森林蓄积水平的提高会对本地森林碳汇量的提高有着推动作用。林业管理水平的提升会提高森林生态环境保护,使森林的生长更加有序,发展更加充分,森林结构更加立体,促进森林碳汇量的增加。森林蓄积水平反应了地区的优势碳汇树种的数量,单位蓄积水平的提高会提升森林碳汇总量。林业产业发展水平和森林灾害程度为负向效应,表明当今林业产业的发展模式在一定程度和森林碳汇量之间存在着不匹配的因素,可能的解释是林业产业政策的侧重点和产业类型的偏好与森林碳汇不能很好的融合发展。而森林灾害与病虫害则直接导致树木生物结构的恶化,遏制林木的正常生长。

二是各变量的间接效应,林业产业发展水平、劳动力要素投入、林业管理水平、森林采伐面积、森林灾害程度和森林蓄积水平的间接效应分别为-0.0784、-0.0259、0.0802、-0.0443、0.0520和0.903。其中,只有林业管理水平和森林蓄积水平通过了1%的显著性检验,剩余变量均未通过10%的显著性检验,即本地区的林业产业发展水平、劳动力要素投入、森林采伐面积和森林灾害程度对相邻地区森林碳汇量的间接影响并不明显。林业管理水平的间接效应系数为正,表明一个地区的林业管理水平的提升会产生示范效应,使得空间相邻地区也会在林业管理水平上向该地区学习,造成相邻省份的林业管理水平也有一定程度的提升,从而提高森林碳汇水平。森林蓄积水平的间接效应系数为正,由于相邻省份的树种在地理位置上的相似性,树种蓄积能力类似,森林碳汇水平有着带动作用。

三是各变量的总效应,林业产业发展水平、劳动力要素投入、林业管理水平、森林采伐面积、森林灾害程度和森林蓄积水平的总效应分别为:-0.169、0.00802、0.101、-0.0344、0.00846、1.758。其中,林业产业发展水平、林业管理水平和森林蓄积水平均通过了1%的显著性检验,剩余变量没有通过10%的显著性检验。说明劳动力要素投入、森林采伐面积和森林灾害程度对森林碳汇量的空间溢出影响并不显著。林业产业发展水平的总效应、直接效应和间接效应系数均为负,但间接效应未通过10%的显著性检验。说明林业产业发展水平对本地区的森林碳汇量有负向影响,但是对相邻省份的影响并不明显。林业管理水平的总效应、直接效应和间接效应系数均为正,且均通过了1%的显著性检验。林业管理水平在直观上对森林结构的丰富度,森林树种的抚育度以及森林的生态保护水平等方面存在正向影响。表明林业管理水平不仅对本地区的森林碳汇量有推动作用,而且对相邻地区产生示范和拉动效应,带动相邻地区森林碳汇水平提升。森林蓄积水平的总效应、直接效应和间接效应系数同样均为正。森林蓄积水平反映了当地树种的综合碳汇能力,且临近地区之间由于地理空间的相似性,树种碳汇能力差异不大,因此,森林蓄积水平对本地区和相邻地区的森林碳汇水平都存在着促进作用。

4 结论与启示

4.1 主要结论

根据我国森林碳汇量的核算及空间外溢效应和影响因素分析结果,本文从森林碳汇量的整体和变化情况、空间外溢效应以及各因素的影响三个方面得到以下结论:

一是我国整体上具有丰富的森林碳汇量。年均碳汇量为1509780.85万t。六次森林资源清查期间,森林碳汇量整体呈上升趋势。由于森林资源禀赋的差异,各地区的森林碳汇量存在较大差异。西藏、云南、四川、黑龙江四个省份的森林碳汇总量较大,占我国碳汇总量的40%以上。天津、上海和宁夏的森林碳汇量处于较低的位置,但是在增长幅度上面,上海市森林碳汇量增速最快。

二是我国各省(市、自治区)森林碳汇量空间关联性分析。我国森林碳汇具有显著的全局空间关联性。1993—2003年期间,森林碳汇Moran′sI指数为倒“V”形的先增加后降低的趋势,在2003年之后,又以较为稳定的趋势上升。Moran′sI指数的变动说明我国各省(市,自治区)的森林碳汇在空间上的集聚程度并不恒定。局部Moran′sI指数反映出高高相关省份为:西藏和云南。低高相关省份为青海和贵州,各个时期变动不大,表明我国森林碳汇已形成高高相关格局,且高高相关的地区主要集中在我国西南地区,由于丰富的森林资源,森林管理和林业重点保护工程的实施,推进了西南地区林业碳汇集聚发展。

三是我国森林碳汇的空间外溢效应及影响因素分析。各省(市、自治区)的森林碳汇量之间存在显著的空间外溢效应。且相邻区域的森林碳汇水平每变化1%,本区域的森林碳汇量正向变化0.238%。其中林业产业发展水平、劳动力要素投入和森林采伐面积对临近省份为负向影响。林业管理水平、森林灾害程度和森林蓄积水平对临近省份为正向影响。对空间外溢效应进行分解:直接效应中,林业管理水平和森林蓄积水平对本地森林碳汇量有正向影响,林业产业发展水平和森林灾害程度为负向效应。间接效应中,本地区的林业产业发展水平、劳动力要素投入、森林采伐面积和森林灾害程度对相邻地区森林碳汇量的间接影响并不明显。林业管理水平和森林蓄积水平对相邻地区森林碳汇量有正向影响。总效应中,林业管理水平和森林蓄积水平对本地和相邻地区森林碳汇量有正向影响,林业产业发展水平对本地区的森林碳汇量有负向影响,但是对相邻省份的影响并不明显。

4.2 政策启示

根据上述研究结论,我国在制定和实施差异化林业碳汇政策时有以下启示:

一是我国森林碳汇量在空间上存在着显著的空间相关性与空间外溢效应,对差异化林业碳汇政策的制定与实施应考虑区位因素,结合相邻地区的森林资源禀赋和林业产业发展水平进行林业发展的空间规划,并综合协调各地区林业政策,在“山水林田湖草沙”的命运共同体理念引领下,合理实现森林碳汇林业绿色发展。

二是利用好森林碳汇的空间外溢效应,林业产业发展水平和森林采伐面积对森林碳汇量的总效应影响为负,说明林业产业绿色发展进程中应进行结构上的调整,要调整林业第一产业占比,优化林业第二、第三产业融合发展,促进森林碳汇水平的提升。森林采伐量应合理制定采伐配额,优化采伐侧重于人工林。发挥西南地区与东北林区的辐射作用,带动相邻地区森林资源的合理利用与规划,充分吸收碳汇林项目实施经验。

三是根据各因素对我国森林碳汇影响程度的不同进行相应变革,林业产业的发展模式和森林碳汇之间存在着不匹配的因素,未能很好的与森林碳汇融合发展,要加快林产品加工业发展规模的扩大和发展方式的转变。在林业管理水平和森林蓄积水平方面,森林环境稳定、树种结构多样化、优势碳汇树种扩大化是林业绿色发展的基础,要保持生态、社会、环境三方的协调,差异化实施林业政策制度,确保森林使用权在流转、抵押和担保上的合规性,提供森林碳汇项目实施的政策保证。严格实施采伐限额制度,提高森林管护能力,科学建立应对森林火灾、病虫害的方案。