羅紅霞， 戴聲佩， 李茂芬， 謝錚輝

(1.中國熱帶農(nóng)業(yè)科學院科技信息研究所/海南省熱帶作物信息技術應用研究重點實驗室，海南儋州 571737；2.農(nóng)業(yè)部農(nóng)業(yè)遙感重點實驗室，北京 101010)

聯(lián)合國政府間氣候變化專門委員會(Intergovernmental Panel on Climate Change，簡稱IPCC)第四次評估報告指出，過去100年地表溫度升高了0.74 ℃，變暖幅度自20世紀90年代以來明顯加速。氣候變暖已經(jīng)成為全世界人們關注的焦點，越來越多的學者加入到研究氣候變化的行列中[1]。海南島地處東亞大陸東南端，屬熱帶季風性氣候，年均氣溫22～26 ℃，年降水量1 600 mm以上，其中8、9月份降水最為充沛，素有“東方夏威夷”之稱。2010年國務院批準建立海南島為國際旅游島，使海南島的經(jīng)濟與文化發(fā)展速度更快，同時也對海南島的氣候環(huán)境提出了更高、更多的要求。了解海南島不同區(qū)域的氣候變化特征，充分利用氣候資源要素，對于指導農(nóng)業(yè)生產(chǎn)、合理開發(fā)島內(nèi)旅游資源等都有相當重要的作用。

近年來，有學者對海南島氣候變化特征開展了系列研究[2-10]，結(jié)果表明，海南島氣溫(平均氣溫、平均最高氣溫、平均最低氣溫)總體呈現(xiàn)上升趨勢；各氣候區(qū)降水量變化大，但總體仍呈增加趨勢，其中南部增長率最高；海南島平均風速、相對濕度、日照時數(shù)和蒸發(fā)量降低趨勢明顯。盡管相關學者已開展過關于海南島不同氣候區(qū)域的氣候變化特征研究，但這些研究較為深入的氣候因子主要是氣溫和降水，雖然對諸如日照時數(shù)、相對濕度因子有部分報道，卻主要集中于研究其變化趨勢，對于其突變情況及未來變化等并未開展研究，因此有必要對海南島各氣候區(qū)多氣候因子進行綜合探討。鑒于以上原因，本研究采用1959—2015年海南島7個氣象站點的平均氣溫、平均最低氣溫、平均最高氣溫、極端最低氣溫、極端最高氣溫、平均相對濕度、日照時數(shù)、降水量等8項氣候要素作為分析指標，應用之前海南島的氣候區(qū)劃成果，研究海南島的氣候變化特征，基于Mann-Kendall法分析各氣候要素的突變性，基于變標度極差分析法(rescaled range analysis，簡稱 R/S 分析)及Hurst指數(shù)試驗法對各氣候要素的未來變化趨勢作出一定的推斷及預測。

1 數(shù)據(jù)來源與方法

1.1 研究區(qū)域氣候區(qū)劃分

考慮到海南島復雜的地形氣候差異性，為了對海南島氣候變化進行區(qū)域差異性分析，筆者將海南島劃分為如圖1所示5個氣候區(qū)[6]：①東北區(qū)，屬于較濕潤的大風氣候區(qū)，主要位于瓊山、文昌沿岸臺地；②西北區(qū)，屬于較干燥的大風氣候區(qū)，主要位于儋州、臨高沿岸臺地；③中部山區(qū)，屬于涼爽濕潤氣候區(qū)，主要位于瓊中、昌江、保亭屯昌等市(縣)；④東南區(qū)，屬于常臺風、多雨、大風氣候區(qū)，主要位于瓊海、萬寧沿岸；⑤西南區(qū)，屬于干熱大風地區(qū)，主要位于東方、樂東、三亞等市(縣)。

1.2 數(shù)據(jù)來源及數(shù)據(jù)分析

數(shù)據(jù)來源于中國氣象數(shù)據(jù)網(wǎng)的海南島1959—2015年的月度數(shù)據(jù)，包括平均氣溫、平均最低氣溫、平均最高氣溫、極端最低氣溫、極端最高氣溫、平均相對濕度、日照時數(shù)、年降水量8項氣候要素。對于出現(xiàn)缺測或異常值的數(shù)據(jù)，采用線性回歸法進行插補。運用平均值方法形成8項氣候要素在1959—2015年時間序列的數(shù)據(jù)集。

1.3 研究方法

應用線性回歸分析、氣候傾向率揭示海南島的氣候變化特征；運用Mann-Kendall和滑動t檢驗法揭示氣候變化規(guī)律；最后采用R/S分析方法對海南島氣候的未來變化趨勢作出科學預測。

1.3.1 氣候傾向率 氣象要素x的趨勢變化用一次回歸方程表示，公式如下：

xi=a+bti。

(1)

式中：xi表示樣本量i的某一氣候變量；ti表示xi對應的時間；a表示回歸常數(shù)；b表示回歸系數(shù)，即線性趨勢項，b>0表示x隨著時間t的增加而增加，反之則隨時間t的增加而減小；b×10為氣候傾向率，表示某氣候變量每10年的變化速率。

1.3.2 Mann-Kendall突變檢驗 本研究采用Mann-Kendall非參數(shù)趨勢檢驗方法對氣候要素進行時間序列突變檢驗。在趨勢檢驗中，原假設H0表示數(shù)據(jù)集X的數(shù)據(jù)樣本獨立同分布，沒有趨勢存在?？蛇x假設H1表示數(shù)據(jù)集X存在1個單調(diào)的趨勢。Mann-Kendall統(tǒng)計檢驗如下[11-12]：

對于具有n個樣本量的時間序列X，構造一秩序列：

(2)

其中：

式中：xi、xj分別表示氣象數(shù)據(jù)集中的相應i、j序列所對應的數(shù)值。

可見秩序列Sk是第i時刻數(shù)值大于第j時刻數(shù)值數(shù)量的累計數(shù)。在時間序列隨機獨立的假定下，定義統(tǒng)計量：

(3)

式中：UF1=0；E(Sk)、Var(Sk)分別是累計數(shù)Sk的均值、方差，在x1，x2，…，xn相互獨立，且有相同連續(xù)分布時，由式(4)計算可得：

(4)

UFk是按時間序列x的順序x1，x2，…，xn計算出的統(tǒng)計量序列，給定顯著性水平α，若|UFk|>Uα，則表明序列存在明顯的趨勢變化。

按時間序列x逆序xn，xn-1，…，x1，再重復上述過程，同時使UBk=-UFK；k=n,n-1,…,1；UB1=0。

若UFk大于0，表明序列呈上升趨勢；若UBk或UFk超過臨界值(顯著水平為0.05)時，表明上升或下降趨勢顯著；將超過臨界線的范圍確定為出現(xiàn)突變的時間區(qū)域。若UFk、UBk2條曲線出現(xiàn)交點，并且在臨界線之間，則交點對應的時刻便是突變開始的時間。

1.3.3 R/S分析法 R/S分析方法由英國科學家Hurst提出，后來經(jīng)Mandelbrot與Wallis對其進一步完善，最終發(fā)展成為研究時間序列的分形理論。其主要原理如下：

給定一時間序列{ζ(t)}，t=1，2，…，n對于任意正整數(shù)τ=1，定義均值序列：

(5)

累積離差：

(6)

極差：

(7)

標準差：

(8)

若存在如下關系：R/S∝τH，則說明所分析的時間序列存在Hurst現(xiàn)象，H值稱為Hurst指數(shù)，其值可在雙對數(shù)坐標系(lnτ，ln(R/S))中用最小二乘法擬合得到。R/S分析法主要依據(jù)Hurst指數(shù)的大小，判斷時間序列是完全隨機的抑或是存在趨勢性成分。若H=0.5，表明該時間序列氣候要素獨立分布，變化完全獨立；若0

表1 Hurst指數(shù)強度分類

2 海南島各氣候要素變化特征

2.1 海南島氣溫變化特征

由圖2可以看出，海南島各氣候區(qū)的平均最高氣溫、平均最低氣溫、平均氣溫總體呈現(xiàn)上升趨勢(除西南區(qū)的三亞站)，且平均氣溫、平均最低氣溫通過α=0.05水平的顯著性檢驗，上升趨勢明顯。從年平均氣溫的總體增幅來看，平均最低氣溫增幅最大，平均最高氣溫增幅最小。在5個氣溫指標中，除三亞的所有氣候區(qū)極差均表現(xiàn)為極端最低氣溫極差>極端最高氣溫極差>平均最低氣溫極差，表明整個海南島極端最低氣溫的波動最大。而西南區(qū)極端最低氣溫、極端最高氣溫、平均最低氣溫、平均最高氣溫、平均氣溫大于其他氣候區(qū)，反映了海南島西南區(qū)氣溫波動幅度較其他氣候區(qū)大。

2.1.1 極端最高氣溫 除海南島西南區(qū)外，年極端最高氣溫在東北區(qū)、西北區(qū)、中部山區(qū)、東南區(qū)均呈現(xiàn)上升趨勢，但只有東南區(qū)的2個站點通過α=0.01、0.001水平的假設檢驗。西南區(qū)的東方、三亞極端最高氣溫在近60年來呈下降趨勢，其氣候傾向率分別為-0.128、-0.188 ℃/10年，但2個站點均未通過α=0.05水平的假設檢驗。

2.1.2 極端最低氣溫 由表2可以看出，海南島各氣候區(qū)年極端最低氣溫均上升，在各氣候區(qū)中上升幅度最大的為西南區(qū)的東方站點，為0.671 ℃/10年，上升幅度最小的瓊海的增長速率也達到了0.373 ℃/10年。其中?？?、儋州、瓊中、瓊海、陵水、東方6個站點均通過α=0.001水平的假設檢驗，三亞站點通過了α=0.05水平的假設檢驗。2007—2008年的拉尼娜事件使得整個海南島出現(xiàn)持續(xù)偏冷的天氣，其中三亞最為明顯。這與何溪澄等報道的在拉尼娜影響下我國南方溫度偏冷的結(jié)果一致[14]。西南氣候區(qū)的極差最大，東方、三亞的極差分別為9.00、10.40 ℃；極差最小的是瓊海，為 6.30 ℃(表2)。

2.1.3 平均最低氣溫 由表2還可看出，在海南島年平均最低氣溫方面，所有站點均呈現(xiàn)出增長態(tài)勢，具體增長率：?？?.308 ℃/10年、儋州0.313 ℃/10年、瓊中0.350 ℃/10年、瓊海0.308 ℃/10年、陵水0.217 ℃/10年、東方0.351 ℃/10年、三亞0.036 ℃/10年，上升幅度最大的是東方、瓊中站點，除三亞站點外其他6個站點的增長趨勢均通過了α=0.001水平的顯著性檢驗。通過計算各氣候區(qū)的極差發(fā)現(xiàn)，西南區(qū)的極差較其他地區(qū)更大，這也表明西南區(qū)的年平均最低氣溫波動更大。

2.1.4 年平均氣溫 海口、儋州、瓊中、瓊海、陵水、東方6個站點的年平均氣溫整體均呈增長趨勢，且全部通過了α=0.001 水平的假設檢驗，增長率分別為0.223、0.235、0.254、0.246、0.202、0.232 ℃/10年，三亞年平均氣溫呈下降趨勢，降幅為0.145 ℃/10年，但這種下降趨勢并未通過顯著性檢驗。三亞年平均氣溫極差達到4.81 ℃，比極差最小的陵水高2.81 ℃(表2)。

2.1.5 平均最高氣溫 ?？?、儋州、瓊中、瓊海、陵水、東方、三亞年平均最低氣溫的氣候傾向率分別為0.117、0.138、0.134、0.201、0.194、0.063、-0.262 ℃/10年，可以看出，除了三亞站點，其他所有氣候區(qū)的年平均最高氣溫均呈現(xiàn)增長趨勢，其中?？?、儋州、瓊中、瓊海、陵水5個站點的增長趨勢通過了α=0.05、0.01或0.001水平的顯著性檢驗。西南區(qū)三亞的極差最大，達到5.62 ℃；最小的是陵水，其極差為2.18 ℃(表2)。

2.2 海南島平均相對濕度變化特征

由圖2還可以看出，在1959—2015年間，海南島?？?、儋州、瓊中、瓊海、陵水、東方的平均相對濕度均呈現(xiàn)出不同程度的下降趨勢，均通過α=0.001水平的顯著性檢驗，其傾向率分別為-0.882%、-0.551%、-0.547%、-0.905%、-0.648%、-0.708%/10年，瓊海地區(qū)相對濕度的下降幅度最大，表明瓊海地區(qū)這些年來變得越來越干。主要是因為隨著城市化進程的加速，下墊面變?yōu)樗苄孕〉乃嗟兀髿饨邓杆佥斎胂滤?，而白天的建筑材料接收輻射能量后快速升溫，導致城市溫度升高。同時，由于城市綠地面積的減少，使植物的蒸騰向大氣輸送的水汽量減少，從而使得相對濕度減少[15]。三亞站點的相對濕度波動幅度較大，總體來看，三亞站點的平均相對濕度在57年來呈增加趨勢，但并未通過α=0.05水平的顯著性檢驗，1959—2008年的曲線較為平緩，但變化幅度不大，自2008年來急劇增加，增長率為21.64%。

2.3 海南島年日照時數(shù)變化特征

總體來看，海南島在1959—2015年間，中部山區(qū)瓊中站點的年日照時數(shù)以61.840 h/10年的速率在增加，且通過α=0.001水平的假設檢驗，表明海南島中部山區(qū)日照時數(shù)增加趨勢明顯。其他氣候區(qū)6個站點均呈現(xiàn)不同趨勢的減少，傾向率絕對值表現(xiàn)為海口(97.990 h/10年)>陵水(86.640 h/10年)>三亞(79.500 h/10年)>瓊海(67.070 h/10年)>儋州(28.340 h/10年)>東方(20.490 h/10年)，6個站點中僅有西南區(qū)的東方站點未通過α=0.01或0.001水平的假設檢驗，這也表明海南島東北區(qū)、西北區(qū)、東南區(qū)日照時數(shù)下降趨勢明顯。主要因為隨著城市化的發(fā)展，大氣污染增加，城市熱島效應導致上升氣流增強，攜帶到大氣中的污染顆粒增多，從而造成城市上空的云量增多，進而引起日照時數(shù)的減少[15-16]。西南區(qū)三亞、中部山區(qū)最大日照時數(shù)與最小日照時數(shù)之間的差值較大，分別為 102.62、88.97 h，西南區(qū)東方的差值最小，為50.54 h，最大值與最小值的比值達到2.03。

2.4 海南島降水量的變化特征

由圖2可以看出，1959—2015年間海南島各氣候區(qū)的降水量除瓊中外，其他地區(qū)均呈現(xiàn)不同程度的增長，具體表現(xiàn)為三亞(84.820 mm/10年)>?？?57.480 mm/10年)>儋州(52.010 mm/10年)>陵水(46.880 mm/10年)>瓊海(38.590 mm/10年)>東方(24.190 mm/10年)。然而各氣候區(qū)年降水量波動較大，年降水量除三亞站點外，均未通過α=0.01水平的假設檢驗。5個氣候區(qū)中中部山區(qū)最大降水量與最小降水量的極差最大，達到196.74 mm。由結(jié)果可以看出，總體降水量的波動性較大，降水量的無規(guī)律性可能是由降水空間分布的隨機性、區(qū)域性和復雜性造成的[6，17]。

3 海南島各氣候要素突變特征

由表3可知，海南島各氣候區(qū)的氣溫突變有一定的差異，極端最低氣溫只有東北區(qū)和東南區(qū)的陵水有突變點，突變年份分別為1987、1994年，東北區(qū)的突變年份早于東南區(qū)；以東北區(qū)例，突變前的極端最低氣溫為6.4 ℃，突變后為9.9 ℃，突變后比突變前增加54.69%，突變前后的極端氣溫差異明顯。極端最高氣溫在5個氣候區(qū)均有突變現(xiàn)象發(fā)生，突變時間最早的是西北區(qū)(1979年)，中部山區(qū)發(fā)生突變的時間最晚(2005年)。平均最低氣溫有東北區(qū)、西北區(qū)、東南區(qū)陵水、西南區(qū)三亞發(fā)生突變，時間分別為1987、1995、1994、1982年。年平均最高氣溫均有突變現(xiàn)象發(fā)生，西南區(qū)的三亞突變發(fā)生時間最早(1977年)，中部山區(qū)的突變發(fā)生時間最晚(2002年)。年平均氣溫有東北區(qū)、東南區(qū)瓊海、西南區(qū)東方、西南區(qū)三亞發(fā)生突變，對應的時間分別為1987、1993、1990、1980年。東北區(qū)的突變時間與陳小麗等研究的海南島全島近42年(1961—2002年)平均氣溫的突變時間吻合[2]。西南區(qū)三亞的突變點比我國年平均氣溫突變點(1978年)推遲了2年[18]，其他地區(qū)的年平均氣溫突變點均在尹云鶴等提出的我國氣候發(fā)生突變的時間(1989年)[19]附近徘徊。這說明在20世紀80年代該區(qū)氣溫出現(xiàn)了突發(fā)性的上升，經(jīng)歷了一個相對偏暖的氣候態(tài)，而這與全國及全球性的氣候變暖有直接聯(lián)系[20-21]。

海南島所有氣候區(qū)的平均相對濕度的UF曲線均超過 0.05 顯著水平的信度線。西北區(qū)、中部山區(qū)、東南區(qū)瓊海、東南區(qū)陵水、西南區(qū)東方、西南區(qū)三亞突變對應的時間分別為1986、1997、2000、1992、1996、1961年，相對濕度突變前后變化幅度分別為5%、1%、-1%、1%、-2%、2%。東北區(qū)相對濕度不存在突變現(xiàn)象。西北區(qū)、中部山區(qū)、東南區(qū)瓊海、東南區(qū)陵水、西南區(qū)東方、西南區(qū)三亞日照時數(shù)發(fā)生突變，突變年份分別為1966、1993、1993、1993、1973、2008年。以上地區(qū)日照時數(shù)突變前分別為181.3、163.7、174.4、173.7、220.3、205.1 h，突變后分別為152.3、154.2、162.6、190.8、215.9、200.5 h。東北區(qū)日照時數(shù)未發(fā)生明顯突變。降水量只有西南區(qū)東方、三亞有突變發(fā)生，突變年份為1995年，表現(xiàn)為突變后降水量增多，降水量在突變前分別為103.7、124.6 mm，在突變后分別為113.6、131.8 mm(表3)。

表2 海南島8項氣候要素的傾向率、決定系數(shù)和極差

注：“*”“**”“***”分別表示通過顯著性水平為0.05、0.01、0.001的檢驗。下表同。

表3 海南島8項氣候要素的突變年份

注：“*”“**”分別表示通過0.05、0.01顯著性水平的檢驗。

4 海南島各氣候要素R/S分析及Hurst指數(shù)試驗

4.1 海南島各氣候要素R/S分析

由表4可以得出，海南島各氣候要素的R/S整體分析結(jié)果：(1)極端最低氣溫的Hurst指數(shù)為-4級的有瓊海，-3級的有海口、儋州、瓊中、陵水、東方和三亞，反向持續(xù)性均較強。從海南島平均相對濕度長時間序列的變化趨勢來看，均呈現(xiàn)下降趨勢。極端最高氣溫的Hurst指數(shù)為1級的有東方，表現(xiàn)出很弱的同向持續(xù)性；-3級的有儋州和三亞，-4級的有海口、瓊海和陵水，-5級的有瓊中，均表現(xiàn)出較強的反向持續(xù)性。年平均氣溫的Hurst指數(shù)為-1級的有三亞，表現(xiàn)出很弱的反向持續(xù)性；-2級的有瓊中，表現(xiàn)出較弱的反向持續(xù)性；-3級的有海口、儋州、瓊海、陵水和東方，均表現(xiàn)出較強的反向持續(xù)性，表明這些站點的增溫趨勢會有較大程度的減弱。年平均最低氣溫的Hurst指數(shù)為-2級的有?？凇傊?、陵水、三亞，表現(xiàn)出較弱的反向持續(xù)性；-3級的有儋州、瓊海，-4級的有東方，表現(xiàn)出較強的反向持續(xù)性。年平均最高氣溫的Hurst指數(shù)為-2級的有陵水、東方和三亞，表現(xiàn)出較弱的反向持續(xù)性；-3級的有?？?、儋州、瓊中、瓊海，表現(xiàn)出較強的反向持續(xù)性。

(2)相對濕度的Hurst指數(shù)為-2級的有瓊中、陵水和三亞，反向持續(xù)性較弱；-3級的有?？凇①僦?、瓊海和東方，反向持續(xù)性較強。由日照時數(shù)R/S分析能夠看出，所有站點的Hurst值均小于0.5，海南島日照時數(shù)的長期相關特征表現(xiàn)為反向持續(xù)性，西南區(qū)三亞站點的Hurst指數(shù)達到-2級，其余氣候區(qū)均為-3級，表現(xiàn)出較強的反向持續(xù)性。海南島所有站點降水量的Hurst值均小于0.5，長期相關特征表現(xiàn)為很強的反向持續(xù)性，其中東方站點的Hurst指數(shù)為-4級，Hurst指數(shù)為0.22，其余站點的Hurst指數(shù)為-3級。

表4 海南島8項氣候要素的Hurst指數(shù)

4.2 海南島各氣候要素Hurst指數(shù)試驗

對海南島氣候要素的R/S進行分析，從整體變化來看，有著幾乎完全一致的反向持續(xù)性，但是關于這種變化趨勢能持續(xù)的時間、持續(xù)性強度會繼續(xù)增強還是減弱等方面的問題，都需要進一步研究。筆者采用馮新靈設計的Hurst指數(shù)試驗分析方案[22]，分析部分臺站(選取?？?、瓊中、三亞，分別代表海南島東北部、中部、西南部)的平均氣溫、平均相對濕度、日照時數(shù)、降水量等時間序列的未來變化趨勢，以期從總體上進行推斷。為了能夠更清晰地反映其變化趨勢，本研究將時間序列長度定為20年。具體的試驗方案如下：從時間序列數(shù)據(jù)開始年起，計算第1個20年序列(1959—1978年)的Hurst指數(shù)，接下來計算第2個20年序列(1960—1979年)的Hurst指數(shù)，再緊接著計算第3個20年序列(1961—1980年)的Hurst指數(shù)，以此類推，連續(xù)計算形成1個Hurst指數(shù)序列，最終形成以20年為序列長度的Hurst指數(shù)序列，部分臺站20年Hurst指數(shù)試驗結(jié)果見圖3。通過對比分析20年Hurst指數(shù)序列與之前研究所得傾向率的線性變化趨勢，預測未來的變化趨勢。

4.2.1 年平均氣溫Hurst指數(shù)試驗結(jié)果分析 海南島東北區(qū)年平均氣溫20年試驗趨勢呈波狀起伏。在1959—1978年的20年試驗中，Hurst指數(shù)開始出現(xiàn)小于0.5的0.347 2，隨后連續(xù)12個20年試驗的Hurst指數(shù)都小于0.5。1979年的年平均氣溫開始降低，而這種下降趨勢一致持續(xù)了12年。在1991—2010年的20年試驗中，Hurst指數(shù)開始大于0.5，連續(xù)6個20年試驗的Hurst指數(shù)都大于0.5，但Hurst指數(shù)一直在0.5附近徘徊，由此可以推測，未來10年海南島東北區(qū)年平均氣溫仍將有微弱的增加趨勢。

海南島中部山區(qū)年平均氣溫20年指數(shù)試驗變化幅度不大，年平均氣溫Hurst指數(shù)總體小于0.5。在1974—1993年的20年指數(shù)試驗中，Hurst指數(shù)小于0.5，而這種現(xiàn)象持續(xù)了21年，盡管1995—2014年的20年指數(shù)試驗出現(xiàn)大于0.5的0.506 4，但并沒有保持穩(wěn)定，在接下來的1個20年Hurst指數(shù)試驗中又躍回原態(tài)，表明中部山區(qū)未來10年依然會按照氣候傾向率0.254 ℃/10年有微弱下降的趨勢。

西南區(qū)三亞年平均氣溫20年指數(shù)試驗曲線呈現(xiàn)波狀起伏，1992—2011年Hurst指數(shù)中出現(xiàn)大于0.5的0.547 0，此后連續(xù)6個20年試驗序列中的Hurst指數(shù)都大于0.5，且Hurst指數(shù)值在不斷增加，于是出現(xiàn)了自2012年起三亞年平均氣溫下降的趨勢，由此可以推測，未來10年西南區(qū)三亞年平均氣溫依氣候傾向率0.145 ℃會有繼續(xù)下降的趨勢。

4.2.2 日照時數(shù)Hurst指數(shù)試驗結(jié)果分析 海南島東北區(qū)日照時數(shù)的20年Hurst指數(shù)試驗曲線呈波狀增加趨勢。在1959—1978年的20年試驗中，Hurst指數(shù)小于0.5，此后連續(xù)18個20年試驗的Hurst指數(shù)都小于0.5。對比分析東北區(qū)日照時數(shù)變化趨勢可知，1979年的日照時數(shù)開始增加，這種上升趨勢十分明顯。在1988開始的20年指數(shù)試驗中，Hurst指數(shù)開始出現(xiàn)大于0.5的0.565 4，此后連續(xù)9個20年試驗中，Hurst指數(shù)均大于0.5，由此可以推測，未來東北區(qū)日照時數(shù)仍將依氣候傾向率有繼續(xù)下降的趨勢。

在海南島中部山區(qū)日照時數(shù)的20年Hurst指數(shù)中，1962—1981年的Hurst指數(shù)開始出現(xiàn)小于0.5的0.418 3，此后連續(xù)22個20年試驗的Hurst指數(shù)都小于0.5，1984—2003年的Hurst指數(shù)開始出現(xiàn)大于0.5的0.510 8，日照時數(shù)開始出現(xiàn)由反向持續(xù)性向持續(xù)性轉(zhuǎn)變。2002年中部山區(qū)日照時數(shù)降低到1 647.8 h，遠遠低于多年平均值。自1988年以后的6個20年試驗中，Hurst指數(shù)一直大于0.5，但卻在不斷減少，到1994—2013年的20年試驗中，出現(xiàn)Hurst指數(shù)小于0.5的情況，由此可以推測，在未來10年內(nèi)，中部山區(qū)日照時數(shù)將會有微弱下降的趨勢。

海南島西南區(qū)三亞日照時數(shù)的20年試驗曲線在波動中有十分微弱的增加趨勢，自1962年以后的18個20年試驗中，Hurst指數(shù)一直小于0.5，1991年日照時數(shù)增加到 2 774.5 h。在1981—2000年的20年指數(shù)試驗中，Hurst指數(shù)開始出現(xiàn)大于0.5的0.525 6，之后連續(xù)16個20年試驗的Hurst指數(shù)都大于0.5，而且一路攀升，推測西南區(qū)三亞日照時數(shù)將依氣候傾向率平均每10年減少79.5 h。

4.2.3 年平均相對濕度Hurst指數(shù)試驗結(jié)果分析 在1988—2007年的20年試驗中，東北區(qū)年平均相對濕度的Hurst指數(shù)開始出現(xiàn)大于0.5的0.521 0，此后連續(xù)9個20年試驗的Hurst指數(shù)都大于0.5，相應的，2008年東北區(qū)的平均相對濕度開始增加，且上升趨勢十分明顯，由此可推測，未來10年東北區(qū)的年平均相對濕度平均每10年下降0.882%。

在1965—1984年中部山區(qū)年平均相對濕度的20年試驗中，Hurst指數(shù)開始出現(xiàn)了大于0.5的0.502 5，此后連續(xù)12個20年試驗的Hurst指數(shù)都大于0.5，平均相對濕度也就出現(xiàn)由反向持續(xù)性到持續(xù)性的轉(zhuǎn)折，在1985年達到最大值 87.25%，之后持續(xù)性下降，這種變化一直持續(xù)了12年。在1988年之后的9個20年試驗中，Hurst指數(shù)一直大于0.5，而且一路攀升，由此可以推測，在未來10年，中部山區(qū)的年平均相對濕度仍將以氣候傾向率平均每10年下降0.547%的速率發(fā)展。

西南區(qū)的三亞臺站年平均相對濕度20年試驗曲線為波狀，“U”形波谷特別明顯，在1966—1985年的20年試驗中，Hurst指數(shù)開始出現(xiàn)大于0.5的0.509 8，在此后連續(xù)8個20年試驗序列中的Hurst指數(shù)都大于0.5。在1976—1995年的20年試驗中，Hurst指數(shù)開始出現(xiàn)持續(xù)14個20年試驗的Hurst指數(shù)小于0.5，1996年的年平均相對濕度開始下降，且下降趨勢越來越明顯，這種下降趨勢一直持續(xù)了14年。在自1990年以后的7個20年試驗中，Hurst指數(shù)又開始轉(zhuǎn)為大于0.5，相應的平均相對濕度變化曲線也開始出現(xiàn)增加的相關性特征，由此可以推測，在未來10年，西南區(qū)的三亞臺站年平均相對濕度仍將按照氣候傾向率以平均每10年增加 0.623%的速率發(fā)展。

4.2.4 降水量Hurst指數(shù)試驗結(jié)果分析 海南島東北區(qū)降水量總體表現(xiàn)增加趨勢，其氣候傾向率為57.48 mm/10年。在1959—1978年的20年試驗中，Hurst指數(shù)小于0.5，在此后32個20年的試驗中，Hurst指數(shù)都小于0.5，直到1990—2009年的20年試驗中，Hurst指數(shù)開始出現(xiàn)大于0.5的0.532 3，降水量的趨勢開始由反向持續(xù)性向持續(xù)性轉(zhuǎn)變，在此后連續(xù)6個20年試驗的Hurst指數(shù)都大于0.5，由此可以推測，東北區(qū)未來降水量仍將有微弱增加趨勢。

海南島中部山區(qū)降水量總體有微弱的下降趨勢，其氣候傾向率為15.30 mm/10年。通過20年試驗可以看出，在1977—1996、1987—2006、1991—2010、1992—2011年這4個20年試驗中出現(xiàn)了Hurst指數(shù)大于0.5的情況，對比分析中部山區(qū)降水量變化曲線可以看出，降水量減少，曲線下滑，但這種變化并沒有保持穩(wěn)定，隨即躍回原態(tài)，總體上看，20年試驗的Hurst指數(shù)小于0.5，表明中部山區(qū)降水量在未來會有一定程度的增加。

西南區(qū)三亞降水量變化曲線呈增加趨勢，氣候傾向率為84.82 mm/10年。20年Hurst指數(shù)試驗的結(jié)果與中部山區(qū)十分類似，在1990—2009年的20年試驗中，Hurst指數(shù)開始出現(xiàn)大于0.5的0.561 6，但持續(xù)時間不長，僅連續(xù)3個20年試驗后，Hurst指數(shù)又開始小于0.5，表明西南區(qū)三亞降水量始終保持著反向持續(xù)性的特征，降水量減少的態(tài)勢是明顯的。

5 結(jié)論

通過對海南島1959—2015年平均氣溫、平均最低氣溫、平均最高氣溫、極端最低氣溫、極端最高氣溫、平均相對濕度、日照時數(shù)和年降水量8項氣候要素的變化特征進行分析，并結(jié)合氣候區(qū)差異，揭示了近57年來海南島氣候變化特點及未來變化趨勢。主要結(jié)論如下：

(1)在全球變暖大趨勢下，海南島各氣候區(qū)的極端最高氣溫、極端最低氣溫、平均最高氣溫、平均最低氣溫、平均氣溫總體均呈上升趨勢(三亞除外)；其中極端最高氣溫、極端最低氣溫、平均最高氣溫、平均最低氣溫、平均氣溫最高增長率分別位于東南區(qū)的陵水、西南區(qū)的東方、東南區(qū)瓊海、西南區(qū)的東方、中部山區(qū)。海南島日照時數(shù)(瓊中除外)、平均相對濕度(三亞除外)總體呈下降趨勢，日照時數(shù)和平均相對濕度下降率較高的位于東北區(qū)；降水量總體呈上升趨勢(瓊中除外)，降水量的最高增長率位于西南區(qū)的三亞。

(2)在近60年的時間內(nèi)，海南島5個氣溫指標的突變特征差異較大。海南島各氣候區(qū)極端最低氣溫只有東北區(qū)、東南區(qū)陵水、西南區(qū)三亞發(fā)生突變，其中東北區(qū)突變發(fā)生時間較早；大部分氣候區(qū)的極端最高氣溫均有突變現(xiàn)象發(fā)生，突變發(fā)生時間最早的是西北區(qū)，發(fā)生時間最晚的是西南區(qū)；年平均氣溫只有東北區(qū)、東南區(qū)、西南區(qū)發(fā)生突變；大部分氣候區(qū)的年平均最低氣溫均發(fā)生突變，但突變年份相差并不大，基本在10年內(nèi)，其中西南區(qū)發(fā)生突變較早；所有氣候區(qū)的年平均最高氣溫均發(fā)生突變，東南區(qū)發(fā)生突變的時間較早，中部山區(qū)發(fā)生突變的時間較晚。

(3)海南島各氣候區(qū)的平均相對濕度除東北區(qū)外，下降趨勢均發(fā)生突變現(xiàn)象，西南區(qū)突變現(xiàn)象發(fā)生得較早，東南區(qū)較晚發(fā)生突變。海南島日照時數(shù)有西北區(qū)、中部山區(qū)、東南區(qū)、西南區(qū)發(fā)生突變，其中西北區(qū)突變現(xiàn)象發(fā)生較早，西南區(qū)發(fā)生突變時間較晚。海南島降水量只有西南區(qū)有突變現(xiàn)象，突變時間均為1995年。

(4)對海南島各氣候要素的R/S分析表明，海南島各氣候區(qū)的8個氣候要素的Hurst指數(shù)均小于0.5(除西南區(qū)東方的極端最高氣溫)，表現(xiàn)出不同程度的反向持續(xù)性，各氣候區(qū)極端最高氣溫、極端最低氣溫、平均最高氣溫、平均最低氣溫、平均氣溫總體上升趨勢將會有不同程度的減弱(三亞除外)，海南島日照時數(shù)(瓊中除外)、平均相對濕度(三亞除外)將會有適當?shù)幕厣Ｍㄟ^對海南島部分臺站的平均氣溫、年平均相對濕度、日照時數(shù)、降水量的20年Hurst指數(shù)試驗可以看出，海南島未來10年日照時數(shù)將按照氣候傾向率持續(xù)下降，這一結(jié)果與陳小敏等的研究結(jié)果[16]一致；而平均氣溫、年平均相對濕度、降水量的未來變化存在不同的區(qū)域性差異。