基于LCZ景观格局的各向异性影响量化框架：以长江三角洲NO2变异为例

《International Journal of Applied Earth Observation and Geoinformation》：Quantifying anisotropic impacts of LCZ landscape patterns on NO 2 variabilities in the Yangtze River Delta, China

　　

本研究针对景观—污染物分析中存在的三个关键局限：对各向异性过程的忽视、对任意固定半径缓冲区的依赖以及标准化分类框架的缺失，构建了一个数据驱动的分析框架，整合了高精度LCZ制图、多尺度空间建模与新型形状自适应影响区。将该框架应用于长江三角洲（YRD）城市群，研究人员首先确定了4000 m圆形区域为NO₂建模的最优尺度，揭示了静态景观格局尺度与动态污染物过程尺度之间的明显解耦。在此基础上，自适应方向加权扇形影响区（ADWSIZ）模型提升了解释力（调整后R²：0.464 → 0.506）并降低了空间自相关。结果表明，在各向异性感知分析单元中开展环境研究具有重要价值。研究发现，工业用地与紧凑高层建筑格局与较高NO₂浓度相关，而破碎化开阔高层与稀疏建成区则充当汇。这些见解可为排放分区与通风廊道保护等规划措施提供依据。总体而言，该框架将影响区的概念化从简单的距离缓冲推进至尺度—形状自适应形式，为研究城市形态与空气质量之间的联系提供了尺度与形状自适应框架，而应用于其他区域时需对尺度和方向权重进行本地化校准。

研究背景方面，NO₂作为全球性污染物，是酸雨和光化学烟雾的前体物，长期暴露可导致呼吸系统和心血管疾病。世界卫生组织（WHO）2021年全球空气质量指南设定年均指导值为10 μg/m³，而长江三角洲（YRD）年均浓度达20 μg/m³，为WHO标准的两倍。人类活动已成为现代NO₂污染的主要驱动因素，化石燃料燃烧是主要人为源。NO₂大气寿命短，近源特征使其地面浓度成为排放强度的强指示指标，并呈现显著时空异质性。现有NO₂分析主要依赖地面监测与卫星遥感，前者精确但成本高且分布不均，后者覆盖广但受云雨影响且柱浓度与地面浓度相关性差。基于土地利用/土地覆盖（LUCC）框架的研究虽已揭示NO₂与景观特征的联系，但存在三大缺口：忽视各向异性、缺乏标准化分类、固定半径缓冲区的局限。为此，研究人员以YRD城市群为研究区，提出整合基于局部气候区（LCZ）的分类、景观指标与形状自适应影响区的分析路径，利用Landsat-8 LCZ图和2024年179个国家监测站NO₂数据开展研究。

研究得出以下结论：LCZ分类总体精度达0.91，Kappa系数为0.90；半变异函数分析显示LCZ格局特征尺度约为1098.22 m，且在各方向上未检测到显著差异；多尺度地理加权岭回归（GWRR）分析表明模型性能在4000 m处达到峰值（调整后R²=0.464），且在2500—6000 m区间保持稳健；ADWSIZ模型相较4000 m圆形缓冲区提升了调整后R²（0.464→0.506），降低了RMSE（2.950→2.834 μg/m³）和AICc（551.614→537.214）。这些发现的重要意义在于：揭示了景观格局尺度与污染物响应尺度之间的解耦现象，证明了各向异性感知分析单元在环境研究中的价值，并为排放分区和通风廊道保护等规划措施提供了科学依据。

在技术方法方面，研究人员主要采用了以下关键技术：基于Landsat-8 OLI影像（30 m分辨率）在Google Earth Engine平台上利用随机森林分类器（100个估计器）进行LCZ分类，训练与验证样本各5100个（每类300个），样本通过百度街景影像目视解译获取；采用7个核心景观格局指数（类别层面：PLAND、LPI、ED；景观层面：CONTAG、CONNECT、SHDI、SHEI），利用FRAGSTATS 4.2软件在500—6000 m（间隔500 m）的圆形影响区内计算；采用多尺度GWRR进行最优尺度识别，通过网格搜索交叉验证优化正则化参数λ；提出ADWSIZ方法，基于GWRR导出的PLAND系数计算八个方向的净影响量并动态调整扇形半径，保持总面积恒定。数据来源为中国国家环境监测中心179个监测站2024年年均NO₂浓度。此外，还构建了主导影响椭圆形影响区作为各向异性基线进行对比。

研究结果方面：（1）LCZ分类结果：通过混淆矩阵评估，大多数LCZ类型分类精度超过90%，仅LCZ 8和LCZ D略低于80%，总体精度0.91，Kappa系数0.90；半变异函数分析表明Matern模型拟合最佳，变程为1098.22 m，方向性半变异函数分析显示在6000 m滞后距离内各方向空间变异性无显著差异（p>0.05）。（2）空间影响区优化：多尺度GWRR分析显示模型性能随尺度非线性变化，从500 m到2500 m调整后R²从0.208快速升至0.452，在4000 m处达到峰值（调整后R²=0.464，AICc=551.614），2500—6000 m区间性能保持稳健；ADWSIZ通过方向加权调整扇形半径，相较圆形缓冲区提升了模型性能（调整后R²：0.464→0.506，RMSE：2.950→2.834 μg/m³，AICc：551.614→537.214），而主导影响椭圆基线未优于圆形缓冲区（调整后R²=0.451）。（3）关键景观指数识别NO₂源汇效应：在4000 m圆形影响区下，工业用地（PLAND_10，96.6%正系数）和紧凑/开阔高层区域（ED_1、PLAND_4、PLAND_5）为主要源景观，开阔高层边缘密度（ED_4，86%负系数）为主导汇景观；在ADWSIZ下，源汇机制更加清晰集中，工业用地正系数比例升至99%，ED_1成为最强源指标，ED_4汇效应增强，稀疏建成区（PLAND_9）新识别为有效汇。

讨论部分总结：研究人员从三个维度展开了讨论。第一，关于尺度依赖的LCZ—NO₂响应，约1000 m的LCZ格局尺度与4000 m的NO₂响应尺度之间的差异可通过景观生态学中格局与过程可能在不同特征尺度上运作的原理来解释，前者反映静态城市景观镶嵌体的空间自相关，后者反映周边源汇景观与观测NO₂变异相关联的更广泛邻域范围；大气过程方面，监测站浓度可反映周边排放环境的传输、稀释和混合，4000 m尺度体现了捕获周边LCZ源汇背景与避免过度空间平均之间的经验平衡。第二，关于LCZ源汇效应的方向异质性，虽然圆形缓冲区在YRD区域尺度上足以表征净影响分布（八个方向无显著差异），但站点尺度上方向异质性仍然明显，ADWSIZ通过允许分析单元适应局部源汇景观配置提供了更局地化的表征；与主导影响椭圆的对比进一步证实简化单轴各向异性影响区不足以捕获多方向LCZ—NO₂交互。第三，关于识别出的景观源汇效应的机制解释，工业用地（PLAND_10）作为最强稳定正效应源，归因于工厂、物流和重型交通设施的密集点源和线源排放；高层边缘密度（ED_1）通过增加地表粗糙度和街谷效应阻碍水平通风而成为源；开阔高层破碎化（ED_4）通过增大建筑间距和天空视角因子促进污染物扩散而成为汇；稀疏建成区（PLAND_9）通过低建筑密度、较高植被覆盖和增强沉积表面而成为汇。研究结论部分翻译如下：本研究针对景观—污染物分析中存在的三个关键局限：对各向异性过程的忽视、对任意固定半径缓冲区的依赖以及标准化分类框架的缺失，构建了一个数据驱动的分析框架，整合了高精度LCZ制图、多尺度空间建模与新型形状自适应影响区。将该框架应用于长江三角洲城市群，研究人员首先确定了4000 m圆形区域为NO₂建模的最优尺度，揭示了静态景观格局尺度与动态污染物过程尺度之间的明显解耦。在此基础上，ADWSIZ模型提升了解释力（调整后R²：0.464→0.506）并降低了空间自相关。这些结果证明了在各向异性感知分析单元中开展环境研究的价值。研究显示，工业用地与紧凑高层建筑格局与较高NO₂浓度相关，而破碎化开阔高层与稀疏建成区则充当汇。这些见解可为排放分区与通风廊道保护等规划措施提供依据。总体而言，该框架将影响区的概念化从简单的距离缓冲推进至尺度—形状自适应形式，为研究城市形态与空气质量之间的联系提供了尺度与形状自适应框架，而应用于其他区域时需对尺度和方向权重进行本地化校准。