《Environmental and Sustainability Indicators》：Distinguishing heatwave occurrence from ozone co-occurrence reveals contrasting autumn phenological associations across Chinese ecosystems
编辑推荐：
复合热浪-臭氧事件代表了气候-空气污染暴露的一种重要形式，可能对植被物候产生影响。利用多源大气、气象、土地覆盖和卫星物候数据集，本研究刻画了2014-2023年中国复合热浪-臭氧事件的时空格局，并检验了其与森林和草地秋季物候的关联。研究人员采用线性固定效应面板
复合热浪-臭氧事件代表了气候-空气污染暴露的一种重要形式，可能对植被物候产生影响。利用多源大气、气象、土地覆盖和卫星物候数据集，本研究刻画了2014-2023年中国复合热浪-臭氧事件的时空格局，并检验了其与森林和草地秋季物候的关联。研究人员采用线性固定效应面板框架，在分离热浪总体发生与热浪内臭氧复合发生的前提下，量化了年内像元尺度的生长季结束期（EOS）关联。复合事件主要集中于华北平原和汾渭平原，并在7月和8月最为频繁。在固定热浪总频率的条件下，热浪内臭氧复合发生频率的增高与森林EOS提前相关联（-0.293天·事件-1），而草地对应的关联为正值但不确定（+0.114天·事件-1）。草地与森林的对比差值为+0.407天·事件-1，表明存在明确的生态系统特异性差异。将热浪总频率与臭氧复合发生分离后，也显著改变了复合事件关联的解释，尤其在草地中。在森林中，该关联还随海拔变化，低海拔处负相关更强，并随海拔升高逐渐减弱。这些发现将热浪事件组成识别为复合暴露的一个重要维度，并证明其物候相关性在不同生态系统和环境背景下存在差异。区分热浪总体发生与臭氧复合发生，为评估植被对气候与空气污染复合极端事件的响应提供了一个更具信息量的框架。
论文解读：复合热浪-臭氧事件与秋季物候的关联——基于中国森林与草地的固定效应面板分析
一、研究背景与科学问题
植被物候是生态系统对环境变化响应的敏感指标，在调节生物圈与大气间碳、水、能量的季节交换中发挥关键作用。除气候变暖外，中国陆地生态系统还暴露于多种大气污染中，这些污染可改变植被功能和季节发育。已有卫星研究表明，PM2.5污染与春季返青延迟相关，而地表臭氧升高可加速秋季衰老并缩短植被生长季。热浪和臭氧升高经常共享有利的气象条件，使其同时发生成为复合气候-空气污染暴露的重要形式。然而，现有关于复合极端事件的植被研究多集中于干旱-热浪或大气-土壤干旱组合及其对生产力、抗性和植被损失的影响，大尺度上热浪-臭氧同时发生与秋季物候关联的证据仍然有限。
复合热浪-臭氧事件的生态后果取决于植物生理调节和生态系统特征。臭氧主要通过气孔进入叶片，可诱发氧化损伤影响光合作用和衰老，而气孔导度对温度、大气水分需求和植物水分状态响应强烈。森林和草地表现出不同的干旱响应和冠层交换策略，为期望生态系统特异的复合胁迫响应提供了生理基础。此外，年度复合事件频率的解释存在挑战：复合热浪-臭氧事件增多可能源于热浪总体更频繁、臭氧在热浪中的占比升高，或两者兼有。这些情形代表不同的暴露配置，但可能产生相似的年度复合事件数量变化。因此，区分热浪总体发生与热浪内臭氧复合发生，为评估各自物候关联提供了更清晰的基础。环境异质性可能进一步调节这种关联，海拔尤其相关，因为植被物候表现出系统但空间异质的海拔梯度，地形条件可改变森林对气候极端事件的响应。
基于上述背景，研究人员整合了高分辨率地表臭氧数据（CHAP）、ERA5气象再分析资料、MCD12Q2卫星物候产品以及土地覆盖信息，对2014-2023年中国森林和草地展开分析，研究两个相关问题：在考虑热浪总体发生变化后，热浪事件内臭氧复合发生是否与秋季物候关联；这种关联是否在主要自然生态系统间和环境梯度上存在差异。研究人员假设在给定热浪事件数量下，臭氧复合发生频率越高与EOS提前相关，且关联强度在森林和草地间存在差异。
二、数据与方法概述
本研究采用多源数据集：ERA5提供日最高气温（水平分辨率0.1°×0.1°），用于热浪识别；CHAP提供1 km分辨率日最大8小时平均臭氧（MDA8）浓度；MCD12Q2提供500 m分辨率年尺度生长季结束期（EOS）；CLCD提供土地覆盖类型；GEBCO提供高程数据。热浪定义为日最高气温连续3天超过1983-2012年基线期第90百分位阈值的事件；臭氧胁迫日为MDA8臭氧浓度超过100 μg m
-3的日子。复合事件定义为热浪期间至少出现一天臭氧胁迫。研究构建了包含24,104个森林像元和24,618个草地像元的平衡面板，采用带像元和年份固定效应的线性面板回归，通过分离热浪总频率（TotalHW = 复合事件频率 + 纯热浪事件频率）与臭氧复合发生，量化像元内EOS的年际关联。同时采用200-km空间块自助法（2000次重复）估计系数不确定性，并用广义加性模型和多种敏感性分析检验稳健性。
三、主要研究结果
3.1 复合胁迫的时空格局
2014-2023年间，复合热浪-臭氧事件在中国呈现明显空间异质性，高发区集中于华北平原和汾渭平原，约21%的陆地区域属于较高复合频率等级。纯热浪事件主要分布于华南和塔里木盆地，超过52.55%的区域纯热浪频率接近零。臭氧单独胁迫日广泛分布在中纬度地区，31%以上陆地区域每年经历超过114个臭氧单独胁迫日。复合事件平均持续5.69天，较长事件（>8.2天）零散分布在西南复杂地形区。时间上，复合影响面积在2022年达到十年最大，季节上呈现"错位-重叠"模式：臭氧污染通常在6月初夏达到峰值，热浪在7-8月仲夏至夏末增强，复合胁迫主要集中于7月和8月，与生长季峰值期重叠。
3.2 森林和草地生态系统的像元内关联
固定效应面板回归显示，在联合估计三个暴露项时，热浪总频率（TotalHW）系数在两个生态系统中均为正（森林+0.237天·事件
-1，草地+0.278天·事件
-1），臭氧单独胁迫日的关联较小（森林+0.007天·日
-1，草地-0.009天·日
-1）。关键的是，复合事件频率（CFREQ）与EOS的关联在生态系统间存在明确差异：森林中CFREQ与EOS显著负相关（-0.293天·事件
-1，95%置信区间：-0.448至-0.133），草地系数为正但统计不确定（+0.114天·事件
-1，95%置信区间：-0.063至0.289）。合并交互模型给出的草地减森林对比为+0.407天·事件
-1（95%置信区间：0.168至0.627），表明两生态系统间的CFREQ组成关联存在显著差异。
3.3 固定热浪总频率下复合事件组成的关联
比较两种热浪计数参数化后发现，在分离计数参数化下（包含CFREQ、HeatFreq和OzoneFreq），森林CFREQ系数弱负（-0.056天·事件
-1），草地为正（+0.392天·事件
-1）；而在固定总频率参数化下，森林系数变为-0.293天·事件
-1，草地系数降至+0.114天·事件
-1且不确定。这表明考虑热浪总频率后，森林的负向组成关联得到加强，而草地的正向关联被大幅减弱，说明草地原始关联中相当部分反映了与热浪总体发生的共享变异。
3.4 CFREQ-EOS关联的海拔依赖性
森林中CFREQ×海拔交互项为正（+0.223天·事件
-1·海拔标准差
-1，95%置信区间：0.121至0.325），相当于每升高1000 m CFREQ系数增加约0.242天·事件
-1。因此，森林负向CFREQ-EOS关联在低海拔处最强，随海拔升高逐渐减弱：在森林第5-95百分位海拔范围（155-3258 m）内，条件CFREQ系数从-0.518升至+0.234天·事件
-1。草地虽然交互项点估计为正（+0.199天·事件
-1·海拔标准差
-1），但置信区间包含零，海拔依赖性未明确解析。
四、讨论与结论
研究识别出复合事件高度集中于华北平原和汾渭平原，这与区域臭氧前体物排放密集和热浪期间持续反气旋异常导致的大气停滞有关；华南虽热浪频繁，但季风相关的云量、降水和较强通风抑制了光化学臭氧累积。季节上7-8月的重叠窗口恰逢中国大部分地区生长季峰值，复合暴露的生态意义尤为突出。
森林和草地的对比关联可能源于水分利用调节和冠层气体交换的差异。臭氧的生理剂量取决于气孔吸收而非仅仅环境浓度，在炎热干旱条件下气孔导度降低可限制臭氧吸收，同时改变臭氧沉降。森林中更强烈的负向关联提示，热浪内更高臭氧复合发生可能对晚季冠层活动施加额外限制。海拔梯度整合了温度、降水、辐射、大气水分需求和植被组成等多重环境变异，低海拔处更强的负向关联可能反映了在更温暖或更需水条件下森林秋季物候对臭氧复合热浪的更高敏感性。
研究存在若干局限：年度事件计数压缩了事件时间、持续时间、强度和前期环境条件等信息；浓度阈值法未直接反映气孔臭氧吸收剂量；卫星EOS代表冠层尺度的物候变化，与生理活动停止可能不完全一致。未来可结合通量型臭氧剂量指标、叶绿素荧光和生态碳通量观测，更直接地解析复合暴露向物候变化传递的路径。
研究结论：本研究揭示，2014-2023年中国复合热浪-臭氧事件主要集中于华北平原和汾渭平原，并在7-8月最为频繁。物候分析显示，在固定热浪总频率下，热浪内臭氧复合发生频率的增高与森林EOS提前显著关联，草地关联不确定，且森林中的关联随海拔升高而减弱。热浪事件组成是联系复合气候-空气污染极端事件与植被物候的重要维度；明确区分热浪总体发生与臭氧复合发生，为解释生态系统对复合极端事件的响应及识别物候敏感性最高的区域和环境条件提供了更清晰的框架。
