淮河上中游流域日照時數時空分布特征

甄自強+王景才+徐蛟+汪雨+郭佳香

摘要：為揭示淮河上中游流域日照時數時空變化特征及突變情況，基于1960-2014年日照時數資料，采用氣候傾向率、MannKendall趨勢與突變檢驗、滑動t檢驗、克里金插值等方法，對流域日照時數的時空變化特征進行剖析。結果表明，淮河上中游流域近55年來日照時數以20世紀80年代為界由多變寡，其中年日照時數大約每10 a減少1022 h，四季日照時數的氣候傾向率排序為夏[-488 h（10a）]、冬[-293 h（10a）]、秋[-197 h（10a）]、春[-54 h（10a）]。流域日照時數的減小趨勢除春季不顯著外，年與夏、秋、冬三季的日照時數均呈現出顯著的減少趨勢。空間上，流域北部的中間地帶（商丘站附近）與西北部（寶豐、鄭州站附近）的日照時數減小趨勢較明顯（Z<-2），減少幅度較大，而東部（盱眙站除外）減少幅度則偏小。年與夏季日照時數的極顯著突變年份為1980年，年日照時數的弱顯著突變年為1984年。秋季、冬季的突變年份分別為1999年、1983年，春季的突變年份在20世紀70年代內。研究成果可為淮河上中游流域的氣候研究、農業生產等提供科學的依據。

關鍵詞：日照時數；趨勢分析；突變檢驗；特征分析；淮河上中游流域

中圖分類號：P466文獻標識碼：A文章編號：

16721683（2017）06003306

Abstract：In order to reveal the spatialtemporal variation characteristics of sunshine duration in the upper and middle reaches of Huaihe River Basin and determine the year of mutation，we adopted such methods as the climate tendency rate，MannKendall trend and mutation test，sliding ttest，and Kriging interpolation，and analyzed the spatialtemporal variation characteristics of sunshine duration based on the data of sunshine duration from 1960 to 2014The results indicated that the sunshine duration in the past 55 years decreased to a lower level in the 1980s；the annual sunshine duration decreased by about 1022 h every 10 years；and the climate tendency rates of sunshine duration in the four seasons were as follows：summer （-488 h（10a）），winter （-293 h（10a）），autumn （-197 h（10a）），and spring（-54 h（10a））The sunshine duration in a year and in all three seasons except spring showed a significant decreasing trendSpatially，the decreasing trend of sunshine duration was more obvious and sharp in the middle of the northern basin （near Shangqiu station） and in the northwest （near Zhengzhou and Baofeng stations） （Z<-2），but the decrease was small in the east （except for Xuyi station）The highly significant mutation year of annual and summer sunshine durations was 1980，and the weakly significant mutation year of annual sunshine duration was 1984The mutation years of autumn and winter were 1999 and 1983 respectively，and the mutation year of spring was in the 1970sThe research results are expected to provide scientific basis for the climate studies and agricultural production in the upper and middle reaches of Huaihe River Basin

Key words：sunshine duration；trend analysis；test of mutation；characteristic analysis；upper and middle reaches of Huaihe River Basin

[HJ19mm]日照時數的變化情況對氣候研究具有重要意義，其數量大小能夠合理地表征太陽能源的豐富程度[1]。具體來講，日照對水文要素具有間接影響；日照時數影響植物的光合作用，其直接決定著植物的生長發育尤其是農作物的產量與質量；日照還與異常氣候現象的發生有關，如近年來出現的厄爾尼諾現象等[2]。對氣候因子時間序列的變化情況進行分析，借助統計學方法進行研究更加有效[3]。近年來，國內很多學者運用各種統計學方法對日照時數時間序列進行了趨勢及突變情況的分析。吳昊旻等[4]通過Morlet連續小波變化、功率譜檢驗及MannKendall檢驗方法，得出浙江省近39年日照時數的時間變化趨勢、周期演變規律和突變點。趙娜等[5]由線性趨勢法與累計距平曲線方法分析出城市化影響下的北京地區日照時數變化趨勢。王釗等[6]利用日照時數距平與5年滑動平均曲線得到西安地區近60年的日照時數變化趨勢，并計算出該城市內部各地區間的日照時數下降百分率的相關系數。張山清等[7]通過線性傾向率與MannKendall突變檢測方法對新疆近50年101個站點的日照時數進行了突變年的討論。王珂清[8]通過線性傾向估計法對近48年淮河流域日照時數進行了研究，得出流域日照時數呈下降趨勢且氣溶膠為其減小的主要影響因素。第15卷 總第93期·南水北調與水利科技·2017年12月endprint

甄自強等·淮河上中游流域日照時數時空分布特征

淮河上中游流域跨越了秦嶺淮河南北分界線，近年來流域的區域性極端氣候現象頻發，日照時數作為流域的關鍵氣候因子，對流域的氣候變化情況具有重要的影響[911]。此外，日照時數對光合作用至關重要，決定著植物的產出效益，直接影響當地農業經濟的發展狀況[1213]。本文利用趨勢與突變檢驗方法對日照時數的時間變化情況進行分析，再通過ArcGis100軟件對研究區內各個站點進行空間插值分析，實現對淮河上中游流域近55年來日照時數的時空變異性比較，以期了解淮河上中游流域日照時數的時間變化特征及空間差異性，從而為研究區內季節及區域性的農業規劃、極端氣候分析等相關研究提供科學的參考。

1材料及方法

11研究區概況

淮河流域位于中國東部地區，流經湖北、河南、安徽、江蘇、山東五大省份。氣候上屬于濕潤氣候與半濕潤氣候的過渡區域，日照充足，農業較為發達，作物主要包含稻作與旱作兩類。研究區選取為淮河上中游流域，其中淮河上游以洪河口為界，淮河中游以洪澤湖出口處的三河閘為界[1415]。研究區多年平均日照時數約為2 080 h，地形西高東低，高程多位于20～140 m之間[16]。本文基于淮河上中游流域內的19個氣象站點（圖1） 1960-2014年的日照時數資料作為計算的基礎數據（圖2）。

12研究方法

121氣候傾向率

利用氣候傾向率對日照時數進行變化趨勢的分[HJ18mm]析，確定趨勢變化的方向及變化幅度的大小[1719]。

[AKx^]t=a1·t+a0（1）

式中：[AKx^]t為日照時數的擬合值（h）；a1×10為日照時數氣候傾向率大小[h（10a）]；a0為回歸常數；t為時間因子。

122MannKendall趨勢檢驗

對日照時數時間變化趨勢的檢驗選用MannKendall趨勢檢驗方法，此法已經被證明為氣象因子時間序列趨勢檢驗中一種比較有效的檢驗途徑。MannKendall趨勢檢驗方法通過計算統計量Z，并得到顯著性水平α下的統計值大小，進而可確定出日照時數的時間變化趨勢。Z值為正代表增加趨勢，Z值為負代表減少趨勢。|Z|的大小表征了變化趨勢的顯著性程度：當|Z|>Z1-α2，認為在顯著性水平為α的條件下拒絕原假設，日照時數有顯著的變化趨勢，否則認為變化趨勢不顯著。經與氣候傾向率進行結合，二者互相驗證，可實現對日照時數變化趨勢情況的綜合性分析[2022]。

123突變檢驗

對氣象因子的時間序列進行突變檢驗，一般可應用累計距平曲線、滑動t檢驗、MannKendall突變檢驗或信噪比方法等。為避免檢驗過程中人為干擾性的存在，研究優先采用MannKendall突變檢驗法。MannKendall突變檢驗方法不受主觀因素的影響，得到正逆統計量兩條曲線，只需曲線交點處的統計量達到一定的顯著度水平即可認為突變結果可信，研究中選取α=005的顯著性水平作為突變點的界限進行判斷[3，5]。通過MannKendall突變檢驗的年份可以確定為突變年，而未通過信度線檢驗的交點年份僅能考慮為可能突變點，此時仍需借助其他方法進行驗證。為此本研究采用滑動t檢驗、累計距平曲線法結合MannKendall方法進行日照時數突變年份的檢驗，即利用滑動t檢驗法驗證MannKendall突變檢驗未通過信度線的交點及累計距平曲線出現拐點處的可能突變點[2324]。

2日照時數變化特征分析

21年際變化特征

近55年淮河上中游流域的日照時數在時間上總體呈下降趨勢。年日照時數最大值發生在1966年，最小值發生在2003年，峰谷差值8392 h，占近55年日照時數的年均值的比例高達404%。由表1可知，近55年研究區年日照時數的氣候傾向率為-1022 h（10a），說明每10 a日照時數減少量大約為1022 h。四季內夏季日照時數的減少幅度較大[-488 h（10a）]，其次為冬季與秋季[-293、[HJ]-197 h（10a）]，春季較小[-54 h（10a）]。說明四季中夏季對年日照時數減少的貢獻量較大，其次為冬季與秋季，而春季對年日照時數減少的貢獻量較小。年、夏季與冬季的Z值均通過了0001顯著性水平檢驗，說明三者的日照時數減小趨勢均表現為極顯著的狀態。秋季Z值通過001顯著性水平檢驗，說明此季節的日照時數減小趨勢為顯著，而春季未通過005顯著性水平檢驗，說明此季的日照時數變化趨勢并不明顯。這與王珂清[8]對近50年來淮河流域日照時數的變化趨勢分析結果基本一致。信度，極顯著。22年代際變化特征

如圖3所示，從近55年內的日照時數年代際變化來看，近55年中，20世紀60年代（年均值）達到最大值，為2 3353 h；20世紀80年代之前的年日照時數為偏多階段（距平為正），20世紀80年代之后一般表現為距平為負的偏少階段；21世紀初出現最小值，為1 8674 h。從20世紀60年代以來，時間上隨著年代的更替，日照時數年代均值表現為依次遞減的趨勢。以近55年內各年代每10 a的年均日照時數變化量作為年代日照時數的變化幅度，各個相鄰年代間的變化幅度存在一定的差異，最大變幅發生于20世紀70年代到80年代的交替期，日照時數明顯減少（減小量為1912 h）；最小變幅為20世紀80年代到90年代的交替期，日照時數幾乎未變動（減小量僅42 h）。總體來看，近55年內淮河上中游流域日照時數年代（年均值）峰谷點的差值高達4679 h，約占年代均值的222%。

23空間變化特征

根據淮河上中游流域各個站點計算得到的日照時數氣候傾向率與MannKendall趨勢檢驗的Z值，通過ArcGis軟件中的空間克里金插值功能可以實現研究區內的日照時數變化趨勢的空間分析（圖4）。endprint

如圖4所示，從MannKendall趨勢檢驗結果來看，流域內19個站點的年日照時數均存在減小趨勢，除蚌埠站（Z=-182）減小趨勢不明顯外，其余站點（Z=-7～-413）均有極顯著的減小趨勢。從氣候傾向率來看，年日照時數氣候傾向率的空間變化范圍為-1717～-282 h（10a），時間上年日照時數減小量在蚌埠、宿州兩站出現較小值[-282、-487 h（10a）]；減小量較大的為商丘站[-1717 h（10a）]，其次為寶豐、阜陽、鄭州、亳州站附近[范圍在-1450～-1255 h（10a）]。

四季中，春季日照時數僅在4個站點（商丘、鄭州、寶豐、亳州站）具有顯著的減小趨勢（Z=-32～-22），其余各站點的變化趨勢均不顯著（Z=-2～2）。其中有6個站點表現為不顯著的增加趨勢，增幅較大的為蚌埠站[-124 h（10a）；另外5站增幅偏小，大約每10 a僅增加0～35 h；其它9個站點出現不顯著的減小趨勢且減小量也偏小（大約每10 a僅增加2～25 h）。夏、秋、冬三季所有站點的日照時數在時間上均為減小趨勢。其中，夏季所有站點日照時數的減小趨勢均為極顯著（Z=-670～-331），減小量大約每10 a減小270～693 h：夏季期間，蚌埠站的減小量較小[-270 h（10a）]，其次為徐州、宿州站[-305、-346 h（10a）]；減小量較大的站為商丘站[-693 h（10a）]，其次為阜陽、鄭州、寶豐、亳州站，這與年日照時數變化趨勢的空間分布情況較為類似。研究區內冬季日照時數的氣候傾向率的變化范圍為-470～-114 h（10a），秋季日照時數的氣候傾向率的變化范圍為-356～-22 h（10a）。冬季與秋季的減小趨勢在宿州、蚌埠站附近不明顯，兩個季節的日照時數每10 a減小量在這兩站附近也較小，而兩個季節的日照時數每10 a減小量的較大值則出現在商丘、鄭州，西華、寶豐站附近。

總體來看，年與四季日照時數的氣候傾向率與Z值的空間大小分布情況較為一致。對日照時數空間變化特征的分析也再次驗證了各個站點對年日照時數的變化貢獻較大的季節為夏季，較小的為春季。研究區內西部區域日照時數一般比東部區域日照時數遞減的趨勢更為明顯。其中東部的蚌埠、宿州等站點附近對研究區日照時數的減小趨勢貢獻程度較小，甚至這些站點春季的日照時數的變化對年日照時數的減少具有削弱作用，但春季日照時數的增加趨勢并不顯著，對年日照時數減小情況的削弱效果也不明顯；貢獻程度較大的位置集中于研究區北部中間地區的商丘、亳州等站附近，其次為西部的鄭州、寶豐站附近。

3突變檢驗

MannKendall突變檢驗結果如圖5所示。對MannKendall突變檢驗結果曲線未通過信度線的交點年份，再利用滑動t檢驗進行驗證，并選擇累計距平曲線拐點進行滑動t檢驗的突變計算，計算結果如表2所示。近55年研究區的年日照時數突變檢驗中，1986年未通過滑動t檢驗，1987年僅通過了顯著性水平為005的滑動t檢驗，說明可將1987年作為研究區年日照時數的弱顯著突變年份。由年日照時數累計距平曲線查到1980年可能為突變年，而此年又通過了顯著性水平為0001的滑動t檢驗，可以確定此年為年日照時數的極顯著突變年份。春季日照時數的UFUB曲線在信度線內存在若干交點，可以認為20世紀70年代期間春季日照時數存在突變。夏季UFUB曲線交點未通過α=005的MannKendall突變檢驗，但根據曲線交點處的滑動t檢驗結果仍能夠確定夏季日照時數的突變年份為1980年，且與年日照時數突變相同，夏季在此年的突變結果也通過了0001信度檢驗。比較t檢驗值絕對值的大小可以得到1980年研究區夏季日照時數較年日照時數相比突然減小的情況更為明顯。進一步說明1980年附近可作為近55年研究區內年與夏季日照時數由偏多期明顯轉變為偏少期的分界點，這與何斌方[25]對近50年來安徽省日照時數的分析結果基本相同。秋季UFUB曲線交點位于信度線以內，其統計量的值為-0650（>-1960），確定突變年份為1999年。冬季UFUB曲線交點統計量的值為-1956（>-1960），故冬季突變年份可確定為1983年。

4結論

通過近55年淮河上中游流域日照時數的變化趨勢及突變年份檢測的計算結果，分析研究區日照時數的時空分布特征，得到的具體結論如下。

（1）1960-2014年內淮河上中游流域日照時數

在20世紀60年代（年均值）較大，并向之后各個年代依次減小。全年、夏季、冬季的日照時數均有極顯著的減小趨勢，減小量分別為每10 a減少1022 h、488 h、2934 h。四季中夏季日照時數對年日照時數減小的貢獻較大，其次為冬季、秋季，春季日照時數對年日照時數減小的貢獻較小。

（2）淮河上中游流域的東部（除盱眙站）對整個研究區日照時數減小的貢獻程度較小，而北部的中間地帶（商丘站附近）與西北部（寶豐、鄭州站附近）對整個淮河上中游流域的日照時數減小的貢獻程度較大。

（3）研究區內全年與夏季的日照時數的時空變化情況不僅具有較好的一致性，突變年份也有類似之處。由累計距平曲線結合滑動t檢驗方法，得到年與夏季日照時數的極顯著突變年份均為1980年；根據MannKendall突變檢驗，可以得到年日照時數的另一個弱顯著突變年份為1984年。秋季、冬季的突變年份為1999年、1983年，而春季的突變年份確定為20世紀70年代。

研究成果期待為淮河上中游流域的氣候研究及農業發展等提供科學的參考。

參考文獻（References）：

[1]董旭光，邱粲，劉煥彬，等山東省日照時數的氣候變化特征及其影響因素[J]中國農業氣象，2013，34（2）：138145（DONG X G，QIU C，LIU H B，et alSunshine duration variation and its influence factors in Shandong Province[J]Chinese Journal of Agrometeorology，2013，34（2）：138145（in Chinese）） DOI：103969jissn10006362201302003endprint

[2]寧和平，敖澤建，趙紅巖甘南白龍江流域光照資源的氣候變化特征[J]自然資源學報，2014，29（7）：12071215（NING H P，AO Z J，ZHAO H YClimatic change of sunshine duration in the Bailong River Valley in Gannan from 1973 to 2012[J]Journal of Natural Resources，2014，29（7）：12071215（in Chinese）） DOI：1011849zrzyxb201407011

[3]符淙斌，王強氣候突變的定義和檢測方法[J]大氣科學，1992，16（4）：482493（FU C B，WANG QThe definition and detection of the abrupt climatic change[J]Chinese Journal of Atmospheric Sciences，1992，16（4）：482493（in Chinese）） DOI：103878jissn1006989519920411

[4]吳昊旻，廖良清浙江省日照時數時空分布特征[J]中國農學通報，2011，27（32）：294298（WU H M，LIAO L QSpatial and temporal distribution of sunshine hours in Zhejiang Province[J]Chinese Agricultural Science Bulletin，2011，27（32）：294298（in Chinese））

[5]趙娜，劉樹華，杜輝，等城市化對北京地區日照時數和云量變化趨勢的影響[J]氣候與環境研究，2012，17（2）：233243（ZHAO N，LIU S H，DU H，et alThe effects of urbanization on sunshine duration and the trend of cloud cover amount variation in Beijing area[J]Climatic and Environmental Research，2012，17（2）：233243（in Chinese）） DOI：103878jissn10069585201110140

[6]王釗，彭艷，白愛娟，等近60年西安日照時數變化特征及其影響因子分析[J]高原氣象，2012，31（1）：185192（WANG Z，PENG Y，BAI A J，et alAnalysis on variation of sunshine duration in Xi′an in recent 60 years and its influencing factor[J]Plateau Meteorology，2012，31（1）：185192（in Chinese））

[7]張山清，普宗朝，李景林近50年新疆日照時數時空變化分析[J]地理學報，2013，68（11）：14811492（ZHANG S Q，PU Z C，LI JinglinThe spatialtemporal variation of sunshine duration in Xinjiang during 19612010[J]Acta Geographica Sinica，2013，68（11）：14811492（in Chinese）） DOI：1011821dlxb201311004

[8]王珂清近五十年淮河流域氣候變化與未來情景預估[D]南京：南京信息工程大學，2013：2132（WANG K QResearch on climatic change of Huaihe Basin in recent years and its projections in the future[D]Nanjing：Nanjing University of Information Science & Technology，2013：2132（in Chinese））

[9]周光裕淮河流域植被的過渡性特點及南北分界線的探討[J]植物生態學與地植物學叢刊，1965，3（1）：131137（ZHOU G YTransitional characteristics of vegetation and boundary between south and north in Huaihe drainage area[J]Acta Phytoecologica et Geobotanica Sinica，1965，3（1）：131137（in Chinese））

[10][ZK（#]趙思雄，孫建華近年來災害天氣機理和預測研究的進展[J]大氣科學，2013，37（2）：297312（ZHAO S X，SUN J HStudy on mechanism and prediction of disastrous weathers during recent years[J]Chinese Journal of Atmospheric Sciences，37（2）：297312（in Chinese）） DOI：103878jissn10069895201212317

[11]高輝，王永光2007年海溫和大氣環流異常及對我國氣候的影響[J]氣象，2008，34（4）：108112（GAO H，WANG Y GSea surface temperature and the general circulation in 2007 and their influences on the climate of China[J]Meteorological Monthly，2008，34（4）：108112（in Chinese）） DOI：107519jissn10000526200804014endprint

[12]許信旺，孫滿英，方宇媛，等安徽省氣候變化對水稻生產的影響及應對[J]農業環境科學學報，2011，30（9）：17551793（XU X W，SUN M Y，FANG Y Y，et alImpact of climatic change on rice production and response strategies in Anhui Province[J]Journal of AgroEnvironment Science，2011，30（9）：17551793（in Chinese））

[13]蘇坤慧，延軍平，白晶，等河南省境內淮河南北氣候變化的小麥適應度比較[J]地理科學進展，2012，31（1）：6371（SU K H，YAN J P，BAI J，et alComparative studies on degree of adaption of wheat under climate change between areas south and north of Huaihe River in Henan Province[J]Progress in Geography，2012，31（1）：6371（in Chinese）） DOI：1011820dlkxjz201201009

[14]汪美華，謝強，王紅亞未來氣候變化對淮河流域徑流深的影響[J]地理研究，2003，22（1）：7988（WANG M H，XIE Q，WANG H YImpact of the climate change on runoff depth of the Huaihe drainage basin[J]Geographical Research，2003，22（1）：7988（in Chinese）） DOI：103321jissn：10000585200301010

[15]楊玉璋，程至杰，李為亞，等淮河上、中游地區史前稻旱混作農業模式的形成、發展與區域差異[J]中國科學：地球科學，2016，46（8）：10371050（YANG Y Z，CHENG Z J，LI W Y，et alThe emergence development and regional differences of mixed farming of rice and millet in the upper and middle Huai River Valley，China[J]Science China Earth Sciences，2016，46（8）：10371050（in Chinese）） DOI：101007s1143001550884

[16]曹琦萍，王傳輝，周建平，等江淮流域夏半年日照時數氣候特征及趨勢分析[J]氣象，2014，40（5）：548553（CAO Q P，WANG C H，ZHOU J P，et alClimate characteristics and trend analysis of sunshine hours in summer half year over the YangtzeHuanhe Basin[J]Meteorological Monthly，2014，40（5）：548553（in Chinese）） DOI：107519jissn10000526201405004

[17]王健，呂新，王江麗，等石河子市近57年積溫變化趨勢特征研究[J]干旱區資源與環境，2013，27（6）：154159（WANG J，LYU X，WANG J L，et alThe variation tendency of accumulated temperature of Shihezi City in recent 57 years[J]Journal of Arid Land Resources and Environment，2013，27（6）：154159（in Chinese）） DOI：1013448jcnkijalre201306032

[18]孫宇輝，韋杰，黃小燕近50年重慶市氣候干濕變化時空特征[J]長江科學院院報，2016，33（5）：1117（SUN Y H，WEI J，HUANG X YSpatialtemporal characteristics of drywet condition in Chongqing in the past five decades[J]Journal of Yangtze River Scientific Research Institute，2016，33（5）：1117（in Chinese）） DOI：1011988ckyyb20150052

[19]劉志娟，楊曉光，王文峰，等氣候變化背景下我國東北三省農業氣候資源變化特征[J]應用生態學報，2009，20（9）：21992206（LIU Z J，YANG X G，WANG W F，et alCharacteristics of agricultural climate resources in three provinces of Northeast China under global climate change[J]Chinese Journal of Applied Ecology，2009，20（9）：21992206（in Chinese）） DOI：1013287jissn1001933220090362

[20]MANN H BNonparametric tests against trend[J]Econometrica：Journal of the Econometric Society，1945，13（3）：245259

[21][ZK（#]康淑媛，張勃，郭玉剛，等1959年至2008年石羊河流域日照時數時空變化特征[J]資源科學，2011，33（3）：483488（KANG S Y，ZHANG B，GUO Y G，et alChanges of sunshine duration over the Shiyang River Basin during the last 50 years[J]Resources Science，2011，33（3）：483488（in Chinese））endprint

[22]楊佳，錢會，高燕燕，等西安市多年降水特征分析及降水量預測[J]南水北調與水利科技，2016，14（3）：3035（YANG J，QIAN H，GAO Y Y，et alMultiyear precipitation characteristics analysis and precipitation forecast of Xi′an city[J]SouthtoNorth Water Transfers and Water Science & Technology，2016，14（3）：3035（in Chinese）） DOI：1013476jcnkinsbdqk201603006

[23]劉義花，汪青春，王振宇，等1971年2007年青海省日照時數的時空分布特征[J]資源科學，2011，33（5）：10101016（LIU Y H，WANG Q C，WANG Z Y，et alDistribution characteristics of sunshine duration in Qinghai Province[J]Resources Science，2011，33（5）：10101016（in Chinese））

[24]張林梅，苗運玲，李健麗，等新疆阿勒泰地區近50 a夏季極端降水事件變化特征[J]冰川凍土，2015，37（5）：11991208（ZHANG L M，MIAO Y L，LI J L，et alVariations of summer extreme precipitation events in Altay Prefecture，Xinjiang Region，from 1961 to 2010[J]Journal of Glaciology and Geocryology，2015，37（5）：11991208（in Chinese）） DOI：107522jisnn1000024020150134

[25]何彬方，馮妍，荀尚培，等安徽省50年日照時數的變化特征及影響因素[J]自然資源學報，2009，24（7）：12751285（HE B F，FENG Y，XUN S P，et alClimatic change of sunshine duration and its influencing factors over Anhui Province during the last 50 years[J]Journal of Natural Resources，2009，24（7）：12751285（in Chinese）） DOI：1011849zrzyxb200907015endprint