楊舒楠 , 孟慶濤 , 周寧芳 , 趙 威

1.國家氣象中心， 北京 100081 2.中國氣象局-河海大學水文氣象研究聯(lián)合實驗室， 北京 100081

0 引 言

國家“一帶一路”倡議的提出和推進，對我國“全球監(jiān)測、全球預報、全球服務”氣象業(yè)務提出了極為迫切的需求。近年來，越來越多的學者開始關注全球的氣候及氣象災害特征研究(張井勇等,2018；金紅梅等,2019；劉甜等,2019)。高溫天氣是影響“一帶一路”沿線國家或地區(qū)的主要氣象災害之一，而部分國家的防災減災能力較弱，且近幾十年來高溫極端天氣事件增多趨勢明顯(Della-Marta et al,2007；Tanarhte et al,2015)，因此，對全球高溫天氣的實時監(jiān)測和準確預報具有重要意義。

針對高溫極端天氣的研究一直備受氣象學家關注。Morak等(2013)利用觀測和氣候模擬數(shù)據(jù)對全球極端高溫頻率進行研究，指出幾乎所有季節(jié)均可看到高溫頻率的增加。Stone等(2000)、Degaetano和Allen(2002)、Rusticucci和Barrucand(2004)、Dash和Mamgain(2011)分別就美國、印度、阿根廷和加拿大等地區(qū)的極端高溫天氣或氣候事件進行了深入研究。國內(nèi)專家針對中國不同區(qū)域的高溫分布、變化特征及區(qū)域差異等的研究工作也很多(李艷等,2014；李縱橫等,2015；陳穎等，2016)，而近年來也逐漸將研究擴展至“一帶一路”沿線國家或地區(qū)(Feng et al,2018；金紅梅等,2019)。

在統(tǒng)計研究中，高溫天氣的定義指標有很多種。如絕對值定義：以35 ℃、37 ℃或40 ℃等作為高溫閾值進行分析(史軍等,2009；Ding et al,2010)，但此方法僅適用于日最高氣溫分布相對均勻的地區(qū)，無法用于全球。為方便不同地區(qū)相互對比，國內(nèi)外更常見的方式是采用某百分位值作為高溫閾值(Joneset al,1999；Yan et al,2002；Zhang et al，2005；李慶祥和黃嘉佑,2011)，一般使用90百分位值作為高溫指標，但具體百分位計算及閾值劃分方式有多種(陳穎等,2016)。Zhang等(2005)對比分析了不同概率分布近似條件下，多種統(tǒng)計時間窗的百分位高溫閾值計算方法。李慶祥和黃嘉佑(2011)利用站點日最高氣溫資料對中國高溫閾值進行探討，指出基于實際樣本頻率分布的高溫閾值計算方法，對中國高溫閾值確定更為合理。

面對全球高溫災害性天氣，中國的全球氣象業(yè)務目前尚無高溫天氣標準。世界氣象組織建議以32 ℃為高溫熱浪標準，但該閾值顯然并不適用于全球大多數(shù)地區(qū)，這致使熱帶地區(qū)高溫日數(shù)過多(張井勇等,2018)。全球大多數(shù)國家都有各自的高溫或高溫熱浪標準，各國閾值差別非常大，如中國以35 ℃作為高溫天氣標準，北歐多國將閾值定在25 ℃，印度的高溫標準則高達40 ℃。由于全球氣溫分布具有不均勻性，無法采用同一數(shù)值作為高溫閾值，因此，亟需建立全球分區(qū)域的高溫監(jiān)測、預報閾值標準，以推動我國全球氣象業(yè)務發(fā)展。

文中利用具有一定時間長度和空間覆蓋率的全球地面觀測數(shù)據(jù)，探討全球高溫閾值計算的方法和結果，旨在為全球分區(qū)域高溫閾值制定和高溫監(jiān)測預報業(yè)務提供理論基礎，也為區(qū)域防災減災提供參考依據(jù)。限于非洲觀測數(shù)據(jù)的質(zhì)量，文中僅涉及陸地范圍(不包含非洲)的高溫閾值探討。

1 資料選取與整理

1.1 數(shù)據(jù)選取

文中所用的日最高氣溫資料來自NOAA的全球地面日值數(shù)據(jù)(Global Surface Summary of the Day，GSOD)，該數(shù)據(jù)由全球地面逐小時觀測數(shù)據(jù)集(The Integrated Surface Hourly dataset，ISH)計算所得，是目前數(shù)據(jù)量較全、較穩(wěn)定的全球站點觀測資料集之一。數(shù)據(jù)集起始時間為1929年，其中1973年以后的數(shù)據(jù)相對較完整，到目前每天有包含大約9 000多站的日值觀測資料，其數(shù)據(jù)計算統(tǒng)計區(qū)間為00：00—23：59 UTC。

GSOD數(shù)據(jù)來自地面氣象觀測數(shù)據(jù)，避免了衛(wèi)星反演數(shù)據(jù)、分析數(shù)據(jù)等在計算和差值過程中產(chǎn)生的誤差。事實上，中國的高溫災害性、極端性天氣監(jiān)測預報業(yè)務中，閾值和標準等的建立均采用氣象站觀測資料作為統(tǒng)計依據(jù)(李慶祥和黃嘉佑，2011)。通過對比GSOD數(shù)據(jù)和CPC、ERA等分析資料的日最高氣溫數(shù)據(jù)發(fā)現(xiàn)，與觀測數(shù)據(jù)相比，分析數(shù)據(jù)雖然平均誤差較小，但其累積概率分布在很多站點均存在較大差異(圖略)，尤其在高溫一側(cè)，這將給高溫閾值的計算帶來誤差。因此，文中采用地面觀測數(shù)據(jù)展開高溫閾值的探討。

1.2 資料整理及高溫數(shù)據(jù)集構建

由于全球不同國家或地區(qū)的氣象業(yè)務發(fā)展水平參差不齊，其參與全球交換的氣象數(shù)據(jù)在時間連續(xù)性和空間覆蓋率等方面均存在較大差異，因此在進行統(tǒng)計前需對資料進行整理。

首先，對1973年以后的日最高氣溫觀測資料序列進行整理，綜合考慮有效觀測站點數(shù)及觀測數(shù)據(jù)的到報率(有效數(shù)據(jù)日數(shù)占統(tǒng)計時段總?cè)諗?shù)的比例)。圖1顯示，GSOD資料序列中日最高氣溫觀測從1973年以后較為完整，有效觀測站點維持在7 000余個，到2012年迅速增至8 000余個。雖然1973—2012年總觀測站點數(shù)變化不大，但觀測數(shù)據(jù)的穩(wěn)定性在不斷提升：20世紀70—80年代，日最高氣溫觀測的平均到報率僅為70%—75%，其中高到報率(≥95%，下同)站點僅3 000余個；進入90年代，數(shù)據(jù)到報率迅速增加，平均到報率從74%增長為80%以上，高到報率站點也穩(wěn)定在4 000個以上；進入21世紀，平均到報率繼續(xù)攀升至90%以上，同時高到報率站點數(shù)迅速增長，可見全球日最高氣溫觀測數(shù)據(jù)越來越穩(wěn)定。近年來，全球80%的觀測站點(超過7 000個)到報率均超過95%。

圖1 GSOD資料序列中日最高氣溫觀測總站點數(shù)、到報率≥95%站點數(shù)和站點平均到報率演變Fig. 1 Evolution of the number of daily maximum temperature observation stations, the number of stations with reported rate ≥95%, and station average reported rate in the GSOD data series

綜合以上分析，自20世紀90年代開始全球日最高氣溫觀測數(shù)據(jù)來源較為穩(wěn)定。在氣候統(tǒng)計分析中，通常以30年作為一個氣候基準，因此文中選用1990—2019年為統(tǒng)計時段。在此期間，GSOD數(shù)據(jù)有日最高氣溫記錄的站點共計12 650個(圖2)。但上述站點在統(tǒng)計時間段內(nèi)的有效記錄覆蓋時間長短差異明顯，大多數(shù)站點在20世紀90年代甚至21世紀初資料缺測現(xiàn)象仍十分嚴重，部分站點在2010年以后才開始有相對較穩(wěn)定的觀測(圖略)。從30 年平均到報率看(圖2)，其空間分布十分不均勻：在北美、歐洲中西部、中國、日韓等經(jīng)濟較發(fā)達地區(qū)，能夠獲得穩(wěn)定資料的站點密集程度要明顯高于全球其他地區(qū)，這些地區(qū)大多數(shù)站點的平均到報率都超過95%；與之相反，中亞、非洲和南美洲的站點到報率總體較低，甚至有的站點的到報率在10%以下。

圖2 全球1990—2019年日最高氣溫觀測站點的空間分布和平均到報率(單位：%)Fig. 2 Spatial distribution and average reported rate (units:%) of daily maximum temperature observation stations from 1990 to 2019

閾值統(tǒng)計需要觀測數(shù)據(jù)時間序列具有一定時間長度和穩(wěn)定的高到報率，因此需對統(tǒng)計時間段內(nèi)的觀測站點進行篩選：要求缺測資料日數(shù)要控制在總資料日數(shù)的5%以內(nèi)(孫軍和張福青,2017)，即到報率95%以上。但由于全球資料缺測嚴重、站點空間分布極不均勻，為保證觀測站點具有一定空間密度和覆蓋率，在統(tǒng)計時對某些地區(qū)適當放寬了要求，具體站點篩選方法如下：

首先，統(tǒng)計1990—2019年所有年份到報率均高于95%的站點(圖3中的“1”)。在此基礎上，增加僅1年到報率低于95%但高于80%(圖3中的“2、3”),或僅2年到報率低于95%但高于90%的站點(圖3中的“4”)。經(jīng)上述篩選，中亞、南亞、中東等地區(qū)的站點仍很少或幾乎沒有，因此需要對上述地區(qū)人工增選觀測站點。所增選的站點在某幾年的到報率低于95%，但其30年平均到報率高于95%，因此可用來進行閾值分析，但在氣候變化等趨勢性分析時需加以剔除，圖3中的“5”)。

圖3 全球篩選日最高氣溫觀測站點分布(總站點數(shù)量1 969個)Fig. 3 Distribution of selected global daily maximum temperature observation stations (the total number of stations is 1 969)

通過上述步驟，共計篩選出站點1 969個(圖3)，建立了較完整連續(xù)的全球1 969個站點30年逐日最高氣溫資料集。該資料集除非洲觀測站點缺乏、南美洲北部及南北極圈內(nèi)的站點較為稀疏外，全球其他地區(qū)的站點均具有較好的空間覆蓋率。

對于北半球高緯度地區(qū)(北美洲北部和俄羅斯北部)常年嚴寒，很多地區(qū)無人居住，氣象觀測站點非常稀少，其中大多數(shù)站點建立于21世紀以后，有些站點甚至是近10年才設立(圖4a)，因此數(shù)據(jù)的有效觀測長度非常短。南美洲北部的數(shù)據(jù)情況與之相仿，南美洲北部的廣袤地區(qū)為熱帶雨林氣候區(qū)，亞馬遜平原遼闊的熱帶森林人口稀少，每平方千米不足1人，氣象觀測發(fā)展緩慢，很多站點建站時間較晚，且觀測設備維護和數(shù)據(jù)傳輸?shù)葐栴}導致數(shù)據(jù)傳輸非常不穩(wěn)定，缺測現(xiàn)象嚴重(圖4b)。

非洲氣象觀測發(fā)展緩慢、觀測儀器缺乏維護、電子化信息化程度低，數(shù)據(jù)缺測程度非常嚴重。統(tǒng)計非洲觀測資料到報率結果顯示，其主要問題在于設備的維護和數(shù)據(jù)的傳輸。非洲觀測站點中，僅有不到0.5%的站點擁有較穩(wěn)定的高到報率(圖略)。在所有站點中約有36%的站點建立時間晚于1990年，其中部分站點建立于2000年之后；而其他63.5%的站點則存在數(shù)據(jù)穩(wěn)定性非常差的問題，甚至經(jīng)常連續(xù)數(shù)年沒有數(shù)據(jù)(圖4c)。從30年平均到報率看，除非洲北部少數(shù)站點外，其他地區(qū)均在50%以下，甚至很多站點低于10%(圖2)。即使近幾年，非洲的地面觀測也僅有33%的站點到報率能達到95%。鑒于以上原因，無法利用地面觀測資料開展高溫閾值統(tǒng)計。針對非洲的高溫閾值，將在其他研究中，依靠全球再分析資料和衛(wèi)星反演資料等進行探討和補充。

圖4 代表站點逐年到報率(a.俄羅斯北部和加拿大北部站點，b.南美洲北部站點，c.非洲站點)

2 高溫閾值計算方法

2.1 實際概率分布百分位閾值統(tǒng)計

文中選用90百分位數(shù)對應的溫度值作為該站點的高溫閾值。在傳統(tǒng)計算中，通常將統(tǒng)計時間段內(nèi)的日最高氣溫序列按從小到大的順序排列(t1，t2，t3，…，tn)，則某一百分位數(shù)對應的溫度值為

t0=(1-a)tm+atm+1

(1)

式中，m為溫度序列的序號；n為總樣本數(shù)；a為權重系數(shù)。實際計算時，m和a的計算有多種方法(Zhang et al,2005;李紅梅等,2008)，一般常取m=[p0(n+1)]，a=p0(n+1)-m，其中p0為所計算百分位值對應的概率，方括號表示數(shù)值取整(黃嘉佑，2016)。

傳統(tǒng)計算方法均對溫度序列的概率分布有所假設，而實際上日最高氣溫序列既不遵從均勻分布，也不完全遵從正態(tài)分布，所以嚴格意義上講，采用傳統(tǒng)的百分位閾值計算方法，會存在2 ℃左右的誤差(李慶祥和黃嘉佑，2011)。因此，為更精確地計算高溫閾值，參照李慶祥和黃嘉佑(2011)的方法，利用最高氣溫數(shù)據(jù)的實際累積概率分布來確定百分位閾值。計算時首先對30年的站點日最高氣溫時間序列進行排序，并確定樣本頻率分布的分組數(shù)量。組數(shù)計算公式為G=1+3.322lgn(黃嘉佑，2016)，其中n為總樣本數(shù)，文中取為30年夏季的總?cè)諗?shù)。因所選月份差異和站點數(shù)據(jù)缺測，不同站點的樣本數(shù)n值為2 660—2 760，該樣本數(shù)的變化對分組數(shù)量的影響不大，根據(jù)計算所得，分組組數(shù)近似取為12。

分組確定后，分別計算每組數(shù)據(jù)的頻率分布，然后利用序列的實際累積概率分布，通過線性差值的方式確定相應百分位的高溫閾值。此方法可較為真實地反映溫度序列的分布情況，具體計算方法：

t0=ati+(1-a)tj

(2)

(3)

當p0(90)百分位數(shù)落入某兩組數(shù)據(jù)(ti，tj)的累積頻率值(pi，pj)之間時，采用式(2)和式(3)線性差值求得相應的氣溫閾值。

2.2 統(tǒng)計時間窗

選取1990—2019年夏季為高溫閾值統(tǒng)計時間窗。在大多數(shù)研究中，北半球夏季一般定義為6—8月，而南半球夏季則定義為12月至次年2月。但實際上，受不同氣候類型、大氣環(huán)流特征、地形及洋流等影響，全球不同站點的氣溫月變化存在較大差異(圖5)。如果僅按南、北半球區(qū)分夏季，所計算的高溫閾值將會存在較大誤差。

圖5 代表站點月平均日最高氣溫分布(a.瑞士日內(nèi)瓦，b.澳大利亞堪培拉，c.日本東京，d.老撾萬象，e.印度新德里，f.菲律賓達沃)Fig. 5 Monthly mean daily maximum temperature distribution of representative stations (a.Geneva, Switzerland; b.Canberra, Australia; c.Tokyo, Japan; d.Vientiane, Laos;e.New Delhi, India; f.Davao, Philippines)

對最熱三個月統(tǒng)計(圖6)顯示，北(南)半球溫帶和寒帶大部分站點以6—8月(12月至次年2月)的日最高氣溫為峰值(如瑞士日內(nèi)瓦和澳大利亞堪培拉，圖5a、b)；但日韓、中國東南沿海、歐美沿海等地區(qū)，則以7—9月為最熱(如日本東京，圖5c)；南亞、東南亞等熱帶地區(qū)，全年氣溫均較高，最熱的月份往往集中在3—5月或4—6月的干季(如老撾萬象和印度新德里，圖5d、e)；而某些熱帶地區(qū)站點，全年日最高氣溫的月較差非常小，最熱三個月甚至并不連續(xù)(如菲律賓達沃，圖5f)。

圖6 全球1 969個氣象觀測站最熱月份(三個月)分布(數(shù)值字代表月份)Fig. 6 The hottest three months of global 1 969 meteorological observation stations (the number represents each month)

考慮全球站點地理位置的特殊性，選取兩種夏季定義方式：第一種直接選取站點真實最熱三個月(圖6)為夏季；第二種對北半球中高緯度地區(qū)夏季定義為6—8月，南半球中高緯度地區(qū)夏季定義為12月至次年2月，而低緯地區(qū)及地形影響較大區(qū)域夏季定義為真實最熱三個月。經(jīng)過計算，此兩種夏季定義方式的計算高溫閾值相差較小，大多數(shù)站點相差不到0.5 ℃，1 969個站的均方差僅0.15 ℃，因此文中主要就第一種方法的計算結果進行討論。

3 計算結果分析

3.1 不同日最高氣溫值的百分位分布

在對全球高溫閾值進行討論前，首先對目前WMO和中國現(xiàn)有高溫標準在氣候統(tǒng)計中的百分位情況進行比較。WMO建議的高溫熱浪標準為日最高氣溫高于32 ℃；中國業(yè)務高溫標準為最高氣溫高于35 ℃，其中全國高溫預警的發(fā)布則以37 ℃高溫的影響范圍為參考。因此，針對上述三種閾值的百分位進行了計算(圖7)。

圖7 32 ℃(a)、35 ℃(b)和37 ℃(c)對應的1990—2019年夏季日最高氣溫氣候統(tǒng)計百分位值(單位：℃)

圖7a顯示，對于歐洲大陸45°N以北的大部分地區(qū)，32 ℃的最高氣溫對應90或以上百分位；45°N以南地區(qū)，32 ℃大多對應85—95百分位，且隨緯度降低而遞減。WMO選用32 ℃作為高溫建議閾值，與其總部日內(nèi)瓦所在地理位置(約46°N)有很大關系，其附近站點32 ℃對應90—95百分位。對于亞洲大陸，情況較為復雜：俄羅斯、蒙古、日本北海道和中國東北地區(qū)及青藏高原高海拔地區(qū)，32 ℃的最高氣溫對應90及以上百分位，但中亞、中東、阿拉伯半島及南亞東南亞等地區(qū)，32 ℃均對應80以下百分位。對于中國35 ℃高溫標準(圖7b)，除高海拔站點外，中國大部分站點35 ℃均對應90以上百分位，有些站點甚至達95百分位。對比整個亞歐大陸，除中亞、阿拉伯半島、中東、印度半島和中南半島外，其他地區(qū)35 ℃均對應90以上百分位。

從WMO和中國高溫標準對應的百分位看，選用夏季90百分位數(shù)作為閾值計算指標是較為合理的，這也與大多數(shù)高溫極端天氣或高溫熱浪研究相符(Zhang et al,2005；李慶祥和黃嘉佑,2011；李艷等,2014)。對于高溫閾值的選取，既要反映大氣的極端性，又要具有一定的代表性，即高溫事件不能太多，也不能太少。對于更高的百分位數(shù)，其能夠體現(xiàn)溫度的更極端狀態(tài)，但以此確定的閾值過高，年高溫日數(shù)將明顯變少。以95百分位為例，在統(tǒng)計上，其年平均高溫日數(shù)僅有4.5 d(90 d×5%)，高溫日數(shù)明顯偏少。比如對于中國，95百分位對應的溫度值僅在江淮和江南地區(qū)之間為35 ℃左右，但在此區(qū)域以北及以南的大部分地區(qū)，溫度值都要超過37 ℃，與中國的現(xiàn)行高溫標準有明顯出入。

對比亞歐區(qū)域不同溫度值的百分位分布(圖7)發(fā)現(xiàn)，歐洲大陸的最高氣溫百分位數(shù)較亞洲略偏大，因此其高溫閾值偏低，閾值確定可參考WMO的高溫標準。亞洲大陸中東部的最高氣溫百分位分布特征與中國相近，閾值確定可參考中國的35 ℃高溫標準；東南亞地區(qū)溫度偏高，37 ℃的最高氣溫對應90百分位(圖7c)；而對于印度半島、中亞、阿拉伯半島等地，37 ℃對應90以下百分位(圖7c)，因此需要更高的高溫閾值。

3.2 90百分位高溫閾值

圖8給出了全球1 969個觀測站夏季(實際最熱三個月)日最高氣溫90百分位對應的高溫閾值計算結果。從全球范圍看，高溫閾值總體呈從赤道向兩極遞減趨勢，閾值最高的區(qū)域位于阿拉伯半島到非洲北部的沙漠地區(qū)(高于40 ℃)，閾值最低的區(qū)域集中在南、北極圈內(nèi)。

圖8 全球1 969個氣象觀測站1990—2019年夏季日最高氣溫90百分位值(單位：℃)Fig. 8 90 percentile of daily maximum temperature (units：℃) at global 1 969 stations during summer from 1990 to 2019

受氣候特征、海陸差異、海拔高度等影響，即使相同緯度的不同地區(qū)閾值差異也較大。例如，同樣位于25°N附近的站點，高溫閾值為27—48 ℃(圖9)：阿拉伯半島(30°—60°E)受熱帶沙漠氣候影響，高溫閾值可達44—48 ℃，站點平均45.5 ℃；印度半島北部(66°—85°E)處于熱帶季風氣候區(qū)，高溫閾值為41—44 ℃；在95°—110°E附近，25°N線貫穿中國云貴高原，此區(qū)域站點海拔普遍高于1 000 m，高溫閾值僅為27—31 ℃，站點平均值為29.8 ℃；110°—120°E之間為中國江南、華南地區(qū)，受季風氣候影響，高溫閾值為35—37 ℃；而西北太平洋島嶼(124°—130°E)受海洋氣候影響，高溫閾值僅為32—33 ℃。

此外，受海陸差異影響，沿海站點的高溫閾值普遍低于內(nèi)陸。例如，北歐挪威海沿岸高溫閾值(19—22 ℃)較同緯度東歐平原北部(26—28 ℃)低6—7 ℃(圖9)；北美阿拉斯加半島沿海高溫閾值較加拿大東北部低5—6 ℃。全球不同區(qū)域地理地形和氣候差異顯著，高溫閾值空間分布極不均勻，以下將分區(qū)域?qū)ζ溥M行探討。

圖9 (24°—26.3°N；30°—135°E)區(qū)域氣象觀測站日最高氣溫90百分位值和海拔Fig. 9 90 percentile high temperature threshold and stations’ altitude of daily maximum temperature observation stations between 24°N and 26.3°N；30°E and 135°E

亞洲地區(qū)，不同地區(qū)或氣候條件下，高溫閾值存在顯著的差異。日本、韓國等東亞國家受東亞季風氣候影響，高溫閾值為30—36 ℃(圖8)。中南半島和印度半島雖然位于近似緯度，但印度大陸面積廣闊，北部有青藏高原阻擋冷空氣南下，地形起伏閉塞，氣溫更易升高，在夏季形成印度熱低壓中心(朱乾根等，2007)，加上印度半島的印度季風爆發(fā)晚于中南半島的南海季風，因此印度半島(尤其北部)高溫閾值明顯高于中南半島。中南半島高溫閾值分布較均勻(37—40 ℃),而印度半島高溫閾值為38—45 ℃，且內(nèi)陸地區(qū)由南向北遞增。印度尼西亞、馬來西亞群島及其他太平洋熱帶島嶼國家雖更接近赤道，但受熱帶雨林氣候影響，高溫閾值普遍為32—35 ℃。地處歐亞大陸腹地的中亞地區(qū)氣候為典型的沙漠和大陸性氣候，雨水稀少、溫度極高，其與阿拉伯半島高溫閾值均較高，基本都在40 ℃以上。其中，阿拉伯半島中部高溫閾值為42—48 ℃，為全球高溫閾值之最。

歐洲地區(qū)，大陸區(qū)域的高溫閾值總體呈從南向北遞減趨勢，沿海地區(qū)閾值低于內(nèi)陸，其閾值大值中心主要位于伊比利亞半島中東部和巴爾干半島中部。這兩個地區(qū)緯度相對較低且均受地中海氣候影響，高溫閾值分別為34—40 ℃和34—37 ℃。此外，歐洲45°N以北地區(qū)高溫閾值普遍低于32 ℃(圖8)，其中，大不列顛島和北歐地區(qū)的閾值為20—25 ℃，其他地區(qū)為27—32 ℃;45°N以南地區(qū)的高溫閾值變化較大，為26—40 ℃。

北美大陸主要受溫帶大陸和溫帶海洋性氣候影響，高溫閾值總體呈由低緯向高緯遞減、沿海低于內(nèi)陸的特征(圖8)。加拿大高溫閾值普遍低于32 ℃，閾值最低的站點位于北極圈內(nèi)的巴芬島和格陵蘭島。北美高溫閾值的最大值(40—44 ℃)出現(xiàn)在美國西南部，此區(qū)域為北美莫哈維沙漠和索諾拉沙漠地區(qū)，受熱帶沙漠氣候影響，溫度顯著高于北美其他地區(qū)。墨西哥到中美洲一帶，氣候特征復雜多樣，高溫閾值為33—40 ℃。

南美洲中北部地區(qū)，高溫閾值隨緯度增加而遞減的趨勢并不明顯，除沿海和西部高海拔地區(qū)外，總體分布較均勻，都集中在34—37 ℃，中部地區(qū)局地可達39—40 ℃(圖8)。南美西部的高海拔和沿海地區(qū)，受高原山地氣候和冷洋流影響，高溫閾值僅為23—27 ℃。

大洋洲氣候分布較為復雜，高溫閾值空間分布極不均勻，閾值跨度較大(圖8)。澳大利亞大陸北部和東部分別受草原氣候和季風氣候控制，氣溫較適宜，高溫閾值為32—35 ℃；而澳大利亞中西部為廣袤的沙漠，在熱帶沙漠氣候影響下，高溫閾值普遍為40—41 ℃。塔斯馬尼亞島和新西蘭為溫帶海洋性氣候，高溫閾值為20—28 ℃。而所羅門群島至庫克群島一帶的太平洋島嶼，受海洋性氣候影響，高溫閾值為31—34 ℃。

3.3 代表站點對比分析

為更直觀對比不同氣候影響下，全球不同區(qū)域的日最高氣溫分布特征及計算高溫閾值，并與現(xiàn)有的高溫標準進行對比，選取全球9個代表站點進行分析，其地理位置、主要影響氣候和高溫閾值如表1所列。

表1 全球不同區(qū)域代表站點信息

分析代表站點日最高氣溫的30年概率密度函數(shù)分布(圖10)發(fā)現(xiàn)，不同區(qū)域、不同氣候背景下，日最高氣溫的分布具有明顯差異，且溫度分布也與典型的正態(tài)分布有一定差別，這是文中選用實際概率分布來計算高溫閾值的原因，也是對不同區(qū)域和氣候類型需要分別考慮不同高溫閾值的原因。

圖10顯示，對于溫帶地區(qū)的4個站點(北京、日內(nèi)瓦、華盛頓、瑞典，圖10a、b、f、i)，其日最高氣溫概率密度函數(shù)大值的分布區(qū)間較寬泛，有的站點會有雙峰或多峰特點。從閾值分布看，溫帶地區(qū)除高海拔及沿海高溫閾值較低、沙漠氣候區(qū)閾值較高之外，其高溫閾值主要受緯度影響較大(圖8)，其中中緯度地區(qū)的溫帶季風、溫帶大陸氣候影響地區(qū)較廣，閾值大多為35—37 ℃(如北京、華盛頓)；越過45°N，高溫閾值迅速降低至32 ℃及以下(如瑞士)，緯度較高地區(qū)，閾值通常為25—27 ℃(如瑞典)。

圖10 1990—2019年全球不同區(qū)域代表站點日最高氣溫概率密度函數(shù)分布(紅線為各站點的相應計算90百分位高溫閾值)(a.54511號站，中國北京； b.6700號站，瑞士日內(nèi)瓦； c.42182號站，印度新德里； d.94926號站，澳大利亞堪培拉； e.40437號站，沙特阿拉伯利雅得； f.72405號站，美國華盛頓； g.83612號站，巴西大坎普； h.48665號站，馬來西亞馬六甲； i.2464號站，瑞典斯德哥爾摩)Fig. 10 Statistical probability distribution of daily maximum temperature at representative stations in different regions of the world from 1990 to 2019 (the red line is the corresponding calculated high temperature threshold of each station) (a.54511, Beijing, China; b.6700, Geneva, Switzerland; c.42182, New Delhi, India; d.94926, Canberra, Australia; e.40437, Riyadh, Saudi Arabia; f.72405, Washington, USA; g.83612, Grand Camp, Brazil; h.48665, Malacca, Malaysia; i.2464, Stockholm, Sweden)

對于熱帶地區(qū)的站點(新德里、利雅得、大坎普、馬六甲，圖10c、e、g、h)，全年氣溫普遍較高，因此其峰值非常集中。尤其熱帶雨林氣候區(qū)(圖10h)，全年高溫多雨，最高氣溫分布集中在29—35 ℃(累積概率達85%以上，圖略)。由于熱帶地區(qū)溫度普遍較高，且概率分布集中，因此該地區(qū)的高溫閾值受緯度影響較小，氣候特征對閾值的影響占主要地位(圖8)。例如，位于熱帶雨林氣候區(qū)的馬六甲高溫閾值僅為34.85 ℃；位于熱帶草原氣候區(qū)的巴西中南部，高溫閾值普遍在37 ℃左右(大坎普36.50 ℃)；位于熱帶季風氣候區(qū)的印度半島高溫閾值超40 ℃(新德里42.90 ℃)；而熱帶沙漠地區(qū)溫度閾值為全球最高(如阿拉伯半島、非洲北部及澳大利亞中西部)，高溫閾值均非常高，其中沙特阿拉伯的利雅得，閾值達45.90 ℃。

澳大利亞堪培拉位于澳大利亞東南部，主要受亞熱帶季風氣候影響，夏熱冬溫，四季分明，其日最高氣溫分布的偏態(tài)分布特征明顯(圖10d)，在較高溫度一側(cè)擁有更高的概率分布。其在緯度上與中國江淮地區(qū)相對應，溫度閾值也較為接近。

計算高溫閾值與國內(nèi)外常用標準對比顯示，北京的計算高溫閾值為35.34 ℃，與中國實行的業(yè)務標準35 ℃較為吻合。瑞士日內(nèi)瓦的計算閾值為31.82 ℃，非常接近WMO建議的高溫熱浪標準32 ℃。印度對高溫熱浪的定義為平原地區(qū)高溫標準40 ℃，瑞典等北歐國家則以25 ℃為高溫標準，文中計算的新德里和斯德哥爾摩高溫閾值較其本國標準略偏高，這是因為選取站點為較大城市，高溫現(xiàn)象更明顯。

分析全年日最高氣溫概率密度函數(shù)分布(圖10)發(fā)現(xiàn)，北京的計算高溫閾值所對應概率密度約為1.4%，計算所得高于該高溫閾值的累積概率約2.9%。與之相對比，其他代表站點的高溫閾值所對應的概率密度及高于該高溫閾值的累積概率分別為1.0%—1.7%和2.2%—3.8%，與北京站的計算值較為接近。而北京站的高溫閾值又與實際業(yè)務標準非常吻合。在此種高溫閾值下，各個站點30年的高溫日數(shù)約為241—416 d，平均每年約為8—14 d，此計算高溫閾值可以較好地反映全年溫度的極端性，且與中國的業(yè)務標準和習慣較為吻合。

4 結論與討論

選用1990—2019年GSOD全球地面觀測數(shù)據(jù)，在對長時間序列日最高氣溫資料分析基礎上，篩選出滿足統(tǒng)計需求的全球1 969個觀測站，建立相對完整連續(xù)的逐日最高氣溫資料集，探討全球(不包括非洲)高溫閾值的統(tǒng)計方法和計算結果，得到：

1) 利用站點實際最熱三個月的日最高氣溫實際累積概率分布，采用90百分位數(shù)作為高溫閾值統(tǒng)計方法，可以較好反映全球高溫天氣，也與大多數(shù)高溫極端天氣或高溫熱浪研究相符。對比國內(nèi)外現(xiàn)有高溫標準以及最高氣溫概率密度函數(shù)分布等，計算最高氣溫閾值能較好地反映全年氣溫的極端性，且與中國的業(yè)務標準和習慣較為吻合。

2) 全球高溫閾值總體從赤道向兩極遞減，最高值位于阿拉伯半島到非洲北部的熱帶沙漠地區(qū)，最低值集中在南、北極圈內(nèi)，同時受氣候特征、海陸差異、海拔高度等影響，同一緯度不同地區(qū)的閾值也會存在顯著差異，海洋附近站點及高海拔地區(qū)閾值較同緯度其他地區(qū)偏低。溫帶地區(qū)高溫閾值受緯度分布影響最明顯，熱帶地區(qū)受氣候類型影響更顯著。

3) 亞洲高溫閾值區(qū)域差異顯著，東亞季風氣候區(qū)閾值為30—36 ℃，中南半島37—40 ℃，印度半島38—45 ℃，中亞及阿拉伯半島均在40 ℃以上，其中阿拉伯半島可達42—48 ℃。歐洲大陸與北美洲的高溫閾值總體從南向北遞減、沿海低于內(nèi)陸。歐洲閾值大值中心位于伊比利亞半島和巴爾干半島(34—40 ℃)，北美則出現(xiàn)在美國西南部的沙漠氣候區(qū)(40—44 ℃)。歐洲45°N以北及加拿大閾值普遍低于32 ℃。南美洲最高氣溫閾值受緯度影響不明顯，除沿海和西部高海拔地區(qū)外(23—27 ℃)，高溫閾值分布總體較均勻(34—37 ℃，局地超過39 ℃)。大洋洲閾值分布極不均勻，澳大利亞北部和東部閾值32—35 ℃，而中西部沙漠地區(qū)閾值為40—41 ℃。此外，亞洲與大洋洲的太平洋熱帶島嶼國家高溫閾值普遍在31—35 ℃。

文中高溫閾值計算基于全球站點實況觀測數(shù)據(jù)，對于數(shù)據(jù)缺測嚴重的非洲大部分地區(qū)，尚缺乏有效的統(tǒng)計站點，因此需在未來的研究中利用再分析資料或衛(wèi)星反演溫度資料等進行討論。此外，為方便全球最高氣溫監(jiān)測預報業(yè)務和服務的需要，需在一定空間范圍內(nèi)實現(xiàn)閾值的統(tǒng)一。在進行閾值區(qū)域統(tǒng)一時，需綜合參考該區(qū)域的計算高溫閾值、氣候類型、國家或地區(qū)的地域范圍及當?shù)刈罡邭鉁貥藴实?，這部分工作有待今后進一步深入研究。