2000—2015年山東省耕地產(chǎn)能的時空格局

王赫彬 王文娟 商令杰

(1.中國科學院 東北地理與農(nóng)業(yè)生態(tài)研究所，長春 130102；2.中國科學院大學，北京 100049；3.山東師范大學 地理與環(huán)境學院，濟南 250014)

新中國成立以來，我國長期致力于糧食增產(chǎn)增收和糧食安全工作，采取一系列改善農(nóng)業(yè)生產(chǎn)條件、提高糧食生產(chǎn)量的措施，取得了很大成效[1]。糧食產(chǎn)量穩(wěn)定是一項關系國計民生的大事，而耕地生產(chǎn)力和耕地產(chǎn)能的研究對估算耕地糧食產(chǎn)量和保障糧食安全具有重要意義。利用耕地生態(tài)系統(tǒng)生產(chǎn)力模擬耕地產(chǎn)能，是一種估算耕地糧食產(chǎn)量的重要手段，并且遙感也是模擬陸地生態(tài)系統(tǒng)生產(chǎn)力的重要途徑。基于遙感觀測的生態(tài)系統(tǒng)NPP研究早已開始應用于實時、連續(xù)監(jiān)測和年際波動和長期變化趨勢的探測[2]。NPP能夠以統(tǒng)一的尺度標準反映不同生態(tài)系統(tǒng)的生產(chǎn)力，是一個非常客觀的農(nóng)田生產(chǎn)力衡量指標，避免了農(nóng)業(yè)結構調整和作物品種變化等因素對耕地生產(chǎn)力測定的干擾[3]。

國外學者從個體、斑塊、景觀、區(qū)域和全球等各個尺度展開NPP研究，1932年丹麥Boysen-Jensen P首次通過公式提出生物生產(chǎn)力的總生產(chǎn)量(Gross produce)和凈生產(chǎn)量(Net produce)概念[4]。20世紀60年代，NPP測定工作在世界范圍內(nèi)展開，1963—1972年，國際生物圈計劃(IBP)研究世界范圍內(nèi)主要生態(tài)系統(tǒng)的物質循環(huán)和生物生產(chǎn)力，基于多種估算模型評估區(qū)域和全球范圍的NPP空間分布及碳平衡；1973年Helmut[5]建立Miami模型，認為NPP是年平均氣溫和年降水量的函數(shù)；20世紀80年代，國際地圈-生物圈計劃(IGBP)致力于NPP研究，1985年Zenbei等[6]建立適用于干旱半干旱區(qū)域的Chikugo模型，1987年Hunt等[7]利用BIOME-BGC模型估算當年全球陸地NPP；Kazuhito等[8]通過BIOME-BGC模型和NCEP氣象數(shù)據(jù)研究熱帶生態(tài)系統(tǒng)碳循環(huán)對氣象因子敏感性。90年代后期遙感方法估測NPP 迅速發(fā)展，David等[9]和Prince等[10]基于植物生理脅迫和呼吸作用建立全球第一個純遙感NPP模型，估算總、凈初級生產(chǎn)力；Christopher等[11]利用CASA模型研究全球陸地NPP分布，Ruimy等[12]修訂CASA模型參數(shù)，提升其實用性，Tate等[13]利用CASA模型估算了新西蘭的NPP。我國陸地生態(tài)系統(tǒng)生產(chǎn)力研究始于20世紀70年代，研究重點包括植被NPP估算模型、全球氣候變化背景下的我國自然植被NPP時空格局。20世紀70年代后期，我國完成國內(nèi)陸地生態(tài)系統(tǒng)生產(chǎn)力的測定[14]；20世紀80年代后期，我國NPP研究正式開始，在全國尺度上，何勇等[15]利用植被與大氣相互作用模式(AVIM)模擬中國陸地生態(tài)系統(tǒng)NPP分布；孫睿等[16]分析植被指數(shù)與植被吸收的光合有效輻射比例之間的線性關系，利用光能利用率研究我國植被NPP分布及其季節(jié)變化；陳利軍等[17]運用遙感數(shù)據(jù)和光能利用模型估算中國陸地植被凈初級生產(chǎn)力；冀詠贊等[18]利用VPM模型分別估算近15年全國陸地生態(tài)系統(tǒng)總初級生產(chǎn)力、耕地與農(nóng)田生產(chǎn)力；區(qū)域尺度上，曾慧卿等[19]應用BIOME-BGC模型估算1993—2004年紅壤丘陵區(qū)濕地松林總第一性生產(chǎn)力和凈第一性生產(chǎn)力；閆慧敏等[20]利用VPM模型模擬并分析黃淮海多熟農(nóng)業(yè)區(qū)的農(nóng)田生產(chǎn)力時空特征。利用“3S”技術定位獲取和反演地表植被信息和相關生物物理學參數(shù)，現(xiàn)已成為陸地生態(tài)系統(tǒng)植被NPP估算的重要手段，這也成為近十幾年來NPP模型研究和估算的最突出特點。

目前關于NPP的研究趨于成熟，呈現(xiàn)多元化趨勢，從最初的特定氣候帶NPP估算和生態(tài)系統(tǒng)碳循環(huán)研究，到目前以不同驅動模型模擬區(qū)域甚至全球尺度的生態(tài)系統(tǒng)生產(chǎn)力。以NPP作為估算的重要手段，這也成為近十幾年來NPP模型研究和估算的最突出特點。CASA、GLO-PEM和VPM模型為代表的光能利用模型逐漸成為NPP模型研究的重點，但目前多數(shù)研究利用單一模型模擬區(qū)域生態(tài)系統(tǒng)生產(chǎn)力，針對某種土地利用類型的NPP研究較少，研究尺度也多集中在全國尺度及東北三省、內(nèi)蒙古自治區(qū)和西北地區(qū)，針對黃淮海地區(qū)的研究較少。山東是我國重要的產(chǎn)糧大省，2017年山東省糧食總產(chǎn)量4 723.2萬t，居全國第三位，研究山東省耕地產(chǎn)能的時空格局和產(chǎn)能變化對保障我國糧食安全至關重要。

本研究通過MODIS遙感、土地利用和其他非遙感數(shù)據(jù)驅動VPM模型模擬山東省耕地生產(chǎn)力，評價山東省耕地NPP的時空分布和動態(tài)變化特征，以期為實時監(jiān)測耕地產(chǎn)能和估算耕地糧食產(chǎn)量提供可靠手段，以及為區(qū)域糧食預測、耕地產(chǎn)能檢測、農(nóng)業(yè)區(qū)劃調整和制定農(nóng)業(yè)政策提供重要參考。

1 研究方法與數(shù)據(jù)來源

1.1 研究方法

1.1.1耕地生產(chǎn)力計算

植被光合模型(VPM，Vegetation Photosynthesis Model)是一個基于衛(wèi)星遙感和通量觀測數(shù)據(jù)發(fā)展起來的光能利用率模型。自21世紀初模型建立以來，已在全國范圍內(nèi)的20多個站點和多個生態(tài)系統(tǒng)上開展模型的應用、模擬、校驗與驗證工作。如VPM模型運行技術路線圖(圖1)所示，將基礎數(shù)據(jù)輸入VPM模型，基礎數(shù)據(jù)又通過MODIS地表反射率產(chǎn)品據(jù)(MOD91A1)計算得LSWI和EVI數(shù)據(jù)，通過AUNSPLINE氣象插值軟件得到的平均氣溫柵格數(shù)據(jù)集，土地利用數(shù)據(jù)，利用MODIS 1B數(shù)據(jù)反演得到PAR數(shù)據(jù)，從文獻資料中查找不同植被類型的三基溫數(shù)據(jù)等。將以上數(shù)據(jù)輸入，驅動VPM模型，模擬估算出陸地生態(tài)系統(tǒng)總生產(chǎn)力，再將生態(tài)系統(tǒng)總生產(chǎn)力減去植物的呼吸消耗的能量，得到生態(tài)系統(tǒng)凈初級生產(chǎn)力。

圖1 VPM模型技術路線圖Fig.1 VPM model technology roadmap

1.1.2耕地產(chǎn)能計算

耕地產(chǎn)能指在一定地域及時期，在一定經(jīng)濟、社會和技術條件下所形成的耕地生產(chǎn)能力，分為耕地理論、可實現(xiàn)和實際產(chǎn)能3個層次。理論產(chǎn)能是指農(nóng)業(yè)生產(chǎn)條件得到充分保證，在經(jīng)濟技術和投入利用水平最優(yōu)的狀況下，僅考慮光、溫、水和土壤等環(huán)境因素農(nóng)作物所能達到的最高產(chǎn)量；可實現(xiàn)產(chǎn)能是指在農(nóng)業(yè)生產(chǎn)條件得到基本保證，光、熱和水土等環(huán)境因素均處于正常狀態(tài)，在現(xiàn)有經(jīng)濟技術和利用水平的最高投入下正常年景農(nóng)作物能夠達到的最高產(chǎn)量；實際產(chǎn)能指目前已經(jīng)實現(xiàn)的生產(chǎn)能力，即某年農(nóng)作物已達到的產(chǎn)量[21]。

通過耕地糧食單產(chǎn)和耕地NPP之間的線性擬合，實現(xiàn)由耕地生產(chǎn)力向耕地產(chǎn)能的轉變。根據(jù)2000—2015年山東省年鑒獲取山東17地市耕地糧食產(chǎn)量，篩選產(chǎn)糧大市作為樣本市，統(tǒng)計產(chǎn)糧大市的糧食單產(chǎn)數(shù)據(jù)與其單位耕地NPP均值擬合，得到產(chǎn)能關系式，實現(xiàn)通過產(chǎn)能關系式估產(chǎn)的效果。選取產(chǎn)糧大市作為樣本市，產(chǎn)糧大市的農(nóng)業(yè)較為發(fā)達，種植結構相對穩(wěn)定，因此作為樣本更為精確。以2010—2015年山東省各市主要農(nóng)作物播種面積統(tǒng)計數(shù)據(jù)為基準，選取糧食播種面積占農(nóng)作物播種總面積比例>70%的產(chǎn)糧大市作為計算樣本市。具體選擇的樣本有濟南、青島、淄博、棗莊、煙臺、濰坊、德州、聊城和濱州市9個地級市，共54個樣本點。各地市耕地NPP均值數(shù)據(jù)是將2010—2015年山東省500 m×500 m分辨率的耕地NPP數(shù)據(jù)按山東省各地市數(shù)據(jù)計算獲得。

1.2 數(shù)據(jù)來源

遙感數(shù)據(jù)為2000—2015年美國國家航空航天局(NASA)地球觀測系統(tǒng)(EOS)Terra衛(wèi)星搭載的中分辨率成像光譜儀傳感器陸地產(chǎn)品第5版數(shù)據(jù)集中的地表反射率產(chǎn)品MOD09A1(http:∥ladsweb.nascom.nasa.gov)，空間分辨率為500 m，時間分辨率為8 d。

溫度數(shù)據(jù)為國家氣象局觀測站記錄的2000—2015年逐日平均氣溫，該通過中國氣象科學數(shù)據(jù)共享服務網(wǎng)(http:∥data.cma.gov.cn/site/index.html)地面氣候資料日值數(shù)據(jù)集下載，為觀測站點數(shù)據(jù)。

土地利用數(shù)據(jù)來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心(http:∥www.resdc.cn)，包含山東省2000、2005、2010和2015年4期土地利用數(shù)據(jù)。

高分辨率EVI數(shù)據(jù)來源于Google Earth Engine云平臺中的Landsat TM/OLI EVI數(shù)據(jù)集，具體數(shù)據(jù)集包括Landsat 5 TM 8-Day EVI Composite和Landsat 8 8-Day EVI Composite，2種數(shù)據(jù)集分別是以Landsat 5 TM和Landsat 8 OLI的大氣表觀反射率(Top of atmosphere reflectance，TOA)數(shù)據(jù)計算得來。

2010—2015年山東省糧食作物類型、播種面積和作物產(chǎn)量等數(shù)據(jù)主要來源于山東省統(tǒng)計局公開的年度數(shù)據(jù)(http:∥www.stats-sd.gov.cn/col/col6279/index.html)，通過2011—2016年山東省統(tǒng)計年鑒獲得。

2 山東省耕地NPP的時空格局分析

2.1 耕地NPP的時間變化特征

從2000—2015年山東省耕地NPP均值和總量變化圖(圖2)觀測，2000年的NPP均值和總量最低，分別為711.67 g/(m2·年)和7.05×1013g/年；2015年NPP均值和總量最高，分別為834.43 g/(m2·年)和8.08×1013g/年。2000—2015年多年耕地NPP均值為766.60 g/(m2·年)，多年耕地NPP總量的均值為7.51×1013g/年。2000—2015年山東省耕地NPP的均值和總量有相似的變化規(guī)律，雖然變化和波動較大，但總體呈波動上升趨勢。通過NPP總量和均值變化折線圖可以看出兩者有著相同的變化趨勢，即表現(xiàn)為同增同減。具體變化為：2000—2004年耕地NPP上升幅度較快，2004—2006年大幅滑落，2006—2007年大幅上升，2007—2013年波動下降，2013年之后穩(wěn)健上升。

圖2 2000—2015年山東省NPP均值和總量變化(E代表10的次方)Fig.2 Changes in the mean and total amount of NPP in Shandong Province from 2000 to 2015 (E is a power of 10)

2.2 耕地NPP的空間格局分布

經(jīng)過VPM模型運算，得到山東省2000—2015年的陸地生態(tài)系統(tǒng)凈生產(chǎn)力分布圖，結合山東省土地利用數(shù)據(jù)得到2000—2015年的山東省耕地NPP空間分布圖(圖3)。總體來看，除個別區(qū)域外，山東省大部分區(qū)域NPP值較高。山東省為暖溫帶季風氣候，耕地植被類型中以冬小麥-夏玉米的兩熟耕地為主，一年兩作，高強度的農(nóng)業(yè)活動促使耕地植被生產(chǎn)力較高。2015年耕地生產(chǎn)力>850 g/(m2·年)的耕地分布在魯西南平原和魯西北平原地區(qū)，該地區(qū)土壤肥沃、土地平坦和灌溉設施較完善，耕地生產(chǎn)力高；耕地生產(chǎn)力為700～850 g/(m2·年)的耕地分布廣泛，集中在魯西北平原北部和山東半島中部地區(qū)；耕地生產(chǎn)力為550～700 g/(m2·年)的耕地廣泛分布于魯中南山地和山東省北部的濱州和東營地區(qū)，其多為山間丘陵地帶，耕地較少且分散，灌溉設施不足，生產(chǎn)力較低；耕地生產(chǎn)力<550 g/(m2·年)的耕地最為集中，位于黃河三角洲的沿海地區(qū)，其地勢低，易受海水倒灌影響，多為鹽堿地，耕地生產(chǎn)力低。從山東省農(nóng)業(yè)區(qū)劃分析，山東省中耕地植被凈生產(chǎn)力均處于較高水平，耕地NPP均值達680 g/(m2·年)以上，雖然每個農(nóng)業(yè)區(qū)劃的耕地面積和NPP總量差異較大，但耕地NPP均值相差不大。其中山東省丘陵水澆地二熟旱坡地花生棉花一熟區(qū)耕地面積最多，達到63 178.25 km2；耕地面積占比為62.59%；耕地生產(chǎn)力均值最低，為680.81 g/(m2·年)；耕地NPP總量為4.30×1013g/(m2·年)，占比為57.85%。黑龍港缺水低平原水澆二熟旱地一熟區(qū)耕地面積和NPP總量最少，分別為806.75 km2和0.06×1013g/(m2·年)，耕地面積占比為0.8%，耕地生產(chǎn)力均值為754.32 g /(m2·年)。黃淮平原南陽盆地旱地水澆地兩熟區(qū)的耕地面積為14 796.50 km2，耕地面積占比14.66%，但是耕地生產(chǎn)力均值最高，為888.56 g/(m2·年)(表1)。

圖3 山東省2000—2015年耕地NPP空間分布Fig.3 Spatial distribution of NPP in cultivated land in Shandong Province from 2000 to 2015

2.3 山東省耕地NPP空間分布狀況及變化原因

對山東省17地級市NPP分布情況進行分析，耕地NPP均值>800 g/(m2·年)的地市有4個，分別是菏澤、聊城、德州和濟寧市，其中菏澤市最高，為910.45 g/(m2·年)；這些地級市主要分布在魯西北平原區(qū)，土地肥沃、灌溉條件好和生產(chǎn)力高。耕地生產(chǎn)力<650 g/(m2·年)的城市有萊蕪、淄博、煙臺和東營市，其中東營市耕地生產(chǎn)力最低，為516.54 g/(m2·年)。究其原因，萊蕪和淄博市位于魯中南山地丘陵區(qū)、煙臺市位于山東半島丘陵區(qū)，受地形地貌和水源條件的限制，這些地區(qū)的耕地土壤貧瘠，灌溉水源條件差，耕地生產(chǎn)力較差；東營市位于黃河入?？诘娜侵摞}堿地平原中，受海水倒灌影響，這些耕地多為鹽堿地，土壤貧瘠，耕地生產(chǎn)力最差(表2)。

表1 山東省各農(nóng)業(yè)區(qū)劃耕地凈生產(chǎn)力Table 1 Net productivity of cultivated land in various agricultural areas in Shandong Province

表2 山東省17地級市耕地凈生產(chǎn)力Table 2 Net productivity of cultivated land in 17 cities of Shandong Province

3 山東省耕地產(chǎn)能的時空格局分析

將山東省2010—2015年產(chǎn)糧大市的耕地單產(chǎn)與耕地NPP均值數(shù)據(jù)利用回歸分析法擬合(圖4)，建立耕地產(chǎn)量與單位耕地面積NPP間的回歸關系，進而得到山東省耕地產(chǎn)能計算公式。

y=7.475x+7 473.2(R2=0.722 3)

式中:y為耕地產(chǎn)能，kg/hm2；x為單位耕地面積NPP，g/(m2·年)。

圖4 回歸分析擬合結果Fig.4 Regression analysis fitting result

NPP均值為500～1 100 g/(m2·年)，耕地產(chǎn)量為10 000～16 000 kg/hm2，僅考慮此區(qū)間內(nèi)擬合的物理意義。R表示NPP均值與產(chǎn)量的相關性)

3.1 耕地產(chǎn)能的時間變化特征

由耕地NPP數(shù)據(jù)與產(chǎn)能計算公式，得到2000—2015年山東省耕地產(chǎn)能空間分布圖5，并根據(jù)耕地產(chǎn)能變化趨勢得到圖6。2000—2015年山東省耕地產(chǎn)能總體升高，其中產(chǎn)能<12 000 kg/hm2、12 000～13 000 kg/hm2和>13 000～14 000 kg/hm2的耕地面積均處于減少狀態(tài)，分別從2000年的19 278.50、33 611.75和33 041.50 km2下降到2015年的13 235.25、19 060.25和19 594.50 km2；而產(chǎn)能>14 000 kg/hm2的耕地面積大幅增加，從2000年 13 034.00 km2增至2015年為45 031.00 km2，增長幅度最大，標志著山東省耕地產(chǎn)能提升明顯。2000年，山東省耕地產(chǎn)能>14 000 kg/hm2的耕地面積小，零散分布于魯西南平原；2005和2010年產(chǎn)能>14 000 kg/hm2的耕地持續(xù)增加，分布于魯西北平原中南部和魯西南平原地區(qū)；2015年，魯西南平原、魯西北平原和山東半島內(nèi)青島和濰坊地區(qū)的大部分耕地的耕地產(chǎn)能均>14 000 kg/hm2；2000—2015年，山東省耕地產(chǎn)能在12 000～14 000 kg/hm2的耕地基本都實現(xiàn)了產(chǎn)能提升，絕大多數(shù)此類耕地產(chǎn)能提升至>14 000 kg/hm2。

圖5 山東省2000—2015年耕地產(chǎn)能空間分布圖Fig.5 Spatial distribution map of cultivated land productivity in Shandong Province from 2000 to 2015

圖6 2000—2015年山東省耕地產(chǎn)能變化趨勢Fig.6 Trends of cultivated land productivity in Shandong Province from 2000 to 2015

3.2 耕地產(chǎn)能的空間分布格局

根據(jù)產(chǎn)能分類標準將耕地產(chǎn)能分類為4級，將耕地產(chǎn)能等級按照從大到小依次編號為Ⅰ～Ⅳ類(圖7)。耕地產(chǎn)能在12 000 kg/hm2以下的耕地主要集中于東營市境內(nèi)黃河三角洲的沿海地區(qū)為Ⅰ類，山東半島丘陵地帶和魯中南山區(qū)也有零散分布，面積較小，此類產(chǎn)量耕地占山東省耕地面積16.61%；耕地產(chǎn)能在12 000～13 000 kg/hm2的耕地主要分布于山東半島丘陵區(qū)、魯中南山區(qū)、濱州市和東營市范圍內(nèi)，耕地分布較為分散，面積占比26.16%為Ⅱ類；耕地產(chǎn)能在>13 000～14 000 kg/hm2的耕地分布最廣，全省各地市均有分布，魯西北平原、山東半島中部地區(qū)相對集中，面積占比為28.92%為Ⅲ類；耕地產(chǎn)能>14 000 kg/hm2的耕地集中分布于魯西北平原的聊城市和魯西南平原的菏澤市和濟寧市，面積占比為28.83%為Ⅳ類。

山東省17地市的耕地產(chǎn)能均值為12 798.83 kg/hm2，耕地產(chǎn)能相對較高。其中耕地產(chǎn)能總量最高的為濰坊市，菏澤市、聊城市、德州市、濟寧市和棗莊市的耕地產(chǎn)能均值>13 000 kg/hm2，這些地市主要分布在魯西北和魯西南平原區(qū)，土地肥沃、灌溉條件好、生產(chǎn)力高；濰坊市、泰安市、青島市、濱州市、臨沂市、日照市、威海市、淄博市和煙臺市耕地產(chǎn)能均值為12 000～13 000 kg/hm2；萊蕪市和東營市耕地產(chǎn)能<12 000 kg/hm2，由于萊蕪市位于魯中南山地丘陵區(qū)，受地形地貌和水源條件的限制，耕地土壤貧瘠，水源主要來源于自然降水，灌溉水源條件差，耕地生產(chǎn)力較差，產(chǎn)能較低；東營市位于黃河入?？诘娜侵摞}堿地平原中，多為鹽堿地，土壤貧瘠，耕地產(chǎn)能最差。

3.3 耕地產(chǎn)能的變化原因

影響耕地產(chǎn)能的原因非常復雜，通常會將這些因素分為自然要素和非自然要素。自然要素有氣候變化、土壤條件和地貌，而非自然要素主要是人為的影響，主要包括農(nóng)用土地的整合。本研究結果表明，2000—2015 年山東省耕地產(chǎn)能年際變化波動明顯，總體呈現(xiàn)增加趨勢，與同時期中國耕地產(chǎn)能總體變化趨勢一致。自2000—2015年，山東省整體平均溫度和降水年際間變化較大，但是在長時間序列上并未呈現(xiàn)較大變化，溫度上升了0.41 ℃，降水下降11 mm，這些自然要素的改變對耕地產(chǎn)能的提升影響并不明顯。

但是近年來山東省土地整合力度非常大，山東省2011—2015年農(nóng)業(yè)綜合開發(fā)項目共建設改造項目區(qū)959個，其中高標準農(nóng)田建設項目區(qū)有455個，中低產(chǎn)田改造項目區(qū)有504個。在項目區(qū)中，產(chǎn)能提升幅度超過1 500 kg/hm2的項目區(qū)僅有一個，提升幅度為1 709.48 kg/hm2；產(chǎn)能提升幅度為1 000～1 500 kg/hm2的項目區(qū)共246個，分布相對集中，主要分布在魯西南平原、魯西北平原和山東半島中部平原的高產(chǎn)田上；產(chǎn)能提升幅度為500～1 000 kg/hm2的項目區(qū)共343個，這些項目區(qū)分布較廣泛，在全省各地區(qū)均有較為零散的分布；產(chǎn)能提升幅度為0～500 kg/hm2的項目區(qū)共369個，數(shù)量最多，占比最大，主要分布在在魯西南和魯西北平原，且在魯中南山地和山東半島丘陵區(qū)中也有較多分布。山東省糧食作物種植項目區(qū)中，沒有發(fā)生產(chǎn)能下降的項目區(qū)。

這些農(nóng)業(yè)綜合開發(fā)項目的建設，極大力度的提升耕地的使用效率，提高了耕地的產(chǎn)能，使得山東省近年來的耕地平均產(chǎn)能出現(xiàn)了明顯提升。

圖7 山東省2015年耕地產(chǎn)能值分類Fig.7 Classification of cultivated land productivity values in Shandong Province in 2015

4 結論與討論

4.1 結論

通過GIS方法和衛(wèi)星遙感手段，實現(xiàn)了省級和農(nóng)業(yè)綜合開發(fā)項目區(qū)尺度的耕地產(chǎn)能估算，并從時間和空間尺度分析2000—2015年山東省耕地生產(chǎn)力和耕地產(chǎn)能的時空分布特征和動態(tài)變化狀況，主要結論如下：

1)從時間看，山東省耕地NPP均值和總量均呈波動上升趨勢，NPP均值從2000年711.67 g/(m2·年)提升到2015年834.43 g/(m2·年)，NPP總量從2000年的7.05×1013g/年提升到2015年的8.07×1013g/年，耕地NPP處于較高水平。從空間上看，NPP較高的地區(qū)分布在魯西北平原區(qū)和魯西南平原，NPP值較低的地區(qū)集中在魯中南山地區(qū)、山東半島丘陵區(qū)和黃河入?？诘娜侵摞}堿地區(qū)，其中黃淮平原南陽盆地旱地水澆地兩熟區(qū)耕地生產(chǎn)力均值在4個農(nóng)業(yè)區(qū)劃中最高，菏澤市、聊城市、德州市和濟寧市的耕地NPP均值在17地市中最高。

2)山東省耕地產(chǎn)能總體升高，耕地產(chǎn)能>14 000 kg/hm2的耕地明顯增多，產(chǎn)能<12 000 kg/hm2的耕地面積逐漸減少，其中魯西南平原、魯西北平原、濰坊市東部、青島市中部地區(qū)的耕地產(chǎn)能較高，>13 000 kg/hm2；山東半島丘陵區(qū)、魯中南山地區(qū)的耕地產(chǎn)能多為12 000～13 000 kg/hm2；在黃河三角洲周圍的沿海地區(qū)耕地產(chǎn)能最低，<12 000 kg/hm2。

3)山東省自然環(huán)境要素，如溫度、降水改變較小，對耕地產(chǎn)能的提升影響并不明顯。但山東省2011—2015年農(nóng)業(yè)綜合開發(fā)項目共建設改造項目區(qū)959個，其中高標準農(nóng)田建設項目區(qū)有455個，中低產(chǎn)田改造項目區(qū)有504個。山東省糧食作物種植項目區(qū)中，耕地產(chǎn)能均有較大提升，沒有發(fā)生產(chǎn)能下降的項目區(qū)。

4.2 討論

本研究通過模型分析、GIS空間分析等方法，利用遙感手段研究了山東省耕地產(chǎn)能的時空格局，但研究僅使用VPM模型反演耕地NPP，未來將嘗試運用多種遙感反演NPP模型，諸如CASA模型等，同時對比分析幾種模型的優(yōu)缺點和適用情境。