基于多元線性回歸模型的華南地區(qū)≥10℃積溫空間模擬研究

戴聲佩+李海亮+羅紅霞+劉海清+曹建華

摘要選用1980～2011年華南地區(qū)111個(gè)及周邊地區(qū)39個(gè)氣象站地面觀測(cè)資料，基于ArcGIS10.0平臺(tái)，利用多元線性回歸模型和IDW、Kriging、Spline等插值方法模擬了華南地區(qū)≥10℃活動(dòng)積溫的空間分布。結(jié)果表明：基于多元線性回歸模型的≥10℃活動(dòng)積溫空間模擬方法與IDW、Kriging、Spline等插值方法相比，解決了積溫空間分布模擬復(fù)雜性的問(wèn)題，并考慮了≥10℃活動(dòng)積溫受海拔、經(jīng)緯度差異等因素的影響，具有較高的模擬精度，適合于較大范圍的積溫?cái)?shù)據(jù)空間化模擬。另外，1980～2011年華南地區(qū)≥10℃活動(dòng)積溫隨海拔的升高而減?。浑S緯度的增大而減??；隨經(jīng)度的變化幅度較大，在研究區(qū)西部云南地區(qū)受地形的影響變化幅度強(qiáng)烈，在研究區(qū)東部廣西、廣東、福建等地變化幅度較小。

關(guān)鍵詞多元線性回歸模型 ；活動(dòng)積溫 ；空間模擬 ；華南地區(qū)

分類(lèi)號(hào)P423

Spatial Simulation of AAT10 (Active Accumulated Temperature≥10℃)

based on Multiple Linear Regression Model

DAI Shengpei1,2)LI Hailiang1,2)LUO Hongxia1,2)LIU Haiqing1,3)CAO Jianhua4)

(1 Institute of Scientific and Technical Information, CATAS, Danzhou, Hainan 571737, China

2 Key Laboratory of Practical Research on Tropical Crops Information Technology in Hainan,

Danzhou, Hainan 571737, China

3 Research Center of Tropical Agricultural Economics，CATAS，Danzhou， Hainan 571737，China;

4 Rubber Research Institute, CATAS, Danzhou, Hainan 571737, China)

AbstractThe spatial distribution of active accumulated temperature ≥10℃ (AAT10) was simulated by using the multiple linear regression model (MLRM) and IDW, Kriging, Spline interpolation method based on the daily meteorological observation data from 111 meteorological stations in Southern China and surrounding areas during 1980 to 2011. The result shows that compared with the conventional methods such as IDW, Kriging and Spline, it solves the complexity of the spatial distribution of AAT10, and considers the elevation, longitude and latitude differences to AAT10. It shows higher simulation accuracy,and is suitable for a wide range of spatialization of accumulated temperature. At the same time, the AAT10 in South China was reduced with increasing altitude from 1980 to 2011; and reduced with increasing latitude. There is a large variation with longitude in the western Yunnan region affected by the impact of changes in topography, and a little variation with longitude in the eastern part of Guangxi,Guangdong, Fujian province.

Keywordsmultiple linear regression model (MLRM) ; active accumulated temperature≥10℃ (AAT10) ; Spatial simulation ; Southern China

政府間氣候變化委員會(huì)(Intergovernmental Panel on Climate Change，IPCC)于2013年9月發(fā)布了第五次全球氣候變化評(píng)估報(bào)告(Fifth Assessment Report，AR5)第一工作組報(bào)告——《氣候變化2013：自然科學(xué)基礎(chǔ)》。報(bào)告指出：氣候變暖是毋庸置疑的，自1880年以來(lái)，地球平均的表面溫度上升了0.85℃，預(yù)計(jì)到2100年全球氣溫將上升2～4.8℃，并且極有可能(95%可信度)是人為活動(dòng)導(dǎo)致全球變暖，同時(shí)伴隨出現(xiàn)大氣和海洋升溫、冰雪量下降、海平面上升、溫室氣體濃度增加、極端天氣頻繁出現(xiàn)等現(xiàn)象[1]。

氣候變化給人類(lèi)的生存和發(fā)展帶來(lái)了一系列重大影響，已經(jīng)危及到農(nóng)業(yè)安全、能源安全、生態(tài)安全、水資源安全、公共衛(wèi)生安全等各方面[2]。農(nóng)業(yè)是受氣候變化影響最敏感的領(lǐng)域之一。積溫是一個(gè)地區(qū)十分重要的熱量資源，是作物生長(zhǎng)的重要因素之一，對(duì)指導(dǎo)農(nóng)業(yè)生產(chǎn)具有重要意義[3-4]，也是進(jìn)行農(nóng)業(yè)氣候區(qū)劃、合理配置農(nóng)作物以及預(yù)報(bào)物候期、病蟲(chóng)害發(fā)生期的重要依據(jù)[5]。了解復(fù)雜地形下各界限溫度下精細(xì)網(wǎng)格的積溫對(duì)農(nóng)、林、牧業(yè)的規(guī)劃、布局是十分重要的[6]。由于受地理?xiàng)l件、維護(hù)條件等因素的限制，氣象站點(diǎn)的布設(shè)很不均勻。如何根據(jù)有限的氣象站資料獲取空間化的氣象要素是近年來(lái)生態(tài)學(xué)、資源科學(xué)和環(huán)境科學(xué)的重要任務(wù)[4]。近年來(lái)，一些學(xué)者開(kāi)展了關(guān)于華南地區(qū)積溫方面的相關(guān)研究[7-11]，因受當(dāng)時(shí)資料、技術(shù)方法等方面的限制，尚存在一些不足，如缺乏對(duì)積溫的精確空間化模擬。

華南地區(qū)屬熱帶、亞熱帶季風(fēng)氣候，氣候溫暖濕潤(rùn)、雨熱同季，區(qū)內(nèi)動(dòng)植物種類(lèi)繁多，同時(shí)擁有廣闊的熱帶海洋。不僅是我國(guó)重要的糧食生產(chǎn)基地，還是我國(guó)主要熱帶作物的生產(chǎn)基地。因此，充分利用該區(qū)豐富的熱量和水分資源，發(fā)展農(nóng)業(yè)生產(chǎn)以及熱帶經(jīng)濟(jì)作物，合理利用和保護(hù)熱帶性植物和動(dòng)物資源，開(kāi)發(fā)熱帶海洋資源等，對(duì)于我國(guó)自然資源的開(kāi)發(fā)利用具有重要意義。因此，在此背景下，研究華南地區(qū)≥10℃積溫的空間模擬、揭示熱量資源分配的新格局，對(duì)現(xiàn)有農(nóng)業(yè)結(jié)構(gòu)和品種布局進(jìn)行重新規(guī)劃和調(diào)整具有重要意義。

endprint

1資料與方法

1.1研究區(qū)域概況

華南地區(qū)位于歐亞大陸南端(3°58′N(xiāo)～26°23′N(xiāo)，104°29′E～117°50′E)，包括廣東、廣西、福建、云南和海南等省(區(qū))，總陸地面積為96.99萬(wàn) km2(圖1)，屬熱帶、亞熱帶季風(fēng)氣候，氣候溫暖濕潤(rùn)，雨熱同季。年均氣溫21.40 ℃，降水量1 900.20 mm。華南地區(qū)大部分農(nóng)作物屬于我國(guó)的晚三熟和熱三熟區(qū)域，農(nóng)業(yè)生產(chǎn)水平相對(duì)較高，主要糧食作物有水稻、旱稻、小麥、番薯、木薯、玉米、高梁等，珠江三角洲、潮汕平原和廈門(mén)、漳州、泉州一帶是我國(guó)雙季稻高產(chǎn)地區(qū)[12]。經(jīng)濟(jì)作物主要有橡膠、甘蔗、麻類(lèi)、花生、芝麻、茶等；該區(qū)的熱帶林木、熱帶水果、熱帶水產(chǎn)在我國(guó)農(nóng)業(yè)生產(chǎn)中占有突出地位。

1.2數(shù)據(jù)來(lái)源

本文選用華南地區(qū)111個(gè)及周邊(四川、貴州、湖南、江西、浙江等省份)39個(gè)氣象站1980～2011年逐日平均氣溫實(shí)測(cè)數(shù)據(jù)及各氣象臺(tái)站地名、經(jīng)度、緯度和海拔高度數(shù)據(jù)，資料來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)(http：//cdc.cma.gov.cn/)和中國(guó)氣象局。部分站點(diǎn)缺測(cè)數(shù)據(jù)，通過(guò)已有數(shù)據(jù)與臨近站點(diǎn)數(shù)據(jù)作一致性檢驗(yàn)之后進(jìn)行插補(bǔ)，以保證數(shù)據(jù)的完整性和連續(xù)性。數(shù)字高程(Digital Elevation Models，DEM)數(shù)據(jù)采用美國(guó)地質(zhì)調(diào)查局(U.S. Geological Survey，USGS)發(fā)布的全球90 m空間分辨率DEM數(shù)據(jù)；研究區(qū)邊界數(shù)據(jù)采用國(guó)家基礎(chǔ)地理信息中心發(fā)布的1∶400萬(wàn)中國(guó)行政區(qū)劃數(shù)據(jù)。

1.3研究方法

1.3.1≥10℃積溫的計(jì)算

本文計(jì)算積溫的界限溫度確定為10℃，采用五日滑動(dòng)平均法確定日平均氣溫穩(wěn)定≥10℃起止日期，其起始日定義為五日滑動(dòng)平均氣溫≥10℃的日期，終止日定義為五日滑動(dòng)平均氣溫<10℃的日期，具體算法見(jiàn)參考文獻(xiàn)[13]。由于在熱帶地區(qū)全年日平均氣溫均大于10℃，所以研究區(qū)內(nèi)部分站點(diǎn)的日平均氣溫穩(wěn)定≥10℃起止日期分別為儒略歷第一日和最后一日。

對(duì)于某一年某一站點(diǎn)而言，在確定界限溫度(10℃)起止日期之后，計(jì)算該起止日期之間的日平均氣溫總和即為日平均氣溫穩(wěn)定≥10℃期間積溫[14]。

1.3.2基于多元線性回歸模型(Multiple Linear Regression Model，MLRM)的積溫空間模擬方法

在多要素的地理系統(tǒng)中，多個(gè)(多于2個(gè))要素之間也存在著相關(guān)影響、相互關(guān)聯(lián)的情況[15-17]。假設(shè)某一因變量y受k個(gè)自變量x1，x2，…，xk的影響，其n組觀測(cè)值為yα，x1α，x2α，…，xkα，α=1，2，…，n。那么，多元線性回歸模型的結(jié)構(gòu)形式為：

ya=β0+β1x1a+β2x2a+…+βkxka+εa(1)

上式中，β0，β1，…，βk為待定參數(shù)，εα為隨機(jī)變量。如果b0，b1，…，bk分別為β0，β1，…，βk的擬合值，則回歸方程為：

y=b0+b1x1+b2x2+…+bkxk(2)

上式中，b0為常數(shù)，b1，b2，…，bk為偏回歸系數(shù)，y為ya的預(yù)測(cè)值，則觀測(cè)值(實(shí)際值)為預(yù)測(cè)值 y與殘差y′的和，公式為：

y=y+y′(3)

在本研究中，考慮到積溫受經(jīng)度、緯度和海拔的影響。因此，建立積溫與海報(bào)、經(jīng)度、緯度之間的多元線性回歸模型，公式為：

y=b0+b1θ+b2φ+b3λ(4)

上式中，θ為經(jīng)度、φ為緯度、λ為海拔；將站點(diǎn)數(shù)的90%(134個(gè))≥10℃積溫及其經(jīng)度、緯度、海拔代入上式進(jìn)行建模，得出回歸方程如下：

yj=22 235.969 8-65.185 5θ-313.936 2φ-1.575 1λ(5)

經(jīng)統(tǒng)計(jì)檢驗(yàn)，模型的置信水平α為0.001，復(fù)相關(guān)系數(shù)為0.941 6，R2為0.886 0，達(dá)到建模的精度驗(yàn)證要求。

1.3.3積溫空間模擬過(guò)程

(1)參數(shù)獲?。豪醚芯繀^(qū)DEM數(shù)據(jù)在ArcGIS10.0平臺(tái)上獲取的研究區(qū)經(jīng)度和緯度90 m分辨率柵格數(shù)據(jù)，獲得模擬所需的三個(gè)參數(shù)(海拔、經(jīng)度、緯度)(圖2)。

(2)預(yù)測(cè)柵格圖：利用獲取的模擬參數(shù)和回歸模型公式(5)求得研究區(qū)1980～2011年≥10℃積溫的預(yù)測(cè)柵格圖。

(3)模擬結(jié)果：利用建模站點(diǎn)回歸模型模擬殘差進(jìn)行反距離權(quán)重(Inverse Distance Weighted，IDW)插值，獲得模擬殘差柵格圖，并根據(jù)公式(3)計(jì)算得出研究區(qū)1980～2011年≥10℃積溫空間分布圖(圖3a)。

(4)精度驗(yàn)證：利用未參與建模的16個(gè)氣象站點(diǎn)的實(shí)測(cè)數(shù)據(jù)作為模擬結(jié)果的驗(yàn)證數(shù)據(jù)，經(jīng)檢驗(yàn)?zāi)M值與實(shí)測(cè)值的線性一致性較強(qiáng)，R2為0.953 8(圖3b)，表明模擬結(jié)果符合精度驗(yàn)證要求，說(shuō)明基于多元線性回歸模型的≥10℃活動(dòng)積溫空間模擬方法可行，適合于較大范圍的積溫?cái)?shù)據(jù)空間化模擬。

(5)模型對(duì)比：利用反距離權(quán)重(IDW)、普通克里格(Kriging)、樣條函數(shù)(Spline)等插值方法對(duì)數(shù)據(jù)進(jìn)行空間插值，獲取不同插值方法模擬的積溫空間分布圖(圖4)，并對(duì)比模擬精度(表1)。

(6)結(jié)果分析：分析研究區(qū)1980～2011年≥10℃積溫的空間分布規(guī)律，并利用ArcGIS10.0三維分析(3D Analyst)工具的剖面分析功能，獲取不同海拔高度和不同經(jīng)緯線(101°E、110°E、117°E、24°N)上的剖面曲線(圖5)，分析≥10℃活動(dòng)積溫隨海拔、經(jīng)緯度的變化情況。

2結(jié)果分析

2.1積溫空間模擬結(jié)果

由圖3a和圖3b可知，基于多元線性回歸模型的≥10℃積溫空間模擬分布圖較好體現(xiàn)了海拔、經(jīng)緯度對(duì)積溫的影響，更真實(shí)地體現(xiàn)了積溫的空間分布狀況。從模型的精度驗(yàn)證也可看出，模擬值與實(shí)測(cè)值的線性一致性較強(qiáng)，R2為0.953 8，表明基于多元線性回歸模型的積溫空間模擬方面適合于較大范圍的積溫?cái)?shù)據(jù)空間化模擬，對(duì)農(nóng)業(yè)結(jié)構(gòu)調(diào)整、空間布局規(guī)劃等方面有積極的指導(dǎo)意義。

研究區(qū)1980～2011年≥10℃積溫小于 3 000℃的區(qū)域分布在云南西北部的德欽、貢山、香格里拉、麗江和東北部的昭通、會(huì)澤等高原地區(qū)；≥10℃積溫大于3 000℃小于5 000℃的區(qū)域分布在云南的沾益、瀘西、玉溪、大理、瀘水、保山和廣西北部融安、桂林以及福建的武夷山、九仙山、漳平等地區(qū)；≥10℃積溫大于5 000℃小于7 000℃的區(qū)域分布在云南中部的騰沖、楚雄、景東和東部的文山、硯山、屏邊，廣西中部的百色、都安、河池、來(lái)賓、蒙山、梧州，廣東北部的連州、廣寧、韶關(guān)、佛崗、南雄、連平、五華、梅縣，福建北部等低山丘陵地區(qū)；≥10℃積溫大于7 000℃的區(qū)域在云南主要分布在怒江、瀾滄江、元江、金沙江等河谷地區(qū)和西雙版納地區(qū)，還分布在廣西、廣東的南部地區(qū)，以及福建南部的漳州、廈門(mén)一帶；其最高值出現(xiàn)在海南島，其值為9 415℃。

2.2積溫隨海拔、經(jīng)度、緯度的變化分析

圖5a曲線為≥10℃積溫隨海拔的變化情況。由該圖可知，研究區(qū)1980～2011年≥10℃積溫隨海拔的升高而減小。圖5b、5c、5d曲線分別為不同經(jīng)線(101°E、110°E、117°E)上的≥10℃積溫變化剖面曲線，可知，≥10℃積溫隨緯度的增大而減小。圖5e曲線為24°N緯線上的≥10℃活動(dòng)積溫變化剖面曲線。由此可知，≥10℃活動(dòng)積溫隨經(jīng)度的變化幅度較大，在研究區(qū)西部云南地區(qū)受地形的影響變化幅度強(qiáng)烈，在研究區(qū)東部廣西、廣東、福建等地變化幅度較小。

endprint

2.3與其它插值方法對(duì)比分析

由圖3、4可知，對(duì)比多元線性回歸模型模擬的積溫空間分布結(jié)果與常用的空間化方法IDW、Kriging、Spline等插值方法結(jié)果，發(fā)現(xiàn)4種方法對(duì)≥10℃積溫總體分布趨勢(shì)表現(xiàn)基本一致，都是隨緯度的增大而減小，呈現(xiàn)出由南向北逐漸降低的趨勢(shì)。但是基于多元線性回歸模型的≥10℃積溫空間模擬方法的精度較優(yōu)，主要表現(xiàn)在以下幾方面。

(1)≥10℃積溫隨海拔的變化趨勢(shì)和精細(xì)度更明顯，體現(xiàn)了海拔對(duì)活動(dòng)積溫的顯著影響。

(2)對(duì)預(yù)留10%(16個(gè))氣象站點(diǎn)的模擬值與觀測(cè)值線性一致性更強(qiáng)，R2為0.953 8，分別大于IDW(0.818 6)、Kriging(0.840 7)、Spline(0.753 2)對(duì)插值方法的模擬結(jié)果(表1)，并且均方根誤差(Root-mean-square Error，RMSE)和相對(duì)均方根誤差(Relative Root-mean-square Error，RMSEr)都小于IDW、Kriging、Spline插值方法(表1)，表明其模擬精度高、誤差小。

3結(jié)論與討論

本文選用華南地區(qū)111個(gè)及周邊地區(qū)39個(gè)氣象站點(diǎn)1980～2011年逐日平均氣溫實(shí)測(cè)數(shù)據(jù)及各氣象臺(tái)站地名、經(jīng)度、緯度和海拔高度數(shù)據(jù)，利用多元線性回歸模型和IDW、Kriging、Spline等插值方法模擬了華南地區(qū)≥10℃積溫的空間分布，得出結(jié)論如下。

(1)通過(guò)分析1980～2011年華南地區(qū)≥10℃積溫的空間分布情況，及其隨海拔、經(jīng)緯度的變化情況，可知研究區(qū)≥10℃積溫最低值出現(xiàn)在云南省北部地區(qū)，最高值出現(xiàn)在海南島；并且≥10℃積溫隨海拔的升高而減小，隨緯度的增大而減小，隨經(jīng)度的變化幅度較大。在研究區(qū)西部云南地區(qū)受地形的影響變化幅度強(qiáng)烈，在研究區(qū)東部廣西、廣東、福建地區(qū)變化幅度較小。

(2)基于多元線性回歸模型的≥10℃積溫空間模擬方法較IDW、Kriging、Spline等插值方法的精度較優(yōu)。其R2為0.953 8，分別大于IDW(0.818 6)、Kriging(0.840 7)、Spline(0.753 2)插值方法的模擬結(jié)果，并且其模擬結(jié)果隨海拔、經(jīng)緯度的變化趨勢(shì)和精細(xì)度更明顯，體現(xiàn)了海拔、經(jīng)緯度對(duì)積溫的顯著影響。

(3)本文在進(jìn)行≥10℃積溫空間模擬時(shí)，雖然考慮了海拔、經(jīng)度和緯度的影響，使模擬結(jié)果的精度和精細(xì)化程度都有所提高。但也存在一些不足，比如未考慮坡度、坡向、地形遮蔽、坡地反射輻射、散射輻射各向異性及天氣過(guò)程等因素對(duì)積溫的影響，同時(shí)，還應(yīng)進(jìn)一步加強(qiáng)遙感衛(wèi)星資料在本研究中的應(yīng)用。

參考文獻(xiàn)

[1] http：//www.ipcc.ch/report/ar5/wg1/[EB/OL].

[2] 柏秦鳳，霍治國(guó)，李世奎，等. 1978年前、后中國(guó)≥10℃年積溫對(duì)比[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2008，19(8)：1 810-1 816.

[3] 廖順寶，李澤輝. 積溫?cái)?shù)據(jù)柵格化方法的實(shí)驗(yàn)[J]. 地理研究，2004，23(5)：633-640.

[4] 袁淑杰，谷曉平，向紅瓊，等. 基于GIS的貴州高原復(fù)雜地形下積溫的精細(xì)空間分布[J]. 資源科學(xué)，2010，32(12)：2427-2432.

[5] 李飛，王春，趙軍，等. 中國(guó)陸地多年平均積溫空間化研究[J]. 自然資源學(xué)報(bào)，2010，25(5)：778-784.

[6] 于貴瑞，何洪林，劉新安，等. 中國(guó)陸地生態(tài)信息空間化技術(shù)研究(Ⅰ)——?dú)庀?氣候信息的空間化技術(shù)途徑[J]. 自然資源學(xué)報(bào)，2004，19(4)：537-544.

[7] 顧駿強(qiáng)，楊軍. 中國(guó)華南地區(qū)氣候和環(huán)境變化特征及其對(duì)策[J]. 資源科學(xué)，2005，27(1)：128-135.

[8] 李勇， 楊曉光， 王文峰， 等. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅰ：華南地區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2010，21(10)：2 605-2 614.

[9] 代姝瑋，楊曉光，趙孟，等. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅱ：西南地區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2011，22(2)：442-452.

[10] 趙錦，楊曉光，劉志娟，等. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅱ：南方地區(qū)氣候要素變化特征及對(duì)種植制度界限可能影響[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2010，43(9)：1 860-1 867.

[11] 李勇，楊曉光，王文峰，等. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅴ：氣候變暖對(duì)中國(guó)熱帶作物種植北界和寒害風(fēng)險(xiǎn)的影響分析[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2010，43(12)：2 477-2 484.

[12] 駱繼賓. 我國(guó)農(nóng)業(yè)氣候資源的開(kāi)發(fā)利用與農(nóng)業(yè)區(qū)域開(kāi)發(fā)[J]. 農(nóng)業(yè)區(qū)劃，1990，11(5)：44-47.

[13] 劉勤，嚴(yán)昌榮，何文清，等. 黃河流域近40a積溫動(dòng)態(tài)變化研究[J]. 自然資源學(xué)報(bào)，2009，24(1)：147-153.

[14] 曲曼麗. 農(nóng)業(yè)氣候?qū)嵙?xí)指導(dǎo). 北京： 北京農(nóng)業(yè)大學(xué)出版社，1991.

[15] 徐建華. 現(xiàn)代地理學(xué)中的數(shù)學(xué)方法(第二版)[M]. 北京：高等教育出版社， 2002.

[16] 劉新安，于貴瑞，范遼生，等. 中國(guó)陸地生態(tài)信息空間化技術(shù)研究(Ⅲ)——溫度、降水等氣候要素[J]. 自然資源學(xué)報(bào)，2004，19(6)：818-825.

[17] 高仲亮，周汝良，李智，等. 積溫分布的空間模擬及三維可視化初探[J]. 廣東農(nóng)業(yè)科學(xué)，2012，16：182-184，237-238.

endprint

2.3與其它插值方法對(duì)比分析

由圖3、4可知，對(duì)比多元線性回歸模型模擬的積溫空間分布結(jié)果與常用的空間化方法IDW、Kriging、Spline等插值方法結(jié)果，發(fā)現(xiàn)4種方法對(duì)≥10℃積溫總體分布趨勢(shì)表現(xiàn)基本一致，都是隨緯度的增大而減小，呈現(xiàn)出由南向北逐漸降低的趨勢(shì)。但是基于多元線性回歸模型的≥10℃積溫空間模擬方法的精度較優(yōu)，主要表現(xiàn)在以下幾方面。

(1)≥10℃積溫隨海拔的變化趨勢(shì)和精細(xì)度更明顯，體現(xiàn)了海拔對(duì)活動(dòng)積溫的顯著影響。

(2)對(duì)預(yù)留10%(16個(gè))氣象站點(diǎn)的模擬值與觀測(cè)值線性一致性更強(qiáng)，R2為0.953 8，分別大于IDW(0.818 6)、Kriging(0.840 7)、Spline(0.753 2)對(duì)插值方法的模擬結(jié)果(表1)，并且均方根誤差(Root-mean-square Error，RMSE)和相對(duì)均方根誤差(Relative Root-mean-square Error，RMSEr)都小于IDW、Kriging、Spline插值方法(表1)，表明其模擬精度高、誤差小。

3結(jié)論與討論

本文選用華南地區(qū)111個(gè)及周邊地區(qū)39個(gè)氣象站點(diǎn)1980～2011年逐日平均氣溫實(shí)測(cè)數(shù)據(jù)及各氣象臺(tái)站地名、經(jīng)度、緯度和海拔高度數(shù)據(jù)，利用多元線性回歸模型和IDW、Kriging、Spline等插值方法模擬了華南地區(qū)≥10℃積溫的空間分布，得出結(jié)論如下。

(1)通過(guò)分析1980～2011年華南地區(qū)≥10℃積溫的空間分布情況，及其隨海拔、經(jīng)緯度的變化情況，可知研究區(qū)≥10℃積溫最低值出現(xiàn)在云南省北部地區(qū)，最高值出現(xiàn)在海南島；并且≥10℃積溫隨海拔的升高而減小，隨緯度的增大而減小，隨經(jīng)度的變化幅度較大。在研究區(qū)西部云南地區(qū)受地形的影響變化幅度強(qiáng)烈，在研究區(qū)東部廣西、廣東、福建地區(qū)變化幅度較小。

(2)基于多元線性回歸模型的≥10℃積溫空間模擬方法較IDW、Kriging、Spline等插值方法的精度較優(yōu)。其R2為0.953 8，分別大于IDW(0.818 6)、Kriging(0.840 7)、Spline(0.753 2)插值方法的模擬結(jié)果，并且其模擬結(jié)果隨海拔、經(jīng)緯度的變化趨勢(shì)和精細(xì)度更明顯，體現(xiàn)了海拔、經(jīng)緯度對(duì)積溫的顯著影響。

(3)本文在進(jìn)行≥10℃積溫空間模擬時(shí)，雖然考慮了海拔、經(jīng)度和緯度的影響，使模擬結(jié)果的精度和精細(xì)化程度都有所提高。但也存在一些不足，比如未考慮坡度、坡向、地形遮蔽、坡地反射輻射、散射輻射各向異性及天氣過(guò)程等因素對(duì)積溫的影響，同時(shí)，還應(yīng)進(jìn)一步加強(qiáng)遙感衛(wèi)星資料在本研究中的應(yīng)用。

參考文獻(xiàn)

[1] http：//www.ipcc.ch/report/ar5/wg1/[EB/OL].

[2] 柏秦鳳，霍治國(guó)，李世奎，等. 1978年前、后中國(guó)≥10℃年積溫對(duì)比[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2008，19(8)：1 810-1 816.

[3] 廖順寶，李澤輝. 積溫?cái)?shù)據(jù)柵格化方法的實(shí)驗(yàn)[J]. 地理研究，2004，23(5)：633-640.

[4] 袁淑杰，谷曉平，向紅瓊，等. 基于GIS的貴州高原復(fù)雜地形下積溫的精細(xì)空間分布[J]. 資源科學(xué)，2010，32(12)：2427-2432.

[5] 李飛，王春，趙軍，等. 中國(guó)陸地多年平均積溫空間化研究[J]. 自然資源學(xué)報(bào)，2010，25(5)：778-784.

[6] 于貴瑞，何洪林，劉新安，等. 中國(guó)陸地生態(tài)信息空間化技術(shù)研究(Ⅰ)——?dú)庀?氣候信息的空間化技術(shù)途徑[J]. 自然資源學(xué)報(bào)，2004，19(4)：537-544.

[7] 顧駿強(qiáng)，楊軍. 中國(guó)華南地區(qū)氣候和環(huán)境變化特征及其對(duì)策[J]. 資源科學(xué)，2005，27(1)：128-135.

[8] 李勇， 楊曉光， 王文峰， 等. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅰ：華南地區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2010，21(10)：2 605-2 614.

[9] 代姝瑋，楊曉光，趙孟，等. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅱ：西南地區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2011，22(2)：442-452.

[10] 趙錦，楊曉光，劉志娟，等. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅱ：南方地區(qū)氣候要素變化特征及對(duì)種植制度界限可能影響[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2010，43(9)：1 860-1 867.

[11] 李勇，楊曉光，王文峰，等. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅴ：氣候變暖對(duì)中國(guó)熱帶作物種植北界和寒害風(fēng)險(xiǎn)的影響分析[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2010，43(12)：2 477-2 484.

[12] 駱繼賓. 我國(guó)農(nóng)業(yè)氣候資源的開(kāi)發(fā)利用與農(nóng)業(yè)區(qū)域開(kāi)發(fā)[J]. 農(nóng)業(yè)區(qū)劃，1990，11(5)：44-47.

[13] 劉勤，嚴(yán)昌榮，何文清，等. 黃河流域近40a積溫動(dòng)態(tài)變化研究[J]. 自然資源學(xué)報(bào)，2009，24(1)：147-153.

[14] 曲曼麗. 農(nóng)業(yè)氣候?qū)嵙?xí)指導(dǎo). 北京： 北京農(nóng)業(yè)大學(xué)出版社，1991.

[15] 徐建華. 現(xiàn)代地理學(xué)中的數(shù)學(xué)方法(第二版)[M]. 北京：高等教育出版社， 2002.

[16] 劉新安，于貴瑞，范遼生，等. 中國(guó)陸地生態(tài)信息空間化技術(shù)研究(Ⅲ)——溫度、降水等氣候要素[J]. 自然資源學(xué)報(bào)，2004，19(6)：818-825.

[17] 高仲亮，周汝良，李智，等. 積溫分布的空間模擬及三維可視化初探[J]. 廣東農(nóng)業(yè)科學(xué)，2012，16：182-184，237-238.

endprint

2.3與其它插值方法對(duì)比分析

由圖3、4可知，對(duì)比多元線性回歸模型模擬的積溫空間分布結(jié)果與常用的空間化方法IDW、Kriging、Spline等插值方法結(jié)果，發(fā)現(xiàn)4種方法對(duì)≥10℃積溫總體分布趨勢(shì)表現(xiàn)基本一致，都是隨緯度的增大而減小，呈現(xiàn)出由南向北逐漸降低的趨勢(shì)。但是基于多元線性回歸模型的≥10℃積溫空間模擬方法的精度較優(yōu)，主要表現(xiàn)在以下幾方面。

(1)≥10℃積溫隨海拔的變化趨勢(shì)和精細(xì)度更明顯，體現(xiàn)了海拔對(duì)活動(dòng)積溫的顯著影響。

(2)對(duì)預(yù)留10%(16個(gè))氣象站點(diǎn)的模擬值與觀測(cè)值線性一致性更強(qiáng)，R2為0.953 8，分別大于IDW(0.818 6)、Kriging(0.840 7)、Spline(0.753 2)對(duì)插值方法的模擬結(jié)果(表1)，并且均方根誤差(Root-mean-square Error，RMSE)和相對(duì)均方根誤差(Relative Root-mean-square Error，RMSEr)都小于IDW、Kriging、Spline插值方法(表1)，表明其模擬精度高、誤差小。

3結(jié)論與討論

本文選用華南地區(qū)111個(gè)及周邊地區(qū)39個(gè)氣象站點(diǎn)1980～2011年逐日平均氣溫實(shí)測(cè)數(shù)據(jù)及各氣象臺(tái)站地名、經(jīng)度、緯度和海拔高度數(shù)據(jù)，利用多元線性回歸模型和IDW、Kriging、Spline等插值方法模擬了華南地區(qū)≥10℃積溫的空間分布，得出結(jié)論如下。

(1)通過(guò)分析1980～2011年華南地區(qū)≥10℃積溫的空間分布情況，及其隨海拔、經(jīng)緯度的變化情況，可知研究區(qū)≥10℃積溫最低值出現(xiàn)在云南省北部地區(qū)，最高值出現(xiàn)在海南島；并且≥10℃積溫隨海拔的升高而減小，隨緯度的增大而減小，隨經(jīng)度的變化幅度較大。在研究區(qū)西部云南地區(qū)受地形的影響變化幅度強(qiáng)烈，在研究區(qū)東部廣西、廣東、福建地區(qū)變化幅度較小。

(2)基于多元線性回歸模型的≥10℃積溫空間模擬方法較IDW、Kriging、Spline等插值方法的精度較優(yōu)。其R2為0.953 8，分別大于IDW(0.818 6)、Kriging(0.840 7)、Spline(0.753 2)插值方法的模擬結(jié)果，并且其模擬結(jié)果隨海拔、經(jīng)緯度的變化趨勢(shì)和精細(xì)度更明顯，體現(xiàn)了海拔、經(jīng)緯度對(duì)積溫的顯著影響。

(3)本文在進(jìn)行≥10℃積溫空間模擬時(shí)，雖然考慮了海拔、經(jīng)度和緯度的影響，使模擬結(jié)果的精度和精細(xì)化程度都有所提高。但也存在一些不足，比如未考慮坡度、坡向、地形遮蔽、坡地反射輻射、散射輻射各向異性及天氣過(guò)程等因素對(duì)積溫的影響，同時(shí)，還應(yīng)進(jìn)一步加強(qiáng)遙感衛(wèi)星資料在本研究中的應(yīng)用。

參考文獻(xiàn)

[1] http：//www.ipcc.ch/report/ar5/wg1/[EB/OL].

[2] 柏秦鳳，霍治國(guó)，李世奎，等. 1978年前、后中國(guó)≥10℃年積溫對(duì)比[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2008，19(8)：1 810-1 816.

[3] 廖順寶，李澤輝. 積溫?cái)?shù)據(jù)柵格化方法的實(shí)驗(yàn)[J]. 地理研究，2004，23(5)：633-640.

[4] 袁淑杰，谷曉平，向紅瓊，等. 基于GIS的貴州高原復(fù)雜地形下積溫的精細(xì)空間分布[J]. 資源科學(xué)，2010，32(12)：2427-2432.

[5] 李飛，王春，趙軍，等. 中國(guó)陸地多年平均積溫空間化研究[J]. 自然資源學(xué)報(bào)，2010，25(5)：778-784.

[6] 于貴瑞，何洪林，劉新安，等. 中國(guó)陸地生態(tài)信息空間化技術(shù)研究(Ⅰ)——?dú)庀?氣候信息的空間化技術(shù)途徑[J]. 自然資源學(xué)報(bào)，2004，19(4)：537-544.

[7] 顧駿強(qiáng)，楊軍. 中國(guó)華南地區(qū)氣候和環(huán)境變化特征及其對(duì)策[J]. 資源科學(xué)，2005，27(1)：128-135.

[8] 李勇， 楊曉光， 王文峰， 等. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅰ：華南地區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2010，21(10)：2 605-2 614.

[9] 代姝瑋，楊曉光，趙孟，等. 氣候變化背景下中國(guó)農(nóng)業(yè)氣候資源變化Ⅱ：西南地區(qū)農(nóng)業(yè)氣候資源時(shí)空變化特征[J]. 應(yīng)用生態(tài)學(xué)報(bào)，2011，22(2)：442-452.

[10] 趙錦，楊曉光，劉志娟，等. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅱ：南方地區(qū)氣候要素變化特征及對(duì)種植制度界限可能影響[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2010，43(9)：1 860-1 867.

[11] 李勇，楊曉光，王文峰，等. 全球氣候變暖對(duì)中國(guó)種植制度可能影響Ⅴ：氣候變暖對(duì)中國(guó)熱帶作物種植北界和寒害風(fēng)險(xiǎn)的影響分析[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2010，43(12)：2 477-2 484.

[12] 駱繼賓. 我國(guó)農(nóng)業(yè)氣候資源的開(kāi)發(fā)利用與農(nóng)業(yè)區(qū)域開(kāi)發(fā)[J]. 農(nóng)業(yè)區(qū)劃，1990，11(5)：44-47.

[13] 劉勤，嚴(yán)昌榮，何文清，等. 黃河流域近40a積溫動(dòng)態(tài)變化研究[J]. 自然資源學(xué)報(bào)，2009，24(1)：147-153.

[14] 曲曼麗. 農(nóng)業(yè)氣候?qū)嵙?xí)指導(dǎo). 北京： 北京農(nóng)業(yè)大學(xué)出版社，1991.

[15] 徐建華. 現(xiàn)代地理學(xué)中的數(shù)學(xué)方法(第二版)[M]. 北京：高等教育出版社， 2002.

[16] 劉新安，于貴瑞，范遼生，等. 中國(guó)陸地生態(tài)信息空間化技術(shù)研究(Ⅲ)——溫度、降水等氣候要素[J]. 自然資源學(xué)報(bào)，2004，19(6)：818-825.

[17] 高仲亮，周汝良，李智，等. 積溫分布的空間模擬及三維可視化初探[J]. 廣東農(nóng)業(yè)科學(xué)，2012，16：182-184，237-238.

endprint