喀什站高空溫度資料均一性檢驗(yàn)與訂正

胡義成，張正陽，王秋香*

（1.中國(guó)氣象局烏魯木齊沙漠氣象研究所，新疆 烏魯木齊830002；2.新疆氣象信息中心，新疆 烏魯木齊830002）

目前存檔的探空資料中蘊(yùn)含了氣候變化的重要信息，但是由于探空資料中人為因素造成的不均一性和氣候因素造成的變化混淆，使得研究者很難利用探空資料準(zhǔn)確地判斷氣候變化趨勢(shì)。為了使探空資料很好地運(yùn)用到氣候研究中，很多學(xué)者致力于檢驗(yàn)和訂正探空資料中時(shí)間序列的不均一性。和地面氣象要素均一性檢驗(yàn)相類似的是[1-9]，對(duì)于探空資料不均一性的檢驗(yàn)需要識(shí)別由于人為因素造成的間斷點(diǎn)，估計(jì)間斷點(diǎn)偏離大小并對(duì)間斷點(diǎn)進(jìn)行訂正[10]。我國(guó)高空觀測(cè)站建站時(shí)間較早，最早出現(xiàn)在1950 年，大部分高空站從1960 年陸續(xù)開始建站，站點(diǎn)遍布全國(guó)。均一性研究在20 世紀(jì)90 年代開始進(jìn)行，但只有少數(shù)臺(tái)站進(jìn)行了試驗(yàn)性研究[10]。近期對(duì)全國(guó)多數(shù)探空溫度序列的均一性分析證實(shí)了輻射訂正、儀器換型等方法的修改造成了序列的非均一性，這種非均一性在各站探空資料中不同程度的存在，均一化過程也存在著相當(dāng)?shù)牟淮_定性[11，12]。

喀什高空氣象觀測(cè)站（以下簡(jiǎn)稱“喀什站”）經(jīng)歷了多次儀器變更、臺(tái)站位置變動(dòng)、觀測(cè)規(guī)范變更、觀測(cè)時(shí)間變化以及輻射訂正方法變化（以下簡(jiǎn)稱：輻射訂正）。這些活動(dòng)是否對(duì)高空資料的均一性產(chǎn)生了影響，高空資料是否可以連續(xù)使用，都是值得去研究和探討的問題。本文利用加拿大環(huán)境部氣候研究中心研發(fā)的RHtest[13，14]均一化檢驗(yàn)方法，對(duì)喀什站1958—2016 年各個(gè)標(biāo)準(zhǔn)等壓面層的年平均溫度資料進(jìn)行均一性檢驗(yàn)與訂正，對(duì)喀什站高空資料的完整性、連續(xù)性和均一性進(jìn)行了分析，有利于幫助臺(tái)站提高單站高空資料的質(zhì)量，也可以為新疆高空氣象觀測(cè)業(yè)務(wù)的運(yùn)行提供參考。

1 資料方法與觀測(cè)環(huán)境

1.1 氣象觀測(cè)環(huán)境變遷及現(xiàn)狀

喀什站始建于1954 年7 月，自1957 年6 月1日承擔(dān)著北京時(shí)間08、20 時(shí)常規(guī)探空和測(cè)風(fēng)的同步觀測(cè)任務(wù)。建站以來，喀什站址多次變動(dòng)。1955 年5月1 日因不明原因發(fā)生了站址位置變動(dòng)（圖1）；1957年7 月1 日，站址由疏附縣浩罕莊遷至喀什市西南郊；1964 年1 月1 日又因不明原因?qū)е抡局方?jīng)緯度信息發(fā)生變化；因2014 年1 月1 日觀測(cè)環(huán)境已遠(yuǎn)遠(yuǎn)不能滿足氣象觀測(cè)的要求。

圖1 歷次遷站過程圖

1.2 資料與方法

使用的基礎(chǔ)資料來自國(guó)家氣象信息中心“中國(guó)高空規(guī)定等壓面定時(shí)值數(shù)據(jù)集（V2.1）”及高空臺(tái)站元數(shù)據(jù)集。資料時(shí)段為喀什站1958—2016 年，包括北京時(shí)間08、20 時(shí)（世界時(shí)00 時(shí)和12 時(shí)）各規(guī)定等壓面位勢(shì)高度、氣溫、露點(diǎn)溫度觀測(cè)數(shù)據(jù)和風(fēng)向風(fēng)速。數(shù)據(jù)集中1951—2010 年數(shù)據(jù)來源于高空觀測(cè)紙質(zhì)月報(bào)表數(shù)字化資料，2010—2016 年12 月的數(shù)據(jù)來源于高空全月觀測(cè)數(shù)據(jù)文件（高空月報(bào)G 文件）。該數(shù)據(jù)集在研制過程中借鑒了《無線電探空資料質(zhì)量控制（QX/T 123-2011）》中相關(guān)質(zhì)量控制方法及參數(shù)，以多源數(shù)據(jù)對(duì)比檢測(cè)為主，對(duì)可疑資料實(shí)施人工核查與訂正。不僅更正了紙質(zhì)資料本身存在的觀測(cè)和抄錄錯(cuò)誤數(shù)據(jù)，還更正了信息化過程導(dǎo)致的資料錯(cuò)誤和無數(shù)據(jù)現(xiàn)象，最大程度地確保了研究數(shù)據(jù)的可靠性。

參考序列是ERA-Interim，ERA-Interim 是歐洲中期天氣預(yù)報(bào)中心（ECMWF）發(fā)布的最新的全球大氣數(shù)值預(yù)報(bào)再分析資料，它是繼早期資料（如ERA-15，ERA-40）后的新產(chǎn)品，目的是對(duì)ERA-40 等資料進(jìn)行完善，從而逐漸取代ERA-40。ERA-Interim所采用的同化系統(tǒng)是ECMWF 集成預(yù)報(bào)系統(tǒng)（IFS，32R2 版本的全球譜模式）。相對(duì)于ERA-40 在很多的關(guān)鍵部分都進(jìn)行了改進(jìn)，包括模式的改進(jìn)、四維空間變量分析的運(yùn)用、對(duì)衛(wèi)星資料變量偏差的修正以及對(duì)觀測(cè)系統(tǒng)中其他數(shù)據(jù)的處理等。

RHtestsV5 軟件包是在RHtestsV4 分位數(shù)匹配訂正（Quantile-Matching adjustments，QM）中增加使用參考序列的基礎(chǔ)上，在序列檢驗(yàn)中提供了一階滯后自相關(guān)，還提供了數(shù)據(jù)序列的均值訂正、QM 訂正和序列曲線圖以及結(jié)果的回歸模擬?？梢詫?duì)日或者小時(shí)數(shù)據(jù)序列進(jìn)行QM 訂正。檢驗(yàn)方法包括兩種，一種是懲罰最大t 檢驗(yàn)（PMT）方法，其檢驗(yàn)過程中需要建立參考序列，待檢序列與參考序列的差值（差值序列）是被檢驗(yàn)的對(duì)象；另一種是懲罰最大F 檢驗(yàn)（PMF），適用于無參考序列的檢驗(yàn)過程。

PMT 方法檢驗(yàn)過程如下：

假設(shè)｛Xt｝（t=1，…，N）代表一個(gè)ⅡD 正態(tài)分布序列，要檢驗(yàn)時(shí)間序列｛Xt｝是否存在斷點(diǎn)。

備擇假設(shè)：

μ1≠μ2，并且“ ｛ Xt｝~‖DN（μ，σ2）”代表｛ Xt｝遵循正態(tài)分布，其均值為μ，方差為σ2。當(dāng)Ha為真時(shí)，t=k 時(shí)的點(diǎn)被稱作間斷點(diǎn)，Δ=被稱作平均突變的大小，最可能的間斷點(diǎn)服從以下分布：

其中：P（κ）是建立的經(jīng)驗(yàn)性的懲罰因子，其建立方法可參考文獻(xiàn)[13]。

通過回歸檢驗(yàn)算法來檢驗(yàn)出多個(gè)間斷點(diǎn)，首先在t∈｛Nmin，Nmin+1，…，（N-Nmin）｝時(shí)段找出最可能的間斷點(diǎn)C0，Nmin指的是時(shí)間序列中被分割片段的最小長(zhǎng)度，也就是說，兩個(gè)相鄰斷點(diǎn)之間子序列的長(zhǎng)度或者第一（最后一）個(gè)斷點(diǎn)到序列起點(diǎn)（終點(diǎn)）之間的長(zhǎng)度，類似的，分別在時(shí)段t={Nmin，…，（C0-Nmin）}和t={（C0+1+Nmin），…，（N-Nmin）}找出次可能的間斷點(diǎn)C01和C03，接下來在C01和C03之間的時(shí)段找出最可能的間斷點(diǎn)C02，然后在C01，C02和C03之間尋找最可能的間斷點(diǎn)，即依次分段找出序列中各段最可能的間斷點(diǎn)，計(jì)算所有間斷點(diǎn)的統(tǒng)計(jì)量PTmax，估計(jì)lag-1 自相關(guān)和p值的大小，找出最大的PTmax值對(duì)應(yīng)的間斷點(diǎn)，如果是顯著的，該突變就被認(rèn)為是找到的第一個(gè)間斷點(diǎn)，該間斷點(diǎn)被記為NC=1。尋找該間斷點(diǎn)位置之后每段最可能的間斷點(diǎn)，估計(jì)其顯著性，找出下一個(gè)可能的間斷點(diǎn)，重復(fù)該過程，將間斷點(diǎn)按照顯著性由大到小排列，形成間斷點(diǎn)列表，判斷最小的間斷點(diǎn)是否顯著，當(dāng)不顯著時(shí)剔除該間斷點(diǎn)，需要再次評(píng)估剩余間斷點(diǎn)的顯著性，最終保留都統(tǒng)計(jì)顯著的間斷點(diǎn)即為序列檢驗(yàn)得到的間斷點(diǎn)，受分段檢驗(yàn)的影響，該方法能夠檢驗(yàn)的最短序列長(zhǎng)度是20，詳細(xì)過程見文獻(xiàn)[14]。

PMF 方法檢驗(yàn)過程如下：

假設(shè)εt代表變量，該變量的均值為0，方差為σ2，對(duì)于存在的線性趨勢(shì)β，其時(shí)間序列為{Xt}，要檢查在t=k 時(shí)刻是否存在一個(gè)平均突變，假設(shè)：

備選假設(shè)：

μ1≠μ2，當(dāng)Ha為真時(shí)，t=k 時(shí)的點(diǎn)被稱作間斷點(diǎn)。 Δ=被稱作平均突變的大小，最可能的間斷點(diǎn)服從以下分布：

其中，P（k）是建立的經(jīng)驗(yàn)性的懲罰因子，其建立方法可參考[14]。

其中μ^0和β0是在μ1=μ2=μ 時(shí)的估計(jì)值。通過回歸檢驗(yàn)算法來檢驗(yàn)出多個(gè)間斷點(diǎn)，首先在t={Nmin，Nmin+1，…，（N-Nmin）}時(shí)段找出最可能的間斷點(diǎn)C0。類似的，分別在時(shí)段t={Nmin，…，（C-Nmin）}和t={（C0+1+Nmin），…，（N-Nmin）}找出次可能的間斷點(diǎn)C01和C03，接下來在C01和C03之間的時(shí)段找出最可能的間斷點(diǎn)C02，然后在C01、C02和C03之間尋找最可能的間斷點(diǎn)，即依次分段找出序列中各段最可能的間斷點(diǎn)，計(jì)算所有間斷點(diǎn)的統(tǒng)計(jì)量PFmax，估計(jì)自相關(guān)和P 的大小，找出最大的PFmax值對(duì)應(yīng)的間斷點(diǎn)，如果是顯著的，該突變就被認(rèn)為是找到的第一個(gè)間斷點(diǎn)，該間斷點(diǎn)被記為NC=1。尋找該間斷點(diǎn)位置之后每段最可能的間斷點(diǎn)，估計(jì)其顯著性，找出下一個(gè)可能的間斷點(diǎn)，重復(fù)該過程，將間斷點(diǎn)按照顯著性由大到小排列，形成間斷點(diǎn)列表。

PMF 方法可以使用元數(shù)據(jù)：當(dāng)利用回歸檢驗(yàn)檢測(cè)出多個(gè)間斷點(diǎn)后，輸入元數(shù)據(jù)，再次判斷輸入的時(shí)間點(diǎn)是否存在間斷點(diǎn)（即依次分段找出序列中各段最可能的間斷點(diǎn)），計(jì)算所有間斷點(diǎn)的統(tǒng)計(jì)量，確定第一個(gè)間斷點(diǎn)；尋找該間斷點(diǎn)位置之后每段最可能的其他間斷點(diǎn)，分別估計(jì)其顯著性，找出下一個(gè)可能的間斷點(diǎn)，重復(fù)該過程，分步找出所有的間斷點(diǎn)。然后再將間斷點(diǎn)按照顯著性由大到小排列，形成新的有元數(shù)據(jù)支持的間斷點(diǎn)列表。判斷最小的間斷點(diǎn)是否顯著，當(dāng)不顯著時(shí)剔除該間斷點(diǎn)，再次評(píng)估剩余間斷點(diǎn)的顯著性，最終保留的統(tǒng)計(jì)顯著的間斷點(diǎn)即為序列檢驗(yàn)得到的間斷點(diǎn)。該方法能夠檢驗(yàn)的最短序列長(zhǎng)度是20[14]。

RHtests 方法已被人們廣泛地應(yīng)用于氣溫、降水、風(fēng)速及相對(duì)濕度等氣候序列的均一性分析中。因喀什站高空各層次日溫度資料有大量缺測(cè)值，且待檢序列（1958—2016 年）與對(duì)應(yīng)的ERA-Interim 參考序列（1979—2013 年）時(shí)間不一致，在對(duì)喀什站高空日、年溫度資料進(jìn)行均一性檢驗(yàn)和訂正的過程中發(fā)現(xiàn)，在利用RHtestsV5 軟件包進(jìn)行訂正時(shí)，大量缺省值的代入造成了訂正序列與實(shí)際觀測(cè)序列的偏差量過大，導(dǎo)致訂正序列不可信。故本文利用RHtestsV5 軟件包提供的PMT 方法，采用對(duì)應(yīng)的ERA-Interim 作為參考序列，對(duì)喀什站高空日溫度資料進(jìn)行均一性檢驗(yàn)和訂正；利用RHtestsV5 軟件包提供的PMF 方法，對(duì)喀什站高空年溫度資料進(jìn)行均一性檢驗(yàn)和訂正。

1.3 月觀測(cè)資料序列完整性

喀什站常規(guī)探空和測(cè)風(fēng)的同步觀測(cè)任務(wù)起始于1957 年6 月1 日，以有整年觀測(cè)數(shù)據(jù)（1958 年1 月）以來的平均溫度序列按照850、700、600、500、400、300、200 hPa 和100 hPa 共8 個(gè)規(guī)定等壓面的月值（共708 個(gè)月）來統(tǒng)計(jì)缺測(cè)率及數(shù)據(jù)缺測(cè)年月分布狀況，統(tǒng)計(jì)結(jié)果見表1 和表2，由表可見，喀什站8 個(gè)規(guī)定層1967 年9—10 月均缺測(cè)。

根據(jù)表1 和表2 給出的缺測(cè)統(tǒng)計(jì)結(jié)果，08 時(shí)8個(gè)層次的月平均溫度除100 hPa 的缺測(cè)率>3%，完整性較差外，其他層次完整性較好；20 時(shí)100 hPa的月平均氣溫完整性較差，200 hPa 的完整性一般，其他6 個(gè)層次的完整性較好。綜合分析認(rèn)為喀什站850、700、600、500、400 hPa 和300 hPa 共6 個(gè)層次的高空溫度觀測(cè)資料序列完整性較好。

表1 喀什站高空溫度序列月值缺測(cè)情況（08 時(shí)）

表2 喀什站高空溫度序列月值缺測(cè)情況（20 時(shí)）

2 觀測(cè)資料序列均一性分析

2.1 日觀測(cè)資料均一性分析

利用RHtestsV5 軟件包提供的PMT 方法，采用對(duì)應(yīng)的ERA-Interim 資料作為參考序列，分別對(duì)喀什站08 時(shí)、20 時(shí)對(duì)應(yīng)的850、700、500、400、300、200 hPa 和100 hPa 共14 條（因搜集的ERA-Interim資料中無對(duì)應(yīng)的600 hPa 序列，故此處沒有對(duì)喀什站600 hPa 的日氣溫序列進(jìn)行斷點(diǎn)檢測(cè)和分析。）溫度序列作置信度為0.95 的PMT 檢驗(yàn)。保留其中的兩類斷點(diǎn)，一類是在年或月序列中檢驗(yàn)出的顯著間斷點(diǎn)，并且有元數(shù)據(jù)的支持，如果該斷點(diǎn)出現(xiàn)的時(shí)間與臺(tái)站元數(shù)據(jù)記錄信息相差一年以內(nèi)，根據(jù)元數(shù)據(jù)記錄時(shí)間替換該斷點(diǎn)位置；另一類是沒有元數(shù)據(jù)支持，但是在日和月（或年）序列中被同時(shí)檢驗(yàn)出的相同年份的顯著間斷點(diǎn)，斷點(diǎn)位置依據(jù)日平均溫度序列被檢驗(yàn)出的斷點(diǎn)時(shí)間。

檢驗(yàn)得到的斷點(diǎn)信息顯示，喀什站14 個(gè)溫度序列被檢驗(yàn)出的斷點(diǎn)共有221 個(gè)（圖2），其中有元數(shù)據(jù)支持的斷點(diǎn)占89.14%。從檢測(cè)結(jié)果和元數(shù)據(jù)分析情況看，喀什站始建于1954 年7 月，而20 世紀(jì)50 年代以來中國(guó)高空觀測(cè)儀器和規(guī)范變化較大。1957 年4 月高空觀測(cè)時(shí)間從北京時(shí)間11 時(shí)和23 時(shí)變更為08 時(shí)和20 時(shí)。1959 年前后、1966 年以及1999 年7 月1 日后又經(jīng)歷了輻射訂正方法的變化，1962 年5 月31 日、1968 年10 月1 日、1970年7 月1 日和2005 年喀什站又陸續(xù)更換探空儀器。綜合分析認(rèn)為，喀什站溫度資料非均一性主要來源于臺(tái)站位置變動(dòng)、觀測(cè)儀器換型、觀測(cè)規(guī)范變更和輻射訂正。

圖2 喀什站溫度序列檢測(cè)出的斷點(diǎn)數(shù)量

通過PMT 方法和元數(shù)據(jù)信息相結(jié)合的方法來綜合判斷資料的斷點(diǎn)，所確定的斷點(diǎn)均通過了0.05的顯著性檢驗(yàn)。其中，有一部分比例的斷點(diǎn)雖然有元數(shù)據(jù)支持，但調(diào)查數(shù)據(jù)時(shí)發(fā)現(xiàn)由于各種原因，在斷點(diǎn)處、前或后的數(shù)據(jù)為缺測(cè)。如，14 條數(shù)據(jù)序列均在1967 年9 月12 日—1967 年11 月10 日缺測(cè)，而檢測(cè)出得到的1967 年9 月12 日和1967 年11 月10日兩個(gè)斷點(diǎn)及其間缺測(cè)的數(shù)據(jù)不進(jìn)行訂正。另外，在對(duì)日原始數(shù)據(jù)序列進(jìn)行訂正時(shí)，僅對(duì)有參考序列對(duì)應(yīng)的時(shí)段進(jìn)行訂正。通過判斷，20 時(shí)700 hPa、08時(shí)500 hPa 和20 時(shí)100 hPa 的非均一性最為明顯，分別檢測(cè)出32 個(gè)、26 個(gè)和25 個(gè)斷點(diǎn)（表3），經(jīng)查詢臺(tái)站元數(shù)據(jù)信息，這3 條序列斷點(diǎn)的元數(shù)據(jù)支持率分別為84.38%、57.69%和88%，綜合考慮以上各種情況，這3 條序列最終確定需要訂正的斷點(diǎn)分別為2 個(gè)、6 個(gè)和2 個(gè)，占所檢測(cè)出的斷點(diǎn)比例僅分別為6.25%、23.08%和8%，做到了絕不無根據(jù)地訂正資料。

表3 2 個(gè)觀測(cè)時(shí)次7 個(gè)標(biāo)準(zhǔn)等壓面上斷點(diǎn)數(shù)

通過檢驗(yàn)出的各規(guī)定層氣溫序列斷點(diǎn)的位置和訂正幅度來看，14 條序列均檢驗(yàn)出了斷點(diǎn)，除了08時(shí)100 hPa 以外，其他13 條序列均檢驗(yàn)出一個(gè)相同位置（1997-12-09）的斷點(diǎn)，經(jīng)查詢臺(tái)站元數(shù)據(jù)文件發(fā)現(xiàn)，在上世紀(jì)90 年代前后，喀什站經(jīng)歷了輻射訂正方法的變化。由此可以認(rèn)為，輻射訂正方法的變化導(dǎo)致13 條序列產(chǎn)生了間斷點(diǎn)。另外，08 時(shí)850~400、100 hPa 和20 時(shí)400 hPa 共7 條序列因臺(tái)站位置的變動(dòng)（1980-05-01 和1993-07-01）而產(chǎn)生了間斷點(diǎn)。在利用RHtestsV5 軟件包進(jìn)行斷點(diǎn)檢測(cè)和訂正時(shí)，序列中的缺測(cè)值需用缺省值“-999.99”來代替。在對(duì)各序列進(jìn)行訂正時(shí)發(fā)現(xiàn)，300 hPa 以上各序列的訂正幅度非常大，其中20 時(shí)100 hPa 的日氣溫序列訂正量高達(dá)-184 ℃，這么大的訂正量是不可信的。通過查找原因，發(fā)現(xiàn)300~100 hPa 這6 條日值序列的缺測(cè)率較高（表4），而08 時(shí)和20 時(shí)各序列訂正量的大小與對(duì)應(yīng)的序列缺測(cè)率之間呈顯著的正相關(guān)關(guān)系，相關(guān)系數(shù)分別高達(dá)0.97 和0.99，由此不難理解，大量缺省值的代入，會(huì)導(dǎo)致訂正序列與觀測(cè)序列的偏差過大。為了保證訂正序列的可信度，本研究?jī)H對(duì)缺測(cè)率小于1%的溫度日值序列進(jìn)行訂正，即僅對(duì)08 時(shí)和20 時(shí)850~400 hPa 共8 條序列進(jìn)行了訂正，具體訂正情況見表5。

2.2 年平均氣溫序列均一性分析

利用RHtestsV5 軟件包中的最大懲罰F 檢驗(yàn)（PMFT）方法對(duì)喀什站1958 年以來觀測(cè)的850~100 hPa 共14 條平均氣溫的年值序列進(jìn)行斷點(diǎn)檢驗(yàn)。該方法經(jīng)驗(yàn)性的考慮了時(shí)間序列的滯后一階自相關(guān)，并嵌入多元線性回歸算法，能夠用于檢驗(yàn)、訂正包含一階自回歸誤差數(shù)據(jù)序列的多個(gè)間斷點(diǎn)（平均突變），而且不需要使用參考序列。

表4 各規(guī)定等壓面日氣溫序列缺測(cè)率

表5 日溫度序列檢驗(yàn)出的斷點(diǎn)位置和訂正量

針對(duì)PMFT 方法檢測(cè)出的年序列中的顯著間斷點(diǎn)，基于喀什站詳實(shí)的高空歷史沿革記錄，挑出儀器變更、臺(tái)站位置變動(dòng)、觀測(cè)規(guī)范變更、觀測(cè)時(shí)間變化以及輻射訂正方法變化等有記錄的時(shí)間點(diǎn)（表6），來判斷在該時(shí)間是否為間斷點(diǎn)。如果該斷點(diǎn)出現(xiàn)的時(shí)間與臺(tái)站元數(shù)據(jù)記錄信息相差一年以內(nèi)，根據(jù)元數(shù)據(jù)記錄時(shí)間替換該斷點(diǎn)位置。若檢測(cè)出的間斷點(diǎn)不顯著，也沒有元數(shù)據(jù)支持，則舍棄該間斷點(diǎn)，也不進(jìn)行訂正，做到絕不無根據(jù)地訂正資料。

表6 喀什站存在變化的時(shí)間點(diǎn)及原因

年平均氣溫資料序列均一性檢驗(yàn)結(jié)果如圖3～8所示。

圖3 喀什站08 時(shí)100 hPa 年平均氣溫檢測(cè)結(jié)果

（不連續(xù)點(diǎn)：1997 年，輻射訂正方法變化）

圖4 喀什站08 時(shí)300 hPa 年平均氣溫檢測(cè)結(jié)果

圖5 喀什站08 時(shí)500 hPa 年平均氣溫檢測(cè)結(jié)果

圖6 喀什站08 時(shí)700 hPa 年平均氣溫檢測(cè)結(jié)果

圖7 喀什站08 時(shí)850 hPa 年平均氣溫檢測(cè)結(jié)果

（不連續(xù)點(diǎn)：1967 年，兩個(gè)月無觀測(cè)；1977 年，海拔高度變化、觀測(cè)手冊(cè)變化；1997 年，輻射訂正方法變化；2005 年，雷達(dá)和探空儀換型、L 波段系統(tǒng)運(yùn)行、數(shù)據(jù)文件換型、觀測(cè)規(guī)范變化）

圖8 喀什站20 時(shí)100 hPa 年平均氣溫檢測(cè)結(jié)果

喀什站14 條溫度序列中有6 條序列即08 時(shí)100a、300、500、700、850 hPa 和20 時(shí)100 hPa 存在不連續(xù)點(diǎn)，其中，08 時(shí)850 hPa 的不連續(xù)點(diǎn)高達(dá)4個(gè)，20 時(shí)除100 hPa 外，其他序列均無不連續(xù)點(diǎn)。結(jié)合元數(shù)據(jù)分析，造成喀什站觀測(cè)資料不連續(xù)的原因主要有：1967 年有近2 個(gè)月停止觀測(cè)，1974 年觀測(cè)規(guī)范變化，1977 年海拔高度變化和觀測(cè)手冊(cè)變化，1997 年輻射訂正方法的變化，2005 年雷達(dá)和探空儀換型、L 波段系統(tǒng)運(yùn)行、數(shù)據(jù)文件換型、觀測(cè)規(guī)范變化。其中，08 時(shí)850 hPa 年平均溫度在2005 年出現(xiàn)的斷點(diǎn)，很可能是由于L 波段與59-701 系統(tǒng)的感應(yīng)元件和測(cè)溫原理不同而存在的器差所致[15]。

2.3 年氣溫訂正幅度及訂正前后趨勢(shì)分析

經(jīng)過均一性分析和趨勢(shì)分析發(fā)現(xiàn)，喀什站14 條高空溫度序列中，08 時(shí)850、700、500、300、100 hPa和20 時(shí)100 hPa 共6 條序列不連續(xù)，需要訂正，其他8 條序列連續(xù)，則不需訂正。通過PMFT 方法和元數(shù)據(jù)信息相結(jié)合的方法確定斷點(diǎn)后，主要采用了RHtestsV5 軟件包中的平均值訂正法對(duì)最終被采納的斷點(diǎn)進(jìn)行了訂正，即以序列的最新資料為標(biāo)準(zhǔn)，把斷點(diǎn)前后的平均值差作為偏差，對(duì)斷點(diǎn)前序列進(jìn)行訂正，訂正結(jié)果見表7。最終采納的斷點(diǎn)數(shù)占PMFT方法檢測(cè)出斷點(diǎn)數(shù)的百分比為76.9%，也就是說有23.1%的斷點(diǎn)沒有被采納。

通過訂正，6 條氣溫序列訂正后與訂正前均呈現(xiàn)出相反的變化趨勢(shì)，其中，08 時(shí)100、300 hPa 和500 hPa 均為副訂正，而且訂正前氣溫均為下降趨勢(shì)，訂正后的氣溫均為上升趨勢(shì)；08 時(shí)700 hPa 和20 時(shí)100 hPa 則為正訂正，訂正后的氣溫均為下降趨勢(shì)。可以很明顯地發(fā)現(xiàn)訂正幅度隨著規(guī)定層高度的變化呈現(xiàn)出兩頭大中間小的現(xiàn)象，其中08 時(shí)和20 時(shí)100 hPa 的訂正幅度分別為-1.13 ℃和-0.80 ℃，而08 時(shí)850 hPa 的訂正幅度則高達(dá)2.56℃。08 時(shí)需要訂正的序列要多于20 時(shí)，08 時(shí)7 條序列中有6 條序列做了訂正，訂正率達(dá)到了85.71%，而20 時(shí)僅有100 hPa 做了訂正。

表7 年溫度序列檢驗(yàn)出的斷點(diǎn)位置和訂正量

3 結(jié)論

（1）喀什站08 時(shí)100 hPa 的高空月平均氣溫缺測(cè)率大于3%，完整性較差，20 時(shí)100 hPa 的完整性最差，缺測(cè)率高達(dá)11.44%；其他層次的完整性較好。

（2）喀什站14 條日溫度觀測(cè)序列共檢驗(yàn)出221個(gè)斷點(diǎn)，14 條年平均氣溫序列共檢驗(yàn)出10 個(gè)斷點(diǎn)。結(jié)合元數(shù)據(jù)分析認(rèn)為，喀什站高空溫度資料的非均一性主要來源于臺(tái)站位置變動(dòng)、觀測(cè)儀器換型、觀測(cè)規(guī)范變更和輻射誤差訂正方法的變化。

（3）利用RHtestsV5 軟件包對(duì)喀什站日平均氣溫序列進(jìn)行斷點(diǎn)檢測(cè)和訂正時(shí)發(fā)現(xiàn)，因缺測(cè)值需用缺省值“-999.99”來代替，大量缺省值的代入，會(huì)導(dǎo)致訂正序列與實(shí)際觀測(cè)序列的偏差過大。為了保證訂正序列的可信度，僅對(duì)缺測(cè)率小于1%的溫度日值序列進(jìn)行了訂正。

（4）喀什站高空年平均氣溫序列的訂正幅度隨著高度的變化呈現(xiàn)出兩頭大中間小的現(xiàn)象，08 時(shí)的訂正率明顯多于20 時(shí)，其中，08 時(shí)的訂正率達(dá)到了85.71%，而20 時(shí)僅有100 hPa 一條序列做了訂正。