韩国空气质量监测站迁移引发的PM10与PM2.5时间序列不连续性：基于贝叶斯结构时间序列的分析

《Atmospheric Environment》：Relocation-associated discontinuities in Korean air-quality monitoring: a Bayesian structural time-series analysis of PM 10 and PM 2.5

　　

： 颗粒物（PM）空气污染是导致多种疾病的主要环境风险因素。空气质量监测网络为空气污染监管与流行病学研究提供定量基础，但监测站迁移可能破坏长期时间序列的连续性。本研究检验两个假设：迁移与PM10和PM2.5时间序列不连续相关，且不连续绝对幅度随水平迁移距离增加而增大。研究人员采用贝叶斯中断时间序列分析（Bayesian Interrupted Time Series Analysis, ITSA），以同期未迁移的供体站（donor stations）作为反事实条件，对韩国2010—2023年间93次经核实的监测站迁移事件进行分析。结果显示，在90天距离控制的对数加一（log1p）主模型中，PM10有14次显著下降、14次显著上升、42次不显著；PM2.5有8次下降、11次上升、27次不显著。显著下降与上升事件的PM10中位效应分别为-16.2%和22.1%，PM2.5分别为-18.3%和19.5%。水平迁移距离与效应幅度无显著关联。研究结论表明，监测站迁移可引入双向、站点特异的时间序列不连续，行政标识的保留并不等同于测量连续性的经验证据，亟需标准化可比性证据与版本化元数据。

一、研究背景与问题

颗粒物（Particulate Matter, PM）空气污染已被广泛认定为多种疾病的重要环境风险因素，大量流行病学与实验研究证实，暴露于细颗粒物（PM_2.5）和粗颗粒物（PM₁₀）会增加心血管疾病、肺部疾病、免疫功能障碍、癌症及不良妊娠结局的发病率和死亡率。世界卫生组织（WHO）近期将PM_2.5年均暴露指导值从10 μg/m³收紧至5 μg/m³，但模型显示全球仍有超过90%的人口暴露于该水平以上。

空气质量监测网络为空气污染监管和流行病学研究提供定量基础，通常被视为固定系统。然而，监测站位置可能因城市发展、土地利用变化、人口密度转移或新污染源出现而改变。大型基础设施项目、租约终止或安全问题均可能导致监测站迁移。监测站迁移后，其周边微环境的变化可能使记录到的污染物浓度出现不连续，而时间序列分析——包括趋势估计和政策效应评估——均假设站点位置和周边环境保持稳定。

不同司法管辖区的迁移治理方式各异，但均未建立国际统一的、经实证验证的连续性阈值。美国将迁移视为站点特定网络决策；加拿大不鼓励移动长期趋势站点并要求可比替代数据；欧盟法律强调空间代表性、连续性和文档记录；日本将站点代码连续性判断交由设立机构；中国对城市评估站点规定了距离、功能区和浓度可比性条件。这些制度存在标准异质性和实施随意性的政策问题。

在韩国，城市和路边监测网络由地方政府在国家指导下管理，地方裁量权可适应站点特定约束，但也导致迁移记录和连续性决策的差异。监测站位置选择可能响应先前的浓度、代表性关切、再开发、公众投诉、租约、安全或监管优先事项，形成一种偏好抽样（preferential sampling）——观测位置可随污染过程演变。行政上保留站点标识并不构成测量连续性的经验证据。

二、研究目的与意义

本研究检验两个假设：其一，迁移与PM₁₀和PM_2.5时间序列不连续相关；其二，不连续的绝对幅度随水平迁移距离增加而增大。为应对共同的时间变化，研究人员修改了原始设计，以同期未迁移的供体站（donor stations）作为反事实条件。该研究证明，政府运营的空气质量监测网络中，站点迁移可在颗粒物时间序列中引入时间不连续，这对长期监测数据的国际比较分析构成关键的可比性问题。

三、主要技术方法

研究人员采用多阶段研究流程：数据获取、灰色文献与公共记录审查、数据清洗与校正、信息公开与反馈、重新校准和影响评估。研究时段为2010年1月1日至2023年12月31日，覆盖韩国11种网络类型。每小时PM₁₀和PM_2.5浓度数据来自AirKorea（韩国环境公团运营），汇总为日平均值。每日气象预测变量（温度、降水、相对湿度、风速、亚洲沙尘指标）取自最近的韩国气象厅观测序列。周末和韩国公共假日指标由日历生成。

迁移记录通过年度环境空气质量报告、官方公告、信息公开回应、地方新闻报道、卫星和航空影像以及街景门户进行三角验证，计算原始与新站点坐标之间的大地测量距离。研究人员采用贝叶斯中断时间序列分析（Bayesian Interrupted Time Series Analysis, ITSA），基于状态空间模型，以90天为平衡前后迁移视界，采用对数加一（log1p）变换处理极端浓度。每个事件选取三个最近的同类型、未迁移且有至少80%覆盖率的供体站作为同期预测变量，使用10,000次马尔可夫链蒙特卡洛（MCMC）迭代，包含局部水平和七天季节状态。

四、研究结果

4.1 韩国空气质量监测网络的演变

韩国空气质量监测网络从2010年的289个站点扩展到2023年的661个，增长129%。2016至2021年间增长最快，与公众对颗粒物污染的关注度提高及监测站进气管高度代表性审查相关。记录的迁移事件至2019年增加后下降。

4.2 迁移站点的空间分布

最终数据集包含2010至2023年间的93次迁移事件。6个站点标识有两次记录迁移，因此有三个连续站点版本。事件发生在88个城市、3个路边和2个郊区监测站点。所有二次迁移均为城市站点。15次事件涉及14个长期趋势监测站点。

4.3 迁移距离模式

迁移距离分布高度偏斜，范围为21.80米至4,694.82米。单次迁移平均距离1,063.39米，中位数745.37米；二次迁移平均1,540.50米，中位数1,273.77米。12个单次迁移站点（14.8%）迁移超过2公里，2例超过4公里。

4.4 迁移相关的PM不连续性

90天距离控制的log1p主规范产生了70个PM₁₀和46个PM_2.5成功拟合。PM₁₀结果包括14次显著下降、14次显著上升和42次不显著事件；PM_2.5结果包括8次下降、11次上升和27次不显著事件。显著下降和上升事件的PM₁₀中位效应分别为-16.2%和22.1%，PM_2.5分别为-18.3%和19.5%。距离假设被拒绝：绝对效应幅度与水平迁移距离无显著线性关联（PM₁₀斜率=-1.99个百分点/公里，p=0.300；PM_2.5斜率=-0.08，p=0.998）。按站点类型分析显示PM₁₀效应倾向于正值，而PM_2.5效应更广泛分布于负值。二次迁移倾向于与PM₁₀浓度增加相关。到最近道路距离的变化与估计效应无显著关联。

五、讨论与结论

迁移与双向不连续相关，研究结果不支持迁移普遍降低观测国家PM水平的一般性主张。迁移方向不能识别动机，行政合法性应区别于统计连续性。水平距离不足以作为测量微环境变化的代理指标，短距离迁移可能跨越高排放源或通风梯度（如街道峡谷），而长距离迁移可能保留相似的源混合和代表性。进气管高度、海拔、遮挡、道路或源接近度、土地覆盖和气流的变化可能主导直线位移。

对于流行病学、趋势和政策应用，后果取决于组合记录的使用方式。当分析将迁移前后观测视为未变化目标的测量时，不连续即成为偏倚。不显著估计并非等价性的肯定证据。连续性决策需要前瞻性定义的容差和为此设计的证据。

研究人员提出韩国应采纳国家最低标准：默认采用版本化站点标识；在可行时进行代表性条件下的校准重叠监测；对趋势站点进行中央审查和迁移后审计；公开发布连续性决策和支持数据。研究提出三个连续性决策标准：当行政和站点可比性要求满足且效应落在预先建立的等效区间内的后验概率高时支持连续性；当效应超过该区间且后验概率高并在合理规范下持续时支持不连续性；当数据不完整、不确定性宽或合理选择改变结论时，连续性不确定，不确定的站点版本应保持单独可识别。

研究存在三个局限性：迁移时间和目的地非随机化，不能排除与迁移同时发生的局部事件或行政动机；历史元数据未一致包含进气管高度、遮挡、局部源强度、交通量或微尺度气象；分析仅覆盖2010至2023年记录的PM₁₀和PM_2.5迁移事件，其他污染物可能有不同响应。

研究结论： 站点迁移可在PM时间序列中引入双向、站点特异的不连续，但方向和分类取决于视界和反事实定义。在主要90天距离控制log1p分析中，大多数可分析事件不显著，显著上升和下降均发生，水平迁移距离不能解释效应幅度。同期供体站通过代表共同时间条件强化反事实，但有意迁移仍非外生，未测量的局部变化仍可能存在。韩国应采纳迁移证据、版本化标识、完整站点元数据、可行时重叠监测和前瞻性使用特定等效决策的国家最低标准，以防止行政标识连续性被误认为已证明的测量连续性。