周建平，孫照渤，倪東鴻，李忠賢

(氣象災(zāi)害教育部重點實驗室(南京信息工程大學(xué))，江蘇南京210044)

0 引言

如果一個測站的氣象記錄序列僅僅是氣候變化的反映，那么這樣的資料就是均一的(屠其璞等，1984)。均一性的氣溫序列是研究氣候變率和變化趨勢的基礎(chǔ)，然而在實際上，除了觀測等不可避免的誤差以外，影響均一性的主要原因還包括臺站遷移、觀測儀器和安裝方法的更新以及觀測時制的改變等。中國的大量臺站都發(fā)生了遷移，從而影響了所測資料的均一性，使得研究結(jié)果不能正確反映氣候的變化趨勢。因此，有必要研究中國臺站遷移對臺站數(shù)據(jù)均一性的影響。

自20世紀80年代以來，國內(nèi)外許多氣候?qū)W家在資料非均一性判斷方法及其訂正方面做了大量研究工作。Hawkins(1977)從統(tǒng)計角度提出了標準正態(tài)檢驗(standard normal homogeneity test，SNHT)方法;Alexandersson(1986)、Alexandersson and Moberg(1997)改進發(fā)展了該方法，并應(yīng)用到氣象序列的均一性檢驗中，提出參考臺站相關(guān)系數(shù)達到0.8時檢驗效果較好，并對瑞典觀測數(shù)據(jù)的均一性檢驗做了許多開創(chuàng)性的工作。屠其璞等(1984)、么枕生和丁裕國(1990)較早開展了序列訂正方法的研究。宋超輝等(1995)、宋超輝和孫安健(1995)研究了氣溫序列非均一性檢驗的3種方法和訂正的4種方法，認為逐步多元線性回歸對序列的訂正效果較好。李慶祥等(2003)系統(tǒng)研究了定點觀測均一性檢驗方法和國外此類工作的開展情況。吳增祥(2005)研究了中國氣象臺站歷史沿革資料，并分析了其對觀測資料均一性的影響，認為臺站遷移對觀測記錄均一性的影響最大。此外，許多氣象學(xué)者對某個地區(qū)、省份或單個臺站遷移對各種要素均一性的影響進行了有益探討(劉小寧，2000;高曉榮等，2008;江志紅等，2008;田紅等，2008;王鈺和黃少平，2008;張弦等，2009;李慶祥等，2010;李祥余等，2010)。

不同于單個臺站的分析，要整體上討論中國大量測站觀測數(shù)據(jù)的均一性是一個非常復(fù)雜的工作。一方面是檢驗的準確性問題;另一方面由于缺乏臺站歷史沿革資料(元數(shù)據(jù))以及臺站數(shù)量眾多，所以分析破壞均一性的因素比較困難。因此，本文主要分析中國臺站的遷移狀況以及年平均溫度、年平均最高溫度和年平均最低溫度序列的均一性，以研究臺站遷移對溫度序列均一性的影響。

1 資料和方法

采用中國國家氣象信息中心整編的中國753個地面觀測站的數(shù)據(jù)，針對1951年1月1日—2010年12月31日的逐日氣溫數(shù)據(jù)，按照氣候統(tǒng)計規(guī)范構(gòu)造每個站點月平均溫度、年平均溫度、平均最高溫度、平均最低溫度序列。在中國由于臺站的改制或撤銷，存在許多業(yè)務(wù)代替的臺站，一般這些代替的臺站距離并不遠，所以將兩站數(shù)據(jù)合并，前后可以作為一次遷移考慮。根據(jù)臺站代替記錄和數(shù)據(jù)的起止時間，取合并后的臺站號和經(jīng)緯度作為這一序列的臺站號和經(jīng)緯度，如54714德州站1995年業(yè)務(wù)由54715陵縣代替，將兩個序列合并后，取臺站號為54715。將753站共計39次臺站代替的數(shù)據(jù)序列進行合并，最后得714站(圖1a)。序列的長度并不一致，超過50 a的臺站有621站，55 a以上的臺站超過一半。20世紀50年代是臺站數(shù)量飛速發(fā)展的時期，由1951年不足148站快速增長到1960年648站(圖1b)，此后隨著部分臺站的改制或撤銷以及一些臺站的新建，逐年略微有所增減。

遷移在中國的臺站中普遍存在，但由于缺乏元數(shù)據(jù)，本文通過原始逐日數(shù)據(jù)——經(jīng)緯度的一致性來確定是否有遷移。在某一日該站的經(jīng)緯度與前一天不一致時，則認為該站在該日有遷移的可能性。為了去除錯誤的經(jīng)緯度記錄和部分短期代替的數(shù)據(jù)，作如下處理:如果在新的經(jīng)緯度上記錄超過一年，則記為一次遷移，并通過經(jīng)緯度換算遷移的空間距離;否則不記為遷移。

圖1 中國714站臺站分布(a)和1951—2010年逐年臺站數(shù)分布(b)Fig．1 (a)Distribution of 714 stations in China and(b)the station number from 1951 to 2010

通過上述標準確定的臺站遷移主要包括以下幾種情況:1)臺站位置有遷移，即從一個地點遷移至另一個地點;2)因觀測任務(wù)變動，將兩個區(qū)站號相互更換;3)前述提到的業(yè)務(wù)代替;4)臺站經(jīng)緯度可能因測量方法或查算方法不同或者錯誤記錄而存在的變動。4種情況中，前3種均可能造成資料的非均一性;而第四種情況是記錄錯誤，并沒有臺站遷移，對觀測資料的均一性沒有影響，設(shè)置遷移應(yīng)長于一年的規(guī)定后，已盡可能地減少了這種情況的發(fā)生。這里確定的遷移有可能在實際情況中臺站并沒有搬遷;為了方便敘述，所有情況均統(tǒng)稱作為臺站遷移。

目前已發(fā)展了多種均一性的檢驗方法。本文研究遷移對臺站溫度序列的影響，由于溫度序列可以通過經(jīng)緯度的變動來確定變化發(fā)生的時間點，所以較容易通過t檢驗來確定遷移前后溫度的變化特點;因為序列長度的不一致，采用完整的t檢驗公式(魏鳳英，2007):

因為t檢驗是單站檢驗，其實質(zhì)上是檢驗一點前后平均值的變化;而溫度序列中也包含了氣候變化趨勢，所以檢驗的結(jié)果可能表明氣候發(fā)生了突變。因此，均一性檢驗僅僅做單站檢驗是不夠的，所以本文采用均一性檢驗中廣泛應(yīng)用的標準正態(tài)檢驗(SNHT)方法(Alexandersson，1986;Alexandersson and Moberg，1997)。檢驗?zāi)昶骄鶞囟鹊贸龅牟贿B續(xù)點，與臺站的遷移點作比較，以評估臺站遷移對平均溫度的影響。該方法的核心是建立一個相對均一的參考序列，然后將它和檢驗序列做比值或者差值對比，以檢驗序列的均一性。由于年平均溫度分布近似遵從正態(tài)分布，所以本文檢驗所有臺站。

本文檢驗溫度，利用差值法(Alexandersson and Moberg，1997)建立序列:

式中:Yi為待檢臺站的溫度序列;Xji(j=1，2，…，5)為5個參考臺站的參考序列;采用相關(guān)系數(shù)ρ的平方作為加權(quán)因子;為了減少序列的差異，建立序列時候采用距平值。對Q進行標準化，得到待檢序列Zi，使得 ˉZ=0，σZ=1，即使其近似服從 N(0，1)分布:

max大于閾值則認為該點為序列的不連續(xù)點，再將序列分段檢驗，直至序列長度小于10 a止。通過SNHT檢驗確定序列的斷點，將其與確定的遷移點及其t檢驗結(jié)果進行核對，以確定中國臺站遷移對中國氣溫序列的影響。

Alexandersson and Moberg(1997)也提供了序列的訂正方法，如果序列在a點斷裂，則:

SNHT檢驗關(guān)鍵是選擇參考臺站建立檢驗序列。本文主要從距離最近的15個包含待檢序列長度的臺站中選擇相關(guān)系數(shù)較高的5個臺站作為參考臺站。由圖1b可以看出，1955年以前的臺站稀少，選擇參考臺站較為困難。因此，本文主要檢驗1955—2010年的臺站遷移對溫度均一性的影響;1955年以前的遷移主要采用t檢驗并根據(jù)遷移前后平均溫度有無顯著差異來判斷序列的均一性。受氣溫年際變化和參考臺站的影響，如果不是變化較明顯，SNHT檢驗不一定能將不連續(xù)點準確地定位在遷移點上，所以如果不連續(xù)點在該站遷移前后的3 a內(nèi)，則認為此點是由遷移造成的，序列在遷移前后是不均一的。

2 臺站遷移對資料的影響

2.1 遷移狀況

由通過經(jīng)緯度差異方法確定臺站遷移的結(jié)果來看:近60 a的714站中，共有950次遷移;完全沒有遷移的臺站只有148站(約占21%)，304站有1次遷移，遷移2次及以上的有262站，部分臺站甚至有5～7次遷移，可見大多數(shù)臺站都有過遷移。由圖2可見，1960年前后和1980年前后發(fā)生的遷移較多，這主要是由于中國臺站建立的完善和行政管理等方面的原因;1990年以后，臺站發(fā)生的遷移較少，這主要因為臺站環(huán)境得到法律法規(guī)保護，臺站遷移受到嚴格控制，但由于城市規(guī)劃或臺站處于城市中間而不再適合作為觀測臺站，所以有許多臺站發(fā)生了遷移，而且大多數(shù)臺站級別沿革、撤銷、新建、業(yè)務(wù)代替等也大多發(fā)生在這一時期，對數(shù)據(jù)的影響也比以前的遷移要大。950次遷移中601次遷移的水平距離在10 km內(nèi)，192次遷移的水平距離超過20 km;950次遷移中有533次遷移的海拔高度沒有變化，713次遷移的高度差低于10 m，93次遷移的高度差超過50 m;249次遷移距離超過20 km或(和)海拔高度相差超過50 m，中國臺站遷移明顯。

圖2 1952—2008年臺站遷移次數(shù)的逐年變化Fig．2 Yearly variation of the station relocation times from 1952 to 2008

2.2 遷移前后溫度的變化特點

t檢驗實質(zhì)上是對遷移前后均值變化的一種檢驗。分別對950次遷移的年平均溫度、年平均最高溫度、年平均最低溫度進行檢驗，t值大于臨界值的臺站數(shù)如表1所示。從臺站數(shù)和次數(shù)上來說，遷移對最低溫度影響最大，其次是平均溫度，對最高溫度的影響最小。遷移前后平均溫度增加約0.56℃，總體上來說遷移后是變暖的，遷移后溫度減少的臺站相對較少。這也說明臺站遷移與增溫評估存在一定關(guān)系。

表1 950次遷移的t檢驗結(jié)果Table 1 The t-test results for 950 times station relocation

圖3為臺站遷移前后年平均溫度、年平均最高氣溫、年平均最低氣溫顯著變化的臺站。可見:遷移后年平均溫度減少的臺站大多分布在中國西南部;遷移后年平均最高氣溫減少的臺站多一些，主要分布在中國西南、華南和西北地區(qū)東部;遷移后年平均最低氣溫未通過t檢驗的臺站數(shù)大大增加，這與日最低溫度的時間分布和空間分布有較大關(guān)系。值得注意的是，臺站遷移對平均氣溫、最高氣溫和最低氣溫的影響存在差異，如漠河站1997年有一次較大的遷移，年平均氣溫遷移后只增加了0.04℃(沒有通過顯著性檢驗)，而年平均最高氣溫和年平均最低氣溫的前后變化卻較顯著，但是它們的變化卻是相反的，遷移后14 a年平均最高氣溫較遷移前平均增加了0.9℃，年平均最低溫度減少了0.85℃。這反映出臺站遷移前后環(huán)境變化對不同時段溫度變化影響的不一致性。有些臺站主要因海拔高度發(fā)生較大改變而產(chǎn)生了較大變化，如五臺山站(53588)1998年從山頂遷移到山腳，海拔高度減少了687.5 m，平均溫度增加了5.74℃，最高溫度增加了6.35℃，最低溫度增加了5.37℃。如果這些臺站的溫度數(shù)據(jù)不經(jīng)過訂正而使用，在氣候趨勢研究中是沒有意義的。

在影響遷移前后溫度變化的因子中，海拔高度差無疑是最重要的。235次遷移的年平均溫度差與海拔高度差的相關(guān)系數(shù)達－0.85，即海拔高度減少后，溫度升高;而與水平距離的相關(guān)系數(shù)只有0.22，如果考慮到南北遷移對溫度升降的差異性后，其相關(guān)系數(shù)可增加到0.38。所以，遷移的海拔高度差是考慮遷移前后年平均溫度變化的一個主要因素，但是遷移還不能僅僅從水平和垂直距離上考慮，有時需要考慮遷移前后臺站周圍的地理環(huán)境變動，是市區(qū)還是郊區(qū)等，因為這些變動僅僅通過空間位置數(shù)據(jù)是不能反映出來的。

對于年平均溫度來說，950次遷移前后溫度平均增加了0.21℃，其中641次遷移后平均溫度是升高的，差異明顯的235次遷移在遷移后平均增溫約0.56℃。如果t值很高，表明遷移前后溫度差值較大，據(jù)此就能直接判斷出溫度序列的不均一性。但是大多數(shù)臺站的t值接近于臨界值，由于氣候序列中包含了氣候變化趨勢，加上臺站周圍環(huán)境如城市化等緩慢變化的作用也包含于序列變化之中，序列均一性的直接判斷就有部分困難，所以要作鄰站對比的檢驗，將氣候趨勢變化剔除出來。

圖3 年平均氣溫(a)、年平均最高氣溫(b)和年平均最低氣溫(c)遷移前后差異顯著(t檢驗)的臺站以及溫度平均差值分布(單位:℃;方框表示遷移后溫度降低，實點表示遷移后溫度升高)Fig．3 Stations with significant temperature differences caused by station relocation and their differences(units:℃;panes(solid points)denote the temperature decreasing(increasing)after station relocation)a．a(chǎn)nnual mean temperature;b．a(chǎn)nnual mean maximum temperature;c．a(chǎn)nnual mean minimum temperature

3 平均溫度的均一性檢驗

通過公式(4)中Tsmax來檢驗每個臺站的溫度序列。以七角井站(51495)為例說明SNHT檢驗過程及訂正結(jié)果分析(圖4)。從經(jīng)緯度上來說，共有3次遷移，分別在1955、1962、1998年，前兩次遷移距離較小，只有1 km左右，而1998年遷移距離約為30 km，海拔高度降低了83 m，該站的年平均溫度序列如圖4a所示，遷移后12 a的平均溫度較遷移前平均增加了約3℃，有一個非常明顯的突變。選定52418、51777、52203、51156、51573 站作為參考臺站(最遠距離約為600 km)，5個參考臺站與七角井站(51495)年平均溫度的相關(guān)系數(shù)均超過0.78。通過差值法建立逐年Q值(圖4c)、Ts值(圖4e)，可以看出，1998年有一個突變點，Tsmax值達到50.8，遠超過閾值8.6。分別檢驗1955—1997年和1998—2010年序列，未發(fā)現(xiàn)不連續(xù)點。七角井站1998年的遷移造成了序列的非均一性，通過公式(5)確定的訂正值約為2.15℃，訂正1955—1997年的溫度，如圖4b所示。經(jīng)檢驗，訂正后該序列不存在不連續(xù)點，表明它是均一的。

經(jīng)過SNHT檢驗，有632站檢驗序列與參考序列相關(guān)系數(shù)超過0.8，相關(guān)度較好。年平均溫度序列完全均一的臺站只有172站，而其余的542站有可疑的不連續(xù)點884個，但是有174個可疑點距離起始點較近(在3 a內(nèi))，203次檢驗的Tsmax與閾值差值小于2。也就是說，如果改變信度的話，會影響均一性判斷的結(jié)果，加上參考臺站等不確定因素，如果不連續(xù)點前后有臺站的遷移，則認為此不連續(xù)點是遷移造成的，是不均一點。經(jīng)過統(tǒng)計，有237個不連續(xù)點是遷移造成的，將1955年以前t檢驗的4次遷移計入結(jié)果內(nèi)，1951—2010年遷移造成的序列不均一點個數(shù)如圖5所示。將臺站遷移時間(圖2)與非均一點所處時間進行逐一比較可以發(fā)現(xiàn)，1960年前后、1980年前后不均一點的兩個峰值與遷移頻繁年份較吻合，這兩段時期的非均一點較多，主要是由遷移較多引起的。而2000年以后臺站遷移雖然相對較少，但是非均一點占遷移次數(shù)的比例卻明顯增高，說明該時期臺站遷移對數(shù)據(jù)均一性影響要大。

圖 4 七角井站訂正前(a，c，e)、后(b，d，f)年平均溫度序列(a，b;單位:℃)、Q 值(c，d)、Ts序列(e，f)Fig．4 (a，b)Annual mean temperature(℃)，(c，d)Q and(e，f)Ts(a，c，e)before and(b，d，f)after adjusting at Qijiaojing station

圖5 1951—2010年平均溫度非均一點個數(shù)的逐年變化Fig．5 Yearly variation of number of inhomogeneity points for annual mean temperature from 1951 to 2010

將SNHT檢驗的非均一點與遷移t檢驗結(jié)果作對比發(fā)現(xiàn)，237次非均一點與t檢驗235次遷移造成的前后平均溫度有顯著差異的次數(shù)差不多，其中有99次遷移造成了遷移前后的溫度差值非常明顯(平均約0.71℃)。也就是說，這99次遷移前后平均溫度差異顯著，而且序列是非均一的，而其他138次遷移，盡管遷移前后的溫度差異并不顯著，但是同樣造成了序列的非均一性。

溫度在不同的時間尺度上都有不同的變化，日最高、最低氣溫反映的是一天的溫度變化范圍，而隨著太陽輻射的季節(jié)變化，月平均溫度也有差異，而遷移對不同時間段的溫度序列均一性影響的差異也值得研究。分別以1、4、7、10月作為代表月份和以年平均最高、最低氣溫代表不同時段的溫度，檢驗均一性對溫度變化的影響(表2)。結(jié)果發(fā)現(xiàn)，各月可疑點和遷移造成的非均一點都要遠遠低于年平均溫度。以1月為例，可疑點不連續(xù)點有449個，約為平均溫度可疑點數(shù)的一半，遷移造成的非均一點有149個，這其中有112個與年平均溫度是一致的。4、7、10月可疑點和遷移造成的非均一點均有所增加，但是與年平均氣溫一致的遷移次數(shù)相差不大，約有70%與年平均氣溫是一致的，即這些遷移在月時間尺度上造成了序列的非均一性，在年尺度上也是非均一的。當(dāng)然，遷移對不同月份溫度的影響存在差異，而年平均溫度是否均一則是12個月綜合作用的結(jié)果。

表2 不同溫度序列的均一性檢驗及其與年平均溫度非均一點的關(guān)系Table 2 The homogeneity test of different temperature series and their relations to inhomogeneity points of annual mean temperature 個

年平均最高氣溫與年平均氣溫相比，遷移對其影響要稍微小一點。年平均最高氣溫有688個可疑的不連續(xù)點，少于年平均氣溫，有196次是遷移造成的，其中114次和年平均溫度一致，是由同一次遷移造成的。遷移對年平均最低氣溫的非均一性影響比年平均氣溫要大。遷移對其造成的非均一點達到最多(266個)，其中與年平均溫度的237個非均一點有161次是同一次遷移造成的。從遷移造成的非均一點也可以看出，遷移對年平均最低溫度的影響最大，其次是年平均氣溫，對年平均最高氣溫影響較小，這與t檢驗的結(jié)果是一致的。

4 非均一年平均溫度訂正

訂正是將一個非均一的序列訂正為盡可能滿足均一條件的序列。年平均溫度序列有204站共計237個非均一點需要訂正，用訂正差值對遷移前的序列加以訂正。圖4b為七角井站訂正后的溫度序列，1998年以前的訂正差值為2.15℃，少于遷移后的溫度差值3℃，這說明該訂正方法考慮了氣候變化的因素。但是統(tǒng)一增加同一個訂正值是否合適，還需檢驗訂正結(jié)果。

237次遷移的平均訂正值小于(大于)0的有152次(85次)，即這些臺站遷移后溫度降低(升高)了，需要將遷移前的序列減去(加上)一定的值(圖4b)。表3為訂正值絕對值的分布?？梢?，約33%的遷移訂正差值超過0.5℃，這些遷移對序列的影響較明顯。

將訂正后的714站結(jié)果用SNHT方法檢驗序列的均一性。經(jīng)檢驗，還存在72站78次因遷移造成的非均一性，其中38次是經(jīng)過訂正還存在的非均一性，還需要進一步訂正，而其余40次是新增的因遷移造成的非均一性。

237次因遷移造成的非均一點已經(jīng)訂正好199次，訂正均一的比例約為84%。雖然還存在78次因遷移造成的非均一性，但是訂正溫度差比上一次的訂正溫度差小了很多，此次的訂正差值絕對值超過0.5℃只有3次，沒有超過1℃的遷移訂正差值，所以均一性訂正效果較好。38次新檢驗出來的非均一性序列主要是由于訂正影響了部分檢驗臺站的參考臺站，這說明參考臺站建立的序列只能近似為均一的，難以達到理論上的均一性，參考臺站數(shù)目越少，受單站的影響越大，增加參考臺站數(shù)量可以減少個別非均一臺站的影響。

表3 年平均氣溫訂正差值大小及其次數(shù)分布Table 3 Adjusted values of annual mean temperature and their corresponding frequency distribution

用遷移前后回歸系數(shù)的差異來表示訂正對序列的影響。由圖6可以看出，訂正差值較大的臺站主要分布于中國西部，這些地區(qū)地形復(fù)雜，遷移對臺站數(shù)據(jù)影響較大，許多臺站從高海拔地區(qū)遷移到低處，訂正后溫度變化趨勢明顯降低。當(dāng)然訂正后序列趨勢變化大小不僅與遷移前后位置變動的大小有關(guān)，還與遷移發(fā)生的時間點以及訂正的序列長度關(guān)系較大，遷移之前或者之后序列長度短，對趨勢的影響都較小。

5 結(jié)論和討論

本文在整理原有753站逐日觀測數(shù)據(jù)的基礎(chǔ)上，分析了氣象臺站的遷移狀況，討論了遷移對中國年平均氣溫序列均一性的影響，并對遷移造成的非均一序列進行訂正，得到如下主要結(jié)論:

1)首先將39個業(yè)務(wù)代替的臺站合并起來，作為一次遷移考慮，得到714站。通過經(jīng)緯度以及海拔高度差異的分析發(fā)現(xiàn)，只有148站沒有遷移，其他566站共有950次遷移，601次遷移的水平距離在10 km內(nèi)，192次遷移超過20 km;950次遷移中有533次遷移的海拔高度沒有變化，713次遷移的海拔高度變化低于10 m，93次高于50 m;249次遷移的水平距離超過20 km或(和)海拔高度變化超過50 m，中國臺站遷移明顯。t檢驗結(jié)果表明，對于年平均溫度，950次遷移共造成201站235次遷移前后的溫度差異顯著。

圖6 訂正前后的年平均溫度趨勢差值分布(單位:℃/a;方框表示趨勢減少，實點表示趨勢增加)Fig．6 Distribution of annual mean temperature trend difference before and after adjusting(units:℃/a;panes(solid points)denote decreasing(increasing)tendency)

2)SNHT檢驗發(fā)現(xiàn)，204站的237個非均一點是由遷移造成的，有79個非均一點溫度訂正差值超過0.5℃。t檢驗結(jié)果和SNHT檢驗結(jié)果表明，有99次遷移前后平均溫度差異顯著，且序列是非均一的，而其他138次遷移，盡管遷移前后的溫度差異不顯著，但同樣序列是非均一的。t檢驗?zāi)芎芎玫嘏c臺站遷移點結(jié)合起來，但難以簡單地通過單站平均溫度的差異判斷序列的均一性，而SNHT檢驗考慮了周圍臺站氣候變化趨勢的影響，能較好地判斷序列的均一性，但是如果沒有元數(shù)據(jù)的支撐，難以確定非均一點的合理性。

3)采用差值訂正方法對序列進行訂正，經(jīng)過訂正，遷移造成的非均一性能夠得到較大的改善，84%由遷移造成的非均一序列能夠調(diào)整得相對均一，其他序列也能得到相應(yīng)改善。比較逐月平均溫度訂正的差異可發(fā)現(xiàn)，不同月份的月平均氣溫訂正值并不一致，對更小尺度序列的訂正需進一步深入研究。通過經(jīng)緯度確定的臺站遷移可能存在部分誤差，其他的臺站沿革對均一性的影響也需要考慮，而均一性檢驗中受參考臺站尤其是西部觀測臺站稀疏等因素的影響，相關(guān)系數(shù)并不能都達到理想狀態(tài)的0.8，參考序列并不能都達到理想狀態(tài)的均一，加上統(tǒng)計檢驗的不確定性以及部分臺站城市化等因素對溫度影響的差異，結(jié)果的應(yīng)用性有待于進一步檢驗，建立準確及時有效的均一性氣候序列還是一項長期工作。

高曉榮，李慶祥，董文杰．2008．五臺山站歷史氣候資料的均一性分析［J］．氣象科技，36(1):112-118．

江志紅，黃群，李慶祥．2008．近50年中國降水序列均一性檢驗與訂正研究［J］．氣候與環(huán)境研究，13(1):67-74．

李慶祥，劉小寧，張洪政，等．2003．定點觀測氣候序列的均一性研究［J］．氣象科技，31(1):3-10．

李慶祥，董文杰，李偉，等．2010．近百年中國氣溫變化中的不確定性估計［J］．科學(xué)通報，55(16):1544-1554．

李祥余，黃少鵬，葉紅，等．2010．廈門站氣溫非均一性訂正及其變化特征對比分析［J］．地理科學(xué)，30(5):798-801．

劉小寧．2000．我國40年年平均風(fēng)速的均一性檢驗［J］．應(yīng)用氣象學(xué)報，11(1):27-34．

宋超輝，孫安健．1995．非均一性氣溫氣候序列訂正方法的研究［J］．高原氣象，14(2):215-220．

宋超輝，劉小寧，李集明．1995．氣溫序列非均一性檢驗方法的研究［J］．應(yīng)用氣象學(xué)報，6(3):289-296．

田紅，王雙五，魯?。?008．江淮流域年降水量和溫度均一性檢驗［J］．氣象科學(xué)，28(2):227-231．

屠其璞，王俊德，丁裕國，等．1984．氣象應(yīng)用概率統(tǒng)計學(xué)［M］．北京:氣象出版社:486-506．

王鈺，黃少平．2008．江西省地面氣象臺站沿革及其對觀測資料均一性的影響［J］．氣象與減災(zāi)研究，31(4):43-47．

魏鳳英．2007．現(xiàn)代氣候統(tǒng)計診斷與預(yù)測技術(shù)［M］．北京:氣象出版社:27-35．

吳增祥．2005．氣象臺站歷史沿革信息及其對觀測資料序列均一性影響的初步分析［J］．應(yīng)用氣象學(xué)報，16(4):461-467．

么枕生，丁裕國．1990．氣候統(tǒng)計［M］．北京:氣象出版社:776-826．

張弦，吳必文，嚴平，等．2009．合肥氣象站遷址對氣溫觀測的影響［J］．安徽農(nóng)學(xué)通報，15(3):67-68．

Alexandersson H．1986．A homogeneity test applied to precipitation data［J］．Int J Climatol，6:661-675．

Alexandersson H，Moberg A．1997．Homogenization of Swedish temperature data．PartⅠ:Homogeneity test for linear trends［J］．Int J Climatol，17:25-34．

Hawkins D M．1977．Testing a sequence of observations for a shift in random location［J］．J Amer Statist Assoc，3:180-185．