当前对于高温热浪事件指标的定义缺乏统一标准，不同学者有不同的算法。Russo（2015）提出了一种每日热浪幅度指数（HWMId，the heat wave magnitude index daily）计算方法，适合对全球的热浪强度进行计算和对比分析。该方法也被wiki百科对热浪指数的介绍中采用（wiki）。 该方法作者提供了相关的R语言包和函数，本文对此进行原理和代码的介绍。

放上wiki上的一张图（Zampieri, M., 2016）

HWMId定义为每年最大热浪事件的幅度，包括热浪事件强度、事件开始时间和持续时间。该指数将高温热浪事件定义为连续3天日高温大于热浪事件温度阈值的一段时期。高温热浪事件阈值根据基准期31天窗口内90%分位数的日高温进行计算。 若将1981-2010年作为基准期（基准期最少为30年），若判断某日d是否发生高温热浪，则高温热浪事件阈值计算公式为： A d = ⋃ y = 1981 2010 ⋃ i = d − 15 d + 15 T y , i A_{d}=\bigcup_{y=1981}^{2010}\bigcup_{i=d-15}^{d+15}T_{y,i} Ad​=⋃y=19812010​⋃i=d−15d+15​Ty,i​

U表示数据集的联合， T y , i T_{y,i} Ty,i​是y年中第i天的日高温。该日的高温热浪幅度（Magnitude）为： M d ( T d ) = { T d − T 30 y 25 p T 30 y 75 p − T 30 y 25 p  if  T d > T 30 y 25 p 0  if  T d ⩽ T 30 y 25 p M_{d}(T_{d})=\begin{cases} \frac{T_{d}-T_{30y25p}}{T_{30y75p}-T_{30y25p}} & \text{ if } T_{d}>T_{30y25p} \\ 0 & \text{ if } T_{d}\leqslant T_{30y25p} \end{cases} Md​(Td​)={T30y75p​−T30y25p​Td​−T30y25p​​0​ if Td​>T30y25p​ if Td​⩽T30y25p​​

T d T_{d} Td​是高温热浪事件发生日d的最高气温， T 30 y 25 p T_{30y25p} T30y25p​和 T 30 y 75 p T_{30y75p} T30y75p​分别是基准期31天窗口内日高温时间序列的第25%分位数和第75%分位数。高温热浪事件往往发生数天，则整个事件的幅度为每日幅度之和。

更详细的介绍请参见：Russo（2015）。

HWMId（RDocumentation）计算可使用R包extRemes（cran，RDocumentation）实现。

该指标计算主要通过以下函数实现： hwmid(yTref, Tref, yTemp, Temp) 输入数据分别为：参考数据的开始年份，参考数据，待计算数据的开始年份，待计算数据； 输出的结果分别为：相应年份的热浪指数，持续天数，开始时间；每天的参考温度；75th和25th分位数数值

yTref和yTemp均是具体的年份值，如1970、2000。Tref和Temp均是一维矩阵数据。

仅介绍“相应年份的热浪指数，持续天数，开始时间”三个常用结果，其他可输出的值可在参考文档中自行查看。 通过原理介绍中列出的公式大家也可看出该指标计算的热浪指数强度表示的是气温相对于一定阈值的相对大小，因此没有单位。持续天数和开始时间表示的均是年内最强一次热浪的持续天数和开始时间。