《Frontiers in Plant Science》：Optimizing the structure of urban woody plantings for a warming climate: a frequency-weighted exposure framework for woody species in Chongqing, China
编辑推荐：
随着气候变暖加剧，城市绿地系统中的木本植物日益受热胁迫（heat stress）。既有城市木本物种气候适应评估多聚焦物种水平热耐受性，而系统水平暴露还取决于各物种在现存种植结构中的使用广度。研究人员以中国亚热带内陆特大城市重庆为案例，提出一种面向城市木本物种的
随着气候变暖加剧，城市绿地系统中的木本植物日益受热胁迫（heat stress）。既有城市木本物种气候适应评估多聚焦物种水平热耐受性，而系统水平暴露还取决于各物种在现存种植结构中的使用广度。研究人员以中国亚热带内陆特大城市重庆为案例，提出一种面向城市木本物种的频权热胁迫暴露（frequency-weighted heat-stress exposure, HSE）框架。通过整合物种特异性热生态位阈值与样地水平出现频率，构建频权HSE指数，量化未来气候条件下单个木本物种对系统水平热暴露的贡献。基于重庆310个样地257个木本物种调查数据，研究发现当前被调查木本植物系统的种植结构中HSE高度集中于少数物种：HSE值最高的前10个物种占总系统HSE近30%，暴露集中度的赫芬达尔-赫希曼指数（Herfindahl–Hirschman Index, HHI）为0.01525。为此研究人员提出在维持各功能群（常绿乔木、落叶乔木等）物种比例不变前提下优化物种使用模式的策略，具体为降低高HSE物种的赋权出现频率，并将释放的频率解析再分配至同功能群内低缺口物种。默认优化情景下，优化配置使总系统HSE降低2.02%、Top10_share降低3.26%、HHI降低0.00254，显示对少数高频高HSE物种的依赖减弱。本研究为物种筛选、苗木采购与更新策略提供量化依据，支撑城市绿地木本植物系统的气候响应管理。
研究背景方面，气候变暖使城市绿地木本植物面临增强的热胁迫，系统水平暴露由物种自身特性与样地水平出现频率共同决定。既有研究多基于热耐受阈值、生态位边界或未来适生性变化识别高风险物种，缺乏将物种水平气候敏感性与种植结构使用广度结合的量化诊断，难以判断当前配置是否将未来热胁迫暴露集中于少数高频物种，亦缺少在维持功能群结构前提下的优化路径。为此研究人员选取屡遭热浪的亚热带内陆特大城市重庆，整合木本植物调查数据、物种特异性热生态位阈值与未来气候情景，建立频权HSE框架并开展结构保持型优化。
主要关键技术方法包括：在重庆主城310个400米2样地调查257个木本物种并记录样地水平出现频率；依据中国种子植物属区系类型系统划分广布、热带、北温带、古地中海与东亚5类起源；取WorldClim 2.1历史气候与SSP370情景2050年预估最暖月最高温（maximum temperature of the warmest month, MTWM）第95百分位定义热生态位阈值T95，最干季降水（precipitation of the driest quarter, PDQ）第5百分位定义干旱生态位阈值P5；以Deficit_MTWM=max(0, MTWM_future?T95)表征热暴露缺口，HSE=freq×Deficit_MTWM构建频权暴露；以Total HSE、Top10_share与HHI刻画集中度；采用Kruskal–Wallis检验起源间阈值差异；优化采取“削峰+安全再分配”，在功能群总频率不变约束下对前20高HSE供体物种按强度α削频，释放频率按1/(Deficit_MTWM+ε)逆权重再分配至同群低缺口候选池（由Deficit_MTWM的q分位界定），设保守（α=0.10,q=0.30）、默认（α=0.20,q=0.20）、上限（α=0.40,q=0.10）三情景及α×q网格搜索。
研究结果如下。3.1节“Niche thresholds did not differ significantly among species from different regions”中，Kruskal–Wallis检验显示T95（H=2.193,df=4,p=0.700,ε2<0.001）与P5（H=5.232,df=4,p=0.264,ε2=0.004）在五类起源间无显著差异，Dunn事后检验均无显著两两差异，表明起源背景不能单独预测未来热风险，系统暴露更依赖物种缺口与使用结构。3.2节“Frequency weighting reveals concentrated HSE and changes species prioritization”中，基线Total HSE=8298.98，Top10_share=0.2992，HHI=0.01525；Deficit_MTWM与HSE排序Spearman ρ=0.226（p<0.001），两top20仅5%重叠（Lagerstroemia indica重合），高频物种在HSE下排名跃升，证明频权改变优先序并揭示暴露集中于少数高频高缺口种。3.3节“Structure-preserving optimization redistributes exposure under functional-group constraints”中，三情景功能群总频不变；保守情景移频4.53%，Total HSE降0.92%、Top10_share降2.32%、HHI降0.175%；默认情景移频9.06%，Total HSE降2.02%、Top10_share降3.26%、HHI降0.25%；上限情景移频18.12%，Total HSE降8.49%、Top10_share降8.87%、HHI降0.51%；α越大削减越强，高α下更小q（更严候选池）效果更强，累积暴露曲线初始段下移，暴露由窄集向宽分布再分配。
讨论部分指出，起源分类不解释阈值变异，系统HSE由缺口与使用该结构放大；仅知Deficit_MTWM不足，须结合样地频率贡献；频权HSE是系统内暴露贡献指标而非独立生理阈值，大缺口低频种需监测、中缺口高频种宜逐步替换。削峰再分配属受功能群与苗木供给约束的渐进情景分析，默认情景Total HSE仅降2.02%但Top10_share与HHI同步下降，意义在于削弱少数高频高HSE种的结构依赖，利于年度更新与采购嵌入。局限包括MTWM超限仅表征潜在暴露非死亡率、T95/P5为现实生态位非生理极限、频率代用广度非生物量或盖度、单情景SSP370单城市；跨城应用宜用相对频率与Top10_share/HHI等份额指标。
结论部分翻译：研究人员建立了一种频权热胁迫暴露框架，将物种水平热生态位缺口与城市木本植物系统样地水平出现频率关联。广义起源类别未显现热或旱生态位阈值的明确梯度，而未来系统水平热暴露更紧密关联物种特异性热暴露缺口与物种使用结构的组合。通过减少高HSE物种使用并将释放频率再分配至同功能群低缺口物种，所提优化途径在维持城市种植广泛功能结构的同时降低了暴露集中度。默认情景下Total HSE降幅有限，但暴露集中度同步下降表明组内渐进调整可在不改变种植广泛功能结构前提下减少对少数高HSE物种的依赖。该框架应作为识别优先物种及指导阶段性物种筛选、苗木采购与绿地更新的决策支持工具，而非树木死亡、冠层损失或生态系统服务衰退的直接预测。论文发表于《Frontiers in Plant Science》。
