气候变化、植被动态与莫桑比克玉米生产：来自矩方法分位数回归和机器学习的证据

《Frontiers in Sustainable Food Systems》：Climate change, vegetation dynamics, and maize production in Mozambique: evidence from method of moments quantile regression and machine learning

　　

本研究利用2002—2023年期间莫桑比克10个省份的平衡面板数据集，考察了气候变率和植被动态对玉米生产的影响。在确认横截面依赖后，采用第二代面板技术。Kao和Westerlund协整检验揭示了变量间存在稳定的长期关系。矩方法分位数回归（MMQREG）的结果表明，气候—玉米关系存在显著的分布异质性。降雨和归一化植被指数（NDVI）对玉米产量产生正向影响，且在低产省份观察到更强的效应；而温度则产生负向且统计显著的影响，尤其在中位数和上分位数。中介分析进一步证实，NDVI部分传递了气候变量对玉米生产的影响，为所提出的气候—植被—生产路径提供了实证支持。交互作用分析表明，充足的降雨和更健康的植被条件有助于缓解温度诱导的胁迫。此外，对莫桑比克北部、中部和南部的区域分析揭示了显著的空间异质性，半干旱省份表现出更高的气候敏感性。补充的机器学习分析确定K-最近邻（KNN）模型为玉米生产最准确的预测器，而温度、NDVI和降雨是最具影响力的预测因子。研究结果强调了在气候—农业关联中考虑分布和空间异质性的重要性，并突出了对气候智慧型、区域特定干预措施及数据驱动预测工具的需求，以增强莫桑比克的玉米生产力和韧性。

论文解读：气候变化、植被动态与莫桑比克玉米生产

一、研究背景与问题
气候变化是21世纪最紧迫的全球挑战之一，对农业生态系统和粮食安全构成深远威胁。莫桑比克因频繁遭受热带气旋、干旱、洪水等灾害，属于全球气候脆弱性最高的国家之一。农业占其国内生产总值的23%～25%，而玉米是该国最主要的粮作物，但生产高度依赖雨养条件，对温度升高和降雨变率极为敏感。然而，现有研究多采用均值回归方法，忽视气候效应的分布异质性；NDVI作为植被健康指标，其作为气候影响产量的中介机制尚不明确；同时，子国家级证据缺乏，且经济计量与机器学习方法的整合有限。因此，该研究旨在利用分布敏感的计量方法和预测模型，全面评估气候变率和植被动态对莫桑比克玉米生产的影响，为制定区域针对性适应政策提供科学依据。

二、研究内容与意义
该研究基于莫桑比克10个省份2002—2023年的平衡面板数据，综合运用第二代面板计量技术、矩方法分位数回归（MMQREG）、中介分析、交互作用、区域异质性分析以及机器学习算法，揭示了气候—玉米关联中的分布与空间异质性，并验证了“气候—植被—生产”传导路径。研究成果发表于《Frontiers in Sustainable Food Systems》。

三、主要技术方法
研究人员主要采用了以下关键技术方法：（1）Pesaran横截面依赖检验、第二代单位根检验及Kao和Westerlund协整检验，以确认数据性质与长期关系；（2）矩方法分位数回归（MMQREG）估计不同产量分位数上的气候效应；（3）逐步回归中介分析检验NDVI的传导作用；（4）交互项模型评估温度—降雨、温度—NDVI、降雨—NDVI的调节效应；（5）按北、中部和南部分组回归以考察空间异质性；（6）Dumitrescu-Hurlin（D-H）面板Granger因果检验；（7）九种机器学习模型（含KNN、SVR、RF等）预测比较，以决定系数（R²）、均方误差（MSE）、均方根误差（RMSE）、平均绝对误差（MAE）及10折交叉验证评估。数据来源包括莫桑比克农业与农村发展部/规划与政策局（MADER/DPP）、中分辨率成像光谱仪（MODIS）NDVI产品、世界气候变化知识（WCCK）数据库和世界发展指标（WDI）。

四、研究结果
4.1 描述性统计与初步检验结果
描述性统计和相关性分析显示，主要解释变量与玉米产量呈正相关，方差膨胀因子（VIF）检验排除多重共线性。数据非正态分布支持采用MMQREG。横截面依赖检验表明省份间存在截面相关，第二代单位根检验确认各变量平稳，Kao和Westerlund协整检验证明变量间存在长期稳定关系。

4.2 矩方法分位数回归（MMQREG）结果
MMQREG估计显示：NDVI在各分位数上正向影响玉米产量，但系数从25分位数的1.371降至75分位数的0.643，说明植被改善对低产省份的边际收益更大；温度在50和75分位数显著为负（系数分别为-0.096和-0.132），表明高产省份对热胁迫更敏感；降雨在25和50分位数显著为正且系数递减，显示水分约束对低产省份更强；耕地面积在所有分位数上显著为正但边际效应递减；农村人口效应不显著。

4.3 中介和交互作用分析
中介模型显示，温度显著降低NDVI，降雨显著提高NDVI；纳入NDVI后，温度系数从2.635降至2.053，降雨系数变为不显著，