基于遙感的新疆人工綠洲擴張中植被凈初級生產力動態變化

張 芳，熊黑鋼，馮 娟，許仲林

?

基于遙感的新疆人工綠洲擴張中植被凈初級生產力動態變化

張 芳1,2，熊黑鋼2,3，馮 娟1,2，許仲林1,2

（1. 新疆大學資源與環境科學學院，烏魯木齊830046；2. 教育部新疆綠洲生態重點實驗室，烏魯木齊830046； 3. 北京聯合大學應用文理學院，北京 100083）

人類活動對全球碳循環過程的影響是全球氣候變化研究中爭論的主要焦點問題之一。不同區域、不同自然和人為因素對碳循環的作用機制十分復雜且差異很大，通過探討人類開墾活動對干旱區綠洲生態系統凈初級生產力（net primary productivity, NPP）分布格局的影響，將有助于增加人類活動對干旱區陸地生態系統及碳循環影響的認識和理解。該文基于長時間序列的Landsat遙感數據，利用CASA（Carnegie-Ames-Stanford Approach）模型估算了1976年6月－2016年6月新疆奇臺新人工綠洲擴張過程中NPP的變化情況，結果表明：1）1976年研究區主要是荒漠草原景觀，農田僅有5 150 hm2，占研究區面積的3.25%，到2016年，農田面積擴大到64 675 hm2，增加11倍以上；2）人工綠洲擴張的同時，泉水溢出帶消失，水庫干涸，灌叢面積從17 871 hm2下降到5 202 hm2；3）土壤鹽堿化發育進程中斷，裸堿地面積減少近70%，并且開始向新生草地演替；4）2016年6月，研究區NPP總量達到了90 602 t，是1976年6月的2.2倍，農田對NPP增量的貢獻超過93%，同期，灌叢則凈損失了8 265 t NPP，沼澤凈損失3 319 t；5）1976年的NPP主要由天然草地和灌叢構成，分別占56.68%和26.61%，2016年的NPP主要由農田和新生草地構成，分別占65.23%和31.85%。

遙感；草地；農田；新人工綠洲擴張；NPP；CASA模型；新疆奇臺綠洲

0 引 言

植被凈初級生產力（net primary productivity，NPP）是指綠色植物在單位面積和單位時間內通過光合作用累積的有機物量，用于測定植物體固定的凈碳量，可反映陸地生態系統與大氣間的碳交換過程[1]。近年來，在“糧食安全”和“全球變暖”2個重大問題驅動下，陸地生態系統NPP的相關研究與應用受到廣泛重視[2-4]，建立于1987年的國際地圈-生物圈計劃（IGBP）、全球變化與陸地生態系統（GCTE）和京都協定（KP）均把植被的NPP研究確定為核心內容之一[5]。NPP測量方式包括直接測定法和間接測定法，直接測定方式中以渦度相關法為公認標準，可以直接測定大氣與植被水汽、CO2交換通量，測量精度高，適用于尺度小、下墊面均一平坦環境[6]；間接測量包括模型估測和遙感反演，其中基于光能利用率的遙感模型目前已成為NPP估算的一種有效手段，使區域及全球尺度的NPP估算成為可能[7-9]。CASA（Carnegie- Ames-Stanford-Approach）模型是目前區域NPP估算中應用較多的光能利用率模型，其優點在于：模型相對簡單、所需參數較少；遙感數據有獲得植被及光合有效輻射信息的優勢，可對植被覆蓋情況分類，通過多時相遙感數據可以實現對NPP進行時間和空間上的動態監測[10]。

新疆地處“一帶一路”核心區，在中國經濟發展和國家安全方面占有十分重要的戰略地位，是中國21世紀重要的資源接替區和經濟增長支點，也是后備耕地資源占全國耕地儲備量1/4的國家糧食安全的重要保障基地[11]。目前，只占新疆土地面積的4%的綠洲，聚集著95%的人口和90%以上的社會財富，灌溉農業則是綠洲地區最廣泛的經濟活動[12-13]。中國從1976年到2016年是經濟發展不斷加速的40 a，如今，整個天山北麓經濟帶僅占全疆面積8.22%，但國內生產總值占全疆近50%，耕地面積和人口均超過全疆的25%[14]。新人工綠洲的大規模擴張必然造成陸地生態系統結構和組成的改變，進而影響全球碳循環過程[15-16]。

本研究選用天山北麓奇臺綠洲平原區1976－2016年的5期同時相的Landsat系列遙感影像（間隔10 a左右），結合地面調查數據，氣象數據等，通過基于光能利用率的CASA模型計算該區域40 a來新綠洲迅速擴張過程中NPP分布格局的變化，以期增加人類活動對干旱區陸地生態系統及碳循環影響的認識，為人類活動在全球氣候變化中的作用提供客觀依據。

1 數據與方法

1.1 研究區概況

依據前人研究結果[17-18]，天山北麓的綠洲分布擴張軌跡表現為：古綠洲（千年前形成）多分布在河流三角洲上；老綠洲（20世紀中葉前形成）多分布于沖洪積扇引水方便的中、上游沖積平原；新綠洲（20世紀中葉后形成）則主要分布于沖積平原的中下部。奇臺綠洲位于天山北麓經濟帶上，農業開發強度很高，是多年的全國農業百強縣，因此，研究區選用奇臺綠洲與古爾班通古特沙漠相毗鄰的交錯帶區域，地理坐標范圍為43°57′N～44°13￠N，89°24′E～90°09′E，年均溫5℃左右，年均降水約176 mm，年均蒸發勢2 141 mm，屬于溫帶大陸性干旱半干旱氣候。

1.2 遙感數據選擇及預處理

Landsat系列數據是目前世界范圍內應用最廣泛的民用對地觀測多光譜遙感數據，其空間分辨率及長時間數據序列非常適合中等區域尺度地表環境演變與監測研究。遙感數據的選擇主要從以下2個方面進行考慮：1）數據質量（如云量、條帶等）；2）時相（植物生長最旺盛季節）。受地理位置和大地形影響，新疆自然植被生長季極不同步，研究區3－5月作物歸一化植被指數（normalized difference vegetation index，NDVI）下降，7－8月自然荒漠植被的NDVI在高溫和不連續降水條件下明顯下降[19]，這與多年對研究區實地調查情況一致：奇臺綠洲平原區主要作物為小麥、玉米和甜菜，玉米、甜菜一般在4月份播種，5月中旬之前其幼苗葉片尚不能蓋滿田塊，7月中/下旬小麥開始收割。本研究旨在從新人工綠洲擴張視角出發，討論人類活動對干旱區植被結構和NPP分布格局的影響，因此，6月為該區域提取植被覆蓋類型變化信息的最佳時間。

通過對上述因素的綜合考慮，選用了1976年6月6日Landsat MSS影像和1989年5月23日、1996年6月28日、2007年6月27日、2016年6月19日的Landsat TM影像共5期衛星數據，間隔約10 a左右。

遙感數據在使用前進行了嚴格的預處理：1）大氣輻射校正；2）幾何精校正，校正誤差控制在0.2個像元以內。預處理通過遙感數據處理專業軟件ENVI5.1完成。預處理后，原灰度值圖像轉換為地表反射率圖像，可用于下一步的圖像分類、運算和模型應用。

1.3 NPP計算

CASA模型中植被NPP主要由植被吸收的光合有效輻射（absorbed photosynthetically active radiation, APAR）與光能轉化率（）2個變量來確定[20-21]。公式如下

NPP=APAR·T1·T2·W·ε（1）

APAR=0.5 SOL·FAPAR （2）

式中SOL是當月研究區太陽總輻射量，MJ/m2；ε為最大光能轉化率，g/MJ；FAPAR（fractional APAR）為植被對光合有效輻射的吸收比例，常數0.5表示植被所能利用的太陽有效輻射（可見光波長0.4～0.7m）占太陽總輻射的比例。T1和T2均為溫度脅迫系數；W為水分脅迫影響系數。

FAPAR由NDVI確定。相關計算基于下列公式完成

NDVI=(NIR–)/(NIR+) （3）

SR=(1+NDVI)/(1–NDVI) （4）

FAPARSR=[(SR–SRmin)/(SRmax–SRmin)]·(FAPARmax–

FAPARmin)+FAPARmin（5）

FAPARNDVI=[(NDVI–NDVImin)/(NDVImax–NDVImin)]·

(FAPARmax–FAPARmin)+FAPARmin（6）

FAPAR=(FAPARSR+FAPARNDVI)/2 （7）

式中NIR為近紅外反射率，為紅波段反射率；SR為簡單比值植被指數，與植被類型有關，SRmin=1.08，SRmax取值在4.14～6.17之間，FAPARmax=0.95，FAPARmin=0.001，這些參數無區域限制性，系通用參數[21]。

光能轉化率是植物通過光合作用固定太陽能，將所截獲/吸收的能量轉化為碳（C）/有機物干物質的效率，用g/MJ（以C計，下同）表示。在此，CASA模型ε取值0.389 g/MJ[22-23]。

T1反映了在低溫和高溫時植物內在的生化作用對光合的限制，計算如下

T1=0.8+0.02T–0.00052（8）

式中T為某區域一年內NDVI值達到最高時月份的平均氣溫，℃。

T2表示環境溫度從最適宜溫度T向高溫和低溫變化時植物光能轉化率變小的趨勢，計算如下：

T2=1.1814({1+e[0.2 (T–10–T)]}/{1+e[0.3 (–T–10+T)]}) （9）

式中為影像數據獲取月份的平均氣溫，℃。

W反映植物所能利用的有效水分條件對光能轉化率的影響。隨著環境有效水分的增加，W逐漸增大，取值為0.5（在極端干旱條件下）～1（非常濕潤條件下）[20]。

W=0.5+0.5(EET/PET) （10）

式中PET為可能蒸散量，mm；EET為實際蒸散量，mm。

PET=58.93BT （11）

式中BT為生物溫度，℃。BT是0～30℃月均溫之和除以12[24]。

EET={·R·[P+R2+·R]}/{(+R)·(2+R2)} （12）

式中為當月降水量，mm；R為當月太陽凈輻射，MJ/m2。

氣象數據來自中國氣象科學數據共享服務網[25]。

經過多年田間走訪調查，研究區為灌溉農業，灌溉過程由生產者根據經驗和土壤墑情決定。傳統長期采用大水漫灌方式，造成水資源大量浪費，2008年后開始逐步推廣節水工程，科學測算各類作物需水量，保障作物充足水分供給。因此，研究區農業生產基本不受降水狀況限制和脅迫，因而農田與沼澤的水分脅迫W取值為1。

1.4 影像分類

目前，研究區主要的地物類型包括農田、草地、裸堿地和灌叢（圖1），樹木主要在部分田邊或路邊以防護林形式成線狀分布，未形成寬度超過1個像元（30 m）的成面積的可提取信息，因此林地未進行單獨分類。農田種植面積最大的作物類型主要有小麥、玉米、甜菜、打瓜等；草地主要指低矮草本植物為主如叉毛蓬、駱駝蓬等；灌叢指有株形較大的灌木或半灌木與草本的混合群落，如檉柳+芨芨草+苦豆子等。研究區土壤為堿性土壤，土壤pH值大于9.5的地方植物難以生長，有些就是光板地，因此劃分為裸堿地[26]。

圖1 研究區典型地物景觀照片

2 結果與分析

2.1 研究區40 a新人工綠洲擴張

研究區地勢東南高，西北低，西北處角地勢低洼處建有一平原水庫，泉水溢出帶沿800 m等高線由西南向東北方向呈弧形分布[27]。根據研究區多年實地定位觀測及訪查資料，采用監督分類法對各期影像進行分類（圖2），5期影像分類精度分別為96.38%、95.65%、92.75%、98.60%、98.45%。從圖2可以看出，1976年時，以泉水溢出帶為水源區形成沼澤、坑塘和向北流淌的泉水河，經實地訪查，沼澤及其周圍生長植物以蘆葦和芨芨草混合灌叢為主，在泉水河兩側，灌叢沿河線狀分布；干旱氣候與地表水分條件的共同作用使鹽堿土廣泛發育；農田稀疏，呈小斑塊狀零星散布，人類的活動強度較弱。1989年，農田在原有田塊周圍迅速增加，沼澤特征減弱，水體水面略有縮小，灌叢大量減少。1996年，農田繼續擴張，分散田塊連接成片，泉水溢出帶消失被開墾為農田，水庫干涸僅留下大壩。2007年，地表覆蓋主要為3種地物類型：農田、草地和裸堿地，農田繼續在原有基礎上向外圍拓展，向沙漠逼近，許多鹽堿地也墾為農田。到2016年，研究區原來的天然綠洲景觀完全被新人工綠洲替代，大面積條田相互連接，成為適合大型農業機械作業的現代化農業綠洲。

圖2 研究區不同時期土地覆蓋圖

表1為根據分類結果提取的各類別覆蓋面積。1976年研究區農田僅有5 150 hm2，占研究區面積3.25%，1989年，農田面積增加到20 376 hm2，年均增加1 171 hm2；之后的1989－1996、1996－2007、2007－2016的3個時間段里，年均增加量分別為2 016、1 355和1 698 hm2。20世紀80年代中期之前，受交通、農業機械化程度、灌溉條件等限制，墾荒速度相對較低；20世紀80年代末到20世紀90年代中期墾荒強度最大，這一階段，隨著聯產承包責任制落實到戶，全縣打井總眼數超過了4 000眼[28]；之后在國家退耕還草等政策調整下，墾荒速度有所放緩，近10 a來墾荒強度又有所增加，2016年耕地面積占比已經超過了研究區面積的40%。

表1 地物類型的面積變化

1976年－2016年裸堿地面積減少約41 056 hm2，縮小近70%，這主要是由于開墾導致研究區水文條件發生改變引起的：研究區地處干旱區沖洪積扇下緣，1976年地表水有水庫、坑塘、泉水河、沼澤等（圖2），說明當時地下水位高，極易造成土壤強烈鹽堿化；由于大量開墾農田，打井灌溉，到1996年地下水位迅速下降超過10 m[29]，沼澤消失，水庫干涸，地表水體僅剩幾個小澇壩。地下水位下降超過一定深度，土壤表層積鹽進程停止，活性鹽堿土變為殘余鹽堿土，在研究區近180 mm的年降水量的淋洗作用下，土壤表層的鹽堿化的程度逐漸減弱，有植被開始生長，鹽堿地向新生草地演替或被開墾改良成農田。灌叢原來多生長在水分條件較好的地方，40 a來，灌叢從17 871 hm2減少到5 202 hm2，已經很少能成片分布。草地面積在40 a里變化不大，但其分布位置明顯變化，原有草地被墾為農田，新生草地在鹽堿地上形成，可知，人類的新綠洲擴張活動使農田大量增加的同時也改變了整個區域的水資源消耗方式和時空分布，進而改變不同植被類型的空間分布格局。

2.2 人工綠洲擴張過程中NPP動態變化

依據分類結果計算各類別中每個像元的NPP值，其中裸堿地和水體由于基本無植物信息不參與統計（表2）。2016年，研究區NPP總量為90 602 t，是1976年同期的2.2倍，其中，僅農田擴張就增加NPP 55 561 t，占增量的93%以上。2016年與1976年相比，灌叢則凈損失的NPP總量為8 265 t，沼澤凈損失3 319 t。1976年農田對NPP的貢獻率僅為8.62%，40 a后，其貢獻率達65.23%。1996年以前，草地面積占比在42.72%～57.18%之間波動，但其NPP總量維持在一個相對平穩水平，近20 a則有明顯增加。1989年與1976年相比，草地總NPP略降，這可能是由于20世紀80年代中后期質量及時相對應最好的Landsat數據為1989年5月23日，與其他4期影像存在約半個月的時間差造成的。

表2 1976－2016 年的植被凈初級生產力（NPP）變化

若不考慮在時間上受到一定影響的1989年數據，1976年6月農田單位面積NPP為68.6 g /m2。隨著農業機械化耕作、化肥農藥的施用、鹽堿化土壤改良等耕作模式與管理手段的不斷提高，農田單位面積NPP逐步提高，2007年6月和2016年6月均已超過90g/m2。6月份的農田NPP一定程度上受到作物類型的影響，例如研究區主要農作物為小麥、玉米，2000年后受市場需求驅動，也開始一定規模的種植甜菜、油葵、打瓜等，這使得農田中不同作物的生長旺盛期不同步。40 a里，6月份的草地單位面積NPP在30～40g/m2范圍之內波動，相對比較穩定，其影響因素主要是降水和春季氣溫回升情況，特別是春季融雪給土壤提供的水量在4－5月可以明顯地影響短生植物的萌發和生長，例如1989年草地面積明顯增大（表1），但單位面積NPP卻較小，就與一定量的短生植物尚處于幼苗期有關。灌叢多為多年生深根植物，其NPP除在一定程度上受到氣象因素影響外，也受到人類較大的干擾，例如1976年地面水量大，灌叢多生長于沼澤和泉水河周圍，形成比較長期穩定的群落結構，隨著新綠洲墾荒強度的增加，灌叢被大量開墾，地下水位下降，新的灌叢多出現于田間空地，株齡短，冠幅小，密度低，NPP也較低。沼澤只在1976年和1989年2期圖像上出現，到1996年時已完全消失，但是1989年沼澤單位面積NPP卻比1976年大42%，這可能是因為1989年沼澤已經開始明顯退化，淹水深度降低，更利于淺水植物生長。

2.3 人工綠洲擴張中地表NPP格局的變化

奇臺綠洲40 a的新人工綠洲擴張活動改變了原來天然綠洲的植被構成和NPP空間分布格局（圖3）。1976年6月單位面積NPP>120 g/m2的區域主要零星分布著農田和高大檉柳+芨芨草等灌叢，介于70～120 g/m2之間的的主要是沼澤和芨芨草+苦豆子+鹽爪爪等混生灌叢群落，≤70 g/m2的多為荒漠草原，1976年6月研究區的NPP主要由天然草地（56.68%）和天然灌叢（26.61%）構成，農田的NPP貢獻率僅占8.62%，沼澤為8.10%（表2），大面積的鹽堿地幾乎無植被生長。經過40a的新人工綠洲大開發過程，2016年6月研究區的NPP主要由農田（65.23%）和新生草地（31.85%）構成，灌叢的NPP貢獻率下降到2.92%，天然綠洲的植被構成和荒漠草原景觀已被完全改變。

圖3 研究區不同時期單位面積植被凈初級生產力（NPP）分布

定量分析40a來研究區NPP分布格局的變化（表3），1976年灌叢中的66.93%、沼澤中的75.58%、草地中的41.43%、鹽堿地中的23.86%和水體中的77.96%構成了現在的農田；1976年的草地中有43.62%在40a后依然是草地，只是基本都被新墾農田分割成碎片散布，然后與由59.38%的裸堿地演替而來的新生草地連接成片，另外1976年灌叢中的21.44%和沼澤中的16.57%退化成為2016年草地中的一部分。1976年的裸堿地中僅剩15.01%在2016年依然存在，其余的裸堿地則來自草地、灌叢、沼澤等的退化。

表3 2016年相較于1976年的地物類別轉移

3 討 論

氣候變化對研究區NPP變化的影響。在無人為干擾條件下，NPP的變化與溫度、水和太陽輻射等因素的變化密切相關。有研究表明，近年來的全球氣候變化緩解了對植物生長的限制，使得北半球中高緯度地區的NPP有所增加[30]。研究區40 a來的氣溫也呈暖化趨勢，但是研究區植物生長季節的溫度變化表現為春季上升0.14 ℃/10 a、夏季氣溫呈反向變化，為-0.14 ℃/10a[31]。通過對研究區40 a來NPP變化研究來看，研究區5個研究年份的生物溫度基本在9.52 ℃上下波動，變化幅度不超過0.25 ℃，可知，溫度對研究區NPP的影響很小。研究區年均降水量約180 mm，近40a的年降水量升率為15.8mm/10a[31]。由于研究區地處干旱區，蒸發強度大，地下水位已普遍降于10 m以下，這樣的降水升率在一定程度上雖然可能有助于草地和短生淺根植物的生物量增加，但受降水增加的影響，研究區近年來的日照時數呈下降趨勢，這對光合作用又會起到一定的削弱，因此，對本研究區而言，自然氣候因子對NPP的影響是比較有限的，這也與前人研究一致，有研究認為天山以北潛在植被NPP受氣候變化影響的脆弱性較低[32]。

人類活動是影響研究區NPP分布格局的決定性因素。1976－2016年是中國經濟發展不斷加速的階段，也是人口持續增加，人民生活水平不斷提高的階段，墾荒活動在干旱區綠洲普遍展開，特別是老綠洲外圍的荒漠交錯帶區域，大面積的天然綠洲被開墾，農田延伸直至沙漠，新人工綠洲的劇烈擴張，強烈改變了水資源的消耗方式和時空分配，引起生態環境顯著變化[17]。由于新墾農田快速增加，奇臺縣農業灌溉用水超過全縣用水量的95%，年超采地下水近1.5億m3，形成的地下水漏斗近800 km2 [33]，地下水超采導致泉水溢出帶消失，水庫干涸，給研究區的生態環境帶來負面效應，因此，水資源的可持續性利用與新人工綠洲擴張過程中生態系統碳庫反饋效應的未來發展趨勢應引起重視和進行更深入的研究。

4 結 論

40 a里，研究區農田面積從5 150 hm2增加到64 675 hm2，增加了11倍以上，裸堿地從58 859 hm2減少到17 803 hm2，減少約70%；地表覆被類型發生了劇烈變化，泉水溢出帶消失，水庫干涸，灌叢大幅減少，土壤鹽堿化發育進程中斷，新生草地在裸堿地上形成，1976年的草地中有41.43%在40 a后轉變成農田，而裸堿地則有23.86%轉變成農田，另外有59.38%轉變成草地。

研究區人工綠洲擴張使地表NPP總量大幅增加，2016年6月地表NPP總量是1976年6月的2.2倍，增量中超過93%的NPP來自于農作物。NPP分布格局也發生了顯著變化，1976年6月的NPP由天然草地、灌叢、農田和沼澤構成，分別占56.68%、26.61%、8.62%和8.10%，到了2016年6月，NPP主要由農田（65.23%）和新生草地（31.85%）構成。

[1] 陳曉玲，曾永年. 亞熱帶山地丘陵區植被NPP時空變化及其與氣候因子的關系：以湖南省為例[J]. 地理學報，2016，71(1)：35－48. Chen Xiaoling, Zeng Yongnian. Spatial and temporal variability of the net primary production (NPP) and its relationship with climate factorsin subtropical mountainous and hilly regions of China:A case study in Hunan province[J]. Acta Geographica Sinica, 2016, 71(1): 35－48. (in Chinese with English abstract)

[2] IPCC. Second assessment report of the Intergovernmental Panel on Climate Change (IPCC)[J]. Environ. Policy Law, 1996, 26(5): 234－235.

[3] 盧玲，李新，Frank Veroustraete. 中國西部地區植被凈初級生產力的時空格局[J]. 生態學報，2005，25(5)：1026－1032. Lu Ling, Li Xin, Frank Veroustraete.Terrestrial net primary productivity and its spatial－temporal variability in western China[J]. Acta Ecologica Sinica, 2005, 25(5): 1026－1032. (in Chinese with English abstract)

[4] Alberti M. The effects of urban patterns on ecosystem function[J]. International Regional Science Review, 2005, 28(2): 168－192.

[5] 楊延征，馬元丹，江洪，等. 基于IBIS模型的1960－2006年中國陸地生態系統碳收支格局研究[J]. 生態學報，2016，36(13)：3911－3922. Yang Yanzheng, Ma Yuandan, Jiang Hong, et al. Evaluating the carbon budget pattern of Chinese terrestrial ecosystem from 1960 to 2006 using Integrated Biosphere Simulator[J]. Acta Ecologica Sinica, 2016, 36(13): 3911－3922. (in Chinese with English abstract)

[6] 李守娟，馬杰，唐立松，等. 梭梭群落碳交換的尺度轉換研究: 從葉片到群落[J].干旱區研究，2016，33(2)：362－370. Li Shoujuan, Ma Jie, Tang Lisong, et al. Scale transformation of carbon exchange over haloxylon ammodendron community: From leaf to community[J]. Arid Zone Research, 2016, 33(2): 362－370. (in Chinese with English abstract)

[7] 姜春，吳志峰，程炯，等. 廣東省土地覆蓋變化對植被凈初級生產力的影響分析[J]. 自然資源學報，2016，31(6)：961－972. Jiang Chun, Wu Zhifeng, Cheng Jiong, et al. Analyzing the effects of land cover change on vegetation net primary productivity in Guangdong Province[J]. Journal of Natural Resources, 2016, 31(6): 961－972. (in Chinese with English abstract)

[8] Fang Jingyun, Piao Shilong, Field C B, et al. Increasing net primary production in China from 1982 to 1999[J]. Frontiers in Ecology and the Environment. 2003, 1(6): 293－297.

[9] 孫睿，朱啟疆. 中國陸地植被凈第一性生產力及季節變化研究[J]. 地理學報，2000，55(1)：36－45. Sun Rui, Zhu Qijiang. Distribution and seasonal change of net primary productivity in China from April 1992 to March 1993[J]. Acta Geographica Sinica, 2000, 55(1): 36－45. (in Chinese with English abstract)

[10] 冉慧. 基于CASA模型的吉林省區域NPP遙感研究[D]. 吉林：吉林大學，2010：9－14. Ran Hui. Remote Sensing Research on NPP of Regions in Jilin Province Based on CASA Model[D]. Jilin: Jilin University, 2010：9－14. (in Chinese with English abstract)

[11] 劉新平. 新疆綠洲土地資源可持續利用的經濟學分析[M]. 北京：中國文史出版社，2005.

[12] 韓德林. 新疆人工綠洲[M]. 北京：中國環境科學出版社，2001.

[13] 中國科學院新疆綜合考察隊. 新疆土壤地理[M].北京：科學出版社，1965.

[14] 程維明，周成虎，李建新. 天山北麓經濟發展與綠洲擴張[J]. 地理學報，2002，57(5)：561－568. Cheng Weiming, Zhou Chenghu, Li Jianxin. Economic development and oasis growth at the Northern Foot of the Tianshan Mountains[J]. Acta Geographica Sinica, 2002, 57(5): 561－568. (in Chinese with English abstract)

[15] Defries R S, Field C B, Fung I, et al. Combining satellite data and biogeochemical models to estimate global effects of human－induced land cover change on carbon emissions and primary productivity[J]. Global Biogeochemical Cycles. 1999, 13(3): 803－815.

[16] 何玲，賈啟建，李超，等. 基于生態系統服務價值和生態安全格局的土地利用格局模擬[J]. 農業工程學報，2016，32(3)：275－284. He Ling, Jia Qijian, Li Chao, et al. Land use pattern simulation based on ecosystem service value and ecological security pattern[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(3): 275－284. (in Chinese with English abstract)

[17] 樊自立，馬英杰，艾力西爾·庫爾班，等. 試論中國荒漠區人工綠洲生態系統的形成演變和可持續發展[J]. 中國沙漠，2004，24(1)：10－16. Fan Zili, Ma Yingjie, Alishir·Kurban, et al. On formation, evolution and sustainable development of artificial oasis ecosystem in desert area in China[J]. Journal of Desert Research, 2004, 24(1): 10－16. (in Chinese with English abstract)

[18] 韓德麟. 關于綠洲若干問題的認識[J]. 干旱區資源與環境，1995，9(3)：13－31. Han Delin. Knowledge of a few issues on oasis[J]. Journal of Arid Land Resources and Environment, 1995,9(3): 13－31. (in Chinese with English abstract)

[19] 杜加強，趙晨曦，賈爾恒·阿哈提，等. 近30 a新疆月NDVI 動態變化及其驅動因子分析. 農業工程學報，2016，32(5)：172－181. Du Jiaqiang, Zhao Chenxi, Jiaerheng·Ahati, et al. Analysis on spatio－temporal trends and drivers in monthly NDVI during recent decades in Xinjiang, China based two datasets[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(5): 172－181. (in Chinese with English abstract)

[20] 張峰，周廣勝，王玉輝. 基于CASA模型的內蒙古典型草原植被凈初級生產力動態模擬[J]. 植物生態學報，2008，32(4)：786－797. Zhang Feng, Zhou Guangsheng, Wang Yuhui. Dynamics simulation of net primary productivity by a satellite data－driven CASA model in Inner Mongolian typical steppe, China[J]. Journal of Plant Ecology, 2008, 32(4): 786－797. (in Chinese with English abstract)

[21] 樸世龍，方精云，郭慶華. 利用CASA模型估算我國植被凈第一性生產力[J]. 植物生態學報，2001，25(5)：603－608. Piao Shilong, Fang Jingyun, Guo Qinghua. Application of CASA model to the estimation of Chinese terrestrial net primary productivity[J]. Acta Phytoecologica Sinica, 2001, 25(5): 603－608. (in Chinese with English abstract)

[22] Field C B, Behrenfeld M J, Randerson J T, et al. Primary production of the biosphere: integrating terrestrial andoceanic components[J]. Science, 1998,281(5374):237－240.

[23] Potter C S, Randerson J T, Field C B, et al. Terrestrial ecosystem production：A process model based on global satellite and surface data[J]. Global Biogeochemical Cycles, 1993, 7(4): 811－841.

[24] 任正超. 基于CASA模型的俄羅斯布里亞特共和國植被NPP變化及其對氣候的響應[D]. 蘭州：甘肅農業大學, 2010：18－21. Ren Zhengchao. Vegetation Net Primary Productivity Variation and Its Response to Climate in Buryatiya Republic Russia Based on CASA Model[D]. Lanzhou: Gansu Agricultural University, 2010：18－21. (in Chinese with English abstract)

[25] 中國氣象科學數據共享服務網.中國地面累年值月值數據集[EB/OL].(2016-09-01)[2016-11-01].http://data.cma.cn/dataService/cdcindex/datacode/SURF_CLI_CHN_MUL_MON/show_value/normal.html.

[26] 張芳，熊黑鋼，安放舟，等. 基于鹽(堿)生植被蓋度的土壤堿化分級[J]. 土壤學報，2012，49(4)：665－672.Zhang Fang, Xiong Heigang, An Fangzhou, et al. Classification of soil alkalization based on halophyte coverage[J]. Acta Pedologica Sinica, 2012, 49(4): 665－672. (in Chinese with English abstract)

[27] 張芳，熊黑鋼，田源，等. 區域尺度地形因素對奇臺綠洲土壤鹽漬化空間分布的影響[J]. 環境科學研究，2011，24(7)：731－739. Zhang Fang, Xiong Heigang, Tian Yuan, et al. Impacts of regional topographic factors on spatial distribution of soil salinization in Qitai Oasis[J]. Research of Environmental Sciences. 2011, 24(7): 731－739. (in Chinese with English abstract)

[28] 李豐德. 論奇臺縣開發利用地下水資源存在的問題及建議[J]. 地下水，2003，25(1)：31－33. Li Fengde. On existed problems and suggestions on development&utilization of underground water resources of Qitai County[J]. Groundwater, 2003, 25(1): 31－33. (in Chinese with English abstract)

[29] 夏倩柔，熊黑鋼，張芳. 綠洲不同尺度地下水時空動態變化特征及成因分析[J]. 新疆農業科學，2013，50(6)：1137－1144. Xia Qianrou, Xiong Heigang, Zhang Fang. Study on the characteristics and causes of temporal－spatial variation of groundwater on different Scales in the oasis of Qitai, Xinjiang[J]. Xinjiang Agricultural Seiences, 2013, 50(6): 1137－1144. (in Chinese with English abstract)

[30] Ramakrishna R Nemani, Charles D Keeling, Hirofumi Hashimoto, et al. Climate－driven increases in global terrestrial net primary production from 1982 to 1999[J]. Science, 2003, 300(5625): 1560－1563.

[31] 宋文娟，熊黑鋼，戚冉. 新疆奇臺縣近42a氣候變化特征分析[J]. 干旱區資源與環境，2009，23(4)：100－106. Song Wenjuan, Xiong Heigang, Qi Ran. Climate change in Qitai County in Xinjiang in recent 42a[J]. Journal of Arid Land Resources and Environment, 2009, 23(4): 100－106. (in Chinese with English abstract)

[32] 苑全治，吳紹洪，戴爾阜，等. 過去50年氣候變化下中國潛在植被NPP的脆弱性評價[J]. 地理學報，2016，71(5)：797－806. Yuan Quanzhi, Wu Shaohong, Dai Erfu, et al. NPP vulnerability of China's potential vegetation to climate change in the past 50 years[J]. Acta Geographica Sinica, 2016, 71(5): 797－806. (in Chinese with English abstract)

[33] 楊建光，艾根生. 奇臺：水危機引發一場“用水革命”.昌吉日報[R]. 2009-05-11[2016-10-01].

Dynamic change of net primary productivity during process of new artificial oasis expansion in Xinjiang based on remote sensing

Zhang Fang1,2, Xiong Heigang2,3, Feng Juan1,2, Xu Zhonglin1,2

(1.830046,; 2.,830046,; 3.100083)

The impact of human activities on the global carbon cycle is one of the main concerns in the global climate change research. However, the mechanisms of the carbon cycle for the different regions, natural and human factors are very complex and different. Exploring the impact of human reclamation activities on the distribution pattern of net primary productivity (NPP) in arid oasis ecosystem will help to increase the understanding of human activities on terrestrial ecosystems and carbon cycle impacts in arid regions. This study aimed to investigate the change of NPP during the process of new artificial oasis expansion in Xinjiang based on the remote sensing data. June was selected as the optimal time for extracting vegetation cover types change information. Based on the long-term Landsat remote sensing data, the Carnegie-Ames-Stanford Approach (CASA) model was used to estimate the change of NPP during the expansion of Xinjiang Qitai artificial oasis from June 1976 to June 2016. The results showed that: 1) In 1976, the study area mainly was dominated by desert steppe, and the farmland was only 5 150 hm2, accounting for 3.25% of the area. By 2016, the amount of farmland expanded to 64 675hm2, which increased by more than 11 times; 2) While the artificial oasis was expanding, spring overflow zone disappeared, the swampland and the reservoir dried up, which resulted in the decrease of shrub area from 17 871 hm2to 5 202 hm2; 3) Reclamation resulted in a significant decline in groundwater level, soil salinity development process was interrupted, and bare alkaline area was reduced from 58 859 hm2in 1976 to 17 803 hm2in 2016, which was reduced by nearly 70% and began to succession to the new grassland; 4) In June 1976, the unit area NPP of farmland was 68.6 g/m2, and the unit area NPP of farmland in June 2007 and June 2016 was over 90 g/m2. In the recent 40 years, the unit area NPP of grassland was fluctuated within the range of 30-40 g/m2, which was relatively stable; 5) In June 1976, the total NPP in the study area was 41 018 t, of which the contribution rate of farmland was only 8.62%. In June 2016, the total NPP of the study area reached to 90 602 t, which was 2.2 times higher than that in June 1976. The contribution rate of farmland to total NPP was 65.23%. In the same period, the net loss of shrub and swamp NPP was 8 265 t and 3 319 t, respectively; 6) Human activities were the decisive factors that influenced the distribution pattern of NPP in the study area. And NPP in 1976 were mainly composed of natural grassland and shrub, accounting for 56.68% and 26.61% respectively, NPP in 2016 was mainly composed of farmland and fresh grassland, accounting for 65.23% and 31.85%, respectively. In 1976, 66.93% of the shrubs, 75.58% of the swamps, 41.43% of the grassland, 23.86 % of the bare alkali land and 77.96% of the water bodies were evolved into farmland in 2016. The 43.62% of the grassland in 1976 was still grassland after 40 years, but was divided into fragments by the newly reclaimed farmland. Only 15.01% of the bare alkaline area in 1976 remained in 2016, while the rest came from degradation of grasslands, shrubs and swamps. In the recent 40 years, the original natural oasis landscape in the study area had been completely changed, and large areas of farmlands connected to each other become modern agriculture oasis that was suitable for large-scale agricultural machinery operations.

remote sensing; grassland; farmland; new artificial oasis expansion; net primary productivity (NPP); Carnegie- Ames-Stan-ford Approach (CASA); Xinjiang Qitai oasis

10.11975/j.issn.1002-6819.2017.12.025

X171;TP 79

A

1002-6819(2017)-12-0194-07

2016-12-13

2017-05-10

國家自然科學基金（41261049、41671198、41361098）；中國博士后基金面上資助（2013M532100）

張 芳，女，山東嘉祥人，博士，副教授，研究方向：干旱區環境演變與人地關系。烏魯木齊 新疆大學資源與環境科學學院，830046。 Email：zhangf602@163.com

張 芳，熊黑鋼，馮 娟，許仲林. 基于遙感的新疆人工綠洲擴張中植被凈初級生產力動態變化[J]. 農業工程學報，2017，33(12)：194－200. doi：10.11975/j.issn.1002-6819.2017.12.025 http://www.tcsae.org

Zhang Fang, Xiong Heigang, Feng Juan, Xu Zhonglin. Dynamic change of net primary productivity during process of new artificial oasis expansion in Xinjiang based on remote sensing[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(12): 194－200. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2017.12.025 http://www.tcsae.org