1961—2013年贵州省地面温度时空变化特征

程清平　王平　谭小爱

摘要：利用貴州省1961-2013年逐日0 cm地面温度观测资料，采用ArcGIS空间分析、线性趋势分析、小波分析、Mannkendall非参数统计检验法等对贵州省地温的空间分布、变化趋势、突变特点、周期变化等进行研究。结果表明：平均地温、平均最低地温年际、季节空间分布大致从黔西北向黔东南随经度增加、海拔和纬度的降低而逐渐升高；平均最高地温年际与秋、冬季空间分布大致从黔东北向黔西南随着经度、纬度的降低和海拔的升高而逐渐升高。平均地温与平均最低地温年际、季节变化均呈增温趋势，尤以平均最低地温年际、季节增温趋势最显著；平均最高地温年际与春、夏、秋季均呈增温趋势，冬季呈缓慢下降趋势。平均地温、平均最高地温、平均最低地温年际、季节突变集中于20世纪90年代以来；且均以春、秋季突变十分显著。平均地温、平均最高地温、平均最低地温年际、季节周期变化主要以长周期为主周期。年平均地温、年平均最低地温与海拔呈显著的负相关，与经度呈显著的正相关，而年平均最高地温则与纬度呈显著的负相关。

关键词：贵州省；地面温度；时空变化；周期变化

中图分类号：S159文献标志码：A文章编号：

16721683（2018）02012210

Abstract：

Using the daily 0 cm soil surface temperature data of Guizhou during 1961 to 2013，adopting the methods of ArcGIS spatial analysis，linear trend，wavelet analysis，and Mann Kendall nonparametric test，we analyzed the spatial distribution，change trend，mutation characteristics，and variation period of the surface temperature in Guizhou.The results showed that the annual and seasonal spatial distribution of average surface temperature and average minimum surface temperature was roughly increasing from northwestern Guizhou to southeastern Guizhou along with the increase of longitude and the decrease of altitude and latitude.The annual，autumn，and winter spatial distribution of average maximum surface temperature was roughly increasing from the northeast to the southwest of Guizhou along with the decrease of latitude and longitude and the increase of altitude.The annual and seasonal changes of average surface temperature and average minimum surface temperature all showed a warming trend，which was especially prominent for the annual and seasonal changes of average minimum surface temperature.The annual，spring，summer，and autumn changes of average maximum surface temperature all showed a warming trend，but the winter change showed a slight downward trend.The annual and seasonal mutations of average surface temperature，average maximum surface temperature，and average minimum surface temperature concentrated in the 1990s，and the mutations were very significant in spring and autumn.The annual and seasonal variation periods of average surface temperature，average maximum surface temperature，and average minimum surface temperature were primarily long periods.The annual average surface temperature and annual average minimum surface temperature were both negatively correlated with the altitude and positively correlated with the longitude，while the annual average maximum surface temperature was negatively correlated with latitude.

Key words：

Guizhou；surface temperature；temporal and spatial variation；periodic variation

地温是表征地表自然生态环境的一个重要物理量，对地表向大气长波能量的辐射能力具有决定作用，也是气候系统变化的一个重要参数[12]。地温是指地表面和以下不同深度处土壤温度的统称，包括地面温度（0 cm）、浅层地温（距离地面5 cm、10 cm、15 cm、20 cm） 和深层地温（距离地面40 cm、80 cm、160 cm、320 cm）[3]。土壤温度是重要的土壤物理性质，其变化特征直接影响着土壤生物的生长发育、[JP+1]繁殖和分布，与农业生产和生态环境紧密相关[4]。农业是对气候变化反应最敏感的领域之一，气候变化会影响农业生产的结构和布局，也会对农作物的生长发育等产生影响。在气候变暖背景下，近地面层气温和各层地温升温速率是否一致，关系到利用气温数据评估气候变化对农作物影响的适用性[5]。此外，各种地温数据对林、牧、工程设施等均具有实际作用[6]。因此，对于地温的研究是十分有价值和有意义的工作。近年来，有不少学者采用回归分析方法、MK非参数检验方法和反距离权重插值等方法对地温的变化趋势、周期变化、空间分布等做了研究，发现大部分地区地温升高趋势和周期变化特点明显，且在季节和空间上存在区域差异[712]。

目前，对于地温的研究主要集中在时间变化趋势特点上，对于空间分布特点研究较少，对于贵州地温的研究仅有杨平等[7]对黔东南州0～20 cm浅层地温的周期变化和突变特点进行过分析，而对贵州全省0 cm地温时空变化特征的研究未见报道。此外，对于空间分析方法上采用反距离权重插值方法未能体现特殊地形、地势影响下地温的分布特点。本文利用贵州省33个国家基准/基本气象站0 cm地温观测资料，通过线性趋势分析、MannKendall非参数检验、小波分析和ArcGIS地统计中的协同克里金插值等方法，分析贵州省地温时间变化趋势、突变特征、周期变化和空间分布特点；可为完善和全面认识贵州省气候变化特点，指导农业生产、护林防火、生态环境保护等提供有关科学依据。

1研究区域概况

贵州省（103°40′-109°40′E，24°37′- 29°24 N）地处中国西南云贵高原东部，地势西高东低，地貌以高原、山地为主，平均海拔约1 100 m，国土面积176 167 km2，喀斯特地貌占全省總面积的73%[13]，属亚热带湿润季风气候，[HJ2.1mm]年均气温15 ℃ 左右，年降水800～1 700 mm，光照适中，雨热同季，无霜期270 d左右[14]。地带性植被以常绿阔叶林为主，土壤以黄壤为主，次为石灰土等。[HJ2.1mm]

2数据来源与研究方法

2.1数据来源

逐日平均地温、平均最高地温、平均最低地温来源于中国气象局国家气象信息中心（http：//www.nmic.gov.cn/），DEM高程数据为90 m分辨率的SRTM，来源于地理空间数据云（http：//www.gscloud.cn/），研究站点分布见图1。各站点数据均经过严格的质量控制，包括极值检验、时间一致性检验，最终选择数据完整的1961年1月1日-2014年2月28日逐日0 cm地温数据，利用泰森多边形方法获得贵州省地温平均值（其中，思南、习水平均最高地温、平均最低地温缺1961年，松桃平均最低地温缺1961-1963年，则该站该年不参与计算）。按春季（3月-5月）、夏季（6月-8月）、秋季（9月-11月）、冬季（12月-翌年2月）统计各季平均地温。

2.2研究方法

2.2.1地温倾向率

xi表示样本量为 n 的地温变量，用ti表示xi所对应的时间，建立xi和ti之间的一元线性回归方程：xi=a+bti ，i=1，2，…，n。式中：a为常数项（截距）；b为地温变量的倾向率，b>0 表示直线递增，b<0表示直线递减，b×10表示每10 a 的变化率，单位为℃。

2.2.2MannKendall非参数统计检验方法

Mannkendall[15]是一种非参数统计检验方法，它具有不需要样本遵从一定的分布，也不受少数异常值干扰等优点，能明确地温的演变趋势是否存在突变现象以及突变开始的时间，并指出突变区域。设此气候序列为x1，x2，…，xn，mi表示第i个样本，xi大于xj（1≤j≤i）累计数。定义一秩序列：

查表获得（取α =005，u0.05=1.96），若|UF|>u0，则表明该气候序列呈显著的变化趋势。按气候序列x的逆序列（xn，xn-1，…，x1）再重复上面的过程，同时使UB[WTB1X]k[WTBX]=-UF[WTB1X]k，k[WTBX]=n，n-1，…，UB1=0，若UF>0，则表明序列呈上升趋势，UF<0呈下降趋势，当超过置信水平时，表明上升（下降）趋势显著。如果UF和UB两条曲线出现交点，且交点位于置信区间内，那么交点对应的时刻即是突变开始的时刻。对于MK法突变点较多的情况，同时结合滑动t和pettitt法确认突变点，取3种方法或2种方法检测一致的年份为突变点。

2.2.3小波分析

Morlet连续复小波[16]基函数 （Complex Morlet Wavelet） 在 matlab 系统中简称为cmor，其表达式为：[HJ1mm]

cmor（x）=[SX（]e2iπxfc×[AKσ-D][SX（]x2[]Fb[SX）][][KF（]πFb[KF）][SX）][JY] （4）

式中：Fb是频带宽度系数；Fc是 wavelet 的中心频率。将小波系数的平方值在时间平移域（b）上积分，可得到小波方差 Var，即：

Var（a）=∫∞-∞|wf（a，b）|2db[JY]（5）

Var随时间周期a的变化分布图， 称为小波方差分布图。它能直观地反映信号波动的能量随周期 a的分布，可用于确定不同振荡周期的相对强度和振荡主周期。

2.2.4空间插值方法

有关学者[1719]对气温、降水、空气相对湿度等要素采用多种空间插值分析方法对比分析得出，协同克里金插值方法由于考虑到经度、纬度、地形、地貌等因子的影响，更能提高空间插值的精度。同时，胡丹桂等[14]也指出，站点数量少插值时比站点数量多插值时均方根误差更小，插值效果更佳。基于以上分析，在此采用协同克里金插值法对地温进行空间插值。

[BT2-*3]3地温时空变化特征分析

3.1年际、季节平均地温

3.1.1年际、季节平均地温空间分布

从贵州全省来看，年平均地温大致从黔西北向黔东南，随经度增加、海拔和纬度的降低而逐渐升高（图2），黔南高于黔北，黔东南高于黔西北，黔东高于黔西，其值变化范围在137～227 ℃之间，1961-2013年贵州平均值为179 ℃。低值中心为黔西北[JP3]威宁站（137 ℃）。高值中心为黔西南罗甸站（227 ℃）。四季平均地温空间分布特征（图略）与年平均地温的分布一致，均呈黔东高于黔西，黔南高于黔北的特点。同时，也存在明显的区域差异。四季和年高值中心均为黔西南罗甸站，春、夏、秋、冬季值分别为235 ℃、303 ℃、235 ℃、134 ℃。四季和年低值中心均出现在黔西北，其中春季、秋季低值中心为威宁站，其值分别为155 ℃和135 ℃；夏季低值中心同为普安站和水城站，其值均为232 ℃；冬季低值中心为习水站，其值为56 ℃。

3.1.2年际、季节平均地温变化趋势及突变特征

从图3（a）可知，贵州年平均地温在海拔较高的黔西地区升温幅度较大，而海拔较低的黔东地区升温幅度较小；同时，在升温趋势明显的黔西地区也存在明显的降温中心，而东部区域降温中心和升温中心均有分布，区域差异明显；升温趋势站点中除5个站升温趋势不显著外，其余24个站均通过005及以上显著性检验；降温趋势站点中除兴义站不显著[CM（22]外，其余3个站点均通过005及以上显著性检验。

从图3（b）可知，贵州年平均地温呈显著的升温趋势，其倾向率为012 ℃/（10a）（α<001）。四季平均地温变化趋势与年一致（图略），均呈升温趋势，其中秋季（0161 ℃/（10a），α<001）升温趋势显著，而春（0112 ℃/（10a），α>005）、夏（0081 ℃/（10a），α>005）、冬季（0128 ℃/（10a），α>005）升温趋势不显著。各站点季节变化与年际一致，升温与降温趋势均有所表现，只是在区域上和显著性水平上表现不同。

年平均地温在顺序列UF和逆序列UB之间于2002年交于一点，突变时间区域为2009-2013年，结合滑动t和pettitt法检验突变点分别为2002年（α<001）和1997年（α<001），可见2002为年平均地温显著突变点（图3（c））。春、夏、秋、冬季地温与年平均地温一致，均呈突变升温（图略），MK法、滑动t法（未检测到显著突变点用0表示）、pettitt法检测突变点分别为2004年（α>005）/1997年（α<005）/1997年（α<005）、2009年（α>005）/2004年（α<005）/2004年（α>01）、1997年（α<005）/1997年（α<005）/1997年（α<001）、1990年（α>005）/0/1985年（α>01）；可见，秋季1997年突变最显著，[HJ1.95mm]其次是春季1997年，夏、冬季突变不显著。在005的显著性水平下各站点年际、季节突变特征与整个贵州一致，主要集中突变于20世纪90年代以来（表1）；3种方法均以MK法和pettitt法检测的突变年份较一致。

313年际、季节平均地温周期分析

从图4（a）可知，年平均地温存在2～5 a左右的周期变化，且贯穿整个研究时段，正负相位交替出现。6～9 a左右的短周期变化主要发生在1970-1980年代；而20～30 a左右的长周期变化，主要发生在1970年代以来。从图4（b）可知，年平均地温在24 a时间尺度上，振荡最强，振幅最大，为年平均地温主周期；此外，还存在4 a和8 a的短周期变化。各站点年际变化主周期中盘县、习水、惠水、兴仁、威宁、兴义、正安站以短周期变化为主周期，其余站点均以长周期变化为主周期（表2）。从不同季节来看（图略），春、冬季主要以7 a短周期变化为主周期；同时，还存在一致的4 a短周期变化。夏、秋季均以[CM（22]长周期变化为主周期，分别为27 a、23 a；同时还存[CM）]

在一致的3 a短周期变化。各站点季节变化主周期中春季有10个站点以长周期变化为主周期，集中在16～22 a之间；23个站点以短周期变化为主周期，集中在2～7 a之间。夏季有19个站点以长周期变化为主周期，集中在14～25 a之間；14个站点以短周期变化为主周期，集中在2～3 a之间。秋季有5个站点以长周期变化为主周期，集中在16～26 a之间；28个站点以短周期变化为主周期，集中在3～8 a之间。冬季有25个站点以长周期变化为主周期，集中在11～25 a之间；8个站点以短周期变化为主周期，集中在3～4 a之间。

3.2年际、季节平均最高地温

3.2.1年际、季节平均最高地温空间分布

贵州年平均最高地温呈现出从黔东北向黔西南随着纬度和经度的降低和海拔的升高而逐渐升高的趋势（图5），黔西南高于黔东北，黔西高于黔东，其值变化范围在269～375 ℃之间，1961-2013年贵州平均值为311 ℃；同时，也表现出明显的区域差异。年平均最高地温最大值出现在黔西南的望谟站（373 ℃），最小值为黔西北的习水站（269 ℃）。春、夏季最小值均出现在黔东南的独山站（春季为279 ℃，夏季为381 ℃），而不是纬度最北端的正安站（春季为303 ℃，夏季为452 ℃），春季最大值为黔西南的兴义站（393 ℃），夏季最大值为黔西南的罗甸站（467 ℃）；秋、冬季最小值与年平均最高地温最小值一致均出现在黔西北的习水站（秋季为260 ℃，冬季为131 ℃），最大值也与年平均最高地温最大值一致为黔西南的望谟站（秋季为359 ℃，冬季为244 ℃）。

3.3.2年际、季节平均最高地温变化趋势及突变特征

从图6（a）可知，年平均最高地温在黔中、东部[HJ1.84mm]降温与升温站点均有分布，以降温站点较升温站点多；而海拔较高的西部以增温趋势为主（图6（a））。同时也存在普安站（-0029 ℃/（10a），α>005）、兴义站（-0006 ℃/（10a），α>005）2个略微呈下降趋势的站点。33个站点中，有12个站点呈降温趋势，且主要位于贵州东部，但仅有天柱站（-0523 ℃/（10a），α<0001）、榕江站（-0399 ℃/（10a），α<001）、贵阳站（-0370 ℃/（10a），α<001）、思南站（-006 ℃/（10a），α<005）降温趋势显著；其余21个站点均呈升温趋势，其中有11个站点通过005及以上显著性检验，10个站点未通过显著性检验。1961-2013年贵州年平均最高地温呈升温趋势（图6（b）），其倾向率为0114 ℃/（10a），未通过005的显著性检验。从季节变化来看，春季（0164 ℃/（10a））、夏季（0088 ℃/（10a））、秋季（0219 ℃/（10a））与年平均最高地温一致呈升温趋势，而冬季（-0009 ℃/（10a））呈略微的降温趋势，但四季均未通过显著性检验。各站点季节变化与年际变化表现特征较一致，升温与降温趋势并存，只是在区域上和显著性水平上有所差别。

年平均最高地温在正序列UF和逆序列UB之间于2008年交于一点，但未突破临界线，因此，MK[CM（22]法未检测出显著突变点（图6（c））；结合滑动t和[CM）]

pettitt法检验一致的突变[HJ1.95mm]点为2000年（α<005），

可见，2000为年平均最高地温突变点。春、夏、秋、冬季MK法、滑动t法、pettitt法检测突变点分别为2005年（α>005）/2000年（α<005）/2000年（α<005）、2010年（α>005）/2002年（α<005）/2002年（α>01）、1998年（α<005）/1998年（α<005）/1997年（α>005）、1965年（α>005）/0/1997年（α>01）；可见，春季突变于2000年、秋季突变于1998年较显著，夏、冬季不显著。在005的显著性水平下各站点年际、季节突变特征与整个贵州表现较一致，主要集中突[HJ1.8mm]变于1990年代以来（表1）；3种方法均以MK法和pettitt法检测的突变年份较一致。

323年际、季节平均最高地温周期分析

从图7（a）年平均最高地温小波变换系数实部得知，年平均最高地温在2～5 a左右处，正负相位交替变换且贯穿整个时段；6～9 a左右的周期变化主要发生在1970-1980年代；此外，还存在20～30 a左右贯穿整个时段的长周期变化。从图7（b）小波方差得知，年平均最高地温在23 a年的长时间尺度上峰值最大，振幅最强为年平均最高地温主周期；此外，还存在3 a、8 a的短周期变化。各站点主周期中除惠水、兴仁、威宁站以短周期变化为主周期外，其余站点均以长周期变化为主周期（表2）。从四季来看，春、夏季以长周期变化为主周期，分别为22 a和26 a；同时，春季还存在一致的4 a短周期变化，夏季还存3 a、5 a的短周期变化。秋、冬季主要以短变化周期为主，分别为5 a和8 a；此外，冬季还存在4 a短周期变化。各站点季节变化主周期中春季有18个站点以长周期变化为主周期，集中在15～26 a之间；15个站点以短周期变化为主周期，集中在2～7 a之间。夏季有29个站点以长周期变化为主周期，集中在11～26 a之间；4个站点以3 a短周期变化为主周期。秋季有17个站点以长周期变化为主周期，集中在11～27 a之间；16个站点以短周期变化为主周期，集中在3～10 a之间。冬季有28个站点以长周期变化为主周期，集中在11～27 a之间；5个站点以短周期变化为主周期，集中在4～7 a之间。

33年平均最低地温

3.3.1年际、季节平均最低地温空间分布

从全省来看，年平均最低地温大致从黔西北向黔东南随经度的增加、海拔和纬度的降低而逐渐升高（图8），与年平均地温分布较一致，其值变化范围在61～159 ℃间，1961-2013年贵州平均值为117 ℃。年平均最低地温低值中心位于黔西的威宁站（61 ℃），高值中心位于黔西南的望谟站（155 ℃）。四季平均最低地温空间分布与年平均最低地温空间分布一致，其最小值均位于黔西的威宁站（春季58 ℃，夏季130 ℃，秋季为72 ℃，冬季为-16 ℃）；最大值与年平均最低地温不同，夏季最大值为榕江站（232 ℃），春、秋、冬季最大值均为罗甸站（春季162 ℃，秋季167 ℃，冬季76 ℃）。

33.2年际、季节平均最低地温变化趋势及突变特征[BT）]

从年平均最低地温变化趋势空间分布（图9（a））来看，所研究站点中除仁怀（-0172 ℃/（10a），α<005）和盘县站（-0004 ℃/（10a），α>005）呈降温趋势和贵阳（0052 ℃/（10a））、三橞（006 ℃/（10a））和兴义（0027 ℃/（10a））升温趋势不明显外。其余28个站点均通过005及以上显著性检验，增温趋势十分显著。从年变化趋势来看，贵州1961-2013年年平均最低地温呈极其显著的增温趋势（图9（b）），其倾向率为024 ℃/（10a）（α<0001）；四季与年平均最低地温一致呈极其显著的增温趋势，以冬季（0302 ℃/（10a），α<0001）增温幅度最大，其次是秋季（022 ℃/（10a），α<0001）和夏季（0202 ℃/（10a），α<0001），春季（0179 ℃/（10a），α<0001）增温幅度最小。各站点季节变化中除秋季均呈升溫趋势外，春季、夏季、冬季与年一致，升温与降温趋势均有所表现，只是在站点和显著性水平上有所差别。

从突变特点来看，年平均最低地温在顺序列UF和逆序列UB之间于1996年和1997年交于一点（图9（c）），[JP+1]并于2002-2013年突破0001的临界线，突变增温趋势极其显著；结合滑动t和pettitt法检验突变点分别为1997年（α<005）和1993年（α<0001）；可见，显著突变点为1997年。春、夏、秋、冬季平均最低地温与年一致，呈突变升温趋势，MK法、滑动t法、pettitt法检测突变点分别为2003年（α<005）/2003年（α<005）/1996年（α<0001）、1990年（α<001）/1977年（α<001）/1977年（α<0001）、2002年（α<001）/2002年（α<005）/2002年（α<005）、1988年（α<001）/0/1975年（α<001）；春季一致突变点为2003年，夏季一致突变点为1977年，秋季一致突变点为2002年，冬季未有一致突变点。在005的显著性水平下，各站点年际、季节突变点与整个贵州一致（表1），突变点主要集中在20世纪90年代以来，3种方法均以MK法和pettitt法检测的突变年份较一致。

323年际、季节平均最低地温周期分析

从图10（a）知，年平均最低地温存在5～10 a左右的短周期变化，正负相位交替出现，且贯穿整个研究时段；此外，还存在25～30 a左右的长周期变化。从小波方差（图10（b））知，年平均最低地温在26 a的长时间尺度上振荡最强，为年平均最低地温主周期，同时还存在8 a短周期变化。各站点年际变化主周期中除安顺水、都匀、兴义站以短周期变化为主周期外，其余站点均以长周期变化为主周期（表2）。春、夏、秋、冬季周期变化与年一致，均以长周期变化为主周期，分别为29 a、28 a、26 a、26 a；此外，春季还存在4 a短周期和14 a长周期变化，夏季还存在4 a、8 a短周期和16 a长周期变化；秋、冬季还存在一致的8 a短周期变化。各站点季节变化主周期中春季有13个站点以长周期变化为主周期，集中在15～26 a之间；20个站点以短周期变化为主周期，集中在2～8 a之间。夏季有19个站点以长周期变化为主周期，集中在15～27 a之间；14个站点以短周期变化为主周期，集中在2～8 a之间。秋季有27个站点以长周期变化为主周期，集中在14～26 a之间；6个站点以短周期变化为主周期，集中在3～9 a之间。冬季有15个站点以长周期变化为主周期，集中在11～25 a之间；18个站点以短周期变化为主周期，集中在2～8 a之间。

4讨论

Shao等[20]指出，气温的变化与地形、海拔、经度、纬度有关，影响程度从大小表现为海拔>纬度>地形>经度的形式。贵州年平均地温与经度、纬度、海拔相关系数分别为0424（α<005）、-0328（α>005）、-0750（α<001），影响大小表现为海拔>经度>纬度的特点；年平均最高地温与经度、纬度、海拔的相关性系数分别为-0247（α>005）、-0601（α<001）、-0064（α>005），影响大小表现为纬度>经度>海拔的特点；年平均最低地温与经度、纬度、海拔的相关性系数分别为0602（α<001）、-0184（α>005）、-0875（α<001），影响大小表现为海拔>经度>纬度的特点；可见经度和海拔对年平均地温和年平均最低地温有显著的影响。而年平均最高地温与纬度密切相关，与经度和海拔相关性不明显。此外，年平均地温、年平均最高地温、年平均最低地温的空间分布格局还与大气环流有密切的关系[2122]，尤其是冬半年昆明准静（滇黔准静止锋）止锋在东、西部控制时间的长短和影响强度的空间差异[2324]。贵州年平均地温、年平均最高地温、年平均最低地温均以21世纪以来增温最为显著，与黔东南州地温的变化趋势一致[7]。土壤温度是重要的土壤物理性质，其变化特征直接影响着土壤生物的生长发育、 繁殖和分布，与农业生产和生态环境紧密相关[4]，因此贵州地温呈显著增温这一变化趋势，可能会对农、林、牧和生态环境等造成不同影响，是十分有必要关注的问题。

5结论

（1）贵州年平均地温、年平均最低地温大致從黔西北向黔东南随着经度的增加、纬度和海拔的降低而逐渐升高，黔南高于黔北，黔东南高于黔西北，黔东高于黔西；年平均最高地温从黔东北向黔西南随着纬度和经度的减少和海拔的升高而逐渐升高。四季平均地温和平均最低地温空间分布与年际分布一致。秋、冬季平均最高地温与年平均最高地温一致，春、夏季有所不同。

（2）贵州年平均地温、年平均最高地温、年平均最低地温整体均呈增温趋势。空间分布上表现为海拔较高的黔西地区年平均地温和年平均最低地温升温幅度大，同时也存在降温区域。海拔较低的黔中、东部地区年平均地温升温幅度小且存在降温区域；而年平均最低地温升温幅度大，同时也存在降温区域。年平均最高地温呈现海拔较高的黔西地区以升温为主，海拔较低的黔中、东部地区以降温为主。年平均地温、年平均最高地温、年平均最低地温四季变化趋势与年际变化趋势较一致，升温与降温趋势并存，只是在区域上有所不同。年平均地温与年平均最高地温均以秋季增温最为显著，年平均最低地温以冬季增温最为显著。各站点季节变化趋势中年平均地温和年平均最高地温均以秋季升温趋势较显著；年平均最低地温四季升温趋势均较显著。

（3）贵州年平均地温、年平均最高地温、年平均最低地温均发生显著的突变增温，突变点分别为2002年、2000年、1997年；年平均地温、年平均最高地温、年平均最低地温均以春、秋季突变十分显著，同时年平均最低地温夏季突变也十分显著。

（4）贵州年平均地温、年平均最高地温、年平均最低地温及四季地温均以长周期变化为主周期，长周期变化集中在10～44 a之间，短周期变化集中在2～10 a之间。各站点年平均地温、年平均最高地温、年平均最低地温年际、季节变化主周期也主要以长周期变化为主，长周期主要集中在10～28 a，短周期主要集中在2～9 a。

参考文献（References）：

[1]李崇银.气候动力学引论[M].北京：气象出版社，1995：290296.（LI C Y.Introduction of Climatic Dynamics[M].Beijing：Meteorological Press，1995：290296.（in Chinese））

[2]WANG M Y，LU D R.Diurnal and seasonal variation of clearsky land surface temperature of several representative land surface types in China retrieved by GMS5[J].Acta Meteorologica Sinica，2006，20（4）：475488.

[3]中国气象局.地面气象观测规范[M].北京：气象出版社，2003：8589.（China Meteorological Administration.Ground meteorological observation standard[M].Beijing：China Meteorological Press，2003：8589.（in Chinese））

[4]周绍毅，劳炜，苏志，等.19612010年广西浅层地温变化特征[J].西南农业学报，2012，25（4）：13721375.（ZHOU S Y，LAO W，SU Z，et al.Variation features of fleet soil temperature in Guangxi from 1961 to 2010[J].Southwest China Journal of Agricultural Sciences[J].Southwest China Journal of Agricultural Sciences，2012，25（4）：13721375.（in Chinese））DOI：10.3969/j.issn.10014829.2012.04.051.

[5]陈超，周广胜.19612010年阿拉善左旗气温和地温的变化特征分析[J].自然资源学报，2014，29（1），91102.（CHEN C，ZHOU G S.Characteristics of air temperature and ground temperature in Alxa Left Banner from 1961 to 2010[J].Journal of natural resource，2014，29（1），91102.（in Chinese））DOI：10.11849/zrzyxb.2014.01.009.

[6]戴加洗.青藏高原气候[M].北京：气象出版社，1990：14215.（DAI J X.The climate of Qinghai Tibet Plateau[M].Beijing：China Meteorological Press，1990：14215.（in Chinese））

[7]杨平，陈少华，顾欣.近50年黔东南浅层地温时空变化特征分析[J].西南师范大学学报：自然科学版，2012，37（11）：5458.（YANG P，CHEN S H.On temporal and spatial characteristics of shallow geothermal in Southeast Guizhou in the past fifty years[J].Journal of Southwest China Normal University（Nattural Science Edition），2012，37（11）：5458.（in Chinese））.DOI：10.3969/j.issn.10005471.2012.11.012.

[8]李帅，王萍，陈莉，等.黑龙江省春季浅层（020 cm） 地温变化特征及预测[J].冰川冻土，2014，36（1）：5562.（LI S，WANG P，CHEN L，et al.Variation features and prediction of the topsoil（0～20cm） ground temperature of spring in Heilongjiang Province[J].Journal of Glaciology and Geocryology，2014，36（1）：5562.（in Chinese））DOI：10.7552/j.issn.10000240.2014.0007.

[9]杨晓玲，丁文魁，马中华，等.石羊河流域中下游浅层地温变化及其对气温变化的响应[J].土壤学报，2015，52（6）：14011411.（YANG X L，DING W K，MA Z H，et al.Change in shallow soil temperature and its response to change in air temperature in middle and lower reaches of Shiyang River Basin[J].Acta Pedologica Sinica，2015，52（6）：14011411.（in Chinese））DOI：10.11766/trxb201410290547.

[10][ZK（#]陳超，周广胜.19612010 年桂林气温和地温的变化特征[J].生态学报，2013，33（7）：20432053.（CHEN C，ZHOU G S.Analysis on variation characteristics of air temperature and ground temperature in Guilin from 1961 to 2010[J].Acta Ecologica Sinica，2013，33（7）：20432053.（in Chinese）） DOI：10.5846/stxb201112312022.

[11]李岳坦，李小雁，崔步礼，等.青海湖流域及周边地区浅层地温对全球变化的响应[J].地球环境学报，2010，1（3）：219225.（LI Y T，LI X Y，CUI B L.Response of shallow geotemperature to global change over Qinghai Lake Basin from 1961 to 2007[J].Journal of Earth Environment，2010，1（3）：219225.（in Chinese））

[12]王佳琳， 潘志華， 韩国琳， 等.19612010年中国 0cm地温变化特征及其与气温变化的关系[J].资源科学，2016，38 （9）： 17331741.（WANG J L，PAN Z H，HAN G L，et al.Variation in ground temperature at a depth of 0 cm and the relationship with airtemperature in China from 1961 to 2010[J].Resources Science，2016，38（9）：17331741（in Chinese））DOI：10.18402/resci.2016.09.11.

[13]张波，谷晓平，易俊莲.19612010年贵州省高温气候的时空分布特征[J].贵州农业科学，2015，43（6）：210213.（ZHANG B，GU X P，YI J L，et al.Spatial and temporal distribution characteristic of high temperature in Guizhou during 19612010[J].Guizhou Agrictural Sciences，2015，43（6）：210213.（in Chinese）） DOI：10.3969/j.issn.10013601.2015.06.052.

[14]陈学凯，徐建新，胡娟萍，等.19612012 年贵州省极端降水时空变化特征[J].水资源与水工程学报，2015，26（4）：5061.（CHEN X K，XU J X，HU H J，et al.Feature of spatial and temporal variation of extreme precipitation in Guizhou Province from 1961 to 2012[J].Journal of Water Resources &Water; Engineering，2015，26（4）：5061.（in Chinese）） DOI：10.11705/j.issn.1672643X.2015.04.10.

[15]魏凤英.现代气候统计诊断预测技术[M].北京：气象出版社，2007：6366.（WEI F Y.Modern technology of statistics diagnosis and prediction for climate[M].Beijing：China Meteorological Press，2007：6366.（in Chinese））

[16]FOUFOULAGEORGIOU E，KUMAR P.Wavelet analysis in geophysics：an introduction[J].Wavelet Analysis and Its Applications，1994，4：143.DOI.org/10.1016/B9780080520872.500074.

[17]胡丹桂，舒红.基于协同克里金空气湿度空间插值研究[J].湖北农业科学，2014，53（9）：20452049.（HU D G，SU H.Air Humidity Based on CoKriging Interpolation[J].Hubei Agricultural Sciences，2014，53（9）：20452049.（in Chinese））DOI：10.3969/j.issn.04398114.2014.09.018.

[18]孙然好，刘清丽，陈利顶.基于地统计学方法的降水空间插值研究[J].水文，2010，30（1）：1417.（SUN R H，LIU Q L，CHEN L D.Study on precipitation interpolation based on the geostatistical analyst method[J].Journal of China Hydrology，2010，30（1）：1413.（in Chinese））DOI：10.3969/j.issn.10000852.2010.01.003.

[19]许民，王雁，周兆叶，等.长江流域逐月气温空间插值方法的探讨[J].长江流域资源与环境，2012，21（3）：327333.（XU M，WANG Y，ZHOU Z Y，et al.Discussion of methods of spatial interpolation for monthly temperature data in the yangtze river basin[J].Resources & Environment in the Yangtze Basin，2012，21（3）：327333.（in Chinese））

[20]SHAO J A，LI Y B，NI J P.The characteristics of temperature variability with terrain，latitude and longitude in SichuanChongqing Region[J].Journal of Geographical Sciences，2012，22（2）：223244.DOI：10.1007/s1144201209234.

[21]吴晓娜，孙照渤.欧亚大陆夏季地表温度的气候特征及与大气环流的联系[J].大气科学学报，2015，38（2）：195204.（WU X N，SUN Z B.Climatic features of summer land surface temperature in Eurasian continent and its relationship with atmospheric circulation[J].Transactions of Atmos Science，2015，38（2）：195204.（in Chinese））DOI：10.13878/j.cnki.dqkxxb.20130525001.

[22]赵红岩，杨瑜峰，张久林，等.夏季西太洋副高位置与中国地温场的关系[J].高原气象，2007，26（5）：11191122.（ZHAO H Y，YANG Y F，ZHANG J L，et al.The relationships between the position of west pacific subtropical high and ground temperature field of China Summer[J].Plateau Meteorology，2007，26（5）：11191122.（in Chinese））

[23]段旭，李英，孫晓东.昆明准静止锋结构[J].高原气象，2002，21（2）：205209.（DUAN X，LI Y，SUN X D.The structure of Kunming quasistationary front[J].Plateau Meteorology，2002，21（2）：205209.（in Chinese））DOI：10.3321/j.issn：10000534.2002.02.014.

[24]杜小玲，蓝伟.两次滇黔准静止锋锋区结构的对比分析[J].高原气象，2010，29（5）：1831195.（DU X L，LAN W.Contrastive analysis on frontal structure of quasistationary front in two precipitation processes of YunnanGuizhou[J].Plateau Meteorology，2010，29（5）：11831195.（in Chinese））