该指数是生理模型的替代品
通用热气候指数(UTCI)把四项气象输入——气温、平均辐射温度、风速和相对湿度——换算成一个等效气温,用来描述人的体感。完整计算要运行 Fiala 多节点人体传热与体温调节模型,并耦合一个服装模型;论文把这一组合称为 Fiala-UTCI。它计算量大,且属于专有付费软件,因此实践中几乎没有人直接运行它。
实际业务系统改用一种六阶多项式近似,它拟合 Fiala-UTCI 的输出。论文称之为 Polynomial-UTCI。RayMan、pythermalcomfort、thermofeel 和官方 UTCI FORTRAN 计算器等常用工具内部都依赖它,0.25° 分辨率的 ERA5-HEAT 全球逐时产品也以它为基础。由于这一近似被如此广泛地嵌入,文献常常把它当作完整模型本身。论文的核心主张是,这种混同掩盖了一个影响重大的误差。
误差主要来自近似,而非预报系统
在 Fiala-UTCI 测试数据上,Neural-UTCI 的 RMSE 为 0.36°C(四分位距 0.35–0.38°C),Polynomial-UTCI 为 2.78°C(IQR 2.53–2.99°C),降低了 87%。误差分布与平均值同样重要:Neural-UTCI 的所有误差都落在 ±2°C 以内,而 Polynomial-UTCI 有 16% 的样本超出这一范围;Neural-UTCI 有 90% 的误差在 ±0.5°C 以内,多项式只有 38%。
在多项式已知失效的地方,差距更大。风速超过 17 m/s 时,Polynomial-UTCI 的 RMSE 为 9.03°C,个别误差超过 20°C,该指数在这一区间通常被宣布为超出适用范围。Neural-UTCI 把高风速下的 RMSE 降到 0.89°C,这也是作者称其可连续用于全球、而非局限于某一范围的原因。
论文还重新界定了总业务误差。Fiala-UTCI 模型自身报告的 RMSE 为 1.04°C,以此为基准,Neural-UTCI 的 0.36°C 近似误差约占总业务 UTCI 误差的 19%,而 Polynomial-UTCI 使预期总误差增加 132%。换言之,用神经网络替代多项式后,近似不再是主要误差来源。
类别边界把小误差变成不同的预警
UTCI 按类别使用:固定阈值把指数划分为从极端冷应激、无热应激到极端热应激的若干区间,这些区间用于天气预警和暴露研究。由于分类基于阈值,误差只要越过边界就会改变标签,因此平均精度会低估边界附近的实际影响。
以 Fiala-UTCI 的类别为基准,Neural-UTCI 的误分类率为 1.7%(IQR 1.4–2.0%),Polynomial-UTCI 为 5.3%(IQR 4.8–5.7%),降低了 68%。Neural-UTCI 在 85% 的测试样本中绝对误差更低,并在 35% 的样本中把误差降低超过 1°C。平均改善最大的是极端冷应激类别,为 2.4°C;最小的是无热应激区间,为 0.2°C。热应激类别中,非常强热应激获益最大,但其平均改善只有 0.4°C。
热应激的数字需要论文明确给出的一个限定:测试数据集中没有极端热应激样本,非常强热应激样本也只有 15 个。关于热的更强证据来自训练数据:Polynomial-UTCI 误分类了 21% 的极端热应激样本和 14% 的非常强热应激样本,Neural-UTCI 则分别为 0.6% 和 2%。
网络重现了多项式缺失的湿度与风速效应
仅凭精度无法说明模型的行为是否与被替代的生理过程一致,因此作者用 SHAP 分解比较两种近似对各个输入的响应。两者都把气温列为 UTCI 偏移量最重要的预测因子,但此后出现分歧。
相对湿度在 Polynomial-UTCI 中是重要性最低的预测因子。高温下多项式几乎不随湿度变化,低温下高湿度会压低其输出。Neural-UTCI 把湿度列为仅次于气温的第二重要因子:低于 50% 参考条件时指数下降,高于该值时指数上升,作者称这与 Fiala 模型的热—湿关系一致。这一关系之所以重要,是因为高湿度会限制出汗带来的蒸发散热。
风速呈现类似模式。两种模型都随风速升高而降低 UTCI,但 Neural-UTCI 让这一效应强烈依赖气温——高温下贡献很小,低温下则大得多;Polynomial-UTCI 则在约 20 m/s 处偏离预期关系,作者将其与高风速下的大误差联系起来。SHAP 比较关注响应的方向和形状;论文指出,两种模型的数值大小不可直接比较。
事件层面的热应激天数明显改变
论文把近似差异换算为两次欧洲事件的暴露指标。在 2003 年 6 月至 8 月欧洲热浪期间的罗马,若用 Neural-UTCI 替代 ERA5-HEAT 等业务产品,非常强热应激天数从 15 天增加到 35 天。作者还考察了 2012 年欧洲寒潮,并报告 Neural-UTCI 能可靠地分类极端冷应激,足以支持连续的全球应用。
这些比较是两种近似在同一事件上分类输出的对比,而不是对健康结果的预报。它们的价值在于诊断:说明近似方法的选择会改变越过预警阈值的天数,而这正是热健康警报所依赖的量。论文进一步报告,相对于 Neural-UTCI,Polynomial-UTCI 在 2024 年全球 UTCI 估计中引入了系统性且具有空间结构的误差,从而改变全球热应激暴露估计。
CPU 上速度相当,GPU 上快得多
替代方案只有在适配现有工作流时才有业务价值。对 CPU 上的单个和批量输入,两种方法的速度差异可以忽略。在 GPU 上,Neural-UTCI 计算批量输入比基于 CPU 的多项式实现快 33 倍,具体取决于批大小和 GPU 规格。该网络以开源 Python 包形式发布,定位为 Polynomial-UTCI 的直接替代品。
输入范围在一个方面也更宽:Neural-UTCI 的经验证风速范围为 0.5 至 30.3 m/s,而 Polynomial-UTCI 的有效范围止于 17 m/s。两者在气温(−50°C 至 +50°C)、平均辐射温度(−80°C 至 +120°C)和相对湿度(5% 至 100%)上共享相同限制,这些限制继承自训练数据。
上限来自生理模型及其训练网格
Neural-UTCI 的精度不可能超过 Fiala-UTCI 本身。该模型代表一个体重 73.4 kg、体脂率 14% 的平均个体,不考虑年龄、健康状况、活动水平、适应状态、文化着装差异,也不考虑此前数日热条件的滞后效应。因此已有研究提出对 UTCI 做本地化调整。对一个存在这些盲区的模型做更忠实的近似,并不能消除这些盲区。
训练数据构成第二重上限。用于拟合两种近似的网格是合成数据,因此无法重现真实大气样本中各输入之间的相关性——例如地球系统模型样本中风速与蒸气压、平均辐射温度和气温之间约 −0.4 的相关性。两种近似都在高风速和高相对湿度下表现最差,因为训练网格在这些区间采样最稀疏。测试数据还来自 1971–80 年的控制模拟,因此不含极端热应激样本;作者认为,需要覆盖当今和新兴极端条件的更新参考模拟,才能知道模型在哪里失效。
优先保证极端条件下的精度,而非平均拟合
论文的论点不是多项式不可用,而是它的误差集中在预警最关键的区间。类别层面改善最大的是极端冷应激,以及多项式被官方视为超出适用范围的高风速区间;罗马的对比显示,非常强热应激天数这一越阈计数增加了一倍以上。随着欧洲热应激天数变得更频繁,固定的误分类率也意味着漏报危险天数的绝对数量在增长,这是作者认为这是一项及时方法学改进的理由。
因此,对考虑是否切换的读者来说,相关的划分依据是条件而非平均表现:Neural-UTCI 的优势在高风速、冷极端和类别边界附近最强,而在底层 Fiala 模型或合成训练网格本身最不可信的地方最弱——即高湿度叠加高风速,以及测试集中缺失的热极端。0.36°C 对 2.78°C 的差距和 1.7% 对 5.3% 的误分类差距描述的是对 Fiala-UTCI 的近似保真度,而不是对实际人群健康结果的保真度。
本文引用的文献
A neural network-based Universal Thermal Climate Index for reliable global thermal-stress classification across extreme weather
一个近似 Fiala-UTCI 模型的紧凑神经网络,相对于业务多项式把近似 RMSE 从 2.78°C 降至 0.36°C,把热应力误分类率从 5.3% 降至 1.7%,同时把有效风速范围扩展到 30.3 m/s。
