馬彩虹

(陜西理工學院 秦巴地理研究所， 陜西 漢中 723000)

0 引 言

作物生產(chǎn)潛力是指在一定條件下能夠持續(xù)生產(chǎn)出滿足人類需要的生物產(chǎn)品的內(nèi)在能力，很大程度上受地形地貌、氣候、土壤等自然環(huán)境條件的影響，同時也受經(jīng)濟、技術等和土地資源開發(fā)密切相關的社會因素的影響。作物生產(chǎn)潛力研究與土地、糧食、人口和發(fā)展相關聯(lián)，是當代農(nóng)業(yè)發(fā)展戰(zhàn)略和資源承載力研究的一個重要課題[1]。作物生產(chǎn)潛力的評估方法基本可分為三大類：趨勢外推法、經(jīng)驗公式法以及潛力衰減法[2-3]。其中，潛力衰減法建立在生理生態(tài)學研究基礎上，被公認為估算作物生產(chǎn)力的最基本方法而得到廣泛應用[4-6]。近年來，隨著計算機技術的日漸成熟、數(shù)學模型的發(fā)展以及對作物產(chǎn)量形成機理認識的加深，作物生產(chǎn)潛力研究也逐漸向科學化和精確化的方向推進[7]。隨著空間插值技術的發(fā)展以及其和GIS結合技術的廣泛應用，使以象元為潛力單位計算農(nóng)用地作物生產(chǎn)潛力成為可能[8-9]。陜西黃土臺塬區(qū)是西北黃土臺塬區(qū)的主體[10]，主要分布在河谷盆地周邊以及黃土梁峁丘陵區(qū)之間，塬面平坦，是黃土高原重要的糧食生產(chǎn)區(qū)。該區(qū)域既承受本區(qū)內(nèi)人口與社會經(jīng)濟發(fā)展對其高強度開發(fā)的要求，又潛在承接生態(tài)脆弱區(qū)和河谷盆地經(jīng)濟高密區(qū)轉(zhuǎn)移而來的壓力，存在一定的潛在生態(tài)風險[11]。本文基于GIS和RS技術支持，從土地利用變化視角切入，探討陜西黃土臺塬區(qū)農(nóng)用地作物生產(chǎn)潛力變化，以期為該區(qū)域土地集約利用提供理論依據(jù)。

1 方法與數(shù)據(jù)

1.1 研究區(qū)概況

陜西黃土臺塬位于34°08′～35°52′N、106°20′～110°36′E之間，呈帶狀展布于渭河谷地南北兩側，總面積1 263 773 hm2。按行政隸屬關系分為寶雞臺塬區(qū)、咸陽臺塬區(qū)、西安臺塬區(qū)、銅川臺塬區(qū)和渭南臺塬區(qū)，包含鳳翔縣、岐山縣、扶風縣、乾縣、禮泉縣、涇陽縣、三原縣、富平縣、耀州區(qū)、蒲城縣、白水縣、澄城縣、合陽縣、韓城市的部分區(qū)域。處于由東北季風區(qū)向西北干旱半干旱氣候的過渡區(qū)，年降水量456.3～735.9 mm，主要集中于夏季。土壤以肥力較高的塿土分布最廣，其次是黃墡土，兩者占到研究區(qū)總面積的80.65%，水資源短缺。耕地是構成研究區(qū)自然景觀的基質(zhì)；城鎮(zhèn)用地團塊狀分布，農(nóng)村居民區(qū)點狀散布；林地和草地主要沿溝壑地帶展布，稀疏河流和水庫散布于基質(zhì)上[7]。社會經(jīng)濟方面，2010年總人口1 011.63萬人，人口密度845.84人/km2，其中非農(nóng)人口250.87萬人，占24.8%；GDP為1 827.11億元，其中，第一產(chǎn)業(yè)占12.52%，第二產(chǎn)業(yè)占50.46%，第三產(chǎn)業(yè)占37.01%。

1.2 作物生產(chǎn)潛力估算模型

光合生產(chǎn)潛力是指假定溫度、水分、土壤肥力、農(nóng)業(yè)技術措施均處于最適宜的狀態(tài)時，僅僅由太陽輻射因子所決定的產(chǎn)量。根據(jù)李世奎(1984)的研究思路，光合生產(chǎn)潛力估算模型如下[12]：

(1)

式中：Y為光合生產(chǎn)潛力；ε為光合有效輻射，取值范圍為0.45～0.55；α為植物葉面反射率；β為漏射率；γ為光飽和限制率；ω為植物呼吸作用損耗率(約為30%)；ρ為植物非光合器官的無效吸收，約占光合有效輻射的10%；ψ為量子效率(約為22.4%)；h為單位干物質(zhì)發(fā)熱量(J/g)；∑Q為各月太陽總輻射；CA為植物體的灰分含量(約為5%)；CM為植物體的含水率(約為14%)。

光溫生產(chǎn)潛力又稱熱量生產(chǎn)潛力，是指在水分、土壤肥力和農(nóng)業(yè)技術措施處于適宜狀態(tài)下，由太陽輻射和溫度因子所決定的作物產(chǎn)量上限。計算公式為：

Yt=Y×f(t)。

(2)

①喜涼作物(小麥、馬鈴薯等)溫度訂正系數(shù)

(3)

②喜溫作物(水稻、谷子、玉米等)溫度訂正系數(shù)

(4)

式中Yt為光溫生產(chǎn)潛力，f(t)為溫度訂正系數(shù)，t為月平均氣溫(℃)。

氣候生產(chǎn)潛力又稱光溫水生產(chǎn)潛力，是指在作物的養(yǎng)分供應處于最適宜狀態(tài)時，由太陽輻射、溫度、水分等主要氣候要素所決定的作物生產(chǎn)潛力。氣候生產(chǎn)潛力是在光溫生產(chǎn)潛力的基礎上，進行水分訂正得出，計算公式為：

Yw=Yt×f(w)，

(5)

(6)

E0=0.001 8×(t+25)2×(100-f)，

(7)

式中Yw為氣候生產(chǎn)潛力，f(w)為水分訂正系數(shù)，R為月降水量，E0為月蒸發(fā)量(mm)，t為月平均氣溫(℃)，f為月相對濕度(%)。

土壤生產(chǎn)潛力也稱為作物的自然生產(chǎn)潛力，即由光、溫、水和土壤肥力4個因子共同決定的作物產(chǎn)量潛力。計算公式為：

Ys=Yw×f(s)，

(8)

f(s)=U×O，

(9)

(10)

(11)

式中Ys為土壤生產(chǎn)潛力，f(s)為土壤肥力的訂正系數(shù)，U為土壤中氮磷鉀等養(yǎng)分含量訂正系數(shù)，O為土壤中有機質(zhì)訂正系數(shù)，N為土壤中的堿解氮含量(mg/kg)，P為土壤總的速效磷含量(mg/kg)，K為土壤中速效鉀含量(mg/kg)，M為土壤中有機質(zhì)含量(g/kg)。

1.3 氣候生產(chǎn)潛力目標下作物增產(chǎn)潛力預測

理論增產(chǎn)潛力是理論單產(chǎn)與最高單產(chǎn)的差。計算公式為：

LTi=YFi-YPi，

(12)

式中LTi為i區(qū)域耕地理論增產(chǎn)潛力，YFi為i區(qū)域的理論單產(chǎn)，YPi為i區(qū)域的最高單產(chǎn)。以氣候生產(chǎn)潛力作為最大理論生產(chǎn)潛力，以2010年研究區(qū)耕地面積和現(xiàn)實人口承載力為基礎，計算臺塬區(qū)各區(qū)縣糧食單產(chǎn)、糧食總產(chǎn)潛力。

1.4 數(shù)據(jù)來源與處理

本文所需土地利用變化數(shù)據(jù)主要來源于Landsat陸地資源衛(wèi)星影像數(shù)據(jù)，并采用ERDAS 9.2軟件進行遙感影像校正、投影轉(zhuǎn)換及解譯，得到研究區(qū)土地利用解譯圖。生產(chǎn)潛力計算數(shù)據(jù)及處理流程如下：先將陜西黃土臺塬區(qū)多個氣象站點各月的輻射、海平面平均氣溫、降水量、蒸發(fā)量近30年的均值進行Kriging內(nèi)插，獲取各氣候因子1—12月30 m×30 m的柵格數(shù)據(jù)；然后利用30 m×30 m DEM對月平均氣溫作高程矯正，獲取實際地形氣溫分布圖。再依據(jù)溫度訂正系數(shù)模型、水分訂正系數(shù)模型、土壤肥力訂正系數(shù)模型獲取陜西黃土臺塬區(qū)農(nóng)用地生產(chǎn)潛力各月的溫度、水分及肥力訂正系數(shù)分布圖。最后依據(jù)光溫階乘模型得出研究區(qū)各月作物各類生產(chǎn)潛力。陜西黃土臺塬區(qū)內(nèi)的氣象站點分別是銅川站點(53947)、耀縣站點(57037)、鳳翔站點(57025)、寶雞站點(57016)、武功站點(57034)、西安站點(57036)、華縣站點(57046)。

2 結果分析

2.1 農(nóng)用地作物生產(chǎn)潛力空間分布特征

計算多年平均的作物生產(chǎn)潛力發(fā)現(xiàn)，黃土臺塬區(qū)農(nóng)用地的光合生產(chǎn)潛力密度空間變化區(qū)間為55.69～65.92 t/(hm2·a)，均值為61.47 t/(hm2·a)；光溫生產(chǎn)潛力密度空間變化區(qū)間為22.29～39.56 t/(hm2·a)，均值為32.21 t/(hm2·a)；氣候生產(chǎn)潛力密度空間變化區(qū)間為8.61～16.34 t/(hm2·a)，均值為11.98 t/(hm2·a)；土壤生產(chǎn)潛力密度空間變化區(qū)間為3.05～15.45 t/(hm2·a)，均值為5.26 t/(hm2·a)?？臻g分布情況為(圖1)，光合生產(chǎn)潛力總體趨勢為由東北向西南遞減；光溫生產(chǎn)潛力表現(xiàn)為寶雞臺塬區(qū)、咸陽臺塬區(qū)北部、銅川臺塬區(qū)、西安臺塬區(qū)以及渭南臺塬區(qū)北部邊緣區(qū)較低，臺塬區(qū)中部的咸陽臺塬區(qū)南部及東部、渭南臺塬區(qū)南部和東部較高；氣候生產(chǎn)潛力分布特點表現(xiàn)為西安臺塬區(qū)、渭南臺塬區(qū)南部、咸陽臺塬區(qū)南部和寶雞臺塬區(qū)東部較高，其它區(qū)域較低；土壤生產(chǎn)潛力總體上黃土臺塬區(qū)南部較高，北部較低，河谷區(qū)的土壤生產(chǎn)潛力數(shù)值顯著高于其它地區(qū)，這與河谷地帶土壤質(zhì)量明顯優(yōu)與其它區(qū)域有關。

(a)光合生產(chǎn)潛力 (b)光溫生產(chǎn)潛力

(c)氣候生產(chǎn)潛力 (d)土壤生產(chǎn)潛力

2.2 土地利用變化對農(nóng)用地作物生產(chǎn)潛力的影響分析

1985—2010年，耕地、草地、林地、水域、未利用地均表現(xiàn)為減少，依次減少1 285.01、11 065.22、1 186.44、748.957 7、35.547 6 hm2，分別占1985年相應地類面積的1.29%、7.18%、4.60%、11.64%、11.24%；建設用地增加，共增加了33 902.8 hm2，增加比例為1985年的47.11%。研究時段的土地利用變化導致研究時段內(nèi)陜西黃土臺塬區(qū)農(nóng)用地生產(chǎn)潛力均有不同程度的流失，光合生產(chǎn)潛力、光溫生產(chǎn)潛力、氣候生產(chǎn)潛力及土壤生產(chǎn)潛力損失量分別為16 670.19、8 579.06、3 185.58、1 147.98 t，各占其對應級別潛力總量的2.37%、2.33%、2.33%、1.92%(表1)。

伴隨著土地利用的變化，農(nóng)用地作物生產(chǎn)潛力也發(fā)生相應變化。研究時段中，耕地的作物生產(chǎn)潛力降低的絕對數(shù)量最大，如光合生產(chǎn)潛力由1985年的598 081.38 t減少為2010年的589 447.66 t，凈減少8 633.72 t，占減少總額度的51.79%，草地潛力降低程度居于第二，凈減少6 377 t，占38.25%，林地光合生產(chǎn)潛力凈減少1 659.42 t，占9.95%；其它各類潛力如光溫生產(chǎn)潛力、氣候生產(chǎn)潛力和土壤生產(chǎn)潛力的損失特征，均表現(xiàn)為耕地的相關潛力損失比最小，其次為草地和林地。由于耕地數(shù)量大，其各類潛力損失的總額最大，對臺塬區(qū)農(nóng)用地作物生產(chǎn)潛力的損失程度貢獻最高。研究時段中光合生產(chǎn)潛力、光溫生產(chǎn)潛力、氣候生產(chǎn)潛力和土壤生產(chǎn)潛力的損失分別占黃土臺塬區(qū)生產(chǎn)潛力損失總額的51.79%、55.04%、51.73%、57.34%；草地的作物生產(chǎn)潛力的損失度居于第二，如草地光合生產(chǎn)潛力、光溫生產(chǎn)潛力、氣候生產(chǎn)潛力和土壤生產(chǎn)潛力的損失分別占黃土臺塬區(qū)生產(chǎn)潛力損失總額的38.25%、36.50%、40.89%、37.41%；林地的作物生產(chǎn)潛力損失對臺塬區(qū)生產(chǎn)潛力的影響最小，如林地光合生產(chǎn)潛力、光溫生產(chǎn)潛力、氣候生產(chǎn)潛力和土壤生產(chǎn)潛力的損失分別占黃土臺塬區(qū)生產(chǎn)潛力損失總額的9.95%、8.46%、7.38%、5.26%。因此，保護臺塬區(qū)農(nóng)用地潛力的重中之重仍然是對耕地的大力保護，其次是對草地的保護。

表1 1985—2010年陜西黃土臺塬區(qū)農(nóng)用地土地潛力變化 (單位：t)

對研究時段內(nèi)土地利用變化造成的農(nóng)用地生產(chǎn)潛力轉(zhuǎn)移進行統(tǒng)計，得出表2。

表2 1985—2010年農(nóng)用地轉(zhuǎn)移方向與生產(chǎn)潛力變化表 (單位：t)

分析1985—2010年黃土臺塬區(qū)土地潛力流向，發(fā)現(xiàn)耕地、林地、草地作物生產(chǎn)潛力的轉(zhuǎn)出均大于轉(zhuǎn)入，表現(xiàn)為凈流出；惟獨建設用地表現(xiàn)為凈流入。農(nóng)用地轉(zhuǎn)為建設用地后，一般很難在較短時間內(nèi)重新轉(zhuǎn)為生產(chǎn)性用途。因此，農(nóng)用地轉(zhuǎn)為建設用地等非農(nóng)業(yè)用地是土地利用轉(zhuǎn)移中最值得警惕的轉(zhuǎn)化方向。1985—2010年，由耕地、草地、林地等的轉(zhuǎn)出造成的光合生產(chǎn)潛力、光溫生產(chǎn)潛力、氣候生產(chǎn)潛力和土壤生產(chǎn)潛力損失值分別為30 719.01、16 297.96、6 096.15、2 549.35 t，其中，建設用地擴張造成的農(nóng)用地生產(chǎn)潛力損失約占50%，可見，建設用地的擴張是陜西黃土高原臺塬區(qū)農(nóng)用地生產(chǎn)潛力損失的主要原因。

2.3 氣候生產(chǎn)潛力目標下作物增產(chǎn)潛力預測

在氣候生產(chǎn)潛力為理論潛力水平下，陜西黃土臺塬區(qū)農(nóng)用地的作物生產(chǎn)潛力具有較大的提升區(qū)間(表3)。其中，渭濱區(qū)、金臺區(qū)、華縣、灞橋區(qū)、澄城縣、白水縣、蒲城縣、臨渭區(qū)提升值在8～10 t/hm2之間；韓城市、合陽縣、王益區(qū)、禮泉縣、陳倉區(qū)、大荔縣、耀州區(qū)、藍田縣、永壽縣、長安區(qū)、興平市、乾縣提升區(qū)間在6～8 t/hm2之間；涇陽縣、富平縣、臨潼區(qū)、鳳翔縣、周至縣、閻良區(qū)、秦都區(qū)、扶風縣、渭城區(qū)、三原縣、武功縣、岐山縣提升區(qū)間在4～6 t/hm2之間。因此，黃土臺塬區(qū)的農(nóng)用地的氣候理論水平作物生產(chǎn)潛力提升區(qū)間仍然較大，集約利用的前景較為廣闊。

表3 陜西黃土臺塬區(qū)農(nóng)用地作物生產(chǎn)潛力提升區(qū)間預測 (單位：t/hm2)

3 結 論

(1)陜西黃土臺塬區(qū)的農(nóng)用地作物生產(chǎn)潛力的空間差異較大。其中，光合生產(chǎn)潛力總體趨勢為由東北向西南遞減，光溫生產(chǎn)潛力表現(xiàn)為咸陽臺塬區(qū)南部及東部、渭南臺塬區(qū)南部和東部高于其它地區(qū)，氣候生產(chǎn)潛力表現(xiàn)為西安臺塬區(qū)、渭南臺塬區(qū)南部、咸陽臺塬區(qū)南部和寶雞臺塬區(qū)東部高于其它地區(qū)；土壤生產(chǎn)潛力總體上黃臺臺塬區(qū)南部較高，北部較低。

(2)1985—2010年，陜西黃土臺塬區(qū)農(nóng)用地作物生產(chǎn)潛力均有不同程度的流失，光合生產(chǎn)潛力、光溫生產(chǎn)潛力、氣候生產(chǎn)潛力及土壤生產(chǎn)潛力損失量分別為16 670.19、8 579.06、3 185.58、1 147.98 t，各占其對應潛力類型總量的2.37%、2.33%、2.33%、1.92%。

(3)建設用地的擴張對區(qū)域農(nóng)用地生產(chǎn)潛力有明顯的負面影響。1985—2010年中，農(nóng)用地作物生產(chǎn)潛力的變化結果為，由耕地、草地、林地等的轉(zhuǎn)出造成的光合生產(chǎn)潛力、光溫生產(chǎn)潛力、氣候生產(chǎn)潛力和土壤生產(chǎn)潛力損失值分別為30 719.01、16 297.96、6 096.15、2 549.35 t，其中，建設用地擴張造成的農(nóng)用地生產(chǎn)潛力損失約占50%。

(4)陜西黃土臺塬區(qū)的農(nóng)用地的氣候理論水平增產(chǎn)潛力仍然較大，集約利用的前景較為廣闊。其中，有8個區(qū)縣的臺塬區(qū)的農(nóng)用地增產(chǎn)潛力在8～10 t/hm2之間，12個區(qū)縣的臺塬區(qū)農(nóng)用地作物增產(chǎn)潛力在4～6 t/hm2之間。

[參考文獻]

[1] 李小兵.李子口流域土地生產(chǎn)潛力及人口承載潛力估算[J].中國人口·資源與環(huán)境,2012,22(5):383-384.

[2] 李相璽,尹忠東,何長高.土地生產(chǎn)潛力研究綜述[J].水土保持學報,2001,15(5):33-36.

[3] HARRIS J M,KENNEDY S.Carrying capacity in agriculture:global and regional issues[J].Ecological Econ-omics,1999(29):443-461.

[4] UNESCO F.Carrying capacity assessment with a pilot study of Kenya:A resource accounting methodology for exploring national options for sustainable development[C].Paris and Rome,1985.

[5] 王立輝,黃進良,杜耘,等.南水北調(diào)中線丹江口庫區(qū)土地生產(chǎn)潛力估算[J].長江流域資源與環(huán)境,2011,20(S1):66-73.

[6] 熊利亞,夏朝宗,劉喜云,等.基于RS和GIS的土地生產(chǎn)力與人口承載量[J].地理研究,2004,23(1):10-18.

[7] 劉新衛(wèi),陳百明.黃土丘陵區(qū)安塞縣縣域糧食生產(chǎn)潛力及其開發(fā)[J].農(nóng)業(yè)工程學報,2004,20(6):286-290.

[8] 朱蕾,黃敬峰.基于GIS和空間插值技術的南方丘陵山區(qū)土地生產(chǎn)潛力估算[J].浙江大學學報:農(nóng)業(yè)與生命科學版,2007, 33(4):458-465.

[9] 石淑芹,陳佑啟,李正國,等.基于空間插值技術和輔助信息的吉林省玉米土地生產(chǎn)潛力研究[J].干旱地區(qū)農(nóng)業(yè)研究,2011, 29(5):210-218.

[10] 齊矗華.黃土高原侵蝕地貌與水土流失關系研究[M].西安:陜西人民教育出版社,1991:83.

[11] 馬彩虹,任志遠,李小燕.黃土臺塬區(qū)土地利用轉(zhuǎn)移流及空間集聚特征分析[J].地理學報,2013,68(2):257-267.

[12] 劉揚,賈樹海,那波.土地生產(chǎn)潛力計算方法研究[J].中國農(nóng)學通報,2005,21(12):377-381.