1. 引言
近年来，以全球变暖为主要特征的气候变化已经成为不争的事实。政府间气候变化委员会第五次综合评估第一组的报告决策者摘要提出，气候变化比原来认识的要更加严重。在全球变暖的背景下，不难发现我国整体的增温趋势与我国所处的北半球情况大致是相似的。但就我国来说，我国的西南地区自50年代以来降温的趋势就一直比较显著，主要的增暖出现在东北、华北和西北等西部地区 [1]。再说我国的其他具体省份，例如青海省，青海省的热量资源就是呈显著增加的趋势，而且大体的增加趋势是由西南向东北逐渐变化的 [2]。
气候变化将会对农业资源产生巨大的影响，而气候变化中的热量资源的时空分布和变化，尤其会对农作物种类、产量以及农作物的种植制度产生深刻的影响。因为只有当植物所处的环境中，气温高于其生长所需的下限气温并且满足一定的积温条件时，植物才能正常生长发育 [3]。一个区域的作物品种类型、种植方式以及栽培方式，灾害预防措施等的确定，热量是其必须考虑的因子之一。而作物产量的高低、品质的好坏在很大程度上受到作物在生育期间所处环境中热量条件的限制。因此，为有效地利用热量资源，合理地分配作物种植的方式与方法，分析热量资源的时间、空间分布对作物产量达到理想的预期与结果有着非常重要的作用。
昆明处于我国的西南地区，属于热带季风气候，四季温暖如春，气候特征可归纳为干湿度分明、气候炎热干燥、气温年较差小、气候垂直变化大和热量资源丰富五个方面。本文就昆明地区≥10℃的热量资源进行了数据处理和分析，以方便找到解决云南地区在农业资源利用方面的问题，能够更好更充分地发挥出气候变化积温增长的优势。
2. 资料与方法
2.1. 资料概况
重点是对昆明地区1960~2019年间的热量资源变化进行的分析，所以本文采用的数据资料是中国气象数据网提供的昆明站1960~2019年逐日平均气温。
2.2. 研究方法
用五日滑动平均法来确定稳定通过某界限气温的起止日期，计算各界限温度对应的积温 [4]。制作一元线性回归方程来描述气候要素的变化趋势 [5]。而对于气候的突变情况，选择使用了M-K检验的方法，此方法也可以用于气温时间序列变化趋势的显著性检验 [6]。
3. 研究结果
3.1. ≥10℃积温年际变化特征
通过1960~2019年昆明地区≥10℃的积温年际变化图(图1)可以看出这60a间昆明地区≥10℃的积温是在不断地增加的，增加变化率为185.5℃/10a，其回归方程的相关系数r = 0.746，对其进行T检验(表1)，其T统计量11.273大于其T临界值2.002，可见该相关系数通过了显著性检验，相关性显著。从图1可见，≥10℃的积温在2019年达到最高，最大值为5645.3℃，在1971年达到最低，最小值为3828.9℃，最大值与最小值相差1816.4℃。
Figure 1. 1960 - 2019 interannual variation trend of accumulated temperature ≥10˚C in Kunming area
图1. 1960~2019年昆明地区≥10℃的积温年际变化趋势
Table 1. Correlation coefficient t test of interannual variation of accumulated temperature in Kunming area ≥10˚C
表1. 昆明地区≥10℃的积温年际变化回归方程相关系数T检验
结合≥10℃的积温M-K突变检验结果(图2)来看，≥10℃的积温在1966~1979年间虽有突破α = 0.05的临界值，但是在1980年最终突破了α = 0.05的临界值并保持上升的趋势。虽然在1992年前后出又现过一次短暂的下降过程，但之后UF的值一直大于0。在1994年正序列曲线与逆序列曲线相交于一点后，其上升趋势变得及其显著，可见≥10℃的积温在1966~1969年间多波动起伏变化，但积温整体变化是有所增加的，在1980年后≥10℃的积温显著上升，特别是在1994年≥10℃的积温有了一个突变性的增加之后，其积温的上升趋势及其显著。突变前平均积温为4385.2℃，突变后平均积温为5004.3℃，突变前后≥10℃的积温增长了619.1℃。
Figure 2. 1960 - 2019 cumulative temperature mutation test at ≥10˚C in Kunming
图2. 1960~2019年昆明地区≥10℃的积温突变检验
3.2. ≥10℃积温的初日变化特征
通过1960~2019年昆明地区≥10℃的积温初日变化图(图3)可以看出这60a间昆明地区≥10℃的积温初日有一定的提前迹象，≥10℃的积温初日最早发生在1999年的1月17日，其次便是2019年的1月22日，最晚发生在1969年的4月7日，其回归方程的相关系数r = 0.440，对其进行T检验(表2)，其T统计量4.478大于其T临界值2.002，可见该相关系数通过了显著性检验。从回归方程的线性变化来看，≥10℃的积温初日是在以4.6 d/10a的速度显著提前的。
Figure 3. The daily trend of accumulated temperature ≥10˚C in Kunming from 1960 to 2019
图3. 1960~2019年昆明地区≥10℃的积温初日变化趋势
Table 2. Correlation coefficient t test of regression equation for accumulated temperature in Kunming area ≥10˚C
表2. 昆明地区≥10℃的积温初日变化回归方程相关系数T检验
Figure 4. 1960 - 2019 detection of the first day mutation of accumulated temperature ≥10˚C in Kunming area
图4. 1960~2019年昆明地区≥10℃的积温初日突变检验
结合M-K突变检验结果(图4)来看，正序列曲线UF与负检验线有多个交点，最后一次发生在2012年，UF彻底突破负临界值，并保持下降趋势，还可以看出有四段时间发生了较大的下降：1960~1966年、1974~1982年、1994~2001年、2005~2019年，正序列曲线与逆序列曲线也有多个交点的存在，最后一次相交是在2005年，可见≥10℃积温初日在21世纪初期很不稳定，并且多次发生突变，在2012年以后≥10℃积温初日出现了显著的提前趋势。
3.3. ≥10℃积温的终日变化特征
通过1960~2019年昆明地区≥10℃的积温终日变化图(图5)可以看出这60a间昆明地区≥10℃的积
Figure 5. 1960 - 2019 the daily trend of accumulated temperature ≥10˚C in Kunming
图5. 1960~2019年昆明地区≥10℃的积温终日变化趋势
温终日有一定的延后推迟的迹象，≥10℃的积温终日最早发生在1992年的10月12日，最晚发生在2018年和2012年，都发生在12月21日，从回归方程的线性变化来看，≥10℃的积温终日结束时间在以3.0 d/10a的速度在向后推后的，其回归方程的相关系数r = 0.396，对其进行T检验(表3)，其T统计量3.881大于其T临界值2.002，可见该相关系数通过了显著性检验，即昆明地区≥10℃的积温终日显著推后。
Table 3. Correlation Coefficient t test of the regression equation for the daily variation of accumulated temperature ≥10˚C in Kunming from 1993 to 2019
表3. 昆明地区1993~2019年间≥10℃的积温终日变化回归方程相关系数T检验
结合M-K突变检验结果(图6)来看，正序列曲线UF在1990~1997年之间与负序列曲线UB存在多个交点，曲线UF在2006年突破α = 0.05的临界值，之后一直持续上升，可见≥10℃积温终日在1990~1997年之间很不稳定，变化幅度较大，多次发生突变，最后从2006年开始≥10℃的积温终日便出现了显著的推后趋势。
Figure 6. 1960 - 2019 detection of a sudden change in accumulated temperature of ≥10˚C in Kunming
图6. 1960~2019年昆明地区≥10℃的积温终日突变检验
3.4. ≥10℃积温的持续日数变化特征
通过1960~2019年昆明地区≥10℃的积温持续日数变化图(图7)可以看出这60a年间昆明地区≥10℃的积温持续日数在不断地增长，增长的波动较大，≥10℃积温持续日数最短为1971年的225 d，最长为2012年的331 d，两者相差106 d，从回归线性方程来看≥10℃的积温持续日数增长速度为7.6 d/10a，其回归方程的相关系数r = 0.558，对其进行T检验(表4)，其T统计量6.392大于其T临界值2.002，可见该相关系数通过了显著性检验，即≥10℃的积温持续日数显著增长。
Figure 7. 1960 - 2019 trend of accumulated temperature persistence days ≥10˚C in Kunming area
图7. 1960~2019年昆明地区≥10℃的积温持续日数变化趋势
Table 4. Correlation Coefficient t test of cumulative temperature duration variation regression equation in Kunming area ≥10˚C
表4. 昆明地区≥10℃的积温持续日数变化回归方程相关系数T检验
结合M-K突变检验结果(图8)来看，其正序列曲线UF与正检验线有6个交点，最早发生再1980年，最后一次发生在1998年，正序列曲线与逆序列曲线也有多次交点，最后一次发生在1995年。可见≥10℃积温持续日数每年变化波动较大，特别是在1980~1998年之间。在1995年≥10℃积温持续日数有一个突变性的增加，在1998年后≥10℃积温持续日数增长趋势显著。
Figure 8. 1960 - 2019 sudden change test of accumulated temperature duration days ≥10˚C in Kunming area
图8. 1960~2019年昆明地区10℃的积温持续日数突变检验
4. 结论
通过对昆明地区1960~2019年的≥10℃积温的年际变化、初日变化、终日变化和持续日数变化特征分析，得到如下结论：≥10℃的积温年际变化均显著增加，增加变化率为185.5℃/10a；在1994年前后发生突变，突变前后≥10℃的积温增长了619.1℃。≥10℃的积温初日在1960~2019年整体上在显著地提前，特别是在2012年以后≥10℃积温初日出现了非常明显的提前趋势。≥10℃的积温终日随着年份显著地推后，特别是在2006年开始≥10℃的积温终日出现了非常明显的推后趋势。≥10℃的积温持续日数随着年份是在显著地增加的；在1995年≥10℃积温持续日数有一个突变性的增加。