黃耀賢 邱建秀
(中山大學(xué)地理科學(xué)與規(guī)劃學(xué)院 廣東廣州 510275)
太陽(yáng)輻射量時(shí)空變化是地表能量收支波動(dòng)的重要指標(biāo)[1],影響區(qū)域氣候資源的利用及區(qū)域自然災(zāi)害。日照時(shí)數(shù)是表征太陽(yáng)輻射強(qiáng)弱的重要?dú)夂蛑笜?biāo)[2]。有研究表明,全球日照時(shí)數(shù)呈現(xiàn)逐年降低的趨勢(shì),并稱這種現(xiàn)象為“全球變暗”[1,3-6]。然而,區(qū)域的日照時(shí)數(shù)變化受到多種因素的影響,例如云量、氣溶膠、水汽壓、降水等。分析區(qū)域內(nèi)的日照時(shí)數(shù)變化規(guī)律,對(duì)農(nóng)作物生長(zhǎng)、經(jīng)濟(jì)林培育、光伏發(fā)電等具有重要參考價(jià)值。
在中國(guó),日照時(shí)數(shù)的時(shí)空變化及驅(qū)動(dòng)因素存在顯著的區(qū)域差異。任國(guó)玉等[7]以全國(guó)1956—2002 年日照時(shí)數(shù)為研究對(duì)象,發(fā)現(xiàn)全國(guó)年平均日照時(shí)數(shù)總體呈顯著下降的趨勢(shì),華東因受人為排放的氣溶膠影響,為日照時(shí)數(shù)減少最顯著的地區(qū)。何彬方等[8]研究認(rèn)為,安徽省1955—2005 年日照時(shí)數(shù)大致呈北高南低分布,日照時(shí)數(shù)的顯著下降與能見(jiàn)度的下降、大氣水汽壓的增加、年降水量和年降雨日數(shù)的增加關(guān)系密切。杜軍等[9]采用西藏25 站1971—2005 年日照時(shí)數(shù)資料分析,認(rèn)為西藏年日照時(shí)數(shù)呈極顯著減少的趨勢(shì),其成因主要是大氣水汽壓增大。也有少數(shù)地區(qū)日照時(shí)數(shù)呈上升趨勢(shì)。劉艷艷等[10]研究認(rèn)為:黑河中上游年平均和季節(jié)平均的日照時(shí)數(shù)都以上升趨勢(shì)為主,其中酒泉、金塔、鼎新一線以北為主要的上升趨勢(shì)區(qū);在變化周期方面,黑河流域日照時(shí)數(shù)變化以12 年周期的振蕩最強(qiáng)。而華維等[11]利用小波分析研究青藏高原近35年日數(shù)時(shí)數(shù)指出,高原東南區(qū)、高原中部區(qū)、青海北部區(qū)和藏西南區(qū)存在顯著的16 年周期,而其他地區(qū)顯著周期及其年代變化差異較大。以往的研究表明,在過(guò)去的幾十年里,我國(guó)的日照時(shí)數(shù)普遍減少,但各地區(qū)的降幅有所差異。
針對(duì)日照時(shí)數(shù)時(shí)空變化的驅(qū)動(dòng)因子(如氣溶膠、云量),在不同地區(qū)也有相應(yīng)的研究。例如,劉玉英等[12]研究表明,1961—2012 年吉林省年日照時(shí)數(shù)呈現(xiàn)出顯著減少的年際變化特征,利用相關(guān)系數(shù)法得出低云量的增加和大氣污染物排放量增加所引起的氣溶膠光學(xué)厚度(Aerosol Optical Depth,AOD)增加是導(dǎo)致吉林省日照時(shí)數(shù)減少的主要原因;薛慧[13]基于MODIS影像得出包頭市AOD與日照時(shí)數(shù)三次多項(xiàng)式擬合模型,二者呈顯著的負(fù)相關(guān);彭艷[14]指出,氣溶膠可以直接反射、散射或吸收太陽(yáng)輻射,1970—2000 年關(guān)中地區(qū)太陽(yáng)輻射和日照時(shí)數(shù)的波動(dòng)下降趨勢(shì),與氣溶膠增加所造成的直接輻射效應(yīng)有關(guān);郭軍等[15]指出,天津地區(qū)能見(jiàn)度下降可能主要是對(duì)流層大氣氣溶膠含量上升的結(jié)果,這是造成天津地區(qū)日照時(shí)數(shù)減少的主要原因;俞兆文等[16]利用南京氣象站近50年的日照時(shí)數(shù)逐日觀測(cè)資料反演了南京地區(qū)相應(yīng)年份的AOD,結(jié)果表明該地區(qū)日照時(shí)數(shù)與氣溶膠存在較強(qiáng)的相關(guān)性。
廣東省是華南地區(qū)日照時(shí)數(shù)和太陽(yáng)能資源較多的地區(qū)之一,探討其在近30 年來(lái)快速的城市化進(jìn)程中,人類活動(dòng)等因素造成的氣溶膠排放增加對(duì)該地區(qū)日照時(shí)長(zhǎng)和光能資源的影響,對(duì)該區(qū)的新能源高效利用、生態(tài)環(huán)境政策制定等具有參考價(jià)值。鑒于此,本文對(duì)廣東省近60 年來(lái)的日照時(shí)數(shù)時(shí)空變化規(guī)律進(jìn)行分析,并探討氣溶膠光學(xué)厚度對(duì)其變化的影響,以期為廣東省氣候變化研究提供參考。
本文研究區(qū)域?yàn)閺V東省,將其劃分為粵北、粵西、粵東和珠三角4個(gè)子區(qū)域。日照時(shí)數(shù)數(shù)據(jù)來(lái)源于中國(guó)氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/),選取位于廣東省內(nèi)的25 個(gè)國(guó)家級(jí)氣象觀測(cè)站1960—2019 年逐月氣象資料,數(shù)據(jù)經(jīng)過(guò)質(zhì)量控制,剔除了缺失數(shù)據(jù)和異常值4個(gè)子區(qū)內(nèi)的氣象觀測(cè)站點(diǎn)分布如圖1 所示。氣溶膠厚度數(shù)據(jù)(MCD19A2 v006) 由 MODIS 藍(lán)光波段 (B3) 0.47 μm 的反射率反演得到 (https://lpdaac.usgs.gov/products/mcd19a2v006/),獲取自Google Earth Engine平臺(tái)。
圖1 廣東省內(nèi)國(guó)家級(jí)氣象站點(diǎn)空間分布圖
1.2.1 時(shí)序變化分析
研究采用Mann-Kendall (MK)方法進(jìn)行突變檢驗(yàn),分析日照時(shí)數(shù)的時(shí)空變化關(guān)鍵節(jié)點(diǎn)[17]。MK是一種非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法,在時(shí)間序列隨機(jī)獨(dú)立的假設(shè)下,構(gòu)建統(tǒng)計(jì)量:
其中,UF1=0,Var(Sk)、是累積量Sk的方差和均值;UFi為標(biāo)準(zhǔn)正態(tài)分布,是按時(shí)間序列x1,x2,…,xn計(jì)算出的統(tǒng)計(jì)量序列,給定顯著性水平α,若UFi>α,則表明序列有明顯的趨勢(shì)變化。計(jì)算順序時(shí)間序列的秩序列Sk,并按方程計(jì)算UFk,計(jì)算逆序時(shí)間序列的秩序列Sk,也按方程計(jì)算出UBk,分析繪出的UFk、UBk曲線,當(dāng)UFk或UBk的值大于0,則表明序列呈上升趨勢(shì),小于0 則表明呈下降趨勢(shì)。當(dāng)UFk或UBk的值超過(guò)臨界線時(shí),表明上升或下降趨勢(shì)顯著。如果UFk、UBk兩條曲線出現(xiàn)交點(diǎn),且交點(diǎn)在臨界線之間,那么交點(diǎn)對(duì)應(yīng)的值便是突變開(kāi)始的時(shí)間。
1.2.2 變化周期分析
利用小波分析開(kāi)展日照時(shí)數(shù)的變化周期分析,小波分析是一種源于Fourier變換的分析時(shí)間序列能量的局部變化的數(shù)學(xué)方法[18]。對(duì)于給定的小波函數(shù)ψ(t),時(shí)間序列f(t)的連續(xù)小波變換為:
式中:a為尺度因子,反映小波的周期長(zhǎng)度;b為時(shí)間因子,反映時(shí)間上的平移;Wf(a,b)稱為小波變換系數(shù)。由于氣象要素時(shí)間序列往往是離散的,即f(t) =f(kΔt)(k=1,2,…,n); Δt為取樣時(shí)間間隔,則f(t) 的離散形式為:
通過(guò)增加或減小伸縮尺度a來(lái)改變信號(hào)窗口大小進(jìn)而得到信號(hào)的低頻或高頻信息,分析信號(hào)的細(xì)節(jié),識(shí)別氣象多時(shí)間尺度演變特征和突變特征。
采用反距離權(quán)重法對(duì)日照時(shí)數(shù)進(jìn)行空間插值,得到全省日照時(shí)數(shù)的空間分布[19-20]。
由年平均日照時(shí)數(shù)M-K突變檢驗(yàn)分析結(jié)果(圖2)可知,珠三角、粵東、粵北、粵西地區(qū)及全省趨勢(shì)相似,存在明顯的突變現(xiàn)象,突變年份分別為1980、1970、1970、1972和1972年。
圖2 1960—2019年廣東省年平均日照時(shí)數(shù)M-K突變檢驗(yàn)
對(duì)廣東省年日照時(shí)數(shù)進(jìn)行趨勢(shì)分析和5年滑動(dòng)平均分析,結(jié)果如圖3 所示,近60 年研究期平均日照時(shí)數(shù)呈極顯著下降趨勢(shì)(p<0.01),下降速率為-46.8 h/10a,平均日照時(shí)數(shù)為1 828.5 h,最大值(2 343.0 h)出現(xiàn)在1963 年,最小值 (1 563.0 h)出現(xiàn)在2016 年。據(jù)5 年滑動(dòng)平均曲線,平均日照時(shí)數(shù)在20世紀(jì)60年代初至80年代初顯著下降,90 年代末短暫回升,之后呈波動(dòng)下降態(tài)勢(shì)。80年代初至21世紀(jì)10年代處于日照時(shí)數(shù)低值期。
表1是各子區(qū)域日照時(shí)數(shù)統(tǒng)計(jì)特征值。1960—2019 年,珠三角地區(qū)日照時(shí)數(shù)以-57.8 h/10a 的向率呈極顯著(p<0.01)的遞減趨勢(shì),且于1980 年發(fā)生了突變,突變后年均日照時(shí)數(shù)較突變前減少了204.9 h。此外,粵東、粵北和粵西3 個(gè)子區(qū)的日照時(shí)數(shù)以極顯著(p<0.01)的遞減趨勢(shì),傾向率分別為-34.4、-48.4 和-40.9 h/10a,其中粵北地區(qū)的年均日照時(shí)數(shù)在突變前后的降幅最大(13.9%)??傮w而言,廣東省在日照時(shí)長(zhǎng)方面“變暗”趨勢(shì)十分顯著。
表1 1960—2019年廣東省日照時(shí)數(shù)突變前、后多年平均值及其變化量
據(jù)廣東省春季日照時(shí)數(shù)變化曲線(圖3-a),1960—2019 年,廣東省春季日照時(shí)數(shù)以-10.5 h/10a 的傾向率呈顯著(p<0.05)的遞減趨勢(shì),近60年呈現(xiàn)上升-下降-上升-下降的波動(dòng)趨勢(shì),60年代初起可觀察到一段持續(xù)20 年左右的下降趨勢(shì),對(duì)研究期內(nèi)廣東省春季“變暗”貢獻(xiàn)明顯。夏季日照時(shí)數(shù)以-13.5 h/10a 的傾向率呈極顯著(p<0.01)的遞減趨勢(shì),平均值為581.7 h(圖3-b),總體呈現(xiàn)波動(dòng)下降后略回升的變化趨勢(shì),波谷位于1996 年前后,沒(méi)有明顯的波峰。秋季日照時(shí)數(shù)以-10.1 h/10 a 的傾向率呈顯著(p<0.05)的遞減趨勢(shì)(圖3-c),在70 年代中期至90 年代初波動(dòng)較大,波峰波谷交替出現(xiàn),隨后沒(méi)有出現(xiàn)上升趨勢(shì)。冬季日照時(shí)數(shù)傾向率為-10.1 h/10a,但沒(méi)有通過(guò)0.05 顯著性水平檢驗(yàn),冬季日照時(shí)數(shù)研究期內(nèi)較平穩(wěn)下降,沒(méi)有明顯的波峰波谷。由此可見(jiàn)夏季日照時(shí)數(shù)下降最快也最為顯著,其次是春季、秋季和冬季。
圖3 廣東省日照時(shí)數(shù)年際變化曲線
表2所示是廣東省年代際平均日照時(shí)數(shù)變化特征。全省日照時(shí)數(shù)由20世紀(jì)60年代至90年代逐代遞減,21 世紀(jì)00 年代有所回升,10 年代降至研究時(shí)段內(nèi)最低水平。20 世紀(jì)60 年代變差系數(shù)為最大值,表明該年代日照時(shí)數(shù)充足,同時(shí)年際波動(dòng)也較明顯。春季日照時(shí)數(shù)的最大值為60 年代的394.8 h,到80 年代遞減至274.0 h,90 年代和21世紀(jì)10 年代又出現(xiàn)小幅增加,10 年代再度減少至低于平均值;變差系數(shù)僅70 年代超過(guò)60 年來(lái)平均值,最小值在80 年代,可知80 年代春季平均日照時(shí)數(shù)較少且波動(dòng)平穩(wěn)。夏季日照時(shí)數(shù)年代際變化總體呈波動(dòng)下降,60 年代平均值最大,90 年代最小但變差系數(shù)最大,可見(jiàn)夏季日照時(shí)數(shù)在90 年代降至較低水平,年際波動(dòng)明顯。秋季年代際日照時(shí)數(shù)最大值也出現(xiàn)在60 年代,最小值為21 世紀(jì)10年代,相差74.3 h,70 年代至90 年代穩(wěn)定少變。冬季年代際日照時(shí)數(shù)最大值仍出現(xiàn)在60 年代,同時(shí)變差系數(shù)也最大,表明此階段日照平均時(shí)數(shù)充足但是年際波動(dòng)也較大。
表2 廣東省平均日照時(shí)數(shù)年代際統(tǒng)計(jì)表
圖4 展示的是年平均日照時(shí)數(shù)的空間分布特征。可見(jiàn)廣東省日照時(shí)數(shù)空間分布大致由南部沿海向北部?jī)?nèi)陸遞減,其中以粵西的徐聞-電白、粵東的汕頭-汕尾-惠來(lái)兩個(gè)區(qū)域?yàn)槿照諘r(shí)數(shù)高值區(qū),年平均值均超過(guò)2 000 h;而日照時(shí)數(shù)低值區(qū)包括粵北連縣-佛崗-韶關(guān)-連平珠三角的廣寧(因其靠近粵北南嶺山脈而與粵北低值區(qū)呈現(xiàn)相似性)和廣州站(年平均值均<1 700 h)。廣東省全省日照時(shí)數(shù)在研究時(shí)段初期(即20 世紀(jì)60 年代)整體處于高值時(shí)期,至90 年代逐年代遞減,21 世紀(jì)初有所回升,與已有研究[21]結(jié)論相符。
圖4 廣東省日照時(shí)數(shù)空間分布
圖5反映出廣東省日照時(shí)數(shù)變化傾向率空間分布,可見(jiàn)全省大部分地區(qū)呈顯著下降趨勢(shì),僅粵東的惠來(lái)站點(diǎn)地區(qū)傾向率為正值,但未通過(guò)0.05顯著性檢驗(yàn)。在下降速率方面,位于珠三角核心片區(qū)的廣州-增城-惠陽(yáng)-深圳為下降速率最快的地區(qū)(下降速率超過(guò)-7 h/a),并通過(guò)0.01 的極顯著檢驗(yàn);粵西的羅定-陽(yáng)江也為日照時(shí)數(shù)降低較明顯的區(qū)域,其下降速率超過(guò)-6 h/a,同樣通過(guò)0.01的極顯著性檢驗(yàn)。
圖5 廣東省1960—2019年日照時(shí)數(shù)變化趨勢(shì)空間分布
廣東省季節(jié)性日照時(shí)數(shù)空間分布(圖6)也展現(xiàn)出由南部沿海向北部?jī)?nèi)陸遞減的規(guī)律,夏、秋兩季為日照時(shí)數(shù)高值季。位于粵西雷州半島的徐聞?wù)军c(diǎn)為高值站點(diǎn),與湛江、電白站點(diǎn)組成雷州半島日照時(shí)數(shù)高值片區(qū);另一高值片區(qū)為汕頭-惠來(lái)-汕尾站點(diǎn)的粵東潮汕片區(qū)。而粵北地區(qū)位于南嶺山脈的連縣-韶關(guān)-南雄站點(diǎn)為日照時(shí)數(shù)低值片區(qū)。在春季,低值片區(qū)擴(kuò)大至廣州、增城站點(diǎn)所在的珠三角腹地。已有的文獻(xiàn)[22-24]表明,觀測(cè)到的這一空間分布規(guī)律可能與城市發(fā)展生產(chǎn)與人類生產(chǎn)排放的氣溶膠有關(guān),氣溶膠增強(qiáng)大氣對(duì)太陽(yáng)光的反射和吸收作用,使到達(dá)地面的太陽(yáng)輻射減少,從而造成日照時(shí)數(shù)減少[25-26]。
圖6 廣東省1960—2019年季節(jié)性日照時(shí)數(shù)空間分布
對(duì)1960—2019 年廣東省年平均日照時(shí)數(shù)進(jìn)行Morlet 小波變換(圖7),得到了小波系數(shù)實(shí)部等值線(圖7-a),以分析不同時(shí)間尺度日照時(shí)數(shù)的周期變化[27]。研究時(shí)段內(nèi)廣東省年平均日照時(shí)數(shù)存在的 25~32 年、11~12 年以及 5~6 年這 3 類尺度的周期變化規(guī)律。由于本研究獲取的時(shí)間序列長(zhǎng)度有限,故25~32 年周期的可信度需要在未來(lái)獲取更長(zhǎng)時(shí)序的日照時(shí)數(shù)觀測(cè)后,進(jìn)一步探討。
由小波系數(shù)模方等值線(圖7-b)可見(jiàn),11~12 年和5~6 年的時(shí)間尺度周期變化不明顯,周期性較弱。小波方差圖(圖7-c)展示3 個(gè)較為明顯的峰值,分別對(duì)應(yīng)28 年、11 年和6 年的時(shí)間尺度。圖7-d、7-e和7-f分別展示廣東省日照時(shí)數(shù)3個(gè)周期的小波系數(shù)。其中,11 年時(shí)間尺度上,平均變化周期為7.1 年左右,經(jīng)歷了8.5 個(gè)高值-低值的交替變化,而6 年尺度的變化周期相位較為雜亂,交替規(guī)律不明顯。
圖7 廣東省年平均日照時(shí)數(shù)Morlet小波變換
我國(guó)氣象臺(tái)站日照時(shí)數(shù)觀測(cè)主要是使用康培爾·斯托克日照儀,這種儀器依靠太陽(yáng)光照到筒內(nèi)涂有感光劑的感光紙上的感光跡線來(lái)計(jì)時(shí),大氣氣溶膠削弱到達(dá)地表的太陽(yáng)直接輻射,從而使得記錄到的日照時(shí)數(shù)減少[15]。為定量分析大氣中氣溶膠含量對(duì)廣東省日照時(shí)數(shù)變化的影響,分析2000—2019 年廣東省氣溶膠厚度(AOD)的空間分布規(guī)律(圖8)。其中粵東和粵西地區(qū)AOD多年平均值較小,而粵北與珠三角地區(qū)AOD多年平均值較大,約為0.58。珠三角地區(qū)由于地勢(shì)平坦、人口密集、城鎮(zhèn)化發(fā)展進(jìn)程較快,人為源氣溶膠顆粒排放量較多,造成該區(qū)AOD均值較高;而粵北南嶺山區(qū)的大氣揮發(fā)性有機(jī)物(Volatile organic compounds,VOCs)的轉(zhuǎn)化程度較高[28],VOCs是生成二次有機(jī)氣溶膠的重要前體物[29-30],是該區(qū)AOD均值較高的可能原因。
圖8 廣東省2000—2019年平均AOD空間分布
對(duì)比日照時(shí)數(shù)變化趨勢(shì)空間分布(圖5)可知,日照時(shí)數(shù)劇烈減少地區(qū)與AOD高值區(qū)具有較好的重合度,如珠三角站點(diǎn)廣寧、廣州、佛崗,粵北的南雄、連縣、東源等站點(diǎn)均落在AOD高值區(qū)內(nèi),同時(shí)AOD低值區(qū)也與日照時(shí)數(shù)增高的站點(diǎn)在空間分布上較為吻合,初步證實(shí)AOD對(duì)日照時(shí)數(shù)的影響。
進(jìn)一步定量分析AOD與日照時(shí)數(shù)在年尺度上的相關(guān)關(guān)系。珠三角、粵東、粵北和粵西地區(qū)的樣本數(shù)分別為160、120、120 和100,在剔除異常值后,擬合模型如圖9 所示。珠三角、粵東、粵北、粵西地區(qū)和全省日照時(shí)數(shù)與AOD均呈負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)分別為-0.23、-0.22、-0.21、-0.42 和-0.39,其中珠三角、粵西地區(qū)和全省通過(guò)0.01極顯著性檢驗(yàn),粵東和粵北地區(qū)通過(guò)0.05 顯著性檢驗(yàn)。
圖9 廣東省2000—2019年AOD與日照時(shí)數(shù)關(guān)系
利用1960—2019 年廣東省25 個(gè)國(guó)家級(jí)氣象觀測(cè)站的逐日日照時(shí)數(shù)資料,基于Mann-Kendall 突變檢驗(yàn)法、趨勢(shì)分析法和小波分析等方法分析了廣東省近60 年來(lái)日照時(shí)數(shù)時(shí)空變化特征,并利用氣溶膠厚度(AOD)數(shù)據(jù),分析日照時(shí)數(shù)與其相關(guān)性,得到如下結(jié)論:
(1)區(qū)域內(nèi)的年際日照時(shí)數(shù)都有顯著縮短的趨勢(shì),1972 年是廣東省日照時(shí)數(shù)縮短的轉(zhuǎn)折點(diǎn),珠三角、粵東、粵北和粵西地區(qū)轉(zhuǎn)折點(diǎn)分別為1980、1970、1970和1972年;春季、夏季、秋季、冬季日照時(shí)數(shù)下降速率分別為-10.5、-13.5、-10.1、-10.1 h/10a,夏季日照在過(guò)去60 年減少最為顯著;年代際變化方面,春季和冬季波動(dòng)較大,夏季和秋季相對(duì)穩(wěn)定。
(2)小波周期分析表明,日照時(shí)數(shù)以11年和6年為周期出現(xiàn)交替變化,11年的周期具有7.5個(gè)高值-低值的交替變化,平均變化周期為7.1年左右;28 年左右的周期有待未來(lái)獲取更長(zhǎng)時(shí)序進(jìn)一步討論。
(3)從空間分布看,日照時(shí)數(shù)呈現(xiàn)由南向北的遞減趨勢(shì),區(qū)內(nèi)最高值在粵西徐聞,最低值在粵北連縣;日照時(shí)數(shù)變化速率方面,有珠三角廣州-增城-惠陽(yáng)-深圳地區(qū)和粵西羅定-陽(yáng)江地區(qū)兩個(gè)減速最快區(qū);夏、秋兩季為日照時(shí)數(shù)高值季,春、冬兩季為低值季,春季廣州站點(diǎn)日照時(shí)數(shù)明顯偏少。
(4)廣東省AOD平均高值區(qū)為珠三角地區(qū)和粵北地區(qū),分別為0.58 和0.57;AOD 高值區(qū)在空間分布上與日照時(shí)數(shù)較低站點(diǎn)如廣寧、廣州、佛崗、連縣和韶關(guān)較吻合。省內(nèi)日照時(shí)數(shù)與AOD總體呈顯著負(fù)相關(guān)關(guān)系,相關(guān)系數(shù)為-0.39,且通過(guò)0.05顯著性檢驗(yàn)。