贝叶斯时空计数模型在伦巴第（Lombardy）对流层臭氧日超标天数中的应用

《Environmetrics》：Bayesian Spatiotemporal Models for Counts of Tropospheric Ozone Daily Exceedances

　　

在意大利北部，空气质量近期已成为主要关注问题。伦巴第地区因地理因素、高人口密度和工业活动的综合影响，被列为欧洲污染最严重区域之一。研究人员聚焦地面臭氧（O3），其受欧盟空气质量阈值监管。在伦巴第，臭氧水平频繁超过这些阈值，引发对健康风险和环境影响的关注，使得控制超标天数至关重要。由于臭氧污染仅在一年中最暖期令人担忧，研究人员仅考虑2000至2023年4月至10月的臭氧超标情况。研究人员关注双变量响应，即每月发生至少一次长期和短期臭氧超标的天数。为研究逐年臭氧超标天数的趋势，研究人员拟合双变量混合效应模型，对各响应采用二项似然（binomial likelihood）。研究人员将地域和气象因子作为固定效应协变量，并加入潜在时空随机效应。特别地，研究人员将各计数数据对应的年份和月份作为随机分组因子。分析表明，逐日至少发生一次长期臭氧暴露阈值超标的概率历年保持不变；另一方面，极端臭氧水平出现的概率历年稳步下降。

研究背景方面，地面臭氧（tropospheric ozone, O₃）由氮氧化物（NO_x）与挥发性有机物（VOCs）在日照下二次生成，具强氧化性，与呼吸道及心血管疾病、生态系统损害密切相关。欧盟指令2008/50/EC设两项阈值：长期暴露以8小时滑动窗口平均浓度120 μg/m³为限，短期警报以小时浓度180 μg/m³为限。伦巴第地处北意，工业、交通、地形与气象耦合致臭氧频发超标，近年区域法令虽削减前体物，但长期暴露阈值超标天数是否改善尚不明晰，且极端峰值与平均超标趋势可能异构，需用可量化时空异质性与年/月结构的模型评估政策成效。该文发表于《Environmetrics》，针对ARPA Lombardia 51个固定站2000–2023年4–10月小时浓度，构建贝叶斯框架回答上述疑问。

关键技术方法上，研究人员以ARPA Lombardia 51站2000–2023年4–10月小时O₃为源，按8小时滑动窗（至少6/8有效）与日有效判定（缺窗<6）聚合成月内超标天数双响应：120 μg/m³滑动窗超标天数与180 μg/m³单小时超标天数，月无效日≤3方入模，少量缺失由imputeTS填补；气象协变量（均温、降水总量、降水时数、最大风速、总太阳辐射）取自Open Meteo并月均，另含海拔及站点类型/周边用地哑变量。建模采用双变量二项混合效应（bivariate binomial mixed-effects model），logit链接下各响应条件独立，线性预测项含标准化连续协变量与分类哑变量固定效应，以及年随机效应（AR(1)双变量高斯先验含跨阈值相关）、月随机效应、站–年时空高斯过程（可分离指数空间核×AR(1)时间核，跨阈值相关ρ近似1）；先验块独立，固定效应N(0,100)，超参半柯西/逆伽马，MCMC由Stan两链各6000预热4000抽样得8000后验；按站保证训练集覆盖切分80%训练、20%测试，以WAIC、LOO-CV、MAE、RMSE、95%区间覆盖率（CVG）选型与验证。

研究结果各节要点如下。数据节（Ozone Air Pollution Data）显示两超标计数呈强正相关与一致夏季峰，120 μg/m³计数无肉眼趋势而180 μg/m³计数显逐年降，气象中温度与太阳辐射经验相关最高、降水负向、风速海拔弱相关；无后期建站偏选（preferential sampling）证据。模型节（A Bivariate Binomial Hierarchical Model）确立月天数Y_ist~Binomial(n_m,p_ist)，p经logit(η)=Xβ+b年+b月+b时空，放弃泊松因有上限31且均值不低，双变量联合随机效应捕获两阈值依存。后验推断节（Posterior Inference）得年效应α_1t（120阈值）历年95%CI含零漂移，α_2t（180阈值）稳降，α_2,2023?α_2,2000后验均值2.357约为α₁差0.375之6倍；固定效应中温度、太阳辐射对两响应皆正且对180更强，风速略负，降水总量与时数符号相反且因共线性净效近零，海拔中度正，用地与站点类型哑变量CI含零不显著；月效应两阈值相似仅4月微异；跨阈值相关ρ_年>0.5、ρ_时空≈1、ρ_月弱，证双变量设定必要；测试集120阈值MAE 2.328/RMSE 3.375/CVG 0.898，180阈值MAE 0.780/RMSE 1.499/CVG 0.976；表1先验敏感与表2随机效应结构比较均支持原文设定最优（LOO 31105.5，WAIC 30971.5）。

讨论节浓缩：研究人员指出双变量二项层次模型兼顾两阈值经验关联与异质趋势，揭示2000–2023年长期暴露（120 μg/m³滑动窗）日超标概率基本恒定仅极微降，而极端单小时（180 μg/m³）超标概率显著下降，提示面向峰值前体物管控（机动车、工业）削峰有效，但背景长期暴露仍需结构性治理；温度与太阳辐射为主驱因子，未来变暖或放大形成条件，气候变量须纳入评估；用地与站类直接效应弱于气象主导。结论译文：本研究借贝叶斯双变量二项层次模型详述伦巴第臭氧超标的时、空与气象决定因素；联合建模120 μg/m³与180 μg/m³阈值月超标天数，在捕获各异趋势同时解释长短期污染强经验关联；后验显示长期阈值日超标概率2000–2023年基本稳定、极端单小时超标概率显著下降，反映削峰政策奏效但长期暴露顽存；温度与太阳辐射为核心驱动，降水风速关系弱而复杂，空间异质适中，地域分类直接效应不强；双变量相关随机效应提升预测精度，后验预测校验稳健；政策上需整合削峰与压低背景水平之策略，且气候变迁预期强化臭氧生成气象条件，主动管控对公共卫生与环境质素至关重要。