1971—2015年無棣縣日照時(shí)數(shù)變化及影響因素分析

徐青文　路淑文　劉樹棣

摘要：根據(jù)無棣國(guó)家一般氣象站1971-2015年近45年的日照數(shù)據(jù)資料，利用線性擬合、相關(guān)系數(shù)、Mann-Kendall非參數(shù)檢驗(yàn)等方法分析了日照時(shí)數(shù)的年、季、月變化規(guī)律，并探討了日照時(shí)數(shù)變化與霧日數(shù)、相對(duì)濕度、云量的關(guān)系。結(jié)果表明，①無棣縣年日照時(shí)數(shù)呈下降趨勢(shì)，傾向率為-12.528 h/年。自2001年開始，年日照時(shí)數(shù)低于近30年平均值。②日照時(shí)數(shù)的季節(jié)變化規(guī)律與年變化規(guī)律一致。其中，春季日照時(shí)數(shù)最多，冬季日照時(shí)數(shù)最少，秋季日照時(shí)數(shù)的減少最顯著。③日照時(shí)數(shù)的各月數(shù)值呈減少趨勢(shì)，9月日照時(shí)數(shù)減少最明顯。④年日照時(shí)數(shù)的突變時(shí)間為2000年。⑤總云量、霧日數(shù)與年日照時(shí)數(shù)、季日照時(shí)數(shù)的變化相關(guān)性較強(qiáng)。而相對(duì)濕度、低云量的相關(guān)性較弱。⑥無棣縣日照時(shí)數(shù)的減少與城市化進(jìn)程的推進(jìn)、高樓建筑增多、熱島效應(yīng)增強(qiáng)、汽車、化工廠污染物排放量增加、大氣氣溶膠增加有密切關(guān)系。

關(guān)鍵詞：日照時(shí)數(shù)；霧日數(shù)；變化特征

中圖分類號(hào)：P468.0+27 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：0439-8114（2017）20-3845-04

DOI：10.14088/j.cnki.issn0439-8114.2017.20.014

Abstract： According to the sunshine data of the 45 years from 1971 to 2015 provided by National Meteorological Station of Wudi county， linear fitting， correlation coefficient， and Mann-Kendall non-parametric test were adopted to explore the change law of sunshine duration monthly， quarterly and annually， and the relation between the changes of sunshine duration and fog days， relative humidity and cloud cover was discussed. The results showed that， ①The annual sunshine duration of Wudi county showed a downward trend， and the inclination rate was -12.528 h/a. Since 2001， the annual sunshine duration had been lower than the average value of the recent 30 years. ②The quarterly change law of sunshine duration was consistent with the annual change law. The sunshine duration in spring was the longest throughout the year while it was the shortest in winter. The sunshine duration in autumn reduced significantly. ③The monthly sunshine duration showed a decreasing trend， and it decreased most significantly in September. ④The jump time of annual sunshine duration was in 2000. ⑤Total cloud cover and fog days had a strong dependency with annual sunshine duration. While relative humidity and low cloud cover had a weak dependency with annual sunshine duration. ⑥The reduction of sunshine duration in Wudi county was mainly caused by the advancement of urbanization， the increase of high-rise buildings， the enhancement of heat island effect， the aggravation of pollutants emission of vehicles and chemical plants， and the increase of atmospheric aerosol.

Key words： sunshine duration； fog days； change characteristics

日照時(shí)數(shù)也稱日照時(shí)間，是指太陽每天在垂直于其光線的平面上的輻射強(qiáng)度超過或等于120 W/m2的時(shí)間長(zhǎng)度。日照是重要的氣象因子，也是可供人類開發(fā)利用的可再生能源。日照的變化對(duì)生態(tài)環(huán)境、社會(huì)經(jīng)濟(jì)發(fā)展、人類生產(chǎn)活動(dòng)和動(dòng)植物生長(zhǎng)發(fā)育有重要影響。

國(guó)內(nèi)外對(duì)日照時(shí)數(shù)變化特征及其影響因素的分析已取得了一定的成效。研究表明，近50年來中國(guó)日照時(shí)數(shù)逐年減少，各地季節(jié)性減幅差異性較大，冬、夏兩季遞減幅明顯大于春、秋兩季，日照時(shí)數(shù)減少主要與大氣清潔度有關(guān)[1-3]；曾昭美等[4]研究發(fā)現(xiàn)云是影響中國(guó)東南部溫度變化的可能因子，云與日照時(shí)數(shù)變化的關(guān)系較為復(fù)雜，日照與氣溫、日較差有極好的線性關(guān)系；石慧蘭等[5]對(duì)魯西北52年日照時(shí)數(shù)變化趨勢(shì)進(jìn)行分析發(fā)現(xiàn)，日照時(shí)數(shù)總體呈減少趨勢(shì)，主要是由低云量增加和人類活動(dòng)特別是工業(yè)污染物排放量增加導(dǎo)致大氣氣溶膠增多所致；覃崢嶸[6]對(duì)南寧日照時(shí)數(shù)減少的原因進(jìn)行分析，發(fā)現(xiàn)夏季日照時(shí)數(shù)的減少與降水量、廢氣排放總量的關(guān)系。endprint

本研究對(duì)山東省無棣縣近45年的日照時(shí)數(shù)變化規(guī)律及其與云量、相對(duì)濕度、霧日數(shù)的關(guān)系進(jìn)行分析。

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

1.1 數(shù)據(jù)來源

選取無棣國(guó)家一般氣象站1971-2015年日照時(shí)數(shù)、云量、霧日數(shù)、相對(duì)濕度等資料作為研究對(duì)象，其中，由于臺(tái)站觀測(cè)項(xiàng)目的調(diào)整，云量觀測(cè)于2014年停止，故云量統(tǒng)計(jì)選取1971-2013年的數(shù)值。對(duì)日照時(shí)數(shù)的年、季、月數(shù)值進(jìn)行統(tǒng)計(jì)分析，并結(jié)合云量、霧日數(shù)、相對(duì)濕度的變化規(guī)律，探討影響日照時(shí)數(shù)變化的有關(guān)因素。考慮到清晨多霧、邊界層多逆溫及人類活動(dòng)的影響，本研究的云量數(shù)據(jù)選取每日14：00的總云量。

1.2 數(shù)據(jù)處理方法

1.2.1 線性擬合 利用線性回歸擬合對(duì)日照時(shí)數(shù)、云量、霧日數(shù)、相對(duì)濕度等氣象要素進(jìn)行趨勢(shì)分析[7-9]。利用最小二乘法對(duì)資料數(shù)據(jù)進(jìn)行線性擬合，得出線性回歸方程，從而得出描述數(shù)據(jù)變化的一元線性方程。

計(jì)算過程中，以時(shí)間年份為自變量，以所研究的氣象要素為因變量，建立線性回歸方程y（t）=ax+b，其中，a為氣象要素的變化趨勢(shì)，即傾向率，若a>0，則表示氣象要素呈上升趨勢(shì)；若a<0，則表示氣象要素呈下降趨勢(shì)。|a|的數(shù)值代表了上升或下降趨勢(shì)變化的速率。

1.2.2 相關(guān)系數(shù) 相關(guān)系數(shù)由卡爾·皮爾森（Karl Pearson）在19世紀(jì)80年代提出，又稱線性相關(guān)系數(shù)、皮氏積距相關(guān)系數(shù)，用于衡量隨機(jī)變量之間的線性相關(guān)程度，通常用r來表示[10-12]。當(dāng)|r|≥0.8時(shí)，可認(rèn)為變量之間高度相關(guān)；當(dāng)0.5≤|r|<0.8時(shí)，可認(rèn)為變量之間中度相關(guān)；當(dāng)0.3≤|r|<0.5時(shí)，可認(rèn)為變量之間低度相關(guān)；當(dāng)|r|<0.3時(shí)，說明兩個(gè)變量之間的相關(guān)程度極弱。

1.2.3 Mann-Kendall非參數(shù)檢驗(yàn) 使用Mann-Kendall非參數(shù)檢驗(yàn)來檢測(cè)氣象要素的突變性。突變指的是相對(duì)穩(wěn)定態(tài)下的不連續(xù)跳躍，氣象學(xué)范疇內(nèi)，在天氣不斷變化的過程中存在的不連續(xù)性即為突變。造成突變的原因有兩種，一種是天氣系統(tǒng)內(nèi)部出現(xiàn)的沒有外界影響的狀態(tài)下的突變，這種突變通常會(huì)有周期性的變化；另一種是由于外力作用影響下的突變，如人為遷站而導(dǎo)致的氣象數(shù)據(jù)突變。Mann-Kendall非參數(shù)檢驗(yàn)方法可以用于檢測(cè)氣象因子的變化趨勢(shì)，也可以對(duì)氣象因子的突變點(diǎn)進(jìn)行分析檢驗(yàn)[13-17]。

通過分析UFk、UBk序列可以得出氣象要素的變化趨勢(shì)，確定突變的時(shí)間點(diǎn)與突變區(qū)域。若UFk數(shù)值大于0，則表明氣象要素為上升趨勢(shì)，反之則為下降趨勢(shì)，當(dāng)二者都超過臨界直線時(shí)，表明上升或下降的趨勢(shì)較為明顯。若二者的曲線存在交點(diǎn)，且交點(diǎn)在臨界直線之間，那么交點(diǎn)所對(duì)應(yīng)的時(shí)間就是突變開始的時(shí)間。

2 日照時(shí)數(shù)變化特征

2.1 年變化特征

由圖1日照時(shí)數(shù)年變化曲線可知，無棣縣年日照時(shí)數(shù)呈波動(dòng)性變化，極高值2 936.6 h出現(xiàn)在1995年。年日照時(shí)數(shù)總體呈下降趨勢(shì)，傾向率為-12.528 h/年，通過了0.01顯著性水平檢驗(yàn)。從2001年開始，日照時(shí)數(shù)明顯少于無棣縣近30年平均值2 633.1 h，極低值2 176.5 h，就出現(xiàn)在這一時(shí)期。

2.2 季節(jié)變化特征

無棣縣各季節(jié)日照時(shí)數(shù)的逐年變化見圖2。從圖2中可以看出，1971-2015年45年里無棣縣日照時(shí)數(shù)的季節(jié)變化趨勢(shì)與年日照時(shí)數(shù)變化一致，呈減少趨勢(shì)。其中，秋季減少最為明顯，傾向率為-4.41 h/年；春季減少程度最弱，傾向率為-1.68 h/年。日照時(shí)數(shù)的季節(jié)分布也并不均勻，春季日照時(shí)數(shù)最多，冬季日照時(shí)數(shù)最少。其中，春季平均日照755.02 h，夏季平均日照707.38 h，秋季平均日照625.62 h，冬季平均日照535.94 h。冬季白晝時(shí)間最短，是冬季平均日照時(shí)數(shù)在四季中最少的主要原因之一。

2.3 月變化特征

從表1中可以看出，1971-2015年無棣縣45年各月平均日照時(shí)數(shù)中，5月日照時(shí)數(shù)最多，為280.27 h。進(jìn)入6月以后，日照時(shí)數(shù)減少，8月以后，日照時(shí)數(shù)增加。9月日照時(shí)數(shù)減少最為明顯，傾向率為 -2.01 h/年，符合無棣縣日照時(shí)數(shù)的季節(jié)分布規(guī)律。且每年各月的日照時(shí)數(shù)也都呈減少趨勢(shì)，與日照時(shí)數(shù)的年變化規(guī)律一致。

由此可以看出，日照時(shí)數(shù)的月變化規(guī)律與太陽高度角各月的變化規(guī)律有密切關(guān)系[18，19]。

2.4 年日照時(shí)數(shù)的突變性分析

利用非參數(shù)統(tǒng)計(jì)檢驗(yàn)方法，對(duì)無棣縣1971-2015年45年的年日照時(shí)數(shù)進(jìn)行突變檢驗(yàn)，給定顯著性水平α=0.05，即μ0.05=±1.96。由圖3中曲線變化趨勢(shì)可以看出，自1971年開始，年日照時(shí)數(shù)呈下降趨勢(shì)，且UF與UB曲線在臨界值之間有一個(gè)交點(diǎn)，對(duì)應(yīng)的時(shí)間為2000年，即年日照時(shí)數(shù)的突變發(fā)生在2000年。根據(jù)統(tǒng)計(jì)資料，以突變時(shí)間點(diǎn)2000年為界限，1971-2000年，平均日照時(shí)數(shù)為2 760.66 h，高于近30年平均日照時(shí)數(shù)2 633.1 h，2001-2015年，平均日照時(shí)數(shù)為2 350.76 h，減少了409.9 h。

3 日照時(shí)數(shù)變化的影響因子分析

3.1 云量、霧、相對(duì)濕度對(duì)年日照時(shí)數(shù)變化的影響

日照時(shí)數(shù)的變化受多種氣象因子的影響，云量、霧、相對(duì)濕度等都對(duì)日照的長(zhǎng)短有所影響。由圖4可知，1971-2010年，無棣縣的霧日數(shù)與總云量呈增長(zhǎng)趨勢(shì)。其中，霧日數(shù)的增長(zhǎng)為0.336 d/年，總云量的增長(zhǎng)為0.047 8成/年。

經(jīng)計(jì)算，總云量的變化與日照時(shí)數(shù)呈顯著負(fù)相關(guān)，相關(guān)系數(shù)為-0.86，通過了0.01顯著性水平檢驗(yàn)。由于自動(dòng)氣象站建立在縣城內(nèi)，隨著城市化進(jìn)程的推進(jìn)，縣城內(nèi)的高層建筑變多，熱島效應(yīng)導(dǎo)致對(duì)流云增加，且由于人類活動(dòng)對(duì)大氣的影響，總云量在不斷增加，使得到達(dá)地面的陽光被削弱，從而降低了日照時(shí)數(shù)。endprint

近幾年，隨著無棣縣城市發(fā)展進(jìn)程的加速，空氣污染比以前更為嚴(yán)重，霧日數(shù)也隨之不斷增加。霧日數(shù)的變化與日照時(shí)數(shù)變化的相關(guān)系數(shù)為-0.43，也呈負(fù)相關(guān)關(guān)系，通過0.01顯著性水平檢驗(yàn)。出現(xiàn)霧時(shí)，空氣中的水分含量增加，大氣透明度降低，陽光的穿透力也隨之下降，日照時(shí)數(shù)減少，霧日數(shù)的增多是日照時(shí)數(shù)減少的主要因素之一。

其中，低云量以及相對(duì)濕度的變化與日照時(shí)數(shù)變化的相關(guān)系數(shù)分別為0.26、0.24，相關(guān)性較弱，在文中暫不做分析。

3.2 云量、霧、相對(duì)濕度對(duì)季節(jié)日照時(shí)數(shù)變化的影響

對(duì)總云量的季變化值進(jìn)行線性擬合發(fā)現(xiàn)，45年中，總云量的季節(jié)變化呈增長(zhǎng)趨勢(shì)，秋季增長(zhǎng)最快，傾向率為0.167 4成/年。且日照時(shí)數(shù)的季節(jié)變化與總云量呈顯著負(fù)相關(guān)，相關(guān)系數(shù)在0.70左右，均通過了0.01顯著性水平檢驗(yàn)，其中夏秋兩季相關(guān)系數(shù)高達(dá)0.85。由此可知，總云量是影響無棣縣日照時(shí)數(shù)季節(jié)變化的主要因素。

通過對(duì)霧日數(shù)進(jìn)行線性擬合可知，霧日數(shù)呈逐年增長(zhǎng)趨勢(shì)，45年中，秋季增長(zhǎng)最為明顯，傾向率為0.167 9 d/年；春季傾向率最小，為0.047 8 d/年，且冬季霧日數(shù)最多，為368 d，春季霧日數(shù)最少，為62 d。這與日照時(shí)數(shù)的季節(jié)變化規(guī)律一致。霧日數(shù)的季節(jié)變化也是影響日照時(shí)數(shù)季節(jié)變化的氣象因素之一，二者呈負(fù)相關(guān)。其中，春、秋兩季相關(guān)系數(shù)為0.50，通過了0.01顯著性水平檢驗(yàn)，夏、冬兩季相關(guān)系數(shù)為0.40，通過了0.05顯著性水平檢驗(yàn)。

其中，低云量、相對(duì)濕度的季節(jié)變化與日照時(shí)數(shù)季節(jié)變化的相關(guān)性較弱，暫不做分析。相對(duì)濕度的大小影響著霧日數(shù)的多少，進(jìn)而會(huì)影響日照時(shí)數(shù)的變化，但相關(guān)性并不明顯。

綜上所述，春季風(fēng)大，相對(duì)濕度較低，為55.84%，為四季最低值，不利于霧的形成，且降水較少，晴天多，云層少且薄，大氣對(duì)太陽輻射的削弱程度較弱，到達(dá)地面的太陽輻射就相對(duì)較強(qiáng)，因此，春季無棣縣的日照時(shí)數(shù)最長(zhǎng)。

夏季，陰雨天氣較多，云層厚且多，總云量、低云量數(shù)值比其他三季都多，分別為20.50、8.76成/年，云對(duì)太陽輻射的削弱作用強(qiáng)，到達(dá)地面的太陽輻射就弱，因而夏季日照較春季少。但是，夏季霧日數(shù)較少，只有106 d，比秋季、冬季霧日數(shù)的一半還要少。再者，夏季白晝較秋、冬季時(shí)間長(zhǎng)，日照時(shí)數(shù)也就比較長(zhǎng)。

秋季，霧日數(shù)、總云量的增加程度是四季之最。20世紀(jì)70年代，秋季霧日數(shù)平均為4 d/年，而2006-2015年，霧日數(shù)平均為9.7 d/年，霧日數(shù)的增長(zhǎng)較為迅速，共有291 d。

冬季，受全球氣候變暖的影響，降雪日數(shù)變少。再加上大風(fēng)天氣較少，不利于大氣污染物擴(kuò)散，霧日數(shù)增多，導(dǎo)致冬季日照時(shí)數(shù)減少。

可見，秋冬季節(jié)，隨著無棣縣城市化進(jìn)程的推進(jìn)，建筑物的增多，城市熱島效應(yīng)增強(qiáng)，且汽車數(shù)量的增加，多個(gè)化工廠的建立，導(dǎo)致污染物排放較多，且污染物擴(kuò)散不利、大氣氣溶膠增加，使得秋冬季節(jié)霧、霾日數(shù)增多，從而導(dǎo)致日照時(shí)數(shù)不斷減少。

3 小結(jié)

1）年日照時(shí)數(shù)呈下降趨勢(shì)，傾向率為-12.528 h/年，極高值為2 936.6 h，出現(xiàn)在1995年。自2001年開始，年日照時(shí)數(shù)低于30年平均值，極低值 2 176.5 h出現(xiàn)在這一時(shí)期。

2）日照時(shí)數(shù)的季節(jié)變化規(guī)律與年變化規(guī)律一致。其中，春季日照時(shí)數(shù)最多，冬季日照時(shí)數(shù)最少。且秋季日照時(shí)數(shù)的減少最顯著。

3）日照時(shí)數(shù)的各月數(shù)值呈減少趨勢(shì)，9月日照時(shí)數(shù)減少最明顯，傾向率為-2.01 h/年。日照時(shí)數(shù)的月變化值與太陽高度角的變化規(guī)律有密切關(guān)系。

4）年日照時(shí)數(shù)的突變時(shí)間為2000年，2000年以前，平均日照時(shí)數(shù)為2 760.66 h，2000年以后，平均日照時(shí)數(shù)減少了409.9 h，變?yōu)? 350.76 h。

5）總云量、霧日數(shù)與年日照時(shí)數(shù)、季日照時(shí)數(shù)的變化相關(guān)性較強(qiáng)。而與相對(duì)濕度、低云量的相關(guān)性較弱。

6）春季，風(fēng)力較大，相對(duì)濕度為55.84%，為四季最低值，不利于霧形成，且降水較少，晴天較多，大氣對(duì)太陽輻射削弱程度低，因此，春季日照時(shí)數(shù)最長(zhǎng)。

7）夏季，陰雨天較多，云層厚，大氣對(duì)太陽輻射的削弱程度大，導(dǎo)致日照較少，但由于夏季霧日數(shù)較少且白晝時(shí)間長(zhǎng)，日照時(shí)數(shù)也就比較多。

8）45年中，秋、冬季節(jié)，霧日數(shù)較多，分別為291、368 d，是夏季霧日數(shù)的一倍多。秋、冬季節(jié)降雨、降雪天氣較少，大風(fēng)天氣也相對(duì)較少，不利于大氣污染物的擴(kuò)散，導(dǎo)致日照時(shí)數(shù)較少。

9）無棣縣日照時(shí)數(shù)的減少與城市化進(jìn)程的推進(jìn)、高樓建筑增多、熱島效應(yīng)增強(qiáng)、汽車、化工廠污染物排放量增加、大氣氣溶膠增加有密切關(guān)系。

參考文獻(xiàn)：

[1] 任國(guó)玉，郭 軍，徐銘志，等.近50年中國(guó)地面氣候變化基本特征[J].氣象學(xué)報(bào)，2005，63（6）：942-956.

[2] 陳少勇，張康林，邢曉賓，等.中國(guó)西北地區(qū)近47a日照時(shí)數(shù)的氣候變化特征[J].自然資源學(xué)報(bào)，2010，25（7）：1142-1152.

[3] YANG Y H，ZHAO N，HAO X H，et al. Decreasing trend of sunshine hours and related driving forces in north China[J].Theoretical and Applied Climatology，2009，97（1/2）：91-98.

[4] 曾昭美，嚴(yán)中偉，章名立.近40年我國(guó)云、日照、溫度及日較差的統(tǒng)計(jì)[J].科學(xué)通報(bào)，1993，38（5）：440-443.

[5] 石慧蘭，王新堂，邵志勇，等.魯西北52年日照變化特征及原因[J].氣象，2007，33（2）：93-97.endprint

[6] 覃崢嶸.南寧日照時(shí)數(shù)減少的氣候特征及成因分析[J].廣西氣象，1999，20（2）：38-40.

[7] 謝仁波，李小龍，呂昌濤.思南縣近55a日照變化特征[J].貴州氣象，2010，34（6）：27-28.

[8] 張秀紅，桂翰林，王春華，等.近40年漠河日照變化特征分析[J].現(xiàn)代化農(nóng)業(yè)，2009（5）：23-25.

[9] 龍昭芪，張萬東，文 成.天柱縣近52a日照時(shí)數(shù)變化特征分析[A].貴州省氣象學(xué)會(huì).貴州省氣象學(xué)會(huì)2013年學(xué)術(shù)年會(huì)論文集[C].2013.272-275.

[10] 虞海燕，劉樹華，趙 娜，等.我國(guó)近59年日照時(shí)數(shù)變化特征及其與溫度、風(fēng)速、降水的關(guān)系[J].氣候與環(huán)境研究，2011，16（3）：389-398.

[11] 安 娟，吳福杰，楊丹寧，等.1953-2010年熊岳日照時(shí)數(shù)變化特征及其影響因素[J].中國(guó)農(nóng)學(xué)通報(bào)，2013，29（32）：355-360.

[12] 蔡冬梅.近60年大連市日照時(shí)數(shù)變化規(guī)律及氣象影響因子分析[J].安徽農(nóng)業(yè)科學(xué)，2011，39（11）：6602-6604.

[13] 王 華，牛清明.阿克蘇市日照時(shí)數(shù)的突變檢測(cè)分析[J].新疆氣象，2002，25（3）：14-15.

[14] 史建橋，白淑英，周 寅，等.近50年甘肅日照時(shí)數(shù)時(shí)空變化特征及突變分析[J].南京信息工程大學(xué)學(xué)報(bào)（自然科學(xué)版），2016， 8（1）：64-70.

[15] 蔣 沖，劉曉磊，程楠楠，等.秦嶺南北日照時(shí)數(shù)時(shí)空變化及突變特征[J].干旱區(qū)地理，2013，36（3）：416-424.

[16] 張明捷，馮 杉，姚 勝.濮陽日照變化趨勢(shì)特征和突變分析[J].安徽農(nóng)學(xué)通報(bào)，2011，17（15）：213-215.

[17] 趙 霞，余衛(wèi)東，范志軍.商丘市近40年日照變化特征及突變分析[J].河南氣象，2005（4）：16-17.

[18] 王國(guó)安，米鴻濤，鄧天宏，等.太陽高度角和日出日落時(shí)刻太陽方位角一年變化范圍的計(jì)算[J].氣象與環(huán)境科學(xué)，2007，30（S1）：161-164.

[19] 張海玉，曹 彥，王永波，等.哈爾濱市有效日照時(shí)數(shù)計(jì)算與分析[J].低溫建筑技術(shù)，2015（2）：20-21.endprint