2005-2014年南方農牧交錯帶凈初級生產力時空分布特征

趙唯茜, 杜華明, 董廷旭, 胡利利

(綿陽師范學院 資源環境工程學院, 四川 綿陽 621000)

農牧交錯帶是我國一條重要的生態安全屏障帶，也是對氣候環境變化和人類活動影響反應最為敏感的區域[1]。其日益惡化的生態環境受到了政府和理論界的重視[2]。南方農牧交錯帶是我國農牧交錯帶的重要組成部分[2]，地處青藏高原與四川盆地、云貴高原的過渡地帶，區域內不同植被生態系統在水平空間和垂直空間交錯分布[3]。既是我國重要的生物多樣性保護和水源涵養區[4]，又是長江中下游地區的重要生態安全屏障。因此區域內生態系統的改變，會對周邊區域產生較大影響，從而對區域內生態系統變化的了解顯得極為重要。凈初級生產力可直接反映出生態系統的變化情況。

凈初級生產力(Net Primary Productivity,NPP)是指植被通過光合作用所產生的有機質總量減去植被自養呼吸后的剩余部分，是生態系統中其他生物成員生存和繁衍的物質基礎[5]。凈初級生產力反映了植物群在自然條件下的生長狀況和生產能力，是一個估算地球支持能力和評價陸地生態系統可持續發展的重要生態指標[6]。近年來我國學者對不同區域、不同尺度下植被NPP的研究工作相繼展開[7-9]，但多數集中于某一行政區域內或單一自然地理單元內[10-14]。對于像農牧交錯帶這樣的地理單元過渡帶內的NPP研究較少。因此，借助遙感和地理信息技術，利用MODIS的遙感數據，對南方農牧交錯帶內整體和不同植被生態系統NPP的空間分布和時間變化進行分析，進而掌握南方農牧交錯帶內生態系統變化情況。以期望為南方農牧交錯帶內生態系統的有效保護和植被資源的有效管理與合理利用提供理論參考。

1 研究區概況

南方農牧交錯帶位于25.55°—34.31°N，98.14°—104.42°E，處于我國第一階梯和第二階梯交界處，是青藏高原與四川盆地和云貴高原的過渡地帶(圖1)。交錯帶內整體地勢高，地勢起伏較大，地勢高差較大。整體趨勢東南低，西北高，平均海拔大于3 000 m。主要地貌表現為高山地貌和高原地貌。氣候類型以高原山地溫帶、寒溫帶季風氣候為主，夏季受東南季風和西南季風影響，降水豐沛，冬季受西北季風影響氣候寒冷。區內河流水系分布較多，大體上呈“川”字形排列，多數江河自西北流向東南方向，主要江河有怒江、瀾滄江、金沙江、雅礱江、大渡河、岷江等。植被主要有云杉、冷杉、紅杉等林地，以及高山灌叢和高山草甸等。行政區劃上包括四川省的阿壩藏族羌族自治州、甘孜藏族自治州和云南省迪慶藏族自治州。交錯帶內地廣人稀，少數民族較多。由于自然、社會等客觀因素的限制影響，致使經濟發展相對緩慢。

2 數據來源以及研究方法

2.1 數據來源

文中所涉及到的遙感數據均來自美國宇航局(NASA)的MODIS產品(https:∥lpdaac.usgs.gov/)。植被NPP數據為MODIS版本6的MOD17A3H數據集，數據利用BIOME-BGC模型和光能利用率模型所建立的NPP估算模型，模擬得到全球生態系統的年NPP值。版本6的MOD17A3H與版本5.5的MOD17A3相比，使用了新的生物屬性調查表(BPLUT)和新版的全球模型與融合室(GMAO)的日氣象數據對NPP數值進行模擬，提高了NPP的估算精度，而且數據的空間分辨率提高到500 m×500 m。目前該數據已廣泛應用于不同空間尺度下的NPP和碳循環研究中。本文中選取了2005—2014年共10 a的NPP數據進行分析。土地覆蓋分類所使用的數據為MCD12Q1數據集，空間分辨率為500 m×500 m。該數據通過對Terra和Aqua兩顆衛星1 a的觀測數據進行處理后，根據國際生物圈計劃分類標準，確定了17個土地覆蓋類別。MCD12Q1數據集包含了5種不同的土地覆蓋分類方案，分別為國際生物圈計劃分類方案、美國馬里蘭大學分類方案、基于MODIS數據葉面積指數/光合有效輻射的分類方案、基于MODIS數據凈初級生產力分類方案、植物功能型方案。本文土地覆蓋類型分類方案選用國際生物圈計劃的分類方案，對研究區域內土地覆蓋分類進行提取。

圖1南方農牧交錯帶區域分布

2.2 研究方法

研究區域跨度較大，涉及MODIS三景影像。因此在對數據提取和分析之前，使用MODIS重投影工具(MRT)對HDF格式的MOD17A3H數據集和MCD12Q1數據集進行波段提取、影像拼接，將影像投影坐標重新定義為WGS_1984_Albers投影，文件格式轉換為TIF格式。應用ArcGIS對重投影后的影像進行研究區域的裁切，并對區域內數據進行提取、分析。

本文所涉及的植被NPP年際變化趨勢分析處理方法主要有簡單差值法和一元線性回歸分析法，在像元的基礎上，模擬研究區域內2005—2014年每個柵格像元植被NPP值的變化趨勢。

簡單差值法是對同一研究區域內，不同年份的影像進行相減，利用差值的大小來衡量在研究時間段內植被NPP值的變化程度和變化量。其計算公式為：

NPPminus=NPPt1-NPPt2

(1)

式中:NPPminus為兩個不同年份帶的NPP的差值；t1,t2為不同的年份。

一元線性回歸分析法是在柵格數據的每個像元基礎上，對研究區域內2005—2014年10年間每個像元NPP值的變化趨勢進行模擬，從而揭示不同時期內NPP變化趨勢的空間特征。其計算公式如下:

(2)

式中：θslope為研究區內NPP值的趨勢線的斜率即線性傾向值；n為研究時間段內總年數，本文中n=10；i為年數變量，即i=1～n；NPPi為第i年的NPP的平均值。隨時間i的增加θslope>0時，NPP處于上升趨勢；反之，隨時間i的增加θslope<0時，NPP值則處于下降趨勢。θslope的大小反映了NPP上升或下降的速率，也就是其趨勢線的傾斜程度。

10年間NPP的變化程度用百分比表示，其公式為：

(3)

式中:NPPchange為NPP值的變化百分比；θslope為由公式(2)計算得出的NPP斜率；NPPmean為10 a的平均NPP值；n為研究時間段內總年數，文中n=10。

3 結果與分析

3.1 年均植被NPP值空間分布特征

2005—2014年均植被NPP分布圖(圖2)可以看出，由于受到經緯度地帶性、植被覆蓋、氣候、地形等因素的綜合影響，2005—2014年植被NPP分布差異較大，整體呈自西北部向南、向東增多的趨勢。大約在31°N線以北地區植被NPP自西向東逐漸增加，介于240～503 gC/(m2·a)之間。大約在29°～31°N，呈西北向東南逐漸增加的趨勢，植被NPP值介于274～405 gC/(m2·a)之間。在29°N以南的地區，植被NPP自北向南大幅度增加，介于289～749 gC/(m2·a)之間。從行政區域劃分來看，云南省瀘水縣、玉龍縣、寧蒗縣植被NPP最大，分別為748.72，712.13，671.84 gC/(m2·a)，3縣均處于低緯高原季風區，受高原和高山峽谷地形的影響較大，降水量豐富，日照時間充足，擁有良好的植被生長條件。四川省德格縣、甘孜縣、理塘縣NPP均值最小，分別為240.40，243.34，245.98 gC/(m2·a)，由于3個縣處于青藏高原東部，受高原氣候影響較大，境內氣溫低，冬長夏短，降水量較少，全年無霜期短，氣候不利于植被生長。區域內河谷地帶植被NPP值較高，沿河流呈帶狀分布。位于南部橫斷山脈的怒江、瀾滄江、金沙江流域內河谷植被NPP高于位于北部大渡河、岷江河谷內植被NPP。

圖2 2005-2014年南方農牧交錯帶年均植被NPP分布

3.2 植被NPP年際變化特征

2005—2014年南方農牧交錯帶年植被NPP變化范圍為344～426 gC/(m2·a)，年均植被NPP為389 gC/(m2·a)。平均變化幅度為-0.48 gC/(m2·a)。植被NPP在2012年時最低，在2006年時植被NPP值最高。通過對2005—2014年南方農牧交錯帶內植被NPP的分級構成比例(圖3)分析發現，南方農牧交錯帶年植被NPP主要分布在200～300 gC/(m2·a)，300～400 gC/(m2·a)范圍內，分別占南方農牧交錯帶總面積的18.87%～26.21%和18.43%～26.83%。其次為100～200 gC/(m2·a)，400～500 gC/(m2·a)范圍內，分別占南方農牧交錯帶總面積的7.12%～15.66%和9.67%～13.69%。再其次為500～600 gC/(m2·a)和>700 gC/(m2·a)區間，分別占南方農牧交錯帶總面積的7.83%～10.78%和4.91%～13.01%。植被NPP在600～700 gC/(m2·a)范圍內所占土地面積最少，占南方農牧交錯帶總面積的5.7%～9.73%。

圖3 2005-2014年南方農牧交錯帶平均植被

利用一元線性回歸法分析得出南方農牧交錯帶平均植被NPP變化率(圖4)。分析發現，南方農牧交錯區10年間植被NPP減少的面積大于增加面積，植被NPP呈減少的面積占南方農牧交錯帶總面積的57.41%，呈增加趨勢的面積占總面積的42.59%。其中減少0～15%的面積占總面積的47.88%，占南方農牧交錯帶土地面積最大，且在交錯帶內廣泛分布。減少15%～30%的面積占總面積的8.56%，主要分布于南部貢山縣、德欽縣、得榮縣、中甸縣、木里縣、稻城縣的河谷地帶和西部高原地區的白玉縣、德格縣、甘孜縣、色達縣和爐霍縣；減少大于30%的面積占總面積的0.97%，主要在南部的河谷地帶分和汶川縣集中分布。增加0～15%的面積占總面積的34.41%，除西部高原的德格縣、色達縣、甘孜縣、爐霍縣、白玉縣分布相對較少，在其他縣域廣泛分布；增加15%～30%的面積占南方農牧交錯帶總面積的6.98%，分布較多的區域主要在東部，主要包括九寨溝縣、松潘縣、黑水縣、理縣、小金縣、金川縣、丹巴縣、道孚縣、雅江縣、康定縣、九龍縣、鹽源縣、寧蒗縣、福貢縣、瀘水縣；增加大于30%的面積占總面積的1.21%，在東部分布較多，主要包括松潘縣、茂縣、理縣、汶川縣、小金縣、丹巴縣、康定縣、九龍縣。

圖4 2005-2014年南方農牧交錯帶植被NPP變化百分率

將本文研究時間段端點時間，即2005年和2014年的年植被NPP影像進行簡單差值計算，得出10 a的植被NPP變化特征分布圖(圖5)，2014年與2005年相比平均植被NPP增加了106 gC/(m2·a)。通過對圖5分析發現，2014年與2005年相比，南方農牧交錯帶內植物NPP呈減少的面積占南方農牧交錯帶總面積的42.05%，植物NPP呈增長的面積占南方農牧交錯帶總面積的57.95%。其中，增幅大于80 gC/(m2·a)的地區占總面積的8.38%，主要分布在東北部的九寨溝縣、松潘縣、黑水縣、茂縣、汶川縣、馬爾康縣、丹巴縣，中部的雅江縣以及南部的鹽源縣、寧蒗縣、貢山縣、維西縣、福貢縣、蘭坪縣、瀘水縣。增幅在40～80 gC/(m2·a)的地區占總面積的15.02%，在九寨溝縣、松潘縣、黑水縣、茂縣、理縣、汶川縣、馬爾康縣、小金縣、金川縣、丹巴縣、道孚縣、雅江縣分布較多；增加在0～40 gC/(m2·a)的地區占總面積的34.56%，廣泛分布于除福貢縣和瀘水縣以外其他區域。減少0～40 gC/(m2·a)的地區占總面積的29.18%，在除九寨溝縣、茂縣、汶川縣、馬爾康縣、金川縣、丹巴縣、雅江縣、維西縣、福貢縣、瀘水縣之外的地區分布較多。減少40～80 gC/(m2·a)的地區占總面積的9.67%，在色達縣、爐霍縣、甘孜縣、德格縣、白玉縣、巴塘縣、鄉城縣、得榮縣、德欽縣、中甸縣、木里縣、鹽源縣、寧蒗縣、麗江縣分布廣泛。減少大于80 gC/(m2·a)的地區占總面積的3.2%，主要分布在汶川縣、色達縣、康定縣、德欽縣、中甸縣、木里縣、麗江縣、寧蒗縣、鹽源縣。

圖5 南方農牧交錯帶2014與2005年植被NPP差值

簡單差值法雖然能夠直接反映10年間南方農牧交錯帶內植被NPP的變化趨勢，由于計算所涉及的數據是研究時間段端點數據之間的差異，因此得到的結果可能會受到極端天氣影響[10]。因此本文選用一元線性回歸和簡單差值法對南方農牧交錯帶內植被NPP的變化趨勢進行分析。

3.3 南方農牧交錯帶不同生態系統類型植被NPP化趨勢

南方農牧交錯帶生態系統分為森林生態系統、草原生態系統、灌叢生態系統、濕地生態系統、農田生態系統。其中森林生態系統主要分布在南部橫斷山脈和河谷谷底。草地生態系統主要分布在西部高原地區，灌叢生態系統、農田生態系統等在區域內零星分布(圖6)。

圖6 南方農牧交錯帶土地利用類型

南方農牧交錯帶不同生態系統類型NPP變化趨勢和變化幅度差異明顯(表1)。森林生態系統年均植被NPP為558 gC/(m2·a),其中2011年森林生態系統植被NPP最高為608 gC/(m2·a)，在2012年植被NPP最低為492 gC/(m2·a)。草地生態系統年均植被NPP為275 gC/(m2·a)，其中2006年草地生態系統植被NPP最高，為305 gC/(m2·a)，2012年草地生態系統植被NPP最低為244 gC/(m2·a)。農田生態系統年均植被NPP為500 gC/(m2·a),其中2006年農田生態系統植被NPP最高為536 gC/(m2·a)，在2012年植被NPP最低為459 gC/(m2·a)。灌叢生態系統年均植被NPP為359 gC/(m2·a)，其中2009年灌叢生態系統植被NPP最高，為395 gC/(m2·a)，在2012年、2014年灌叢生態系統植被NPP最低為325 gC/(m2·a)。濕地年均植被NPP為326 gC/(m2·a)，其中2006年濕地生態系統植被NPP最高為357 gC/(m2·a)，在2010年濕地植被NPP最低為298 gC/(m2·a)。不同生態系統植被NPP均值表現為：森林生態系統>農田生態系統>灌叢生態系統>濕地生態系統>草地生態系統。

從圖7中可以看出，南方農牧交錯帶內不同生態系統植被NPP變化相似，呈波動減少變化。其中灌叢植被的NPP下降趨勢最大，平均變化幅度為-3.34 gC/(m2·a)。森林植被NPP下降趨勢為-0.86 gC/(m2·a)，農田植被NPP下降趨勢為-0.52 gC/(m2·a)，濕地植被NPP下降趨勢為-0.32 gC/(m2·a)，草地植被NPP下降趨勢最小，為-0.23 gC/(m2·a)。該地區不同植被類型NPP在2009年之后NPP波動變化幅度較大，可能與2009—2012年開始的川滇旱災有關。

表1 2005-2014年南方農牧交錯帶不同生態系統植被NPP特征 gC/(m2·a)

圖7 2004-2015年南方農牧交錯帶內

4 結 論

(1) 南方農牧交錯帶NPP總體上呈南部高北部底，東北部高，西部低的分布趨勢。河谷地帶分布較高。南部橫斷山脈是南方農牧交錯帶NPP最高的地區，植被NPP大于600 gC/(m2·a)。東北部有部分區域年均植被NPP也超過600 gC/(m2·a)。西北部高原地區年均植被NPP小于300 gC/(m2·a)，是該區域年均NPP最低的區域。而西北部是區域的高海拔分布區，東北部和南部是低海拔區域。

(2) 2005—2014年，南方農牧交錯帶植被NPP年變化范圍為344～426 gC/(m2·a)，年平均植被NPP為470 gC/(m2·a)。南方農牧交錯帶年NPP在200～300 gC/(m2·a)區間分布最多，占南方農牧交錯帶總面積的26%。其次為300～400 gC/(m2·a)占23%。

(3) 2014年與2005年相比57.95%的區域為增長，平均植被NPP增加了106 gC/(m2·a)。但是，2005—2014年，南方農牧交錯帶57%的區域植被NPP趨勢呈略微減少。在10年間，48%的植被覆蓋區域植被NPP減少了0～15%，為植被NPP變化分布最多的區間。34%植被覆蓋區域植被NPP增加了0～15%。

(4) 南方農牧交錯帶不同生態系統平均植被NPP差異明顯，森林生態系統植被NPP最高，變化范圍為492～608 gC/(m2·a)；其次是農田生態系統，變化范圍為459～536 gC/(m2·a)；灌叢生態系統變化范圍為325～395 gC/(m2·a)；濕地生態系統變化范圍為298～357 gC/(m2·a)；草地生態系統植被NPP最低，變化范圍為275～305 gC/(m2·a)。