生态恢复与气候变率对中国黄土高原空气质量改善影响的解耦：一项十年观测评估（2015–2024）

《Ecological Indicators》：Decoupling the effects of ecological restoration and climate variability on air quality improvement over China's loess plateau: A decade-long observational assessment (2015–2024)

　　

解耦生态恢复和气候变率对区域空气质量改善的贡献仍然是基于自然的解决方案评估的关键挑战。在此，研究人员整合了卫星NDVI（归一化植被指数）、地面空气质量监测（PM₁₀、SO₂、NO₂、CO、O₃）和气象数据，构建了中国黄土高原跨十年的（2015–2024）站点-年份面板数据。利用随机森林（RF）气象调整、双向固定效应面板模型、结构方程模型（SEM）和阈值回归，研究人员提供了四项新贡献。首先，研究人员分离了植被恢复的独立净效应：在去除气象波动和时空固定效应后，在半干旱区，每增加0.1单位NDVI可使PM₁₀降低9.06 μg·m^?3（95%置信区间CI: ?11.42至?6.70），在半湿润区降低6.88 μg·m^?3（95% CI: ?9.23至?4.53），而在干旱区无显著效应。其次，研究人员确定了一个关键NDVI阈值（≈0.45）：低于该水平，PM₁₀降低可忽略不计；高于该水平，净化效应急剧加速。第三，研究人员量化了中介路径：植被通过增加湿度、降低风速和辐射间接降低PM₁₀（中介比例：34–47%）。第四，研究人员诊断了O₃权衡：绿化仅在与高NO_x（氮氧化物）区域相关的东南部、高BVOC（生物源挥发性有机化合物）物种（刺槐*Robinia pseudoacacia*、杨属*Populus* spp.）占主导的地区与显著的O₃增加相关；该效应在统计上由温度和辐射增加所中介。研究人员的发现提供了首个区域尺度因果证据，表明在充足水分条件下生态恢复独立改善空气质量，建立了最大化协同效益的定量覆盖目标（NDVI ≥0.45），并证明O₃结果具有物种依赖性——要求在污染区域进行BVOC敏感的造林。该分析框架可推广至全球其他水分受限的恢复区域。

在气候变化与城市化交织的背景下，区域空气污染威胁人类健康与生态系统安全，各国侧重“减排”策略，但大规模植被恢复可能通过改变地表过程调节空气质量，这种基于自然的解决方案（NbS）备受关注。中国黄土高原作为世界最严重侵蚀区之一，自1990年代末实施了退耕还林和三北防护林等大规模生态工程，但其空气质量协同效益尚不明确。现有归因研究几乎只聚焦人为减排，将植被视为被动变量；气象归一化常采用线性或加性模型，忽视非线性交互；植被影响空气质量的机制——直接沉降与间接气象中介——仍属推测，而植被引发O₃权衡的条件尚不清晰。为此，研究人员开展了本研究，旨在提供区域尺度因果证据，量化生态恢复对空气质量的独立净效应，识别关键阈值，并诊断O₃权衡。论文发表在《Ecological Indicators》。研究整合了卫星NDVI、地面空气质量监测（PM₁₀、SO₂、NO₂、CO、O₃）和气象数据，构建了2015–2024年黄土高原155个站点的站点-年份面板数据。采用随机森林气象调整（去除气象波动）、双向固定效应面板模型（控制站点和年份固定效应）、结构方程模型（检验中介路径）以及分段回归与广义可加模型（识别非线性阈值），系统评估了植被恢复的独立作用。

结果部分：
**3.1 植被恢复与空气污染物的时空演变特征**：通过Theil-Sen趋势分析和Mann-Kendall检验，研究人员揭示了NDVI在十年间显著上升（p<0.01），且2018年后加速；PM₁₀、SO₂、NO₂、CO呈显著下降趋势，而O₃略有上升，尤其在2019–2021年间。空间上，PM₁₀和SO₂降幅在东南部大于西北部，O₃上升集中在东南部（如陕西关中平原）。NDVI与污染物趋势的象限分析（同现/权衡关系）表明，大多数站点位于象限II（NDVI上升且PM₁₀、SO₂、NO₂下降），呈现协同；而O₃的象限I（NDVI上升且O₃上升）站点集中在东南部，表明权衡。

**3.2 构建“气象调整”的污染物浓度序列**：利用随机森林模型（以7个气象因子为预测变量）拟合污染物浓度，计算气象异常（观测值减预测值）。模型验证显示测试集R²为0.481–0.646，RMSE合理（如PM₁₀为14.19 μg·m^?3）。调整后的浓度序列趋势更清晰，证实长期改善并非由有利气象驱动。稳健性检验中，部分依赖分析表明NDVI对PM₁₀的边际效应（每0.1 NDVI降低7.8 μg·m^?3）与面板固定效应结果一致，且SO₂的负对照检验（NDVI–SO₂关系不显著）排除了混杂。

**3.3 NDVI对空气污染指标的净效应及异质性分析**：双向固定效应面板模型表明，NDVI仅在半干旱和半湿润区对PM₁₀有显著负向净效应（每0.1 NDVI分别降低9.06和6.88 μg·m^?3），干旱区不显著；对O₃无显著整体效应。NO₂在半干旱区有显著负效应，而SO₂和CO效应不显著。气候分区显著调节NDVI–PM₁₀关系，呈水分依赖性梯度。

**3.4 植被恢复与气象变率对空气污染物的驱动机制**：冗余分析（RDA）确认气象因子解释了污染物变异的绝大部分。结构方程模型（SEM）定量解析了中介路径：对于PM₁₀，NDVI通过降低风速产生正效应（不利），但通过增加湿度与降低辐射产生更强负效应，总间接效应占38%，直接效应（如干沉降）占62%。对于O₃，NDVI通过增加温度和辐射产生正间接效应。阈值分析通过分段回归和广义可加模型（GAM）确定PM₁₀还原的NDVI阈值约为0.45，低于该值效应微弱，高于后急剧加速；O₃也存在非线性关系，阈值前后方向相反。

讨论部分：研究首次在区域尺度提供因果证据，证明生态恢复在半干旱/半湿润区独立降低PM₁₀，这源于植被的生物地球物理调节（如增加蒸散发、提高湿度）。气候分区水分可用性（降水>250 mm）是关键调节因子，干旱区植被功能受限。O₃权衡在高NO_x、高BVOC物种占主导的东南部出现，与温度和辐射增加的中介效应相关，区别于欧美城市研究。NDVI阈值（≈0.45）为生态工程设定了定量目标，强调覆盖度而非面积扩张。局限性包括：无法完全排除反向因果、人为排放时序混杂、年度数据聚合可能掩盖短期动态，且阈值具有区域特异性。

结论部分（翻译论文结论）：本项研究提供了生态恢复对黄土高原空气质量改善贡献的因果评估。四大关键发现包括：（1）植被仅在降水超过250 mm的区域独立降低PM₁₀，效应集中于250–550 mm过渡带；（2）存在关键NDVI阈值（≈0.45），超过后PM₁₀净化加速，强调需提升密度而非仅扩大面积；（3）O₃响应具有物种依赖性，在高NO_x区域且高BVOC物种占主导时，NDVI–O₃呈正关联，需采用BVOC敏感物种造林；（4）分析框架可推广至全球其他水分受限的恢复区域。这些发现支持基于证据、空间靶向的恢复政策，在最大化空气质量协同效益的同时减轻非预期权衡。