1960−2015年淮河流域初终霜日时空变化分析1

马尚谦，张 勃**，唐 敏，马 彬，杨 梅，梁晶晶，张佳琦，加力戈



1960−2015年淮河流域初终霜日时空变化分析1

马尚谦1，张 勃1**，唐 敏1，马 彬1，杨 梅2，梁晶晶1，张佳琦1，加力戈1

（1. 西北师范大学地理与环境科学学院，兰州 730070；2. 北京师范大学社会发展与公共政策学院，北京 100000）

基于淮河流域1960−2015年61个气象站点地面0cm日最低气温资料，采用线性倾向估计、反距离加权、Mann-Kendall突变检验、滑动T检验方法，分析近56a淮河流域初霜日、终霜日和无霜期的时空变化特征及突变年份。结果显示：（1）研究期内，淮河流域平均初霜日期、终霜日期和无霜期分别以2.15、2.49、4.38d·10a−1的速率呈推迟、提前、延长的趋势（P＜0.01），其中，在20世纪90年代的变化最为显著（P＜0.01），速率分别为16.38、5.34、20.6d·10a−1。（2）平均初霜日期在空间上呈北早、南迟，山区早、平原迟的分布；86.9%的站点初霜日期呈显著推迟趋势（P＜0.05）。终霜日期呈西南早、东北迟，平原早、山区迟的分布；83.6%的站点通过0.05水平的显著性检验，以3.44～5.92d·10a−1的速率呈提早趋势。无霜期随纬度和海拔升高而缩短；93.4%的站点通过0.05水平的显著性检验，变化率为3.56～7.59d·10a−1，无霜期整体延长。（3）11月8日线、4月1日线、220d等值线位置较其它各气候基准期和各年代分别偏北约1个和2个纬距，在32°N和34°N附近的偏北趋势最为明显，佐证了该区初霜日期整体推迟、终霜日期整体提前、无霜期整体延长的趋势。（4）初霜日期、终霜日期和无霜期分别在2002年、1995年和1998年发生突变。

淮河流域；初霜日；终霜日；无霜期

霜冻作为一种用温度来表征的农业气象灾害，历来是气象学科、农业学科及其它相关学科研究的重要课题[1]。霜日以地面最低温度降到0℃为指标，霜日是否导致霜冻害，与作物品种、发育期等因素有关。作为衡量温度变化的重要指标，分析初、终霜日期及无霜期的时空变化趋势有利于减少极端低温对农业生产带来的危害，也有利于因地制宜地选择作物品种，及时调整种植制度等，目前已有较多研究在初、终霜日出现早晚、无霜期长短等方面取得重要进展。

Erlat等[2]以日最低温度≤0℃为基础温度，分析发现1950−2010年在土耳其大部分地区霜日均有显著的减少趋势。Terando等[3]观察分析表明，1981−2010年北美中部部分平原区霜冻日数表现出增加趋势，其余地区则表现为显著的减少趋势。宁晓菊等[4]研究发现，近60a中国初霜日期推迟、终霜日期提前和无霜期延长的区域占整个中国的4/5，北方的变化幅度大于南方，东部大于西部，20世纪90年代无霜期发生突变。杨晓玲等[5]揭示了近54a河西走廊东部地区的初霜日显著推迟，终霜日显著提前，无霜期显著延长。张山清等[6]以日最低气温≤0℃为霜冻指标研究发现，近50a新疆地区初霜日以2.2d·10a−1的速率推迟，终霜日以1.4d·10a−1的速率提前，无霜期以3.59d·10a−1的速率延长。杜军等[7]证实，20世纪90年代之后西藏主要农业区部分站点初霜日推迟，终霜日提前，无霜期延长，所有站点的霜冻日数均表现为显著的减小趋势。李芬等发现[8]，山西正常出现初霜冻的概率为62%～82%，平均初霜日发生明显突变的年份为1989年。张霞等[9]指出，太原市平均初霜日为10月14日，平均终霜日为4月12日，无霜期为182d，重度终霜日提前趋势最为明显。学者对于异常霜冻日期的研究也较为深入[10]。不同区域尺度的研究结果，总体上表现为中国多数地区初霜日期推迟、终霜日期提前、无霜期延长。

目前，针对淮河流域的霜冻日研究大多集中在局部地区，只有部分研究对无霜期突变进行了分析，将初、终霜日与无霜期的变化结合起来分析三者的变化趋势及突变的研究报道尚且不多。淮河流域地处中国南北气候过渡带，属暖温带半湿润季风气候区，其特点是冬春干旱少雨，夏秋闷热多雨，冷暖和旱涝转变急剧，是气候变化的脆弱区与敏感区，是中国冬小麦的主产区。在目前农业研究中大多以最低气温作为农作物霜冻指标。研究表明，以地面0cm日最低温度作为基础温度得出的初霜日期与实测数据更为接近。本研究以地面0cm日最低温度作为基础数据，参照《作物霜冻害等级》标准和中国科学技术蓝皮书第5号《气候》[11]，基于不同基准气候期和不同年代，探索淮河流域初、终霜日、无霜期的变化规律，旨在通过系统分析该区域近56a来的霜日变化特征，为合理利用农业气候资源和有效防御霜冻提供一定的科学依据。

1 资料与方法

1.1 研究区概况

淮河流域位于111°55'−121°25'E、30°55'−36°36'N，面积2.7×105km2。西起桐柏山、伏牛山，东临黄海，南以大别山、江淮丘陵、通扬运河及如泰运河南堤与长江分界，北以黄河南堤和泰山为界而与黄河流域毗邻。淮河干流发源于桐柏山太白顶北麓，依次流经河南省、湖北省、安徽省和江苏省。年平均气温11～16℃，由北向南、由沿海向内陆递增，最高月平均气温25℃左右，出现在7月；最低月平均气温在0℃，出现在1月。极端最高气温达44.5℃，极端最低气温达−24.1℃。蒸发量南小北大，年平均水面蒸发量为900～1500mm，无霜期200～240d。多年平均降水量约920mm，由南向北递减，山区多于平原，沿海大于内陆。

1.2 霜日指标计算

根据韩荣青等[12]的分析，北方地区初霜日期基本在8月1日后，因此，以8月1日作为统计起始点，将8月1日定义为1，用日序表示初、终霜日，如1960年9月20日为初霜日期，则该年初霜日就以51计，依次类推，建立所有站点初、终霜日的数据序列。初霜日期为8月1日以后第一次出现地面0cm日最低温度≤0℃的日期，终霜日期为5月31日之前最后一次出现地面0cm日最低温度≤0℃的日期，无霜期为一年中终霜日至初霜日之间的日数。

1.3 数据来源

选用1960年1月1日−2015年5月31日淮河流域内41个站点以及周边20个站点共61个站点（安徽15个、河南15个、湖北5个、江苏18个、山东8个）的地面0cm日最低温度数据，数据来自中国国家气象局气象信息中心（http://www.nmic.gov.cn/）的“中国地面气候资料日值数据V3.0”。数据经过严格的质量控制，主要包括极值控制、缺测站点时间一致性检验。通过查阅气象台站简史，61个站点位置均无较大变化，7个缺测站点插补后与相邻站点数据序列间的一致性指标（决定系数R2）介于0.89～0.97（表1），说明插补之后的数据质量和完整性均较好。

表1 缺测站点与相邻站点地面0cm日最低温度数据序列的一致性检验结果（1960−01−01—2015−05−31）

1.4 突变分析

采用Mann-Kendall突变分析法和滑动T检验法[13]，分析近56a霜期的突变状况。滑动T检验是通过考察两组样本平均值的差异是否显著来检验突变。该方法把同一气候序列中的两段子序列均值的显著差异性看作来自两个总体均值显著差异性的问题来检验。如果两段子序列的均值差异超过了一定的显著水平，可以认为均值发生了突变。对于具有n个样本量的时间序列x，设置某一时刻为基准点，连续设置基准点是从第10年直至倒数第10年，基准点前后两端子序列x1和x2的样本分别为n1和n2，两端子序列平均值分别为`x1和`x2，方差分别为s12和s22，则

其中

式中，t为统计量，s为序列方差。本文滑动T检验选取两个子序列长度为n1=n2=10a，显著性水平α=0.01，按t分布自由度v=n1+n2−2=18，t0.01=±2.23。

1.5 气候基准期

20世纪80年代中后期全球气候态的转变影响中欧湖泊水质和生态系统[14-15]，导致北半球海面温度升高[16]，引起大气温度整体上升[17]，甚至对地下水也造成了重大影响[18]。因此，在气候变化研究中，要更加注重对比气候态转型前后发生气象灾害事件规律的不同点。以往研究大多使用多年的气候平均状态，容易忽视气候转型前后气候所发生的变化。澳大利亚气候监测局主任Karl Braganza表示，由于长期气候变暖趋势，1961−1990年的气候正常值对于目前的气候条件来说不具有代表性。澳大利亚气象局已经在使用1981−2010年新的参考标准值。鉴于此，以30a为一个时间段，选取具有代表性的气候基准期进行研究，分别为1961−1990、1971−2000、1981−2010，90年代之前的年代变化参照1961−1990基准期，90年代之后参照1981−2000年基准期。

2 结果与分析

2.1 淮河流域初霜日的时空变化特征

2.1.1 年际及年代际变化

由图1可见，近56a来全区平均初霜日期发生在10月22日−11月24日，最早在1981年，最迟在2011年；从各年代平均值来看，20世纪70年代初霜发生日最早，在11月4日，21世纪00年代初霜发生日最晚，在11月13日。整个分析期内，全区平均初霜日期以2.15d·10a−1（P＜0.01）的速率呈极显著推迟趋势。20世纪60、70、80、90年代和21世纪00年代的初霜日平均值分别为11月7、6、8、8和5日。多年平均初霜日期具有明显的年代际特征，整体以“V”型分布。90年代初霜日期以16.38d·10a−1（P＜0.05）的速率呈较为明显的推迟趋势，其它年代均没有经过显著性检验，表明淮河流域农作物遭受初霜冻危害的可能性将降低[19]。

图1 淮河流域平均初霜日期的年际及年代际变化

注：r90代表20世纪90年代。

Note：r90stands for simple correlation coefficient of 1990s.

2.1.2 空间变化

依据所有站点（61个）1961−2015年平均初霜日，采用反距离加权模型进行插值，绘制淮河流域初霜日及其线性变化趋势的空间分布，结果见图2。由图2a可见，整个区域内各地初霜发生日期有明显差异，北部沂源、沂蒙山、西部伏牛山受高海拔影响，初霜发生最早，为10月23日，淮河流域东北部的莒县初霜发生在10月下旬，淮河流域中部初霜发生在11月上、中旬，而淮河流域南部固始、信阳初霜发生最晚，在11月24日。最早与最晚之间相差32d。近56a该区域内初霜日期呈“阶梯式”分布，北早、南迟，平原迟、山区早，究其原因主要是受地形、地理位置和冷空气路径的影响，淮河流域大部以平原为主，北方冷空气直驱南下，北部初霜冻必早于南部。由图2b可见，近56a淮河流域及其周边共61个站点中，86.9%的站点其初霜日期均呈显著推迟趋势（P＜0.05），主要分布在34°N带附近。进一步基于不同气候基准期、不同年代数据，选取与农业生产关系密切的多年平均初霜日（11月8日）分布线进行对比，结果见图3。由图3a可见，以1961−1990年为气候基准期计算的11月8日线与以1971−2000年为气候基准期计算的11月8日线基本一致，均沿33°N纬线，即驻马店—阜阳—蒙城—蚌埠—盱眙—阜宁分布，而以1981−2010 年为气候基准期计算的11月8日线较前两个基准期均偏北约1个纬距，到达34°N 线，即宝丰—西华—睢宁—沭阳—灌云一线。图3b显示，20世纪60、70、80和90年代的11月8日等值线在33°N线附近，即驻马店—阜阳—蚌埠—盱眙—射阳一线，2000年以后11月8日等值线沿35°N即定陶—郯城—赣榆—日照分布，比20世纪4个年代偏北2个纬距。可见，气候基准期不同计算区域初霜日期存在明显区别，而11月8日等值线随气候基准期和年代的更替逐渐北移说明该区初霜日期整体有推迟的趋势。

图2 淮河流域平均初霜日期（a）及其线性变化趋势显著站点（b）的空间分布（1961−2015年）

图3 基于不同基准气候期（a）和不同年代（b）计算的淮河流域初霜日期为11月8日地带的位置变化

2.1.3 突变分析

由Mann-Kendall突变分析（图4a）可知，淮河流域地区初霜日期UF和UB曲线交于信度线±1.96（α=0.05）之间，交点为2002年。而根据图4b可知，该区初霜日期在自1970年以来，t统计量有两处超过信度线±2.23（α=0.01），1996−2005年可能存在突变点的年份为1998年、2002年、2004年。结合图4a和图4b认为，淮河流域的初霜日期在2002年发生突变。

2.2 淮河流域终霜日的时空变化特征

2.2.1 年际及年代际变化

图5显示了淮河流域近56a平均终霜日期的年际及年代际变化，平均终霜日期在3月11日−4月15日，最早发生在2007年，最迟发生在1961年，平均终霜日期以2.49d·a−1（P＜0.01）的速率呈极显著提早趋势。从年代际上来看，整体上呈“W”型分布。20世纪60、70、80和90年代终霜日的平均值分别为4月3日、3月31日、3月29日和3月21日。90年代以5.34d·10a−1(P＜0.01)的速率呈极显著提早趋势，但21世纪00年代以11.03d·10−1(P＜0.01)的速率呈极显著推迟趋势，其它年代变化趋势均未通过显著性检验。终霜冻日越早，农作物遭受终晚霜冻害的可能性越低。

2.2.2 空间变化

图6为淮河流域平均终霜日期及变化趋势的空间分布。图6a显示，该区域终霜日期在3月18 日−4月15日，淮河流域西南部终霜冻最早结束，伏牛山、沂蒙山、沂源以及日照附近最迟结束。在空间分布上终霜日呈现为自西南向东北逐渐推迟，平原结束早、山区结束迟的规律。图6b为终霜日期变化趋势的空间分布，83.6%的站点通过P＜0.05水平的显著性检验，均以3.44～5.92d.10a−1(P<0.05)的速率呈显著提早趋势，且32°N和34°N附近的站点提早趋势较明显。

同样选取不同气候基准期和不同年代际的多年平均终霜日4月1日等值线来分析该区域终霜日期的空间变化。图7a为不同气候基准期下的4月1日等值线变化，由图可见，1981−2010年4月1日变化线较1971−2000年和1961−1990年均偏北，分别约偏北0.5°和1°，并沿开封—永城—徐州—睢宁—阜宁一线分布。1981−2010年和1961−1990年的4月1日变化线基本沿1971−2000年变化线呈对称分布，即沿宝丰—西华—蒙城—盱眙一线呈对称分布。就年代际而言，20世纪60、70、80年代3个年代的4月1日线呈波动式往东北方向移动约1°，并由驻马店—阜阳—寿县移动至郑州—亳州—大丰一线，而根据Reid等[14]的研究发现，20世纪80年代后期全球出现规模性的气候转变，导致气候稳定性变差，这可能是导致80年代较70年代相对偏西南的原因。而21世纪00年代较20世纪4个年代均偏北，其中较60年代偏北约2°，分布在兖州—邳县—灌云一线。4月1日的终霜日等值线逐渐向北移动，这也佐证了该区域终霜日期提前的趋势。

图4 淮河流域初霜日期的突变分析

图5 淮河流域平均终霜日期年际及年代际变化

注：r90和r00分别代表20世纪90年代和21世纪00年代。下同。

Note：r90, r00stands for simple correlation coefficient of 1990s and 2000s, respectively. The same as below.

图6 淮河流域平均终霜日期（a）及其线性变化趋势显著站点（b）的空间分布（1961−2015年）

图7 不同基准气候期（a）和不同年代（b）计算的终霜日期4月1日地带的位置变化

2.2.3 突变分析

图8为利用M-K突变检验以及滑动T检验对终霜日期进行的突变分析。由图8a可知，UF线与UB 线在1995年存在一个明显的交点，且通过0.05的显著性水平检验。图8b显示了滑动统计量的变化情况，图中有2处超过0.05显著性水平，即1995、2002年均为突变点。因此，近56a淮河流域终霜日期在1995年发生突变。

图8 淮河流域终霜日期的突变分析

2.3 无霜期的变化特征

2.3.1 年际及年代际变化

分析图9可知，淮河流域近56a年均无霜期在198～247d，平均220d，无霜期最长年份在2008年，最短发生在1962年。平均无霜期整体呈明显的延长趋势，速率达到4.38d·10a−1（P＜0.01）。年代际上，整体呈现“阶梯”上升型分布，各年代平均值分别为218、218、221、222、237d。20世纪90年代的上升速率达到最大，为20.6d·10a−1(P＜0.01)，21世纪00年代以3.69d·10a−1(P＜0.05)的速率呈上升趋势。

图9 淮河流域平均无霜期年际及年代际变化

2.3.2 空间变化

图10a为近55a淮河流域平均无霜期的空间分布，由图可见，各地平均无霜期在201～258d，最高值出现在淮河流域南部信阳、固始，最低值出现在伏牛山区和沂蒙山区以及沂源、莒县，整体分布上仍然呈现自西南向东北逐渐减少的分布规律，平原长，山区短。图10b给出了61个站点无霜期变化趋势分布，由图可见，该区域61个站点中93.4%的站点均通过P<0.05水平的显著性检验，变化率为3.56～7.59d·10a−1，且在32°N和34°N纬度带无霜期延长趋势最大。

淮河流域平均无霜期为220d，以其位置的变化来分析近55a来无霜期的不同时段空间变化状况。图11a为不同气候基准期下的220d等值线的变化，1961−1990年和1971−2000年的220d等值线基本保持一致，并沿宝丰—阜阳—盱眙一线分布。而1981−2010年的220d等值线南移至33.5°附近，偏南了0.5°，分布在南阳—驻马店—阜阳—蚌埠—大丰一线。就年代际而言，220d等值线尤以70年代的变化最为明显。60、80、90年代基本在33°N附近，沿南阳—驻马店—阜阳—蚌埠—盱眙—射阳一线波动式变化；而21世纪00年代北移至34.5°N带一线，较前3个年代北移了近1.5个纬距，分布在宝丰—蒙城—阜宁一线。无霜期220d等值线空间位置、初霜日期11月8日等值线、终霜日期4月1日等值线均是在气候突变之后的21世纪00年代发生较大变化，无霜期的变化受初霜日和终霜日变化的共同影响[4]。

2.3.3 突变分析

如图12a所示，淮河流域近56a平均无霜期UF和UB曲线在1998年有明显交点，并通过α=0.05水平显著性检验。图12b显示了滑动统计量的变化，由图可知，平均无霜期在1994−2005年均存在突变点。结合两种方法的结果，确定淮河流域无霜期在1998年发生突变。

图10 淮河流域平均无霜期（a）及其线性变化倾向率（b）的空间分布（1961−2015年）

图11 基于不同基准气候期（a）和不同年代（b）计算的无霜期220d地带的位置变化

图12 淮河流域无霜期的突变分析

3 结论与讨论

3.1 结论

（1）近56a淮河流域平均初霜日期以2.15d·10a−1（P＜0.01）的速率呈极显著推迟趋势，平均终霜日期以2.49d·10a−1（P＜0.01）的速率呈极显著提前趋势，平均无霜期以4.38d·10a−1（P＜0.01）的速率呈极显著延长趋势。平均初霜日期整体呈“V”型分布，20世纪90年代初霜日期以16.38d·10a−1（P＜0.01）的速率呈极显著推迟趋势；平均终霜日期整体呈“W”型分布，90年代以5.34d·10a−1(P＜0.01)的速率呈极显著提前趋势，21世纪00年代则与之相反，以11.03d·10−1(P＜0.01)的速率呈极显著推迟趋势。平均无霜期整体则呈“阶梯”上升型分布，20世纪90年代的上升速率达到最大，为20.6d·10a−1(P＜0.01)，21世纪00年代以3.69d·10a−1(P＜0.05)的速率呈显著上升趋势。

（2）初霜日期呈“阶梯式”分布，北早南迟，山区早平原迟；86.9%的站点初霜日期呈显著推迟趋势（P＜0.05），在34°N 附近提前趋势最明显；1981−2010年和21世纪00年代的11月8日线位置分别较其它各气候基准期和各年代偏北约1个和2个纬距，表明初霜日期整体呈推迟的趋势。终霜日期呈西南早、东北迟，平原早，山区迟的规律；83.6%的站点通过P＜0.05水平的显著性检验，以3.44～5.92d·10a−1的速率呈提早趋势，在32°N和34°N附近提前趋势最为明显；4 月1 日线在1981−2010年和21世纪00年代位置最为偏北，较其它时间段分别偏北约1个和2个纬距，表明终霜日期整体有提前趋势。无霜期从西南至东北逐渐减少，平原长，山区短；93.4%的站点通过P＜0.05水平的显著性检验，变化率在3.56～7.59d·10a−1，且在32°N和34°N附近延长趋势最为显著；1971−2000年和20世纪70年代的220d等值线位置最为偏北，较其它时间段分别偏北约0.5个和2个纬距，表明无霜期整体延长。

（3）淮河流域初霜日期在2002年发生突变，终霜日期在1995年、无霜期在1998年发生突变。

3.2 讨论

根据IPCC第五次评估报告[20]，北半球过去30a（1983−2012年）可能是过去1400a以来最热的30a。初霜日期的推迟、终霜期的提前和无霜期的延长均在气候突变发生之后表现最为明显（1981−2010年和21世纪00年代），这与高超等研究的结果保持一致。高超等[21]认为，1958−2007年淮河流域年平均气温在1990年以前以降温为主，1990年以后增温显著，可见，11月8日等值线位置北移主要是该区域气温增加所致。气候基准期不同，计算的区域初霜日期有明显区别，而11月8日等值线随气候基准期和年代的更替逐渐北移说明该区初霜日期整体有推迟的趋势，对该区冬小麦的生长和成熟产生一定的影响。首先，冬小麦的“条链病”发生率将会增加；其次，整地播种的进度会明显加快；最后，作物霜冻害是以植株受到伤害为标准的，初霜日推迟导致冬小麦的生长季延长，造成冬小麦生长加速，在抽穗期容易发生霜冻[22]。最新研究发现，气候变暖总体上减少了每年的霜冻天数，20世纪80−90年代霜冻增加最多，90年代−21世纪00年代增加减缓。而霜冻天数增加明显的区域其生长季延长也更明显，无霜期的延长可以使复种指数上升，生长季延长增加了植物对霜冻的暴露，但植物生长对霜冻的敏感性因物种、生长条件、生长阶段而异[23]，霜冻天数增加不一定会对植物带来损害[24]。春秋季节植物物候学的基本机制仍有很大一部分未知，因此，在无霜期延长[25-26]，异常霜冻发生频率下降[27]的气候背景下，迫切需要进行实地观测和野外实验，以提高对植物物候与霜冻损害二者关系的理解。

[1]李芬,于文金,张建新,等.山西近50a初霜冻的气候变化特征及其周期分析[J].大气科学学报,2012,35(6):754-761.

Li F,Yu W J,Zhang J X,et al.Characteristics and period analysis of the climatic changes of first frost in the past 50 years in Shanxi Province[J].Transactions of Atmospheric Sciences,2012,35(6):754-761.(in Chinese)

[2]Erlat E,Turkes M.Analysis of observed variability and trends in numbers of frost days in Turkey for the period 1950-2010[J].International Journal of Climatology,2012,32(12): 1889-1898.

[3]Terando A,Easterling W E,Keller K,et al.Observed and modeled twentieth-century spatial and temporal patterns of selected agro-climate indices in North America[J].Journal of Climate,2012,25(2):473-490.

[4]宁晓菊,张丽君,杨群涛,等.1951年以来中国无霜期的变化趋势[J].地理学报,2015,70(11):1811-1822.

Ning X J,Zhang L J,Yang Q T,et al.Trends of frost-free period in China since 1951[J].Acta Geographica Sinica,2015,70(11): 1811-1822.(in Chinese)

[5]杨晓玲,胡津革,张宇林.河西走廊东部初、终霜冻日和无霜期变化特征分析[J].中国农学通报,2016,32(11):149-155.

Yang X L,Hu J G,Zhang Y L.The analysis of first and last frost date and frost free period variation in the east of Hexi Corridor[J].Chinese Agricultural Science Bulletin,2016, 32(11):149-155.(in Chinese)

[6]张山清,普宗朝,李景林,等.气候变暖背景下新疆无霜冻期时空变化分析[J].资源科学,2013,35(9):1908-1916.

Zhang S Q,Pu Z C,Li J L,et al.Analysis of spatio-temporal variation of frost free period in Xinjiang under the background of climate warming[J].Resources Science,2013, 35(9):1908-1916.(in Chinese)

[7]杜军,石磊,袁雷.近50年西藏主要农区霜冻指标的变化特征[J].中国农业气象,2013,34(3):264-271.

Du J,Shi L,Yuan L.Variation characteristics of frost index in major agricultural areas of Tibet in recent 50 years[J].Chinese Journal of Agrometeorology,2013,34(3):264-271.(in Chinese)

[8]李芬,张建新,闫永刚,等.山西近50年初霜冻的时空分布及其突变特征[J].中国农业气象,2012,33(3):448-456.

Li F,Zhang J X,Yan Y G,et al.Spatial and temporal distribution and mutation characteristics of frost in Shanxi Province in the past 50 years[J].Chinese Journal of Agrometeorology,2012, 33(3):448-456.(in Chinese)

[9]张霞,钱锦霞.气候变暖背景下太原市霜冻发生特征及其对农业的影响[J].中国农业气象,2010,31(1):111-114.

Zhang X,Qian J X.Occurrence of frost in Taiyuan under climate warming and its impact on agriculture[J].Chinese Journal of Agrometeorology,2010,31(1):111-114.(in Chinese)

[10]梁进秋,贾利芳,何正梅,等.大同市近50年霜冻及其异常事件变化特征分析[J].中国农业气象,2010,31(S1):124-128.

Liang J Q,Jia L F,He Z M,et al.Analysis of the variation of frost and its anomalous events in Datong in recent 50 years[J].Chinese Journal of Agrometeorology,2010,31(Suppl1): 124-128.(in Chinese)

[11]国家科学技术委员会.中国科学技术蓝皮书第5号《气候》[M].北京:科学技术文献出版社,1990,25-26.

National Science and Technology Commission.China science and technology blue book No.5《Climate》[M].Beijing:Science and Technology Literature Publishing House,1990,25-26.(in Chinese)

[12]韩荣青,李维京,艾婉秀,等.中国北方初霜冻日期变化及其对农业的影响[J].地理学报,2010,65(5):525-532.

Han R Q,Li W J,Ai W X,et al.Changes in the date of first frost in the north of China and its influence on agriculture[J].Acta Geographica Sinica,2010,65(5):525-532. (in Chinese)

[13]符淙斌,王强.气候突变的定义和检测方法[J].大气科学,1992,16(4):482-493.

Fu Z B,Wang Q.The definition of climate change and detection methods[J].Atmospheric Science,1992,16(4):482-493.(in Chinese)

[14]Reid P C,Hari R E,Beaugrand G,et al.Global impacts of the 1980s regime shift[J].Global Change Biology,2016,22(2):682.

[15]Woolway R I,Dokulil M T,Marszelewski W,et al.Warming of central European lakes and their response to the 1980s climate regime shift[J].Climatic Change,2017,142(3-4):505-520.

[16]Yasunaka S,Hanawa K.Regime shifts in the Northern Hemisphere SST Field:revisited in relation to tropical variations[J].Journal of the Meteorological Society of Japan,2003,81(2):415-424.

[17]Xiao D,Li J P,Zhao P.Four-dimensional structures and physical process of the decadal abrupt changes of the northern extratropical ocean–atmosphere system in 1980s[J]. International Journal of Climatology,2012,32(7):983-994.

[18]Figura S,Livingstone D M,Hoehn E,et al.Regime shift in groundwater temperature triggered by the Arctic Oscillation[J]. Geophysical Research Letters,2011,38(23):23401.

[19]李芬.山西霜冻发生规律及其区域特征研究[D].南京:南京信息工程大学,2012.

Li F.Study on the regularity and regional characteristics of frost in Shanxi[D].Nanjing:Nanjing University of Information Science and Technology,2012.(in Chinese)

[20]IPCC.Working Group I Contribution to the IPCC Fifth Assessment Report(AR5).Climate change 2013:the physical science basis. Final Draft Underlying Scientific-Technical Assessment[R/OL].IPCC,2013-10-30.

[21]高超,姜彤,翟建青.过去(1958-2007)和未来(2011-2060)50年淮河流域气候变化趋势分析[J].中国农业气象,2012,33 (1):8-17.

Gao C,Jiang T,Zhai J Q.Analysis and prediction of climate change in the Huaihe River Basin[J].Chinese Journal of Agrometeorology,2012, 33(1):8-17. (in Chinese)

[22]Zheng B,Chapman S C,Christopher J T,et al.Frost trends and their estimated impact on yield in the Australian wheat belt[J].Journal of Experimental Botany,2015,66(12):3611.

[23]Augspurger C K.Reconstructing patterns of temperature, phenology,and frost damage over 124 years:spring damage risk is increasing[J].Ecology,2013,94(1):41-50.

[24]Liu Q,Piao S,Janssens I A,et al.Extension of the growing season increases vegetation exposure to frost[J].Nature Communication,2018.doi:10.1038/s41467-017-02690-y.

[25]王媛媛, 张勃. 1971年至2010年陇东地区霜冻与积温变化特征[J].资源科学,2012,34(11):2181-2188.

Wang Y Y,Zhang B.Characteristics of frost days and accumulated temperature in eastern Gansu over the last 40 years[J]. Resources Science, 2012, 34(11):2181-2188. (in Chinese)

[26]马琼. 1980-2014年黄土高原霜冻时空变化及冬小麦晚霜冻危险性研究[D]. 兰州：西北师范大学,2016.

Ma Q.Spatio-temporal change of frost and loess risk of winter wheat from 1980 to 2014[D].Lanzhou: Northwest Normal University, 2016.(in Chinese)

[27]马彬,张勃,贾艳青,等.1961-2014年中国内陆农业区异常初、终霜日时空变化及其与环流因子的关系[J].气象学报,2017, 75(4):661-671.

Ma B,Zhang B,Jia Y Q,et al.Temporal and spatial variations of the first and last frost dates in China’s inland agricultural region from 1961 to 2014 and their relationships with circulation factors[J]. Acta Meteorologica Sinica,2017,75(4): 661-671.(in Chinese)

Analysis on the Temporal and Spatial Changes of Frost Date in the Huaihe River Basin from 1960 to 2015

MA Shang-qian1, ZHANG Bo1, TANG Min1, MA Bin1, YANG Mei2, LIANG Jing-jing1, ZHANG Jia-qi1, JIA Li-ge1

(1. College of Geography and Environment Science of Northwest Normal University, Lanzhou 730070, China; 2. The School of Social Development and Public Policy, Beijing Normal University, Beijing 100000)

Ground 0cm daily minimum temperature collected at 61 meteorological stations were used to identify the first and last frost dates. Linear trend estimation was used to describe the trends in the first frost date, last frost date and the length of the frost-free period. The Mann-Kendall mutation test, the sliding T-test and the IDW method were used to analyze the mutation year and spatial variations of the first frost date, last frost date and the frost-free period in the Huaihe River Basin. The results showed that: (1) with the speed of 2.15, 2.49, 4.38d·10y−1(P＜0.01)of the first frost date, the last frost date and the frost-free period, respectively, showing the trends in delaying, advance and extension. The rates of change in the three groups were the most significant in the 1990s, at 16.38, 5.34 and 20.6d·10y−1(P＜0.01), respectively. (2) The first frost date was spatially late in the south, early in the north, early in the mountain and late in the plain. There were 86.9% of the stations postponing significantly (P＜0.05). The last frost date was early in the southwest, late in northeast, early in the plain and late in the mountain, 83.6% stations showed an early trend with the rate of 3.44−5.92d·10y−1(P＜0.05). The frost-free period decreased with the increase of latitude and altitude. 93.4% stations passed the 0.05 significant level test, and the extension rate of change was 3.56−7.59d·10y−1(P＜0.05). (3) The locations of the contour lines of 8 November, 1 April and 220 days were slightly northerly about 1° and 2° than the other climate reference periods and the decadal, with the most obvious trend being around 32°N and 34°N, which corroborated the 2ndconclusion. (4) Mutations of first frost date, last frost date and the frost-free period occurred in 2002, 1995 and 1998, respectively.

Huaihe River Basin; First frost date; Last frost date; Frost-free period

10.3969/j.issn.1000-6362.2018.07.005

马尚谦,张勃,唐敏,等.1960−2015年淮河流域初终霜日时空变化分析[J].中国农业气象,2018,39(7):468-478

2017−12−16

。E-mail：zhangbo@nwnu.edu.cn

国家自然科学基金项目（41561024）；高校博士学科点专项科研基金项目（20136203110002）

马尚谦（1991−），硕士生，主要从事气候变化与农业生态研究。E-mail: msqnwnuedu@163.com