徐吉宏，柳小妮，張德罡，徐吉偉

(1.甘肅省社會(huì)科學(xué)院，甘肅 蘭州 730070； 2.甘肅農(nóng)業(yè)大學(xué)草業(yè)學(xué)院，甘肅 蘭州 730070； 3.蘭州大學(xué)草地農(nóng)業(yè)科技學(xué)院，甘肅 蘭州 730020)

目前，通用的草地分類方法主要以氣候?yàn)闆Q定因素，植被為主體，結(jié)合生長地形、土壤特征，依據(jù)草地類型學(xué)原理，探求草地形成要素和分類指標(biāo)關(guān)系。任繼周和胡自治提出的草地綜合順序分類系統(tǒng)(COSGC)方法，很好地反映了草地植被與氣候的相互關(guān)系，其中的草地分類檢索圖在國際上被稱為任-胡氏檢索圖。

空間插值的實(shí)質(zhì)是通過已知樣點(diǎn)的數(shù)據(jù)來估算未知點(diǎn)的數(shù)據(jù)。長期以來，國內(nèi)外學(xué)者對(duì)空間離散數(shù)據(jù)的插值進(jìn)行了大量的研究，并提出了多種空間插值方法。常用于氣象要素空間插值的方法有多項(xiàng)式插值法(interpolating polynomials)、趨勢(shì)面法(trend surface)、反距離權(quán)重法(inverse distance weighted)、樣條函數(shù)法(spline)、克里金法(kriging)、自然鄰域插值法(natural neighbor)等[1-3]?？臻g插值方法各有利弊，不同的方法插值結(jié)果差異很大。如何根據(jù)數(shù)據(jù)與區(qū)域的特征選擇最優(yōu)的方法一直是研究的一個(gè)難題[4]。為此，根據(jù)可以獲取的氣象數(shù)據(jù)進(jìn)行空間插值方法的選擇成為一個(gè)重要的研究熱點(diǎn)。

本研究采用ArcGIS 9.1技術(shù)[5-7]，比較不同插值方法對(duì)天?？h草地進(jìn)行類的劃分的效果，并驗(yàn)證COSGC分類的適用性，為其進(jìn)一步推廣應(yīng)用奠定基礎(chǔ)。

1 研究區(qū)概況與研究方法

1.1研究區(qū)概況 天祝藏族自治縣(以下簡稱天?？h)位于甘肅中部，地理位置為102°07′～103°46′ E， 36°31′～ 37°55′ N。全縣東西寬124.6 km，南北長158.4 km，土地面積為6 900余km2，地勢(shì)西部高峻，向東南逐漸變低，境內(nèi)地形以山地為主，分布了祁連山系東緣的一系列山脈，地形復(fù)雜。高寒地區(qū)海拔高，氣溫、降水量差異很大，氣候變化明顯，氣溫低，年較差小，日較差大，降水分布不均勻，年際變化明顯，水熱條件較為嚴(yán)酷，易形成較為頻繁的災(zāi)害性天氣(風(fēng)、沙、旱、雪、霜)。日照充足，太陽輻射強(qiáng)烈，具有大陸性高原季風(fēng)氣候特征[8-12]。

1.2研究方法

1.2.1數(shù)據(jù)收集及處理

氣象數(shù)據(jù)與處理：收集了天?？h及周邊氣象站1961―1990年30年57個(gè)點(diǎn)的氣象數(shù)據(jù)資料(來源于天祝縣氣象局，甘肅省氣象局等相關(guān)部門)。依據(jù)COSGC模型，計(jì)算得到>0℃年積溫(∑θ)及濕潤度(K)。

由于天?？h境內(nèi)地形以山地為主，分布了祁連山系東緣的一系列山脈，海拔明顯高于周邊地區(qū)，單純利用天?？h及周邊地區(qū)的∑θ和K，通過空間插值所得到的分類指標(biāo)與實(shí)際差異很大。馬興祥等[12]研究認(rèn)為，積溫隨海拔增高而降低，二者呈線性相關(guān)關(guān)系。因此，結(jié)合天?？h數(shù)字高程模型(DEM)，對(duì)∑θ進(jìn)行修正(式1、2)，然后重新計(jì)算K，以提高∑θ和K指標(biāo)的精度，滿足天?？h草地分類的要求。

T=14.23-0.047h

(1)

∑T0=5 547.23-1.043 7h

(2)

式中，T為年平均氣溫，h為海拔高度,∑T0為≥0℃年積溫。

基礎(chǔ)圖庫：收集天?？h1986年100 m×100 m的DEM圖，1984年1∶250 000的行政區(qū)劃圖，天?？h遙感影像植被分布圖和土壤類型圖，以及2001年1∶100 000植被類型圖等紙質(zhì)圖[13-14]，并對(duì)其進(jìn)行掃描數(shù)字化，不同地圖統(tǒng)一定義為Krasovsky_1940_Albers投影。

天?？h行政區(qū)劃圖、天?？h遙感影像植被分布圖、土壤類型圖和植被類型圖等紙質(zhì)圖來源于甘肅省基礎(chǔ)地理信息中心和甘肅省天?？h農(nóng)業(yè)區(qū)劃委員會(huì)；DEM圖由甘肅省草原生態(tài)研究所遙感協(xié)作組和天祝縣草原工作站提供。

1.2.2草地綜合順序分類圖的生成 以地理信息系統(tǒng)和COSGC論為指導(dǎo)，實(shí)現(xiàn)∑θ和K數(shù)據(jù)的空間化，進(jìn)行空間數(shù)據(jù)疊置分析和重分類，生成COSGC第1級(jí)分類單位類圖。

點(diǎn)數(shù)據(jù)的空間化：結(jié)合天?？hDEM，整理并修正∑θ和K數(shù)據(jù)為jw.dbf和K.dbf格式，導(dǎo)入ArcGIS 9.1，利用Add XY data模塊，實(shí)現(xiàn)點(diǎn)數(shù)據(jù)的空間化表達(dá)。

分類指標(biāo)數(shù)據(jù)的空間化：將空間化的jw.dbf和K.dbf數(shù)據(jù)，導(dǎo)入ArcGIS 9.1，利用ArcToolBox-Spatial Analyst Tool-Interpolation和Geostatistial Analyst--Geostatistial Wizard模塊，分別用IDW、Spline、 Kriging， Natural Neighbor進(jìn)行空間插值，以期獲得連續(xù)有序的空間數(shù)據(jù)。

分類指標(biāo)空間化數(shù)據(jù)的重分類：按照COSGC的劃分體系，以天?？h的∑θ和K為指標(biāo)，進(jìn)行天?？h草地綜合順序分類。

空間疊置分析：采用矢量數(shù)據(jù)和柵格數(shù)據(jù)相結(jié)合的疊置分析方法。

1.2.3草地綜合順序分類圖的可信度分析 采用分區(qū)抽樣，選擇1984年中國南、北草地資源調(diào)查辦公室制作的天?？h遙感影像植被分布圖[15]，均勻選取350個(gè)點(diǎn)，與本研究生成的天祝草地綜合順序分類圖的同域點(diǎn)，通過兩圖同域點(diǎn)原生植被景觀類型異同的比較，判斷天祝縣草地綜合順序分類圖的準(zhǔn)確度[16]。

2 結(jié)果與分析

2.1空間插值 在ArcGIS 9.1軟件下，分別以IDW、Spline、Kriging、Natural Neighbor四種方法對(duì)∑θ和K數(shù)據(jù)進(jìn)行插值[17-21]，結(jié)果如圖1、2所示。

2.2天?？h草地綜合順序分類圖 將∑θ和K的不同空間插值結(jié)果進(jìn)行空間疊置，生成9個(gè)不同的草地綜合順序分類圖(圖3、表1)。

圖1 ∑θ插值結(jié)果

圖2 K插值結(jié)果

圖3 天?？h草地綜合順序分類

表1 草地綜合順序類序號(hào)說明

序號(hào)類1寒冷干旱半荒漠、高山半荒漠2寒冷微干極地干凍原、高山草地3寒冷微潤少雨凍原、高山草甸草地4寒冷濕潤凍原、高山草甸5寒冷潮濕多雨凍原、高山草甸6寒溫干旱山地半荒漠7寒溫微干山地草地8寒溫微潤山地草甸草地9寒溫濕潤山地草甸10寒溫潮濕寒溫性針葉林11微溫干旱溫帶半荒漠12微溫微干溫帶典型草地13微溫微潤草甸草地14微溫濕潤森林草地15微溫潮濕針葉闊葉混交林

2.3草地綜合順序分類圖的可信度分析 天?？h9個(gè)草地綜合順序分類圖(圖3)與天?？h遙感影像植被分布圖的同域點(diǎn)屬性對(duì)照結(jié)果如表2所示。

在350個(gè)對(duì)照點(diǎn)中，圖F與天祝縣遙感影像植被分布圖的植被景觀相同點(diǎn)最多，共有256個(gè)，占73.14%，不同點(diǎn)94個(gè)，僅占26.86%。因此，在ArcGIS 9.1環(huán)境下，∑θ使用Natural Neighbor，K使用Kriging插值方法生成的天祝縣草地綜合順序分類圖具有較高的可信度。

表2 草地綜合順序分類圖的可信度分析結(jié)果

3 結(jié)論

在GIS軟件支持下，以COSGC為指導(dǎo)，利用天?？h1961-1990年30年的降水和氣溫資料，對(duì)>0℃年積溫(∑θ)和濕潤度(K)分別采用IDW、Spline、Kriging和Natural Neighbor四種方法進(jìn)行空間插值，生成天?？h草地綜合順序分類圖。

利用疊加分析生成9類草地綜合順序分類圖與天?？h遙感影像植被分布圖進(jìn)行屬性對(duì)照[16]，結(jié)果表明Natural Neighbor最適合天祝縣∑θ插值，Kriging最適合天?？hK插值。

天?？h草地可劃分為12個(gè)草地綜合順序類，分別是寒冷微干極地干凍原、高山草地；寒溫微干山地草地；微溫微干溫帶典型草地；寒冷微潤少雨凍原、高山草甸草地；寒溫微潤山地草甸草地；微溫微潤草甸草地；寒冷濕潤凍原、高山草甸；寒溫濕潤山地草甸；微溫濕潤森林草地；寒冷潮濕多雨凍原、高山草甸；寒溫潮濕寒溫性針葉林；微溫潮濕針葉闊葉混交林草地。

[1] Evan J.England Spatial Simulation：Environmental Application,Environmental Modeling with GIS[M]．New York：Oxford University Press,1993:1-66.

[2] Price D T,Mckenney D W,Nald I A,etal.A comparison of two statistical methods for interpolation of Canadian monthly mean climate data[J].Agriculture for Meteorology,2000,101:81-94.

[3] 鄔倫，劉瑜，張晶，等，地里信息系統(tǒng)[M]．北京：北京科學(xué)出版社，2003:178-191.

[4] 劉登偉，封志明，楊艷昭．海河流域降水空間插值方法的選取[J]．地球信息科學(xué)，2006,8(4):75-79.

[5] 馬軒龍，李文娟，陳全功．基于GIS與草原綜合順序分類法對(duì)甘肅省草地類型的劃分初探[J]．草業(yè)科學(xué)，2009,26(5):7-13.

[6] 李紅梅，馬玉壽．基于GIS技術(shù)的青海省草地類型分類研究[J]．草業(yè)科學(xué)，2009,26(12):24-29.

[7] 李紅梅，馬玉壽．改進(jìn)的綜合順序分類法在青海草地分類中的應(yīng)用[J]．草業(yè)學(xué)報(bào)，2009,18(2):76-82.

[8] 王樹青．甘肅省天?？h草地資源及其開發(fā)利用[J]．中國草地，1999,16(6):65-67.

[9] 廖國藩，賈幼陵．中國草地資源[M]．北京：中國科學(xué)技術(shù)出版社，1996:562-563.

[10] 王樹青,張起榮，馬蒼．天?？h天然草地退化原因及治理對(duì)策[J]．草業(yè)科學(xué)，2003,20(6):7-8.

[11] 蔡有壽，張國強(qiáng)．淺談天祝自然條件與農(nóng)業(yè)生產(chǎn)[J]．甘肅農(nóng)業(yè)，2002(5):20-22.

[12] 馬興祥，劉明春，張子廉，等．天祝縣氣候變化對(duì)畜牧業(yè)可持續(xù)發(fā)展的影響及對(duì)策[J]．草業(yè)科學(xué)，2000，17(5):25-30.

[13] 甘肅省地圖院． 武威市地圖集[M]．西安：西安地圖出版社，2006.

[14] 天?？h區(qū)劃辦公室．甘肅省天祝藏族自治縣農(nóng)牧業(yè)區(qū)劃資料匯編[M]．蘭州：甘肅科學(xué)技術(shù)出版社，1986：26.

[15] 甘肅草地生態(tài)研究所遙感協(xié)作組，甘肅省天祝藏族自治縣農(nóng)業(yè)區(qū)劃委員會(huì)．天祝藏族自治縣農(nóng)牧業(yè)自然條件遙感影象系列圖[Z]．蘭州：中國草原與牧草、蘭州大學(xué)學(xué)報(bào)，1984.

[16] 陳全功，任繼周，王珈誼．中國草業(yè)開發(fā)與生態(tài)建設(shè)專家系統(tǒng)．北京：電子工業(yè)出版社，2006.

[17] 潘耀忠，龔道溢，鄧?yán)?等．基于DEM的中國陸地多年平均溫度插值方法[J]．地理學(xué)報(bào)，2004,59(3):336-374.

[18] 張洪亮，倪紹祥，鄧自旺，等．基于DEM的山區(qū)氣溫空間模擬方法[J]．山地學(xué)報(bào)，2002,20(3):360-364.

[19] 黃興元，馬勁松，湯勤．地理信息系統(tǒng)概論[M]．北京：科學(xué)出版社，2001.

[20] 徐吉宏，柳小妮，張德罡．天?？h草地綜合順序分類及生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)價(jià)[J]．中國草地學(xué)報(bào)，2009,31(5):23-29.

[21] 聶欣，吳良林，趙玲．利用DEM數(shù)據(jù)分析山區(qū)氣溫[J]．安徽農(nóng)業(yè)科學(xué)，2009,37(5):2187-2188.