本研究由复旦大学公共卫生学院的孟夏、浙江中医药大学公共卫生学院的张丽娜、牛津大学纳菲尔德人口健康系的陈家康（Ka Hung Chan）、北京大学公共卫生学院的吕筠等多位学者共同完成，通讯作者为复旦大学阚海东、牛津大学陈铮鸣和北京大学李立明。论文发表于 eco-environment & health 期刊2025年第4卷，文章编号100148，于2025年4月24日在线发表，属于原创研究论文。

该研究基于中国慢性病前瞻性队列研究（China Kadoorie Biobank，CKB），纳入了来自中国10个地区（5个城市地区和5个农村地区）的512,691名30–79岁成年人，旨在探讨居住环境绿色度（greenness）与心血管疾病（cardiovascular disease，CVD）及其主要亚型发病风险之间的前瞻性关联，并估算若达到世界卫生组织（WHO）推荐的绿色空间水平，城市地区可预防的心血管疾病发病比例。研究背景在于，心血管疾病是全球首位死因，中低收入国家的疾病负担尤为沉重，而绿色植被暴露作为一种潜在的可改变环境因素，其对心血管疾病发病（incidence）的前瞻性证据在中低收入国家仍十分有限。既往研究多关注心血管疾病死亡，且大多在高收入国家开展，样本量或地域代表性有限。因此，该研究利用中国大规模多中心前瞻性队列，填补了中低收入国家绿色度与心血管疾病发病关联的证据空白。

研究的工作流程包括研究对象招募与基线调查、绿色度暴露评估、心血管疾病发病结局随访、统计分析以及健康影响评估五个主要环节。首先，CKB队列在2004–2008年间招募了来自中国10个城乡地区的512,691名参与者，通过笔记本电脑问卷收集人口学、社会经济状况、生活方式及环境暴露信息，并进行身高、体重等体格测量。其次，绿色度暴露采用卫星遥感归一化植被指数（Normalized Difference Vegetation Index，NDVI）进行估计。数据来源于Terra卫星中分辨率成像光谱仪（MODIS）MOD13Q1产品，时间跨度为2005年1月至2017年12月，时间分辨率为16天，空间分辨率为250 m×250 m。研究计算了参与者居住地周边500 m和1000 m缓冲区内的年最大NDVI（NDVImax）值，以代表年度最大植被覆盖水平，同时排除负值以避免水体干扰。第三，心血管疾病发病结局通过与中国国家医疗保险数据库及死因和疾病登记系统的电子链接获得，随访期间为2005–2017年，心血管疾病事件依据国际疾病分类第十版（ICD-10）编码，包括总心血管疾病（I00–I99）、缺血性心脏病（IHD，I20–I25）、急性心肌梗死（AMI，I21）、总卒中（I60–I61、I63–I64和I69）、缺血性卒中（IS，I63）及出血性卒中（HS，I61）。研究仅纳入首次住院事件，随访至发病、死亡、失访或2017年12月31日。第四，统计分析采用时间依赖性Cox比例风险回归模型，NDVImax作为时间变化变量进入模型，调整因素包括地区与诊所分层以及年龄、性别、吸烟、二手烟暴露、饮酒、教育、收入、体质指数、体力活动、水果、蔬菜、鱼类、肉类摄入、烹饪燃料、取暖燃料和自评健康等。暴露-反应关系使用自然样条函数在完全调整模型中分析。研究还进行了按年龄、性别、教育、收入、体力活动、体质指数、自评健康、PM2.5和O3分层分析，并进行PM2.5和O3的中介分析。敏感性分析包括使用年平均值NDVI（NDVIann_mean）替代最大NDVI、排除基线自报心血管疾病患者、排除随访不足两年者、按中国人群超重标准重新定义体质指数分类以及进一步调整PM2.5和O3。第五，健康影响评估部分基于WHO关于普遍可及绿色空间的建议（至少0.5公顷且线性距离300 m以内），建立NDVI与绿色面积百分比（%GA）的广义加性模型，确定各研究地区对应的NDVI目标值（NDVIgoal），然后比较2020年人口加权NDVI水平与该目标值，估算达到WHO推荐水平时可预防的心血管疾病发病比例。

研究结果显示，经过508万人年的随访，共记录148,032例新发心血管疾病事件，其中缺血性心脏病50,323例，急性心肌梗死4,604例，总卒中57,224例，缺血性卒中50,175例，出血性卒中7,685例。全部参与者NDVImax平均值为0.543。研究发现绿色度与总心血管疾病及其亚型发病风险呈显著负相关。具体而言，在500 m缓冲区内，NDVImax每增加0.1，总心血管疾病发病风险降低2.4%（风险比HR=0.976，95%置信区间CI：0.958–0.994），缺血性心脏病风险降低4.3%（HR=0.957，95% CI：0.934–0.980），急性心肌梗死风险降低6.5%（HR=0.935，95% CI：0.895–0.977），总卒中风险降低4.9%（HR=0.951，95% CI：0.933–0.969），缺血性卒中风险降低5.9%（HR=0.941，95% CI：0.922–0.961），但与出血性卒中无显著关联（HR=1.002，95% CI：0.977–1.027）。1000 m缓冲区内的关联略强且一致。暴露-反应曲线显示除出血性卒中外，各结局风险随绿色度增加呈下降趋势。分层分析显示绿色度保护效应在65岁以下人群中显著（HR=0.957，95% CI：0.934–0.981），而在65岁及以上人群中不显著（HR=0.994，95% CI：0.971–1.019）；在PM2.5高浓度地区保护效应显著（HR=0.966，95% CI：0.941–0.993），低浓度地区则不显著；O3分层结果相反。此外，低教育和低收入人群的保护效应显著，高教育和高收入人群则不显著。中介分析表明，O3显著介导了绿色度与心血管疾病及其亚型发病关联的13.1%–52.8%，而PM2.5的中介作用不显著或极小。敏感性分析结果基本稳健。健康影响评估显示，五个城市地区约60%的居民生活在绿色度低于WHO推荐水平的区域，而农村地区仅24%未达标。若五个城市地区达到WHO推荐绿色度水平，平均可预防3.81%（范围2.24%–7.39%）的心血管疾病发病，其中哈尔滨地区健康收益最高（7.39%），海口地区最低（2.24%）。

该研究的结论是，较高的绿色度暴露与中国成年人较低的心血管疾病发病风险显著相关，且这一保护效应在年轻人群、低社会经济地位人群和高空气污染暴露人群中更为明显。这是中低收入国家首个评估绿色度与心血管疾病发病关联并量化健康收益的大规模前瞻性队列研究。其科学价值在于提供了高质量的前瞻性证据，证实绿色空间对心血管疾病一级预防具有潜在作用；应用价值在于为城市规划和公共卫生政策提供了定量依据，表明提升城市绿色空间可达性可在人群层面带来可观的健康收益，特别是在快速城市化地区。研究的亮点包括：样本量大、地域覆盖广泛、随访时间长、结局事件数量多；首次在中低收入国家同时评估总心血管疾病及多种亚型；采用时间依赖性暴露评估和多种敏感性分析增强结果稳健性；结合WHO推荐值进行健康影响评估，将流行病学效应估计转化为政策可操作的预防比例。研究也指出其局限性包括NDVI无法区分绿色空间类型和质量、缺乏参与者实际活动时间信息、未考虑随访期间迁移以及无法精确获得发病时间等。