， ，

(1.山東科技大學(xué) 測(cè)繪科學(xué)與工程學(xué)院，山東 青島 266590； 2.國家海洋信息中心，天津 300171)

草地不僅是畜牧業(yè)經(jīng)濟(jì)的飼料來源，也是大部分野生動(dòng)物生活的場(chǎng)所，對(duì)于全球碳循環(huán)和氣候調(diào)節(jié)具有重要作用，是重要的國土資源[1]。三江源自然保護(hù)區(qū)是青藏高原腹地，是我國生態(tài)環(huán)境保護(hù)和可持續(xù)發(fā)展的重要保證，三江源強(qiáng)大的生態(tài)功能是其生物多樣性以及植被豐富度的重要保證[2]。

我國對(duì)草地的遙感研究始于20世紀(jì)80年代初，至今已取得了多方面研究成果。周鑫等[3]針對(duì)新疆伊犁地區(qū)的植被指數(shù)建立了回歸模型，分析了地上生物量的空間分布；趙強(qiáng)[4]基于遙感技術(shù)反演了呼倫貝爾草原草產(chǎn)量；董建軍[5]以內(nèi)蒙古錫林河流域的典型草原草場(chǎng)為對(duì)象，通過建模與比較，探討了多源遙感數(shù)據(jù)在典型草原植被估產(chǎn)方面的潛力；方金[6]利用回歸分析模擬出甘南草地生草場(chǎng)物量遙感模型，反演出甘南地區(qū)的草地生物量；Kawamura等[7]研究了中分辨率成像光譜儀(moderate-resolution imaging spectroradiometer，MODIS)歸一化植被指數(shù)與內(nèi)蒙古草地生物量之間的關(guān)系，發(fā)現(xiàn)二者存在顯著的相關(guān)性；Field等[8]利用卡內(nèi)基艾姆期坦福方法(Carnegie-Ames-Stanford Approach，CASA)模型估算出全球陸地年均凈初級(jí)生產(chǎn)力(net primary productivity，NPP)值為56.4 Pg/a，其中草地占4.3%；Nayak等[9]通過CASA模型估算印度的草地NPP為267 gC/m2·a，占總NPP的1.2%；Peterson等[10]利用Landsat/TM數(shù)據(jù)區(qū)別道格拉斯草地類型，得出7月遙感數(shù)據(jù)的評(píng)價(jià)精度高達(dá)81.8%。以上模型大都應(yīng)用SPSS等軟件進(jìn)行實(shí)測(cè)數(shù)據(jù)和遙感數(shù)據(jù)回歸分析，利用最優(yōu)模型反演研究區(qū)草地生物量，泛化性較差，同時(shí)沒有對(duì)海拔高、高差大的草原地區(qū)草產(chǎn)量遙感估算研究。因此，本文利用三江源地區(qū)2000—2010年MODIS數(shù)據(jù)，借助氣象、太陽輻射數(shù)據(jù)和數(shù)字高程模型數(shù)據(jù)，基于改進(jìn)的CASA模型反演草地畝產(chǎn)量，為高海拔、高差大的三江源地區(qū)草地生物量遙感監(jiān)測(cè)、草場(chǎng)資源的可持續(xù)利用和畜牧業(yè)發(fā)展等相關(guān)研究提供借鑒。

1 研究基礎(chǔ)

1.1 研究區(qū)概況

表1 數(shù)據(jù)列表

三江源自然保護(hù)區(qū)位于我國西北部，平均海拔在4 800 m左右。地理位置為31°39′～36°12′N，89°45′～102°23′E[11]。屬于大陸性氣候，由于海拔高，地形復(fù)雜，氣候變化多樣，年際溫差小，降雨年際變化大，太陽輻射強(qiáng)，造成三江源地區(qū)東西植被分布差別大。三江源地區(qū)是世界上海拔最高，分布植被類型最多的典型高寒地帶，是我國的高寒生物自然種質(zhì)資源庫[12]。三江源地區(qū)的代表草種類型有草甸、草原、沼澤及水生草地、墊狀和稀疏草地等。三江源地區(qū)草地草場(chǎng)生長期集中在每年4—9月份，因此，本文選擇4—9月為草場(chǎng)生物量累積時(shí)間。

1.2 數(shù)據(jù)源獲取

采用遙感數(shù)據(jù)和非遙感數(shù)據(jù)相結(jié)合的方法獲取數(shù)據(jù)，然后將反演結(jié)果與平臺(tái)提供數(shù)據(jù)進(jìn)行驗(yàn)證(表1、表2)。

注：西藏草場(chǎng)屬性數(shù)據(jù)參考自國家地球系統(tǒng)科學(xué)數(shù)據(jù)共享平臺(tái)

1.2.1 MODIS凈初級(jí)生產(chǎn)力數(shù)據(jù)——MOD13A1

MOD13A1遙感影像是MODIS數(shù)據(jù)的L3級(jí)產(chǎn)品，數(shù)據(jù)類型為柵格，其中包括歸一化植被指數(shù)(normadized differential，NDVI)和增強(qiáng)型植被指數(shù)(enhanced vegetation index，EVI)，提供時(shí)間分辨率為16天、幾何分辨率為500 m的3級(jí)產(chǎn)品。

1.2.2 氣象數(shù)據(jù)

氣象數(shù)據(jù)主要包括青海省SOL(太陽總凈輻射量，total net radiation of solar)，2000—2010年4—9月份平均溫度和降雨量數(shù)據(jù)。根據(jù)CASA模型需要，提取植被最適宜溫度和水分脅迫因子。三江源地區(qū)溫度和降雨分布成反比，海拔高的地區(qū)溫度偏低，降雨偏多；太陽輻射量分布從東到西呈逐漸增多的趨勢(shì)，與海拔呈正比。

1.2.3 高程數(shù)據(jù)

因?yàn)槿吹貐^(qū)以山地地形為主，地勢(shì)復(fù)雜，數(shù)字高程模型數(shù)據(jù)顯得尤為重要，尤其在土地覆蓋類型劃分過程中。本文運(yùn)用航天飛機(jī)雷達(dá)地形測(cè)繪使命-全球數(shù)字高程模型(shuttle radar topography mission-global digital elevation model，SRTM3-GDEM)數(shù)據(jù)，對(duì)應(yīng)的分辨率精度為90 m。

2 NPP模型的建立與計(jì)算

改進(jìn)的CASA模型[13]中植被凈第一性生產(chǎn)力NPP由牧草的光合有效輻射(APAR，absorbed photosynthetic active radiation)和實(shí)際光能利用率(ε)來決定，其估算公式如下：

NPP(x,t)=APAR(x,t)×ε(x,t)。

(1)

APAR是指在單位時(shí)間內(nèi)牧草吸收的有效輻射能量，gC·m-2·month-1；ε表示在單位時(shí)間內(nèi)實(shí)際的光能利用率，gC·MJ-1。

2.1 光和有效輻射APAR的估算

牧草的太陽有效輻射和對(duì)光合有效輻射的吸收系數(shù)決定APAR的大小，進(jìn)而決定著NPP的變化。

APAR(x,t)=SOL(x,t)×FPAR(x,t)×0.5。

(2)

其中，SOL表示牧草的太陽總輻射量，gC·m-2·month-1；FPAR(fraction of absorbed photosynthetically active radiation，光合有效輻射系數(shù))表示牧草對(duì)入射光合有效輻射的吸收比例，式中的0.5表示這一類型牧草所能利用的太陽有效輻射占太陽在這一區(qū)域總輻射的比例。

2.2 光合有效輻射系數(shù)FPAR的估算

光合有效輻射系數(shù)的變化和植被歸一化指數(shù)NDVI之間存在如下關(guān)系：

(3)

其中，NDVImax和NDVImin分別為三江源內(nèi)不同草場(chǎng)類型的NDVI最大值和最小值。

光合有效輻射系數(shù)(FPAR)和比值植被指數(shù)(simple ratio index，SR)也存在著一定關(guān)系：

(4)

其中，F(xiàn)PARmax和FPARmin為常數(shù)，與草場(chǎng)類型無關(guān)，分別為0.001和0.95[14]；SRmax和SRmin分別對(duì)應(yīng)三江源地區(qū)牧草的NDVI的95%和5%下側(cè)百分位數(shù)，SR可以通過NDVI得到：

(5)

因?yàn)镹DVI和SR對(duì)牧草的FPAR估算均有其優(yōu)缺點(diǎn)，為了彌補(bǔ)兩種指數(shù)的不足，取其加權(quán)平均值作為FPAR的估算值：

FPAR(x,t)=αFPARNDVI+(1-α)FPARSR。

(6)

一般取α為0.5[14]。

2.3 光能利用率ε的估算

光能利用率是指在單位時(shí)間單位面積上生產(chǎn)的生物量與單位時(shí)間內(nèi)照到該面積上的光合有效輻射能量的比值。

ε(x,t)=Tε1(x,t)×Tε2(x,t)×Wε(x,t)×εmax，

(7)

其中，Tε1和Tε2為溫度脅迫系數(shù)，表示溫度變化對(duì)光能利用率的影響；Wε為水分脅迫系數(shù)，表示環(huán)境的濕度變化對(duì)光能利用率的影響；εmax是牧草生長最適宜條件下的最大光能利用率。

2.4 溫度脅迫因子的估算

Tε1的估算：只有在適宜溫度下植被才能最大化進(jìn)行光合作用，而溫度過高或者過低都會(huì)抑制光合作用的進(jìn)行，溫度脅迫因子Tε1的估算：

Tε1(x,t)=0.8+0.02×Topt(x)-0.000 5×[Topt(x)]2，

(8)

其中，Topt是三江源草地牧草生長的最適宜溫度，是指一年內(nèi)NDVI達(dá)到最高時(shí)月份的平均氣溫。溫度脅迫因子Tε2的估算：表示牧草在溫度變化的影響，溫度向過高或者過低變化，光能利用率都會(huì)下降，因此對(duì)偏離最適宜溫度的影響因素進(jìn)行估算：

(9)

當(dāng)某一月平均溫度T(x,t)比最適溫度Topt(x)高10 ℃或低13 ℃時(shí)，該月的Tε2(x,t)值等于月平均溫度T(x,t)為最適溫度Topt(x)時(shí)Tε2(x,t)值的一半。

2.5 水分脅迫因子的估算

水分脅迫系數(shù)Wε表示周圍環(huán)境濕度變化對(duì)牧草生長的影響，進(jìn)而對(duì)光能利用率的影響，取值范圍為0.5～1，0.5代表最干旱的情況，1為最濕潤的情況。

用全年降水量P和>0 ℃年積溫∑θ之比表示濕潤度指標(biāo)K，

(10)

為適應(yīng)改進(jìn)的CASA模型，保證牧草反演的精度，采用月降雨量和>0 ℃月積溫之比估算濕潤度指標(biāo)K，

(11)

根據(jù)上面式子計(jì)算水分脅迫因子W(x,t):

(12)

2.6 CASA模型的優(yōu)化

最大光能利用率εmax定義為0.389 gC·MJ-1[15]，取值受到溫度、水分、土壤和植物個(gè)體發(fā)育等因素的顯著影響，將其作為一個(gè)常數(shù)在全球范圍內(nèi)使用會(huì)引起很大的誤差，εmax還和草場(chǎng)類型有關(guān)，不同草場(chǎng)的優(yōu)勢(shì)草種也不同。

三江源地區(qū)植被光能利用率與植被分布和地形有關(guān)，其分布特征具有明顯地帶性，呈由西北向東南逐漸遞增的規(guī)律，隨季節(jié)變化明顯，本次研究基于三江源氣象因素(溫度、降水、太陽輻射等數(shù)據(jù))和地形因素，同時(shí)參考衛(wèi)亞星等[16]對(duì)三江源地區(qū)草場(chǎng)植被的研究，利用加權(quán)平均值計(jì)算平均最大光能利用率為0.096 gC·MJ-1。

根據(jù)得到的全年NPP和三江源植被覆蓋類型計(jì)算每種類型草場(chǎng)的畝產(chǎn)量分布圖，與國家地球系統(tǒng)科學(xué)數(shù)據(jù)共享平臺(tái)提供的牧草畝產(chǎn)量數(shù)據(jù)進(jìn)行精度驗(yàn)證，并做進(jìn)一步分析。

3 NPP時(shí)空分析

基于改進(jìn)的CASA模型，應(yīng)用MOD13A1遙感影像和氣象數(shù)據(jù)計(jì)算得到2000—2010年三江源地區(qū)NPP分布，從11年NPP分布可以看出，中南部玉樹藏族自治州、果洛藏族自治州和黃南藏族自治州等地區(qū)NPP值普遍在150 gC/m2以上，西部海西蒙古藏族自治州等地區(qū)普遍在50 gC/m2。圖1為2008年NPP平均值分布圖。

圖1 2008年NPP平均值分布圖

3.1 NPP均值空間分布

三江源地區(qū)草地草場(chǎng)NPP平均值空間變化呈東南—西北遞減規(guī)律，高寒牧草地NPP值較低，河密集、海拔較低的林地NPP值相對(duì)較高。其中，研究區(qū)西部地區(qū)NPP值偏低，如治多縣和格爾木市等地區(qū)，NPP值50～70 gC/m2；研究區(qū)中部70～120 gC/m2；研究區(qū)東南部的囊謙縣、玉樹縣，以及東部的達(dá)日縣、甘德縣、久治縣等地區(qū)NPP值大于150 gC/m2。從空間分布特征得出，NPP分布與地形、海拔、河流及植被類型等自然因素相一致。

3.2 NPP均值時(shí)間分布

圖2展示了三江源地區(qū)NPP年平均值變化，2000—2010年的變化區(qū)間為133.03～191.91 gC/m2，整體趨勢(shì)為逐年增加，尤其在2007—2010年，增加速度很快。從2000—2001年三江源地區(qū)NPP增速較快，約為18 gC/m2；從2001—2007年相對(duì)較穩(wěn)定，為140～150 gC/m2；從2007—2010年增速大幅上升，尤其在2009—2010年，增速達(dá)40 gC/m2。這種改善得益于三江源地區(qū)在2000年之后實(shí)施的各種生態(tài)工程項(xiàng)目，尤其在2007年之后效果明顯。

3.3 NPP均值變化趨勢(shì)

根據(jù)NPP均值年際變化折線圖回歸分析，分別得到線性、二次多項(xiàng)式和三次多項(xiàng)式模型(表3)，其中線性和二次多項(xiàng)式誤差較大，其R2分別為0.44和0.58，三次多項(xiàng)式模型很好地模擬出變化趨勢(shì)，R2達(dá)0.91。可以看出，從2010年往后的趨勢(shì)為持續(xù)增加。

圖2 2000—2010年NPP均值變化圖

模型R2線性模型 Y=2.91X+133.570.44二次多項(xiàng)式模型Y=0.6X2-4.25X+149.080.58三次多項(xiàng)式模型Y=0.335X3-5.44X2+25.99X+112.470.91

3.4 草場(chǎng)類型劃分

三江源草地草場(chǎng)劃分主要依據(jù)國家科技資源共享服務(wù)平臺(tái)數(shù)據(jù)實(shí)體草場(chǎng)植被覆蓋度[17]要素和每種草場(chǎng)優(yōu)勢(shì)草種的生長海拔兩個(gè)基本條件，結(jié)合優(yōu)勢(shì)草種生長的溫度和降雨要素，利用決策樹對(duì)三江源草地進(jìn)行分類，得到圖3三江源草地草場(chǎng)分布圖。

畝產(chǎn)量最高的湖盆河灘草甸草場(chǎng)分布在中南部和東部地區(qū)；山地稀疏森林草場(chǎng)、高原寬谷草原草場(chǎng)和山地灌叢草場(chǎng)畝產(chǎn)量中等分布在中部地區(qū)；高山草甸草場(chǎng)畝產(chǎn)量中等偏少，分布在西南地區(qū)；高原寬谷荒漠草場(chǎng)、山地荒漠草場(chǎng)和山地草原草場(chǎng)畝產(chǎn)量最低，分布在西部地區(qū)。

圖3 三江源地區(qū)草場(chǎng)類型分布圖

3.5 草地畝產(chǎn)量驗(yàn)證與分析

如圖4所示，分別顯示了8種草場(chǎng)2000—2010年畝產(chǎn)量的變化，可以看出各類草場(chǎng)的畝產(chǎn)量相對(duì)穩(wěn)定，變化不大。

表4顯示了遙感反演畝產(chǎn)量與國家生態(tài)系統(tǒng)觀測(cè)平臺(tái)的實(shí)際監(jiān)測(cè)畝產(chǎn)量數(shù)據(jù)，從表4反演準(zhǔn)確性看，反演畝產(chǎn)量與實(shí)際畝產(chǎn)量相比：高山草甸草場(chǎng)、山地稀疏森林草場(chǎng)和山地灌叢草場(chǎng)非常吻合；山地草原草場(chǎng)基本吻合；湖盆河灘草甸草場(chǎng)產(chǎn)量偏低，原因是NDVI對(duì)于植被覆蓋度太高的牧草產(chǎn)生了過飽和現(xiàn)象，對(duì)湖盆河灘草甸草場(chǎng)畝產(chǎn)量造成過低估計(jì)；高原寬谷荒漠草場(chǎng)和山地荒漠草場(chǎng)產(chǎn)量偏高，原因是NDVI對(duì)于低覆蓋植被出現(xiàn)過高估計(jì)現(xiàn)象。

圖4 2000—2010年三江源草地畝產(chǎn)量變化圖

草場(chǎng)類型反演畝產(chǎn)量/kg實(shí)際畝產(chǎn)量/kg吻合程度變化趨勢(shì)高山草甸草場(chǎng)110～13035～125非常吻合相對(duì)穩(wěn)定湖盆河灘草甸草場(chǎng)160～180200～300產(chǎn)量偏低略有上升高原寬谷草原草場(chǎng)70～80相對(duì)穩(wěn)定山地草原草場(chǎng)48～565～50基本吻合很穩(wěn)定高原寬谷荒漠草場(chǎng)39～4510產(chǎn)量偏高很穩(wěn)定山地稀疏森林草場(chǎng)130～160>125非常吻合略有下降山地荒漠草場(chǎng)40～5020產(chǎn)量偏高很穩(wěn)定山地灌叢草場(chǎng)128～160>125非常吻合略有下降

注：實(shí)際畝產(chǎn)量數(shù)據(jù)參考自國家地球系統(tǒng)科學(xué)數(shù)據(jù)共享平臺(tái)

圖5 2000—2010年三江源草地總產(chǎn)量變化圖

從畝產(chǎn)量變化趨勢(shì)上看，山地草原草場(chǎng)、山地荒漠草場(chǎng)和高原寬谷荒漠草場(chǎng)等低產(chǎn)量草場(chǎng)畝產(chǎn)量很穩(wěn)定；高山草甸草場(chǎng)和高原寬谷草原草場(chǎng)屬中等產(chǎn)量草場(chǎng)相對(duì)穩(wěn)定；湖盆河灘草甸草場(chǎng)畝產(chǎn)量逐年增加；山地稀疏森林草場(chǎng)和山地灌叢草場(chǎng)略有下降。

由圖5可以看出，從2000—2007年三江源草地總產(chǎn)量先增加后穩(wěn)定又下降，從2007—2010年快速增長，這一改善與三江源地區(qū)管理部門對(duì)三江源保護(hù)區(qū)的治理密不可分。自從三江源草地出現(xiàn)草地退化現(xiàn)象，管理部門就開始大力投資退牧換草、圍欄放牧等措施，2007年開始效果逐步體現(xiàn)。

4 結(jié)論

基于改進(jìn)的CASA模型，利用2000—2010年三江源地區(qū)MOD13A1植被指數(shù)3級(jí)數(shù)據(jù)以及對(duì)應(yīng)的氣象數(shù)據(jù)，計(jì)算得到三江源地區(qū)不同類型草場(chǎng)的畝產(chǎn)量分布，得出以下結(jié)論：

1) 空間上，草場(chǎng)畝產(chǎn)量分布特征呈東南—西北逐漸降低，西部地區(qū)畝產(chǎn)量偏低，40～60 kg/畝，中部地區(qū)60～120 kg/畝，東部地區(qū)畝產(chǎn)量一般大于120 kg/畝。高值出現(xiàn)在東南地區(qū)海拔偏低、地勢(shì)平坦的濕潤地區(qū)，低值一般出現(xiàn)在西北地區(qū)海拔高的干旱高寒地區(qū)。進(jìn)一步說明了三江源地區(qū)NPP的分布東西差異大，受海拔、地勢(shì)及環(huán)境濕潤程度影響大。

2) 時(shí)間上，不同類型草地畝產(chǎn)量整體變化不大，高產(chǎn)量草場(chǎng)畝產(chǎn)量一般在120～180 kg/畝，低產(chǎn)量草場(chǎng)畝產(chǎn)量一般在40～60 kg/畝，2000—2010年草場(chǎng)整體狀況比較穩(wěn)定，草場(chǎng)總產(chǎn)量整體呈上升趨勢(shì)，說明對(duì)草場(chǎng)的保護(hù)措施明顯改善。

3) 三江源草場(chǎng)草產(chǎn)量變化的主要驅(qū)動(dòng)力是人類活動(dòng)，管理部門采取退牧換草、圍欄放牧等措施大大提高了草產(chǎn)量，改善了三江源草場(chǎng)生態(tài)環(huán)境。

4) 基于改進(jìn)的CASA模型對(duì)三江源地區(qū)2000—2010年草地畝產(chǎn)量的遙感反演結(jié)果非常理想，不同類型草場(chǎng)的畝產(chǎn)量基本在給定的實(shí)測(cè)范圍內(nèi)，有望在以后代替落后人工實(shí)測(cè)，節(jié)省大量的人力、物力和財(cái)力，實(shí)現(xiàn)快速有效的產(chǎn)量評(píng)估。