摘要:
为了解寿光市“数九”时节温度变化规律，基于寿光市1961~2018年“数九”时节逐日气温观测资料，利用现代气候诊断方法，分析了寿光“数九”时节气温变化特征。结果表明，1) 寿光近58年“数九”平均气温、最低气温总体呈均上升趋势。冬季变暖极端最低气温比平均气温贡献率大，“五九”、“七九”贡献相对大，“四九”贡献最小。2) ≤−10.0℃日数“五九”最多达2.3天/年。寒冷时段出现在“三九”至“五九”，“五九”中期最冷，且“五九”期间温度起伏波动大。“六九”开始升温，“七九”平均气温升至0℃以上。3) 寿光“数九”平均气温年代际变暖趋势明显，年代最低值出现在60年代，最高值出现在21世纪初10年。“五九”年代升温最快，“四九”年代升温最慢，两者相差达0.514℃/10a。4) “数九”平均气温1988年发生突变，突变前后相差2.2℃。
Abstract:
The paper aims to reveal the variation regularity of “shu jiu” in Shouguang, based on the daily temperature observation data of Shouguang from 1961 to 2018, and using modern climate diagnosis method, the characteristics of temperature change in Shouguang were analyzed. The results showed that: 1) In recent 58 years, the average temperature and the lowest temperature of “shu jiu” in Shouguang have increased. The contribution rate of winter warming extreme minimum temperature is higher than the average temperature, the contribution rate of “wu jiu” and “qi jiu” is relatively large, and the contribution rate of “si jiu” is the smallest. 2) The maximum number of days at −10.0˚C is 2.3 days/year. The cold period is from “si jiu” to “wu jiu”, and the middle period of “wu jiu” is the coldest, and the temperature fluctuates greatly during “wu jiu”. “Liu jiu” began to heat up, and the average temperature of “qi jiu” “rose above 0˚C. 3) In Shouguang, the trend of decadal warming of the average temperature is obvious, with the lowest value in the 1960s and the highest value in the early 10 years of the 21st century. The “wu jiu” and the “si jiu” years are the fastest warming, with a difference of 0.514˚C/10a. 4) The average temperature of shu jiu changed abruptly in 1988, with a difference of 2.2˚C.
1. 引言
中国农历中“九”是习惯用的民间节气，从冬至算起每九天为一个九，一直数到“九九”八十一天，被称为“数九” [1] [2]。俗语“九九消寒歌”曰：一九二九不出手，三九四九冰上走，五九六九沿河看柳，七九河开，八九燕来，九九加一九，耕牛遍地走！农谚有“冷在三九”，“春打五九尾，要饭跑断腿，春打六九头，吃穿不用愁”。而这些俗语都与工农业生产和人们生活息息相关。因此，研究“数九”时节气候变化特征还是十分必要的。
近年来，随着全球气候变暖 [3]，康丽华等 [4] [5] [6] 从不同角度分析了我国冬季气温特征，于群等 [7] 研究了山东冬季气温变化的气候特征和分型，张磊等 [8] 研究了51年山东临沂市冷暖冬气候特征及其成因分析，徐风霞等 [9] 分析了寿光市1961~2018年负积温变化特征，结果显示冬季负积温总体呈减少、负积温最长持续日数总体呈减少趋势，张姝丽等 [10] [11] [12] 分别研究了北京、邢台、西安“数九”时节气温变化特征。山东寿光市地处黄河中下游，属暖温带季风大陆性气候，四季分明。目前对寿光“数九”时节气温变化特征未见报道。本文利用寿光市观测资料，分析寿光近58年“数九”时节气候变化特征及影响，以期更好的为当地农业结构调整和人们生活提供气象服务。
2. 资料与方法
日最低气温≤−10.0℃指低温严寒 [13]；负积温的多少，在一定程度上反映低温强度和持续时间的综合影响，是低温灾害的指标之一 [14]。本文整理了寿光市1960~2018年“数九”时节逐日平均气温、最低气温、极端最低气温≤−10.0℃日数、“数九”负积温资料。 “数九”≤−10.0℃日数从极端最低气温挑取；“某九”负积温是该时段日平均气温小于0℃的累积值。年代以10位数字为年代数字值，2000~2009年为21世纪初10年，2010~2018年为近9年。
使用线性倾向估计法、累积距平分析变化趋势，采用Mann-Kendall、滑动t-检验法进行突变检验 [15]，使用Excel 2003制图。
3. 结果与分析
“数九”变化特征分析
由表1可见，平均气温“一九”逐渐降低，“四九”略升，“五九”又略降，“六九”开始升温，从“七九”开始平均气温升至0℃以上，平均气温最低出现在“三九”为−2.8℃，平均气温“一九”到“九九”年代均为升温趋势，“七九”年代升温最快，达0.931℃/10a，“四九”年代升温最慢，仅0.275℃/10a；相关系数除“四九”没有通过显著性检验外，“一九”、“二九”“六九”均通过0.05显著性检验，“三九”、“五九”、“七九”、“八九”、“九九”均通过0.01显著性检验；极端最低气温−22.3℃，出现在五九，极端最低气温“一九”到“九九”年代均为升温趋势，“五九”年代升温最快，达0.952℃/10a，“七九”达0.938℃/10a，“四九”年代升温最慢，达0.483℃/10a；相关系数除“四九”通过0.05显著性检验外，其余均通过0.01显著性检验。从气候倾向率看，冬季变暖极端最低气温比平均气温贡献率大，“五九”、“七九”贡献相对大。≤−10.0℃日数从“一九”到“五九”逐渐增多，“五九”最多达2.3天/年，从“六九”后依次减少。
Table 1. Average temperature, extreme minimum temperature, days ≤ −10.0˚C in shu jiu period from 1961 to 2018
表1. 寿光1961~2018年“数九”时节平均气温、极端最低气温、≤−10.0℃日数
注：*表示通过0.1显著性检验，**表示通过0.05显著性检验，***表示通过0.01显著性检验。
由图1可以看出，平均气温曲线“一九”开始下降至“三九”达最低点，“三九”至“五九”上下波动不大，“六九”至“九九”依次上升；极端最低气温曲线 “一九”至“四九”依次下降，“五九” 猛降，“六九”至“九九”依次上升；负积温曲线“一九”至“三九”持续下降，“三九”至“五九”略有上下波动，“六九”至“九九”依次上升。
由图2可以看出，“数九”平均气温与极端最低气温变化趋势线走向基本一致。“一九”第1天气温开始下降至“三九”第6天达低点，平均气温、极端最低气温分别为−3.3℃、−7.8℃，然后呈波浪上升至“五九”第一天，平均气温、极端最低气温分别上升为−2.0℃、−6.4℃，又开始下降至 “五九”第5天达最低点，平均气温、极端最低气温分别为−3.5℃、−8.1℃，后呈上升趋势，至“六九”第8天平均气温上升至0℃以上。
Figure 1. From 1961 to 2018, the average temperature, extreme minimum temperature and negative cumulative temperature of Shujiu in Shouguang
图1. 寿光1961~2018年“数九”平均气温、极端最低气温、负积温变化趋势
Figure 2. The 58-year average daily average temperature and extreme minimum temperature in Shouguang from 1961 to 2018
图2. 寿光1961~2018年“数九”日平均气温、极端最低气温58年平均值
Figure 3. The trend of average temperature and extreme minimum temperature in shu jiu from 1961 to 2018 in 58 years
图3. 寿光1961~2018年“数九”58年平均气温、极端最低气温年代变化趋势
由图3可以看出，寿光平均气温、最低气温年代变化总体呈升高趋势，气候倾向率分别为0.65℃/10a、0.72℃/10a，相关系数r分别为0.65、0.72，均通过了α = 0.001的显著性检验。20世纪60~80年代平均气温、最低气温年代值在平均值以下，进入90年代跳变到平均值以上。60年代平均值最低，21世纪初10年最高。由此可以看出，近58年来，寿光“数九”年代际变暖趋势明显，进入20世纪90年代以来这种趋势尤为明显，21世纪初10年最高，近9年略下降。
由图4“数九”每一天年代平均气温累积距平变化曲线判断年代由冷变暖拐点，可以看出，20世纪60年至21世纪近9年，回暖日期均在“六九”，90年代回暖最早在六九第一天，近9年回暖最晚在六九第9天，年代最早与最晚相差8天。
Figure 4. Chronological average daily temperature anomaly of “shu jiu” from 1961 to 2018 in Shouguang
图4. 寿光1961~2018年“数九”日平均气温年代累积距平
4. “数九”温度突变分析
由图5中UF曲线可见，“数九”平均气温1976年以来有明显的上升趋势，1991年这种上升趋势超过0.05临界线，升温十分显著，根据UF和UB曲线交点，1988年发生突变，突变前后相差2.2℃。
Figure 5. The Mann-Kendall statistical curve of the Mean air temperature in Shouguang from 1961 to 2018
图5. 寿光市1961~2018年“数九”平均气温Mann-Kendall统计量曲线
由图6中UF曲线可见，“数九”最低气温1973年以来有明显的升温趋势，1991年这种升温趋势超过0.05临界线，升温十分显著，根据UF和UB曲线交于0.05临界线上，通过滑动t-检验，没有发生突变。
Figure 6. The Mann-Kendall statistical curve of the minimum temperature in Shouguang from 1961 to 2018
图6. 寿光市1961~2018年“数九”最低气温Mann-Kendall统计量曲线
5. 结论与讨论
1) 寿光近58年“数九”平均气温、最低气温总体呈均上升趋势。冬季变暖极端最低气温比平均气温贡献率大，“五九”、“七九”贡献相对大。≤−10.0℃日数“五九”最多达2.3天/年。寒冷时段出现在“三九”至“五九”，“五九”中期最冷，且“五九”期间温度起伏波动大，与农谚“冷在三九”之说有差别。
2) 寿光“数九”平均气温年代际变暖趋势明显，年代最低值出现在60年代，最高值出现在21世纪初10年。“五九”年代升温最快，“四九”年代升温最慢，两者相差达0.514℃/10a；“六九”开始升温，从“七九”开始平均气温升至0℃以上。
3) “数九”平均气温1988年发生突变，突变前后相差2.2℃。
4) 专家对“数九”变化趋势论证结论与本研究结论基本一致，不一致的地方就是寿光“五九”中期最冷，与农谚“冷在三九”之说有差别。本研究的创新点填补寿光区域“数九”气温统计资料空白；对农谚“冷在三九”进行了验证。
5) 讨论
本研究只从气温方面进行研究，今后工作中，加强寿光气候资源变化全面研究，为当地气候区划和工农业生产提供数据支撑。
基金项目
潍坊市气象局气象科学技术研究项目“潍坊地区冷暖冬气候特征研究”(2019wfqxkt12)。
NOTES
*第一作者。