謝毅文, 李 娟，陳偉榮，羅世豪

(1. 水資源與水電工程科學(xué)國家重點實驗室，湖北 武漢430072；2. 東莞理工學(xué)院 化學(xué)與環(huán)境工程學(xué)院，廣東 東莞523808；3. 北京師范大學(xué)水科學(xué)研究院，北京100875；4. 中國環(huán)境科學(xué)研究院，北京100012；5. 東莞市海匯環(huán)?？萍加邢薰?，廣東 東莞523000)

1959-2013年珠江流域平均氣溫時空變化特征*

謝毅文1,2, 李 娟3,4，陳偉榮2，羅世豪5

(1. 水資源與水電工程科學(xué)國家重點實驗室，湖北 武漢430072；2. 東莞理工學(xué)院 化學(xué)與環(huán)境工程學(xué)院，廣東 東莞523808；3. 北京師范大學(xué)水科學(xué)研究院，北京100875；4. 中國環(huán)境科學(xué)研究院，北京100012；5. 東莞市海匯環(huán)保科技有限公司，廣東 東莞523000)

根據(jù)1959-2013年珠江流域43個氣象站點月平均氣溫數(shù)據(jù)，采用線性傾向估計、M-K突變分析和小波分析等方法分析了珠江流域平均氣溫的時空變化特征。結(jié)果表明：珠江流域在統(tǒng)計時間段的平均氣溫均呈上升趨勢，且絕大部分區(qū)域上升顯著(P<0.05)；年平均氣溫序列于1993年存在顯著突變，其余時間尺度突變主要發(fā)生在20世紀(jì)90年代初期；年平均氣溫序列主周期為3.6 a，1963-1970年期間存在2.4～3.4 a的變化周期；1983-2008年期間存在2.8～4.4 a的變化周期。流域內(nèi)年平均氣溫范圍為10.68～22.60 ℃，平均值為19.50 ℃，從西部到東部總體呈增加趨勢；絕大多數(shù)氣象站點在不同時間尺度內(nèi)的平均氣溫均呈上升趨勢，其中珠江三角洲地區(qū)和東江流域受城市化影響平均氣溫上升尤為顯著。

平均氣溫；時空變化；升溫趨勢；全球變暖；珠江流域

近年來，因全球變暖引起的洪澇和干旱災(zāi)害不斷增多，嚴(yán)重影響人類生態(tài)環(huán)境和社會經(jīng)濟發(fā)展[1]。全球氣溫升高是眾多氣象災(zāi)害事件發(fā)生的重要驅(qū)動因素之一，已引起世界范圍內(nèi)的廣泛關(guān)注[2-3]。據(jù)IPCC第五次評估報告指出，全球氣候變化程度有可能比目前人類所認(rèn)識到的還更嚴(yán)重，20世紀(jì)50年代以來的變化劇烈程度是近千年以來所罕見的[4]。從有詳細(xì)氣象記錄以來的19世紀(jì)50 年代開始，在剛剛過去的3個10 a中，每一個10 a都刷新了平均氣溫最高的紀(jì)錄； 1980-2012年期間則可能是北半球自1400年以來最熱的30 a。1880-2012年，全球海陸表面平均溫度呈線性上升趨勢，約升高了0.85 ℃；2003-2012年平均溫度比1850-1900年平均溫度上升了0.78 ℃，因此，全球氣溫變暖已是不爭的事實[5]。全球變暖將影響地球的水循環(huán)系統(tǒng)，導(dǎo)致地球干濕更分明、海水升溫、冰川消融、北半球春季積雪減少等，并由此引發(fā)出一系列災(zāi)害。

在全球氣溫升高的背景下，了解區(qū)域氣溫的變化規(guī)律對及時預(yù)測預(yù)警災(zāi)害性天氣顯得尤為重要。中國為此也開展了一系列探索區(qū)域氣溫變化規(guī)律的研究，以期為氣候災(zāi)害預(yù)報提供有效參考依據(jù)，如屠其璞等[6]發(fā)現(xiàn)中國平均氣溫的升溫趨勢為每百年0.76 ℃，與北半球陸地平均氣溫的相關(guān)系數(shù)為0.78，且中國平均氣溫和北半球平均氣溫目前都處在最高溫期；曾小凡等[7]分析了長江流域年平均氣溫空間變化特征，發(fā)現(xiàn)自1991年以來全流域平均氣溫均為上升趨勢，其中中下游地區(qū)和金沙江流域平均氣溫上升幅度最大；王瑩等[8]發(fā)現(xiàn)在全球變暖背景下，近50 a廣西全區(qū)年平均氣溫和年極端最高、最低平均氣溫呈現(xiàn)總體上升趨勢，其中年平均氣溫和年極端最低平均氣溫的上升趨勢特別顯著；張雪芹等[9]發(fā)現(xiàn)中國干旱區(qū)1961-2007年間的年平均氣溫呈明顯上升趨勢，約上升了1.8 ℃，平均每10 a氣溫升溫速率為0.39 ℃，增暖幅度與波動幅度均高于同期全球以及中國區(qū)域平均水平。上述研究表明，我國大部分地區(qū)平均氣溫與全球變暖的大背景相一致，均呈現(xiàn)上升趨勢。

珠江是中國7大江河之一，流域下游珠江三角洲地區(qū)經(jīng)濟發(fā)達(dá)，是中國最為重要的經(jīng)濟增長極之一[10]。珠江流域是中國南方水循環(huán)最活躍的流域之一，然而，隨著全球氣候變暖引起的流域氣象災(zāi)害也逐漸頻繁，主要表現(xiàn)在兩方面：① 上游地區(qū)降雨偏少，導(dǎo)致這些地區(qū)經(jīng)常遭受干旱災(zāi)害威脅[11]；② 下游地區(qū)水循環(huán)變得更加活躍，導(dǎo)致這些地區(qū)經(jīng)常產(chǎn)生極端降雨并形成洪澇威脅[12-13]。可見，全球氣溫升高給珠江流域社會經(jīng)濟發(fā)展帶來了嚴(yán)峻的挑戰(zhàn)。為了積極應(yīng)對因氣溫升高導(dǎo)致氣候異常而帶來的一系列問題，探討珠江流域氣溫變化規(guī)律顯得非常重要。鑒于此，本文采用1959-2013年珠江流域43個氣象站提供的月平均氣溫數(shù)據(jù)，分析流域平均氣溫的時空變化特征并探討其發(fā)展規(guī)律，以期為流域災(zāi)害預(yù)警預(yù)報、災(zāi)害管理、防災(zāi)減災(zāi)等工作提供科學(xué)的參考依據(jù)。

1 流域概況與數(shù)據(jù)來源

珠江流域(102°14′E-115°53′E; 21°31′N-26°49′N) 面積約47.6 萬km2，由西江、北江、東江、珠江三角洲諸河組成。珠江流域地處熱帶和亞熱帶季風(fēng)氣候區(qū)，多年平均氣溫在14～22 ℃，多年平均降水量1 525 mm，年內(nèi)降水多集中在4-9月，約占全年水量的80%。

站點數(shù)據(jù)源于中國氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http:∥cdc.cma.gov.cn/home.do)。珠江流域共有國家級氣象站點48個，本文選擇其中數(shù)據(jù)系列完整并通過R軟件檢測的站點共43個(如圖1)，對1959-2013年逐月平均氣溫數(shù)據(jù)，通過R軟件中的RClimDex 程序處理，以達(dá)到檢測和提高數(shù)據(jù)質(zhì)量的目的。檢驗處理主要包括以下幾個方面錯誤記錄檢查內(nèi)容：① 月最低溫>月最高溫；②記錄值嚴(yán)重偏離本地區(qū)氣象實際情況，即超出3 倍標(biāo)準(zhǔn)差的值定義為出界值。通過人工檢查同相鄰站點的記錄進(jìn)行比對，對合理的數(shù)據(jù)進(jìn)行保留，不合理的則按缺測進(jìn)行處理，而缺測數(shù)據(jù)采用3次樣條函數(shù)內(nèi)插補齊。為了方便分析，根據(jù)珠江流域內(nèi)各獨立的水系特征及地形地貌相似的原則，把流域劃分為7個子區(qū)域：區(qū)域1包含南盤江和北盤江流域，區(qū)域 2 包括紅水河及柳江流域，區(qū)域 3 包含左江、右江及郁江流域，區(qū)域 4 包含賀桂江、黔潯江及西江干流流域，區(qū)域 5 為北江流域，區(qū)域 6 為東江流域，區(qū)域 7 為珠江三角洲，具體如圖1所示。

圖1 珠江流域地形概況及其氣象站點分布圖Fig.1 Terrain situation and meteorological stations of Pearl River basin

2 研究方法

珠江流域的年平均氣溫在年內(nèi)不同時間尺度呈現(xiàn)明顯的差異特征，因此本文選取年、四季、汛期與非汛期共7個時間尺度，并對于時間尺度作如下劃分：春季(3-5月)，夏季(6-8月)，秋季(9-11月)，冬季(12月-翌年2月)，汛期劃分為4-9月，非汛期劃分為10月-翌年3月。對于年均氣溫的時間變化主要采取趨勢分析、突變分析和周期分析，空間變化則采取Kriging插值法。

根據(jù)月平均氣溫統(tǒng)計出各年平均氣溫后，利用線性傾向估計方法構(gòu)建年平均氣溫值與年份之間的線性關(guān)系，以分析年平均氣溫的時間趨勢變化，計算步驟詳見文獻(xiàn)[14]。突變分析采取Mann-Kendall法(簡稱M-K法)分析空間模態(tài)時間變化特征，該方法又稱無分布檢驗，其優(yōu)點是不需要樣本遵從一定的分布，也不受少數(shù)異常值的干擾，可檢驗?zāi)昶骄鶜鉁亻L時間序列的變化趨勢，定量反映變化趨勢的顯著性，計算步驟詳見文獻(xiàn)[14]。周期分析則采取小波分析，它是建立在經(jīng)典的傅里葉分析和窗口傅里葉變換的基礎(chǔ)上發(fā)展起來的，可以分析出時間序列周期變化的局部特性，具體計算方法詳見文獻(xiàn)[15-16]。

在年平均氣溫空間格局研究中，由于基于地統(tǒng)計學(xué)的Kriging插值方法有較好的擬合效果[17]，利用ArcGIS的地統(tǒng)計模塊將43個站點的多年平均氣溫和變化速率進(jìn)行空間內(nèi)插，得到空間連續(xù)分布的珠江流域多年年平均氣溫及其變化速率圖。

3 研究結(jié)果分析

3.1 多時間尺度變化特征

運用線性傾向估計方法，分析珠江流域43個站點1959-2013年之間55 a平均氣溫數(shù)據(jù)，并計算其及七個子區(qū)域多時間尺度趨勢線斜率slope值和相關(guān)系數(shù)R值。對于年平均氣溫、春季平均氣溫、夏季平均氣溫、秋季平均氣溫和汛期平均氣溫，時間序列長達(dá)55 a，以0.05顯著水平進(jìn)行檢驗查得R0.05,55=0.265 61；冬季以及非汛期平均氣溫時間序列長度為54 a，以0.05顯著水平進(jìn)行檢驗查得R0.05,54=0.268 09。流域及其子區(qū)域slope值以和信度水平檢驗結(jié)果如表1所示?？芍?，珠江流域及其子區(qū)域在大部分時間尺度內(nèi)的平均氣溫均通過了通過95%的信度水平顯著性檢驗，表明在對應(yīng)的時間尺度內(nèi)平均氣溫上升顯著。

表1 子區(qū)域及珠江流域平均氣溫趨勢線斜率1)

1)“*”為通過95%的信度水平顯著性檢驗

圖1為7個子區(qū)域和珠江流域年平均氣溫趨勢變化圖?？芍髯訁^(qū)域和流域均呈上升趨勢，且均通過了0.05置信水平檢驗，表明在1959-2013年55 a間年平均氣溫顯著地上升，其中子區(qū)域7上升速率最高(每10 a為0.19 ℃)，而區(qū)域2上升速率最低(每10 a為0.111℃)。珠江流域最大年均氣溫發(fā)生在1998年的20.39為℃，最小發(fā)生在1984年的18.69 ℃。從平均氣溫年紀(jì)變化來看，珠江流域在1960年代中期到1980年代中期平均氣溫下降到一個較低水平，此后分別在1980年代中期至1990年代中期、1990年代中期至2010年代初期平均氣溫回升到較高水平。

珠江流域四季平均氣溫變化趨勢如圖3(a-d)所示，可知珠江流域四季溫度均呈上升趨勢，其中春季平均氣溫變化沒有通過0.05置信水平檢驗，上升趨勢不顯著；而夏秋冬平均氣溫均通過置信水平檢驗，上升趨勢顯著(P<0.05)。各個季節(jié)中，冬季平均氣溫上升速率最大(每10 a為0.28 ℃)，秋季次之(每10 a為0.24 ℃)。對于子區(qū)域，代表南盤江和北盤江流域的區(qū)域1以及珠江三角洲所在的區(qū)域7同樣是春季平均氣溫呈不顯著上升趨勢，但夏秋冬平均氣溫上升趨勢顯著(P<0.05)，且冬季上升速率最高；其他子區(qū)域春季和冬季平均氣溫上升不顯著，夏季和冬季上升趨勢顯著(P<0.05)。

珠江流域汛期與非汛期平均氣溫變化趨勢如圖3(e-f)所示，整個珠江流域在汛期和非汛期平均氣溫均呈現(xiàn)上升趨勢，且均通過0.05信度檢驗，而非汛期上升速率較大(每10 a為0.17 ℃)。7個子區(qū)域在汛期平均氣溫均呈顯著上升趨勢(P<0.05)，其中代表珠江三角洲的區(qū)域1上升速率最高(每10 a為0.15 ℃)；7個子區(qū)域非汛期平均氣溫也呈上升趨勢，除區(qū)域2和區(qū)域3未通過0.05置信水平檢驗外其余均通過檢驗，而非汛期平均氣溫上升速率最高發(fā)生在代表珠江三角洲的區(qū)域1。

由圖2和圖3可知，在1959-2013年期間珠江流域平均氣溫呈上升趨勢，平均每10 a上升速率為0.134 ℃。本文的上升趨勢與王兆禮等的結(jié)果基本一致，但后者的每10 a平均上升速率為0.10 ℃，顯然，計算時間序列較短(1961-2000年)是導(dǎo)致二者差異的主要原因[18]。

外界條件以及氣象環(huán)流特征是引起珠江流域平均氣溫上升的主要原因。文獻(xiàn)[19-20]研究發(fā)現(xiàn)500 hPa高度場在1970年代末后呈現(xiàn)波動式上升，中緯度環(huán)流由以前的經(jīng)向環(huán)流為主逐漸演變成緯向環(huán)流為主，而此時副熱帶高壓同樣也進(jìn)入偏強期，這種環(huán)流的演變使得進(jìn)入及影響我國的冷空氣次數(shù)偏少、強度偏弱，從而導(dǎo)致珠江流域的氣溫明顯偏高，尤其是近10 a增加趨勢更加明顯。劉燕等[21]發(fā)現(xiàn)冬季西風(fēng)帶系統(tǒng)活動減弱，西風(fēng)槽經(jīng)向度減小使得冷空氣活動減弱，導(dǎo)致對珠江流域的影響也相應(yīng)減弱，從而有利于冬季溫度升高；同時華南低層偏南風(fēng)異常(冬季風(fēng)減弱)，珠江流域上空水汽含量增加而削弱向外長波輻射，不利于夜間輻射降溫而導(dǎo)致夜間最低溫度的顯著增暖。另外副熱帶高壓加強西伸使得華南受副高控制機會增加，副高控制下的下沉氣流利于珠江流域溫度顯著升高；南海夏季風(fēng)強度減弱，導(dǎo)致大氣中水汽含量減少、比熱降低，因而白天吸收短波輻射時更容易導(dǎo)致升溫。

3.2 突變分析

為分析年均氣溫序列的突變情況，采用M-K法對區(qū)域1-7及整個珠江流域的平均氣溫序列進(jìn)行突變分析，發(fā)現(xiàn)7個子區(qū)域和珠江流域的年均氣溫序列均沒有發(fā)生顯著的突變。以珠江流域為例，圖4顯示了流域的年均氣溫突變情況：整個珠江流域UF曲線表現(xiàn)出先降后穩(wěn)定上升的趨勢，UF和UB兩曲線在臨界值±1.96(α=0.05)之間僅在1993年處有交點，表明存在突變點；UF曲線在2004年以后超出臨界線，表明珠江流域在2004年后年平均氣溫上升趨勢顯著(P<0.05)。再統(tǒng)計、分析珠江流域四季、汛期與非汛期尺度平均氣溫突變情況如表2所示，可知珠江流域除了夏季突變發(fā)生時間較早(1983)、春季沒有發(fā)生突變外，其余時間尺度的平均氣溫突變主要發(fā)生在20世紀(jì)90年代初期，且這些突變呈上升趨勢。

圖2 珠江流域及各區(qū)域年平均氣溫變化趨Fig.2 Variation trend of annual average temperature in Pearl River basin and the sub-regions

3.3 周期分析

為分析年均氣溫周期變化情況，對1959-2013年區(qū)域1-7及整個珠江流域的年均氣溫進(jìn)行小波分析。限于篇幅，僅以全流域的年均氣溫為例分析序列的周期特性。圖5為年平均氣溫小波分析圖，左側(cè)為連續(xù)小波變換圖，各包絡(luò)線為能量密度等值線；右側(cè)為全局小波功率圖，點線表示95%置信水平的紅噪聲檢驗曲線，若實線的峰值超過點線則表示對應(yīng)的周期顯著。1959-2013年珠江流域平均氣溫序列存在3.6，8和22.7 a的周期，其中3.6 a的周期通過了95% 的紅噪聲檢驗，故3.6 a為序列的主周期。具體來看，在 COI(Cone of influence，小波影響錐)區(qū)域內(nèi)小波功率譜能量密度的峰值主要集中在1963-1970年，存在2.4～3.4 a的變化周期，并通過顯著性檢驗；1983-2008年存在2.8～4.4 a的變化周期，通過了顯著性檢驗；在其他時間尺度，小波功率譜能量都比較小，均沒有通過95%置信水平的顯著性檢驗。

表2 珠江流域多時間尺度平均氣溫突變分析

圖3 珠江流域各時間尺度平均氣溫變化趨勢Fig.3 Average temperature variation trend of different time scales over Pearl River basin

圖4 珠江流域年平均氣溫突變分析圖Fig.4 Mutation analysis diagram of annual average temperature in Pearl River basin

再對珠江流域不同時間尺度的周期分析結(jié)果匯總于表3，由表可知，流域各時間尺度的主周期在2.6～3.6 a之間。

圖5 珠江流域年平均氣溫小波分析圖Fig.5 Wavelet analysis of annual average temperature over the Pearl River basin

表3 珠江流域多時間尺度平均氣溫周期分析

Table 3 Period analysis of average temperature with multi-time scale

序列序列年份周期/a主周期/a顯著年份全年1959-20132 8～4 4/2 4～3 43 61963-1970/1983-2008春季1959-20131 7～5 23 61963-2004夏季1959-20131 8～3 42 61964-1979/1990-2010秋季1959-20131 8～3/6 5～9 22 61962-2009/1974-2002冬季1959-20121 8～4 8/5 7～8 63 61964-1992/1971-1988汛期1959-20131 8～3 4/1 8～6 22 81973-1987/1993-2006非汛期1959-20122～4/2 8～9 53 11964-1981/1969-1996

3.4 年平均氣溫空間分布及其變化特征

計算各站點多年平均氣溫并通過Kriging插值法對全流域進(jìn)行插值，可得珠江流域年均氣溫空間分布圖(圖6a)。流域內(nèi)年平均氣溫范圍為10.68～22.60 ℃，平均值為19.50 ℃，最大值為深圳站的22.62 ℃，最小值為威寧站的10.64 ℃。東部、中部區(qū)域年平均氣溫較高而西部較低，從西部到東部呈總體增加趨勢。

利用1959-2013年珠江流域43個站點不同時間尺度的平均氣溫與時間序列建立線性相關(guān)關(guān)系，得到各站點對應(yīng)的相關(guān)系數(shù)。給定顯著性水平α= 0.05(置信度95%)，分析平均氣溫時間序列的趨勢變化情況和顯著程度。表4為各站點平均氣溫趨勢統(tǒng)計表，可知流域大多數(shù)站點的在不同時間尺度均呈上升趨勢。

圖6 珠江流域年均氣溫空間分布圖Fig.6 Spatial distribution characteristics of annual average temperature in Pearl River basin

圖7 珠江流域不同時間尺度平均氣溫變化速率空間分布圖Fig.7 Spatial distribution maps of average temperature change rates with different time scales in the Pearl River basin

表4 珠江流域各站點平均氣溫變化趨勢特征統(tǒng)計

Table 4 Trend characteristics of average temperature variation of the stations in Pearl River basin

時間尺度趨勢檢驗通過率/%上升個數(shù)總數(shù)顯著不顯著下降個數(shù)總數(shù)顯著不顯著全年86 042375101春季25 6391029413夏季88 441374211秋季83 741365202冬季39 5431726000汛期79 141338211非汛期65 1432815000

計算各個站點平均氣溫變化速率并利用Kriging插值法插值得到珠江流域平均氣溫變化趨勢空間分布特征，如圖7所示。對于全年尺度，共有42個站點(97.67%)呈上升趨勢，其中37個站點(86.05%)呈顯著上升(P<0.05)，表明珠江流域年平均氣溫呈總體上升趨勢，而流域東部及西南部上升尤為顯著。對于四季、汛期和非汛期時間尺度，大部分站點平均氣溫呈上升趨勢，僅有極個別站點呈下降趨勢。東部地區(qū)的區(qū)域6與區(qū)域7平均氣溫上升尤為顯著，主要是因為近幾十年來東江流域和珠江三角洲地區(qū)城市化進(jìn)程不斷推進(jìn)使得城市的熱島效應(yīng)越來越明顯，從而導(dǎo)致區(qū)域升溫趨勢顯著。

4 結(jié) 論

根據(jù)珠江流域43個氣象站1959-2013年逐月平均氣溫數(shù)據(jù)，采用線性傾向估計、M-K趨突變檢驗、小波周期分析和Krijing插值法等方法分析了珠江流域及其子區(qū)域平均氣溫時空變化特征，得出以下主要結(jié)論：

1)珠江流域及其子區(qū)域在1959-2013年間年平均氣溫呈總體上升趨勢，且絕大部分區(qū)域上升顯著(P<0.05)，其中珠江三角洲地區(qū)的上升速率最高，每10 a達(dá)到了0.19 ℃；四季、汛期與非汛期時間尺度的平均氣溫也有不同程度的升高，且大部分地區(qū)上升顯著(P<0.05)。

2)珠江流域及其子區(qū)域的平均氣溫均存在顯著上升突變，全流域突變點出現(xiàn)在1993年，且在2004年后年平均氣溫上升趨勢顯著(P<0.05)；四季、汛期與非汛期時間尺度平均氣溫突變主要發(fā)生在20世紀(jì)90年代初期，這些突變均呈上升趨勢。

3)珠江流域年平均氣溫序列的主周期為3.6 a。1963-1970年存在2.4～3.4 a的變化周期，1983 -2008年存在2.8～4.4 a的變化周期；其他時間尺度存在2.6～3.6 a的主周期。

4)流域內(nèi)年平均氣溫范圍為10.68～22.60 ℃，平均值為19.50 ℃，從西部到東部呈總體上升趨勢。

5)流域絕大多數(shù)氣象站點在不同時間尺度的平均氣溫呈上升趨勢，其中珠江三角洲地區(qū)和東江流域受城市化熱島效應(yīng)影響上升尤為顯著。

[1] 夏軍，邱冰，潘興瑤，等. 氣候變化影響下水資源脆弱性評估方法及其應(yīng)用[J].地球科學(xué)進(jìn)展，2012，27(4): 443-451.

[2] TOLLEFSON J. Global-warming limit of 2 degrees C hangs in the balance[J]. Nature, 2015,520(7545):14-15.

[3] WALLACE J M. Global warming and winter weather[J]. Science, 2014,343(6174):969.

[4] IPCC. Climate change 2013: the physical science basis [M]. Cambridge: Cambridg University Press,2014.

[5] 曾婷，楊東，郭佩佩，等. 1960-2012 年安徽省降水和氣溫變化特征及其與 ENSO 的相關(guān)性分析[J]．熱帶地理，2014，34(6): 783-793.

[6] 屠其璞，鄧自旺，周曉蘭. 中國近 117年年平均氣溫變化的區(qū)域特征研究[J]. 應(yīng)用氣象學(xué)報，1999,10(增刊):34-42.

[7] 曾小凡，翟建青，蘇布達(dá)，等. 長江流域年平均平均氣溫的時空變化特征[J]. 長江流域資源與環(huán)境，2009，18(5):427-431.

[8] 王瑩，蘇永秀，李政.1961-2010年廣西平均氣溫變化對全球變暖的響應(yīng)[J].自然資源學(xué)報，2013，28(10):1707-1717.

[9] 張雪芹，孫楊，毛煒嶧，等. 中國干旱區(qū)氣溫變化對全球變暖的區(qū)域響應(yīng)[J]. 干旱區(qū)研究，2010,27(4): 592-599.

[10] 賴成光，陳曉宏，王兆禮，等.珠江流域1960-2012年降雨侵蝕力時空變化特征[J].農(nóng)業(yè)工程學(xué)報，2015，31(8):159-167.

[11] 姚玉璧，張強，王勁松，等. 中國西南干旱對氣候變暖的響應(yīng)特征[J]. 生態(tài)環(huán)境學(xué)報 2014, 23(9): 1409-1417.

[12] 黃強，陳子燊. 全球變暖背景下珠江流域極端氣溫與降水事件時空變化的區(qū)域研究[J]. 自然科學(xué)進(jìn)展，2010,29(8): 956-967.

[13] 賴成光，陳曉宏，趙仕威，等. 基于隨機森林的洪災(zāi)風(fēng)險評價模型及其應(yīng)用[J]. 水利學(xué)報，2015,46(1):58-66.

[14] HUANG J, ZHANG J C, ZHANG Z X, et al. Spatial and temporal variations in rainfall erosivity during 1960-2005 in the Yangtze River basin [J]. Stoch Environ Res Risk Assess,2013, 27(2):337-351.

[15] GRINSTED A, MOORE J C, JEVREJEVA S. Application of the cross wavelet transform and wavelet coherence to geophysical time series [J]. Geophysics, 2004, 11(5-6): 561-566.

[16] TORRENCE C, COMPO G P. A practical guide to wavelet analysis [J]. Bull Am Meteorol Soc 1998, 79 (1):62-78.

[17] 杜國明，汪光松，吳超羽，等. 克里金在珠江河道地形空間數(shù)據(jù)內(nèi)插中的應(yīng)用[J]. 中山大學(xué)學(xué)報(自然科學(xué)版)，2007,46(1)：119-122.

[18] 王兆禮，陳曉宏，黃國如. 近40年來珠江流域平均氣溫時空演變特征[J]. 熱帶地理，2007,27(4):289-293+322.

[19] 慕巧珍，王紹武，朱錦紅，等. 近百年夏季西太平洋副熱帶高壓的變化[J].大氣科學(xué), 2001, 25 (6): 787-797.

[20] 林振敏，施能. 北半球冬季大氣環(huán)流遙相關(guān)型特征與我國區(qū)域氣候[J].氣象科技, 2004, 32(5): 333- 342.

[21] 劉燕，程正泉，葉萌. 廣州市氣溫的氣候變化特征及其成因分析[J]. 氣象，2008，32(2): 52-60.

Spatio-temporal variation of average temperature over the Pearl River basin during 1959-2013

XIEYiwen1,2,LIJuan3,4,CHENWeirong2,LUOShihao5

(1. State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan 430072, China; 2. College of Chemical and Environmental Engineering, Dongguan University of Technology, Dongguan 523808, China;3. College of Water Sciences, Beijing Normal University, Beijing 100875, China;4. Chinese Research Academy of Environmental Sciences, Beijing 100012, China;5. Dongguan Haihui Environmental Protection Science & Technology Limited Company,Dongguan 523000, China)

Based on the monthly average temperature data from 43 meteorological stations during 1959 to 2013, the spatiotemporal variation in the Pearl River basin was analyzed using the methods of linear regression, Mann-Kendall and wavelet analysis. The results show that: the average temperature showed an increasing trend in different time scales and the average temperature in most of the areas increased significantly (P<0.05). The annual average temperature series of the Pearl River basin had a mutation in 1993 but the other time scales had mutations mainly in the early 1990s. The primary period of the annual average temperature series was 3.6a and there were 2.4～3.4 years and 2.8～4.4 years change periods in 1963-1970 and 1983-2008 respectively. The value range of the average temperature of the basin was 10.68～22.60℃ with the average value of 19.50℃, which showed an increasing trend totally form the western to eastern. The average temperature of most meteorological stations of different time scales showed an increasing trend, which increased significantly in the Pearl River Delta and the Dongjiang River basin due to the influence of urbanization.

average temperature; spatiotemporal variation; warming trend; global warming; the Pearl River basin

10.13471/j.cnki.acta.snus.2016.03.005

2015-05-26

環(huán)保公益性行業(yè)科研專項基金資助項目(201309003)；水資源與水電工程科學(xué)國家重點實驗室開放基金資助項目(2013B106)；華南地區(qū)水循環(huán)與水安全廣東省普通高校重點實驗室開放基金資助項目(KLB09002)；廣東高?；で鍧嵣a(chǎn)與綠色化學(xué)品工程技術(shù)開發(fā)中心開放基金資助項目(201403)

謝毅文(1980年生)，男;研究方向：水資源水環(huán)境保護；通訊作者：李娟 ;E-mail:lijuan@craes.org.cn

P468.1

A

0529-6579(2016)03-0030-09