The Origins and Control of Forest Fires in the Tropics 基本信息️ 一、 引言部分简要总结 (Introduction Summary) 二、 研究方法详细介绍与分析 (Methodology Framework)图片名Figure 1图例说明Figure 1 | 基于 MODIS 逐日热点连续追踪单次火灾起火点与空间蔓延的算法示意图。黑色边框方块代表第 1 天起火的初始起火像元Ignition area红色方块代表第 2 天及后续连续日期在 1 像元缓冲区内蔓延燃烧的受灾区域Spread area蓝色矩形框 A、B、C、D 分别标识出 4 起独立火灾事件总过火面积定义为黑框起火区与红色蔓延区的空间并集。图片名Figure 4图例说明Figure 4 | 目标像元 6 km 半径缓冲区内地块权属构成与空间分类示意图。中心红色方块为进行点火决策的目标像元红色圆形为 6 km 搜索半径包含 137 个网格像元根据各像元质心所处位置精确识别出橙色区域本企业特许权林地、蓝色区域未出租国有生产林地以及深绿色区域法定自然保护林用以量化潜在火势蔓延波及不同产权主体的面积比例。 三、 主要研究结果与充实数据分析 (Main Results Quantified Data)1. 人类活动与清地工序构成了热带林火的核心根源对印尼主要岛屿卫星火点的空间叠置分析表明林火的发生具有极强的土地利用偏好性。详见 Figure 2。图片名Figure 2 图片位置留空 图例说明Figure 2 | 印度尼西亚苏门答腊岛廖内省Riau起火点与林区权属分布的微观空间映射图。底图展示了 7,700 km² 范围内检测到的逾 3,400 个起火点红色方块黄色虚线与橙色稀疏虚线分别标出木材纤维与油棕特许权边界深绿为保护区占布鲁德国家公园几乎无火浅绿为普通林业辖区白色为辖区外城镇火点极少清晰揭示起火点高度依附于种植园清地作业并刻意避开保护区与聚落。图片名Figure 3 图片位置留空 图例说明Figure 3 | 廖内省局部 2012 年毁林砍伐与 2013 年火灾起火点的高空间分辨率时空重叠图。橙色斑块代表基于 30 m Landsat 影像解译的 2012 年实际森林损毁区域红色方块为 2013 年 MODIS 监测到的起火像元直观佐证绝大多数人工点火高度精准地发生在前一年刚遭砍伐的迹地内部。2. 火灾蔓延的极端长尾风险与气象因子的可预测性对单场火灾物理特征的统计显示87% 的火灾仅燃烧 1 天89% 未蔓延出起火像元但火灾外溢具有极强的长尾毁灭性持续时间最长的火灾达 36 天最大单体火灾蔓延面积达 656 km²蔓延至初始起火面积的 466 倍。特许权内多日火灾蔓延产生的过火面积中有 28% 直接烧到了本特许权边界之外。 气象敏感性分析Table 3表明火灾蔓延面积高度可预测3. 土地权属差异引发的显著外部性与企业策略性避险在全要素固定效应模型中Table 4 Panel B 第 6 列气象蔓延指数与周边土地权属构成的交互项揭示了企业的精明策略4. 私人谈判机制检验科斯定理在低交易成本下的成立为探究私人谈判能否解决林场间的外部性研究考察了周边仅存在自身林地与单一邻近私营特许权的交界样本Table 55. 公共规制路径检验庇古式威慑的现实印证与治理边界对政府执法行为的经验检验Table 9表明行政处罚的预期风险塑造了企业的行为边界 四、 讨论部分的不足与未来展望 (Discussion Limitations Perspectives)消除“公地悲剧”清晰界定边缘林带管护权责 研究表明企业过度用火的外部性主要宣泄在缺乏排他性产权的未分配国有生产林上。未来政策不应单纯在法律层面一味发布不可执行的“全面禁火令”而应将此类管理真空的公有林地尽快明确租赁权属、划入严格管护的自然保护地网络或确权赋予周边原住民社区进行共管借助明确的排他性产权消除无主之地的监管套利空间。完善面向空间外溢的庇古式差别化定额罚则 印尼政府当前的执法表现出对人口密集区受灾的高度威慑但在自然资源毁损上执法羸弱。鉴于火灾蔓延具有高度气象可预测性生态环境法制需推行“高火险气象窗口严惩制度”在干燥大风天气下大幅提高在林区边界点火的起刑点与环境损害连带赔偿倍数重塑企业在气象危险期的风险预期。提供机械清地技术补贴以对冲成本落差 由于火烧清地比机械清地节约 44% 至 70% 的直接资金支出只要重型机械租赁与推土成本居高不下非法纵火便难以从根本上杜绝。国际气候减排资金如 REDD应设立针对热带小农与种植园绿色翻耕的专项补贴机制降低无火机械化整地的经济门槛从技术供给侧消除企业依赖高风险点火技术的经济动力。未来研究展望与政策干预方向中分辨率遥感在亚像元与单日蔓延上的物理局限MODIS 1 km 像元无法精细识别 1 天之内在数百米尺度上的局地微蔓延将起火首日连通的所有火点一律归为起火区在技术上可能略微低估了超小规模火灾的即时外溢扩散比例。特许权边界数据时间维度的静态简化受印尼官方地籍档案更新限制本研究所用特许权矢量面文件锚定在样本中期的 2008 年未能完全捕捉 15 年间个别林场租赁权属转让、抵押或行政边界微调等跨期动态。空间外部性范围的设定约束论文侧重于地表物理蔓延这一直接邻近的局部负外部性Local externality以 6 km 作为关键影响域未将跨越数百公里乃至跨国界的大气烟尘气溶胶扩散与人群呼吸系统长远健康损害进行直接因果联动建模。自然点火信号的潜在微弱残留虽然通过与伐木时序联动和横向土地用途比对证实了绝大多数火灾的人为属性但在极端干燥年份雷击等自然火源对极偏远高海拔针叶林带的点火概率仍无法彻底置零排除。研究不足与不确定性来源庇古式精准惩治的经验证据分析 2015 年大火后被政府立案调查的涉案企业发现火势蔓延波及的绝对人口规模是诱发政府调查最主要的决定性因子系数 0.0007, P 0.01相反火灾烧毁多少未出租国有林并不会显著增加企业被立案追责的几率系数 -0.017, P 0.10。这在实证上完美呼应了企业的点火心理企业之所以极力避开居民区用火正是因为预料到了针对人口聚落受损的“庇古式惩罚税罚金与刑责”具有极高的执行概率。官员反腐与国有产权的制度限度利用印尼反贪会查处 26 起林业贪腐判决的交叠双重差分检验Table 10显示地方主官落马并未显著降低后续辖区内的点火水平或削弱外部性转嫁国有林业企业Table 11虽然整体点火率比私企低约 40%系数 -0.4903, P 0.01但在火灾高风险期向无主公有地转嫁蔓延风险的边际行为上与私营企业同样懈怠。自身与邻里地块的对称待遇在双边垄断的边界上无论是在截面回归还是气象火险交互回归中自身林地与邻居特许权林地的系数均无统计学显著差异Table 5 Panel B 第 6 列交互项系数为 -0.0012, P 0.10。企业对待相邻私有企业的财产风险表现出如同对待自己土地一样的防范水准。林木资产价值对谈判的精细约束结合 Hansen 毁林数据进一步区分林木状态发现Table 6企业对待两类地块表现出惊人一致的资产保护倾向蔓延风险烧向自身未砍伐林地系数 -0.0043, P 0.01与烧向邻近企业未砍伐林地系数 -0.0041, P 0.01的负向抑制强度几乎完全一致而对于双方已经皆伐砍空的废弃地块双方均不在意火势蔓延。这提供了科斯定理在现实中发挥作用的罕见清晰证据当产权清晰且谈判主体单一交易成本极低时私营主体间可以通过默契或私人协调内生化火灾风险。声誉与自愿认证机制的失灵相反大企业规模拥有多个特许权或大面积林场Table 7虽然平均点火率略低但在高风险期并没有额外抑制向公有地转嫁火灾外部性的行为加入可持续棕榈油圆桌倡议RSPO认证的林场向周边公共林地转嫁火灾的负外部性倾向亦未得到任何统计意义上的改善。对“公地”的放任侵害过度用火当周边土地为政府未出租生产林缺乏排他性产权与监管管护时交互项系数为显著的正值 0.0047标准误 0.0013, P 0.01。表明在风大、干燥的天气下企业明知火势极易失控依然会倾向于点火将火势外溢的潜在烧毁风险转嫁给缺乏维权主体的公共山林。对邻近人口聚落的极力规避与此形成极鲜明对照的是当火势可能蔓延至林业局辖区外的村镇私有地时交互项系数为高度显著的负值 -0.0057标准误 0.0019, P 0.01。在蔓延高风险月份企业会高度克制在居民点附近动火以规避可能招致的严重民事纠纷或官方严打。外部性总规模的反事实推演若企业对待所有周边土地都像对待村镇私有地那样谨小慎微印尼木材与油棕林场内的点火总量将锐减 55% 至 57%以此换算若将此类外部性完全内生化仅 2015 年单年即可挽回相当于印尼当年 GDP 0.2% 的直接经济损失减少 0.7 Gt 碳排放占全球化石燃料排放的 7.1%并避免多达 14,630 名成年人与 4,226 名幼童因空气污染而过早死亡。地表风速每增加 1 km/h火灾蔓延面积显著扩大 15.1%P 0.01月均温每升高 1°C蔓延面积激增 57.0% 至 77.7%P 0.01降水量每增加 1 mm蔓延面积则受抑收缩 0.66%P 0.01。特许权利用类型间的极化分布在所有林业特许权内部火灾中需要彻底皆伐并整地的油棕特许权占特许权总面积的 28%和木材纤维人工林占总面积的 22%分别贡献了 46% 和 42% 的起火点而采取可持续选择性采伐的传统伐木特许权占总面积的 51%仅占起火总数的 12%。单位面积计算人工种植园地块的起火发生率是天然林采伐区的整整 8 倍毁林与点火之间的强因果时间序列像元级固定效应泊松面板回归Table 2显示滞后 1 年的森林毁伐对次年点火具有压倒性的正向驱动系数达 1.3472标准误 0.1321, P 0.01。这意味着一个完全发生毁林的 1 km² 像元其在次年发生起火点燃的概率激增 285%更为关键的是滞后 2 年系数 -0.3081, P 0.05与滞后 3 年系数 -0.3492, P 0.05的效应均显著转负。这一先正后负的时序动态排除了“毁林后林地变得容易自燃”的自然假说确凿证明火灾是企业紧随采伐之后、在种植幼苗前清除残余杂木植被的专用作业技术。详见 Figure 3。科斯检验将样本严格限制在仅与单一私营林场接壤的像元子集消解多方谈判交易成本并细分邻地是否包含未砍伐高价值原生立木。庇古检验构建二值 Logit 概率模型以 2015 年印尼政府立案调查企业清单为因变量回归解释火灾蔓延至各权属地块及居民区人口数对政府实质执法立案的边际贡献。科斯谈判与庇古威慑的实证检验结构毁林对点火驱动的泊松面板模型 检验林火是否属于清地工序构建像元-月份-年份层面的泊松固定效应模型 E[Ignitions_imt] exp( beta_1 * Forestloss_it-1 beta_2 * Forestloss_it-2 beta_3 * Forestloss_it-3 gamma_i delta_m delta_t ) 式中加入像元固定效应 gamma_i吸收固定地形与土壤易燃性差异、月份固定效应 delta_m吸收季节性规律与年份固定效应 delta_t控制全域年际干旱冲击标准误在 50 km × 50 km 空间网格层面上进行聚类调整。气象蔓延指数Weather Spread Risk构建 在控制像元固定效应下利用泊松模型拟合实际蔓延像元数与月均风速Windspeed、总降水Precipitation及月均温Temperature的响应弹性估计出客观表征火灾潜在失控蔓延概率的标准化气象火险综合指数 Weather Spread Risk。基于气象火险交互项的外部性因果识别模型 在每个像元周围建立 6 km 半径空间缓冲区包含 137 个像元对应多日火灾最大蔓延距离的第 90 百分位数见 Figure 4。将气象火险指数与周边各权属像元数NeighbourLandType_i^j构建交互项 E[Ignitions_imt] exp( beta_1 * WeatherSpreadRisk_imt ∑ beta_2^j * NeighbourLandType_i^j * WeatherSpreadRisk_imt beta_3 * AvgPopnDens_i * WeatherSpreadRisk_imt beta_4 * X_i * WeatherSpreadRisk_imt gamma_i delta_m delta_t ) 基准遗漏类别为“自身特许权林地”。交互项系数 beta_2^j 的正负与量级直接测度了在天气极易助长火势蔓延的危险月份企业是否会因为火灾将烧向特定外部土地而表现出“更加肆无忌惮地用火正系数存在外部性侵害”或“更加谨慎地克制用火负系数存在有效威慑”。计量经济学模型与识别策略设计火灾起火点与多日蔓延轨迹重构算法 论文首创了基于卫星热点的火灾轨迹重构算法详见 Figure 1 以第 X 天连通的火点像元群 A_X 为基础向外构建 1 像元空间缓冲区若第 X1 天有火点落入该缓冲区内则将其判定为同一场火灾在时空上的连续蔓延随后围绕第 X1 天新火点继续迭代建立缓冲区直至火点熄灭遇到多云遮挡像元时自动顺延至下一无云日。该算法成功将离散的火点聚合为具备明确首日起火位置与多日蔓延边界的独立火灾实体。MODIS 每日热异常火点MOD14A1 Level 3提取全球 1 km 网格每日 4 次过境的卫星热异常数据记录像元是否处于燃烧状态。土地权属与特许权边界数据由全球森林观察Global Forest Watch及政府地籍矢量图层提取精确划分出五大权属空间本企业特许权地、其他私营特许权地木材纤维、油棕或伐木特许权、法定保护林国家公园与水源保护林、林业辖区外私有地村镇与小农集中分布区以及未出租国有生产林Unleased productive forest即缺乏实质管护的“无主林地”。Hansen 全球高分辨率毁林数据利用 30 m 分辨率 Landsat 逐年森林覆盖变化数据计算每个 1 km 像元内部历年被砍伐清地的面积比例。气象再分析驱动场整合 NCEP-DOE Reanalysis 2 逐日地面风速、MSWEP v2.8 逐月降水量11 km 空间网格面积加权聚合及 CRU TS v4.06 逐月气温反距离权重插值至 1 km。政府行政处罚与司法反腐案卷手动匹配印尼环境与林业部在 2015 年大火后对涉案企业立案调查的官方名单与特许权企业名录收集印尼反贪污委员会KPK及最高法院数据库中涉及林火领域的地方官员刑事判决卷宗。技术路线与数据空间范围界定 研究聚焦印度尼西亚主要多林岛屿苏门答腊、加里曼丹、苏拉威西、巴布亚等排除爪哇与小巽他群岛的国家林业辖区Forest estate内部时间跨度为 2000 年 10 月至 2016 年 1 月。全区共追踪识别出 107,334 起森林火灾核心计量分析重点针对发生在油棕与木材纤维特许权地块内部的 39,189 起火灾展开分析样本覆盖超过 22 万个 1 km² 像元。核心数据源提取与数据清洗本文试图解决的核心科学问题 热带林火究竟是气象自然引发还是企业清地驱动的人为产物周边土地的产权属性如何影响企业点火决策的外部性权衡在多大程度上科斯式私人谈判能自主解决邻近林场间的火灾蔓延问题而政府的庇古式惩罚又在多大程度上重塑了企业的点火威慑微观实证证据的严重匮乏虽然科斯的私人产权谈判理论与庇古的政府矫正性惩罚理论在形式体系上已十分成熟但实证研究极少直接观测并量化企业在做出破坏性决策时究竟在多大程度上依据“环境损害是否属于外部性”来调整自身行为。人为点火与自然火险的混淆以往宏观统计往往将卫星火点直接归因于干旱等气象因素缺乏将单次火灾从最初起火点Ignition到后续多日空间蔓延Spread进行全生命周期高分辨率轨迹重构的技术手段无法有效解耦企业的主动用火决策与被动气象驱动。空间外部性在产权交界处的传导机制不明热带林区通常由复杂的土地权属拼接而成自身租约地、邻近私营租约地、法定保护区、未分配国有林及居民区私有地各权属所受法律保护程度不同企业点火行为是否随周边产权配置与气象外溢风险动态调整缺乏大样本空间计量分析。研究的现实背景与科学背景 自马歇尔、帕累托与庇古开创环境经济学以来未受补偿的外部损害外部性始终被视作市场失灵与环境恶化的核心根源。在热带低收入国家利用火烧清理林地以种植油棕或纸浆木材纤维是企业普遍采取的低成本手段其作业成本比推土机等机械化作业便宜 44% 至 70%。然而失控的林火蔓延会引发巨大的跨境雾霾、生态灭失与碳排放例如 2015 年印尼特大火灾排放了约 4 亿吨 CO2 当量并在东南亚多国造成逾 10 万人超额死亡。前人研究的局限性与核心瓶颈中文题目热带森林火灾的成因与管控基于外部性视角的环境经济学实证研究英文题目The Origins and Control of Forest Fires in the Tropics期刊名称Review of Economic Studies发表时间2024年12月接收2025年10月在线出版2026年正式刊出Vol. 93, No. 4, pp. 2348–2389作者信息Clare Balboni伦敦政治经济学院、Robin Burgess伦敦政治经济学院、Benjamin A. Olken麻省理工学院DOIhttps://doi.org/10.1093/restud/rdaf088核心观点 环境外部性是驱动气候变化、毁林与生态退化的根本动因但在私人微观决策层面外部性究竟如何在多大程度上驱动企业行为尚缺乏大规模微观实证证据。本研究利用印度尼西亚 2000–2016 年逾 15 年的高分辨率卫星每日火点数据覆盖超过 107,000 起林火事件结合气象蔓延风险与周围土地产权权属的空间外生变异系统检验了企业在清地用火中的外部性考量。研究发现当火势外溢风险主要由产权薄弱的未出租国有生产林承担时企业表现出明显的过度用火倾向相反当火势可能蔓延至林区外人口密集村镇时企业则显著抑制点火行为。若将所有土地外溢风险完全视同居民区内生化全区点火规模将缩减 55% 至 57%。进一步检验表明在仅有两方私营林场的边界上企业对待邻近私有林地与自身林地的火灾蔓延风险并无差异符合科斯定理Coase Theorem关于低交易成本下私人协商消除外部性的预期政府在 2015 年特大火灾后的行政处罚高度集中于波及人口聚居区的火情印证了企业规避居民区火灾契合庇古Pigouvian威慑效应。