城市骑行者的体力活动获益与PM2.5暴露风险模拟评估：基于162次真实通勤行程分析

《Toxics》：Modelled Physical Activity Benefits and PM2.5 Risks in Two Urban Cyclists

　　

本研究利用波兰格利维采（Gliwice）两名通勤者162次重复骑行行程的既往发表监测数据，评估了体力活动获益与PM_2.5相关死亡风险之间的模拟平衡。每次行程均定义了一个假设的长期模式，即每周重复五天，并与同时段静止参考浓度下休息相同时间进行比较。将观测到的浓度和持续时间与假定通气量及已发表的剂量-反应函数相结合。研究人员将流行病学系数不确定性和替代通气假设传播到净相对危险度（RR_NET）中，并检验了示例性传感器灵敏度和行程特异性盈亏平衡浓度。在中心固定通气假设下，RR_NET中位数为0.840，对应模拟的相对危险度降低16.0%。当流行病学系数变化时，模拟数据集的中位数为0.839（95%不确定性区间0.792–0.892）。替代通气情景也产生了有利的条件平衡。个人浓度变化±10%改变了12次行程的剂量等效增量的符号，但未产生中心RR_NET ≥ 1。在固定其他假设的情况下，模拟的盈亏平衡浓度中位数为196.7 μg/m³。这些估计值并非观察到的健康效应、人群效应或普遍污染阈值。未校准的光学测量和两人样本限制了解释；需要进行更大规模的研究，包括传感器验证和通气量测量。

城市骑行被广泛推广为可持续交通方式和融入日常身体活动的实用途径。既往健康影响评估大多基于背景PM_2.5浓度和假设的骑行暴露比，较少使用真实世界个人暴露测量，因此难以捕捉行程间的巨大变异性。本研究利用Mainka et al.先前发表的波兰格利维采通勤数据集中的自行车子集，旨在模拟评估体力活动获益与PM_2.5相关死亡风险之间的平衡，并定量传播流行病学系数和通气假设的不确定性。研究人员将162次骑行行程中的每一次定义为每周重复五天的假设长期模式，并与同时段静止参考浓度下休息相同时间进行比较。在中心假设下，中位净相对危险度（RR_NET）为0.840，对应模拟的16.0%相对风险降低，表明在观测到的真实骑行条件下，体力活动获益大幅超过PM_2.5相关风险。该结果并非观察到的健康效应或人群水平结论，但为可复现的暴露-风险模拟框架提供了实例。

研究人员将162次骑行行程中的每一次定义为每周重复五天的假设长期模式，并与同时段静止参考浓度下休息相同时间进行比较。个人PM_2.5浓度由Flow 2便携式光学传感器测量，固定参考浓度来自移动实验室BAM1020监测器。采用固定通气率模型（Model A）和两个替代通气场景（B1、B2），结合体力活动剂量-反应函数（0.50次方模型）和PM_2.5死亡剂量-反应函数（Orellano et al.），通过10,000次蒙特卡洛模拟传播流行病学系数和通气量不确定性，并进行了±10%个人浓度传感器敏感性及盈亏平衡浓度分析。该数据集包含78次加热季和84次非加热季行程，样本来自两名通勤者。

**3.1 城市骑行中的真实世界PM_2.5暴露（Real-World PM_2.5 Exposure During Urban Cycling）**

162次骑行行程中，加热季78次，非加热季84次。个人PM_2.5浓度呈右偏分布，加热季中位数4.13 μg/m³（IQR 2.28–12.44），非加热季3.06 μg/m³（IQR 2.11–10.30），均值显著高于中位数，说明少数高暴露行程拉高了平均值。参与者1的加热季中位浓度高于非加热季，而参与者2反之。这些比较仅描述监测记录，不能推断人群或季节效应。

**3.2 从暴露到吸入剂量（From Exposure to Inhaled Dose）**

吸入剂量是浓度、持续时间和通气量的函数。在固定通气率模型A下，加热季平均吸入剂量为23.13 μg/次，非加热季为20.37 μg/次；替代场景B1和B2产生不同估计，例如加热季中位数分别为3.88和6.47 μg/次。剂量为模型计算量，不同通气假设显著影响剂量大小，但属于灵敏度分析而非实测结果。

**3.3 骑行健康获益（Health Benefits of Cycling）**

假设每周骑行5天，观察到的中位骑行时长为33.0分钟/天，对应的体力活动剂量为18.70 MET-h/week，中位体力活动相关相对危险度（RR_PA）为0.836，即相对危险度降低16.4%。标准化每日骑行15–90分钟时，RR_PA从0.886降至0.743，表明体力活动剂量增加带来递减的边际获益。

**3.4 PM_2.5相关健康风险（PM_2.5-Related Health Risk）**

使用Orellano et al.的剂量-反应函数，中位额外等效PM_2.5浓度（ΔC_eq）为0.2 μg/m³，对应的中位PM_2.5相关相对危险度（RR_PM）为1.00182，即风险增加0.182%。加热季中位RR_PM（1.00289）高于非加热季（1.00086）。42次行程的ΔC_eq为负值，表示骑行剂量低于静止参考剂量，但不应解释为PM_2.5的保护效应。采用Tainio et al.历史系数RR=1.07时，中位RR_PM为1.00135。

**3.5 中心净相对危险度（Central Net Relative Risks）**

将体力活动获益与PM_2.5风险相乘，得到中位RR_NET为0.840（IQR 0.815–0.879），对应模拟的净相对死亡风险降低16.0%。所有162次行程在中心假设下RR_NET均<1。加热季中位RR_NET为0.831，非加热季为0.855；参与者1和2分别为0.824和0.861。用历史PM_2.5系数时中位RR_NET约为0.8393，变化很小。

**3.6 联合不确定性、传感器敏感性与盈亏平衡（Joint Uncertainty Sensor Sensitivity and Break-Even）**

在基线传感器因子下，固定通气模型A联合系数不确定性给出的数据集高中位RR_NET为0.8395（95%不确定性区间0.7918–0.8920），B1和B2分别为0.8348和0.8363；在所有模型和传感器场景中，