何英彬，姚艷敏，李建平，黃 青，陳佑啟，陳仲新，楊 鵬，于士凱

(1.農(nóng)業(yè)部資源遙感與數(shù)字農(nóng)業(yè)重點開放實驗室，北京 100081;2.呼倫貝爾草原生態(tài)系統(tǒng)國家野外科學(xué)觀測研究站，北京 100081;3.中國農(nóng)業(yè)科學(xué)院農(nóng)業(yè)資源與農(nóng)業(yè)區(qū)劃研究所，北京 100081)

1 引言

土地適宜性評價，尤其是糧食作物種植用地適宜性評價已經(jīng)成為開發(fā)與保護、人口與資源關(guān)系的重要內(nèi)容［1］。相關(guān)領(lǐng)域研究從大空間尺度到中、小空間尺度過渡的特點仍在延續(xù)［2，3］。鄭宇等對山東省土地適宜性做了空間分析及其優(yōu)化配置研究［4］，邢世和等針對福建省農(nóng)用地資源進行適宜性評價分析［5］，謝樹春等在進行湘中紫色土丘陵地區(qū)土地適宜性評價時，將村級單位作為試點研究［6］。近年來，環(huán)境生態(tài)脆弱區(qū)土地適宜性評價明顯增加:一些研究者在內(nèi)蒙古地區(qū)［7，8］、喀斯特地區(qū)［9］、黃土高原地區(qū)［10］、黃土丘陵溝壑區(qū)［11］等生態(tài)脆弱地區(qū)都開展了土地適宜性評價。而單一用途土地適宜性評價一個較為明顯的特點是基本都圍繞經(jīng)濟作物展開［12-17］。國外對這一領(lǐng)域研究也較為深入，Allen等建議土地利用適宜性應(yīng)從地形、土壤、植被及其內(nèi)部關(guān)系著手研究［18］。Steiner等基于生態(tài)調(diào)查的方法設(shè)計了一個土地適宜性研究框架，框架內(nèi)生態(tài)保護的限制作用及積極因素得以充分體現(xiàn)［19］。還有很多國外學(xué)者致力于應(yīng)用GIS軟件深入評價土地適宜性、開發(fā)土地適宜性管理系統(tǒng)及實現(xiàn)可視化功能［20-23］。

目前，作物種植氣候適宜性研究較為集中深入［24-27］。Cutforth和Shay研發(fā)了作物反應(yīng)函數(shù)并分析了玉米種植的溫度適宜性［28］。趙峰等運用模糊數(shù)學(xué)方法，結(jié)合河南省冬小麥生理生態(tài)特征建立了旬溫度適宜度函數(shù)、旬光照適宜度函數(shù)和旬降水適宜度函數(shù)［29］。何英彬等應(yīng)用綜合作物反應(yīng)函數(shù)計算了中國北方11個省市大豆種植氣候適宜性［30］。FAO發(fā)布了“21世紀(jì)全球農(nóng)業(yè)生態(tài)評價:方法與結(jié)果”，報告中系統(tǒng)的應(yīng)用AEZ法研究了灌溉作物與旱作作物的種植適宜性［31］。以往作物適宜性研究雖然比較有意義，但通常沒有建立體現(xiàn)研究區(qū)各方面綜合情況的指標(biāo)體系，尤其是缺乏社會經(jīng)濟因素的定量與空間表達;而對于氣候要素，也通常使用的是全生育期的氣候數(shù)據(jù)，極少有分生育期研究氣候要素對于種植適宜性影響的案例。

大豆是重要的糧油兼用作物，我國是世界大豆第四生產(chǎn)國，僅次于美國、巴西和阿根廷，我國東北三省大豆播種面積約為全國總播種面積的一半［32］。在過去30年年間，東北三省的大豆播種面積經(jīng)歷了幾次波動，波動原因極為復(fù)雜，但種植適宜性是其中很重要的原因。綜上所述，該文以東北三省為例進行了大豆種植適宜性的研究，建立了一套涵蓋自然和經(jīng)濟影響因子的大豆適宜性評價綜合指標(biāo)體系，指標(biāo)體系中不但綜合考慮自然與社會經(jīng)濟影響因素，而且還應(yīng)用作物反應(yīng)函數(shù)計算溫度、降水和日照氣候3個要素在不同生長期對于大豆種植適宜性的影響;在獲取東北三省大豆適宜性等級圖基礎(chǔ)上，將其與東北三省實際大豆種植區(qū)域進行對比，以便為農(nóng)業(yè)及相關(guān)部門適當(dāng)調(diào)整大豆種植的空間結(jié)構(gòu)而服務(wù)，該文也將為其他相關(guān)研究提供參考。

2 方法與數(shù)據(jù)處理

2.1 研究區(qū)域

東北三省 (包括遼寧、吉林和黑龍江3省)介于北緯38°43＇～53°24＇、東經(jīng)115°20＇～135°之間，總面積約為79.18萬km2，人口1.05億人 (圖1)。東北三省地貌類型多樣，地表結(jié)構(gòu)略呈三面環(huán)山、平原中開的盆地輪廓，地勢西北、東北高，中南部低，分布具有一定規(guī)律;東北三省平原沿著主要山體的走向延伸，主要有三江平原 (包括興凱湖平原)、東北三省中部平原 (以松嫩平原為主)以及周圍的山前臺地平原 (含松遼分水嶺)。氣候?qū)贉貛Т箨懶约撅L(fēng)氣候，冬季干燥、嚴(yán)寒漫長，夏季濕潤、溫暖短促，≥10℃積溫除黑龍江北部外，其余地區(qū)均在2 000℃以上;東北三省年降水量為350～1 200mm，大部分地區(qū)在450～850mm之間，由東南向西北遞減。東北三省氣候資源比較豐富，具有太陽輻射較強，夏季氣溫高，降水比較適中，雨熱同期，晝夜溫差大等特點，對農(nóng)作物生長有利［33］。大豆是東北三省主要作物之一，一年一熟制耕作，雖然該區(qū)有很多培育的品種，但總體可以歸結(jié)為春大豆系列［34］。

2.2 研究方法

2.2.1 建立適宜性指標(biāo)體系

評價指標(biāo)的選取應(yīng)遵循代表性和全面性的原則，為避免指標(biāo)高度相關(guān)，只需從諸多因子中選取少數(shù)幾個能夠真實全面反映土地耕作適宜性的評價指標(biāo)即可;評價指標(biāo)的選取并非固定不變，而是因地制宜靈活選取［35］。參考以往的研究結(jié)果結(jié)合大豆生長的生理物候特征，將影響東北三省大豆種植適宜性的因素分為自然與人文兩大類，包含7個方面:自然因子有氣候因子、地形因子、土壤因子、自然災(zāi)害影響因子，而社會經(jīng)濟影響因子包括基礎(chǔ)設(shè)施因子、科技進步因子及農(nóng)民意愿因子［35-38］。各種因子通過可定量化及空間化的指標(biāo)進行大豆種植適宜性的分析:氣候因子量化為氣候適宜度指標(biāo)，地形因子量化為坡度和坡向，土壤因子量化為pH值、有機質(zhì)含量、堿解氮含量、速效磷含量、速效鉀含量和表土含鹽量，自然災(zāi)害因子量化為災(zāi)害率，基礎(chǔ)設(shè)施因子量化為灌溉保證率和農(nóng)村用電量，科學(xué)技術(shù)進步因子量化為地膜使用量、化肥使用量、機械化率，農(nóng)民的意愿量化為耕作半徑長度 (表1)。其中，自然災(zāi)害率、灌溉保證率、農(nóng)村用電量、地膜使用量、化肥使用量、機械化率等指標(biāo)是以社會經(jīng)濟統(tǒng)計數(shù)據(jù)為基礎(chǔ)進行量化及空間化計算。

圖1 東北三省地理位置

2.2.2 氣候適宜度的計算

東北大豆雖然品種眾多，但大致可歸結(jié)為早熟、中熟和晚熟3大種類，大豆物候期也因不同種類而有所變化［34］，經(jīng)研究，中熟是東北三省比較適宜種植也是實際種植最為廣泛的品種，因此本研究以中熟種類物候期為基準(zhǔn)進行研究區(qū)大豆氣候適宜性計算［30］。為了反應(yīng)氣候因子對于大豆不同生育期的影響，可以將東北大豆分為5個基本的生育階段:種子萌發(fā)和出苗期、幼苗生長期、花芽分化期、開花結(jié)莢期和鼓粒成熟期［39］。而氣候因子由溫度、降水和日照3要素組成，與以往通常采用作物全生育期的數(shù)據(jù)作為計算依據(jù)不同，該文針對大豆5個不同生育期，應(yīng)用大豆綜合作物反應(yīng)函數(shù)方法計算氣候3要素不同生育期的適宜度，再通過加權(quán)獲得東北大豆全生育期的氣候適宜度，這樣可以反映不同的氣候要素對于大豆生長不同階段的影響，提高評價精度［26］。

2.2.3 權(quán)值的獲取

指標(biāo)權(quán)重是指標(biāo)對研究對象影響程度的體現(xiàn)，基于專家打分并通過德爾菲法求可獲取各指標(biāo)的權(quán)值［33］。選擇土壤學(xué)、農(nóng)學(xué)、農(nóng)業(yè)經(jīng)濟管理、地理學(xué)、區(qū)域經(jīng)濟學(xué)、生態(tài)學(xué)、環(huán)境科學(xué)等領(lǐng)域以及東北三省省級專家23人，采取獨立填表選取權(quán)數(shù)的形式，請各專家對指標(biāo)的重要性進行打分;然后將他們各自選取的權(quán)數(shù)進行整理和統(tǒng)計分析，最后確定出各指標(biāo)的權(quán)數(shù)。結(jié)合東北的實際情況，對指標(biāo)的內(nèi)容和權(quán)重進行了調(diào)整，將原始海拔、降水及積溫指標(biāo)統(tǒng)一為氣候適宜度而3個指標(biāo)的權(quán)重統(tǒng)一為一個權(quán)重，見表1。

表1 東北三省大豆種植適宜性評價指標(biāo)體系

2.2.4 大豆種植適宜度等級圖獲取及其與實際種植區(qū)域?qū)Ρ确治?/p>

將所有量化指標(biāo)空間化為1km×1km的柵格，并將數(shù)據(jù)標(biāo)準(zhǔn)化至0～1范圍內(nèi)［37］，結(jié)合表1各指標(biāo)的權(quán)值，進行加權(quán)，可以得到東北三省大豆種植適宜度空間分布圖，應(yīng)用自然斷點法可以將其分為適宜種植區(qū)、中間等級區(qū)及非適宜種植區(qū)，生成大豆種植適宜度等級圖。應(yīng)用MODIS影像獲得的東北三省大豆種植空間分布圖與等級圖進行對比，可以獲取大豆種植適宜性的空間特征。

2.3 數(shù)據(jù)的收集與整理

該文應(yīng)用的DEM高程、土壤和自然災(zāi)害等自然環(huán)境數(shù)據(jù)和社會經(jīng)濟統(tǒng)計數(shù)據(jù)均來源于科技部科研院所社會公益研究專項“中國區(qū)域性耕地資源變化影響評價及其預(yù)警”所建立的“東北三省資源環(huán)境基礎(chǔ)數(shù)據(jù)庫”［33］，此數(shù)據(jù)庫在相關(guān)文獻有詳細(xì)敘述，在此不再贅述。而氣象數(shù)據(jù)來自科技部科技基礎(chǔ)性工作專項“中國農(nóng)業(yè)氣候資源數(shù)字化圖集編制——大豆、牧草光溫數(shù)據(jù)數(shù)字化制圖實施方案項目”的數(shù)據(jù)集，數(shù)據(jù)為1980～2010年日均氣溫、日降雨量和日照時數(shù)的30年日平均值［30］。應(yīng)用“東北三省資源環(huán)境基礎(chǔ)數(shù)據(jù)庫”中的灌溉面積、糧播面積、鄉(xiāng)村勞力、化肥純量、地膜使用量、農(nóng)村用電量和農(nóng)機總動力等社會經(jīng)濟統(tǒng)計數(shù)據(jù)，計算出灌溉保證程度 (灌溉面積/糧播面積)、單位面積均化肥施用量 (化肥純量/糧播面積)、機械化程度 (鄉(xiāng)村勞力/農(nóng)機總動力)、單位面積均地膜使用量、單位面積均農(nóng)村用電量等指標(biāo)值;應(yīng)用GIS軟件將這些屬性數(shù)據(jù)與縣級行政區(qū)劃數(shù)據(jù)連接形成矢量數(shù)據(jù)，并將其轉(zhuǎn)換成1km×1km的柵格數(shù)據(jù)。其中，農(nóng)村用電量、機械化程度、單位面積均化肥施用量、地膜使用量等取多年的平均值，以便降低各指標(biāo)年度波動的不穩(wěn)定性。根據(jù)野外調(diào)查記錄的各土壤采樣點社會經(jīng)濟投入數(shù)據(jù)、應(yīng)用地貌類型和土壤類型資料修正縣內(nèi)社會經(jīng)濟柵格化的平均數(shù)據(jù)［33］。應(yīng)用DEM數(shù)據(jù)衍生出1km×1km的坡度圖和坡向圖。由此，所有評價要素數(shù)據(jù)統(tǒng)一為1km×1km的柵格數(shù)據(jù)，標(biāo)準(zhǔn)化處理后可用于評價。

3 結(jié)果與分析

3.1 東北三省大豆適宜性分析

東北三省大豆適宜性等級圖如圖2所示 (森林地區(qū)被剔除)，適宜度在大于等于0.41的區(qū)域為大豆適宜種植區(qū)、大于等于0.36小于0.41之間的區(qū)域為中間等級區(qū)，而非適宜種植區(qū)為適宜度小于0.36的區(qū)域。東北三省總體平均的適宜度為0.389，屬于中上水平。從圖2可知，評價結(jié)果并為受到指標(biāo)數(shù)據(jù)空間化過程中一些以行政單元為基準(zhǔn)進行量化指標(biāo)的限制，而呈現(xiàn)出明顯的行政邊界痕跡。從圖2也可以明顯分辨出:大豆種植適宜區(qū)多分布在東北中東部的松嫰平原、遼寧中東部的遼河平原及三江平原東部邊緣地區(qū)，總面積約為9.09萬km2，約占東北三省全境總面積的11.5%。吉林省和遼寧省分別有3.32萬km2和3.77萬km2的適宜區(qū)，分別占適宜區(qū)總面積的36.5%和41.5%;而黑龍江省所擁有的適宜區(qū)面積僅占適宜區(qū)總面積的22%。

東北三省大豆種植非適宜區(qū)總面積約為7.99萬km2，約占東北三省全境面積的10.11%，與適宜區(qū)相反，非適宜區(qū)主要集中于黑龍江省，該省有6.47萬km2非適宜區(qū)，約占非適宜區(qū)總面積的81%，且非適宜區(qū)集中于三江平原，大小興安嶺山前臺地。而吉林省和遼寧省僅僅有0.54萬km2和0.98萬km2的非種植適宜區(qū)，分別占非適宜種植區(qū)總面積的6.76% 和12.24%，主要集中于長白山山前臺地地區(qū)。

東北三省共有11.45萬km2的中間等級區(qū)域，約占東北三省全境的14.49%。中間等級區(qū)域主要集中于黑龍江省，面積約為8.48萬km2，約為中間等級區(qū)域總面積的74.08%;而吉林和遼寧省擁有的中間等級區(qū)域各占總面積的12.55% 和13.37%。中間等級區(qū)域分布空間特征很明顯，主要是沿河流分布，雨水較多的年份，河流沿岸土地內(nèi)澇，使得大豆種植受到影響，這也是這一地帶并非種植大豆適宜區(qū)的一個主要原因。

3.2 東北三省大豆適宜性等級區(qū)域與實際種植區(qū)域?qū)Ρ确治?/h3>

為了進行東北三省大豆適宜性等級區(qū)域與實際種植區(qū)域?qū)Ρ妊芯?，作者?yīng)用了東北三省大豆種植矢量圖 (圖3)，該圖是由空間分辨率為0.25km×0.25km的MODIS圖像解譯而成，解譯精度為85%以上，符合研究要求［40］。將結(jié)果重采樣為1km×1km柵格，以便與東北三省大豆適宜性等級圖進行疊加，結(jié)果如圖4。

由圖4可知，大豆種植適宜區(qū)與大豆實際種植區(qū)相符合區(qū)域的面積為1.14萬km2，且符合區(qū)主要分布在松遼平原;在相符合的區(qū)域中0.28萬km2分布在吉林省，0.19萬km2分布在遼寧省，0.67萬km2分布在黑龍江省。中間等級區(qū)域與大豆實際種植區(qū)域相符合的區(qū)域面積為2.13萬km2，主要分布在山脈與平原之間的臺地和三江平原。大豆不適宜區(qū)與大豆實際種植區(qū)的符合面積為1.33萬km2，主要分布在吉林省、遼寧省山前臺地及黑龍江的三江平原;符合面積中只有500km2分布在吉林省，289km2分布在遼寧省，而1.25萬km2分布在黑龍江省。對于東北三省研究區(qū)整體而言，大豆種植適宜區(qū)、中間等級區(qū)域及非適宜區(qū)域與大豆實際種植區(qū)域相符合的面積分別占大豆實際種植面積的24.78%，46.30%和28.92%;如果以省為研究單元，吉林省為73.31%，16.53%，10.16%，遼寧省為62.21%，28.41%，9.38%，黑龍江省為15.39%，51.70%，32.91%;不難看出，相對比于吉林省和遼寧省，黑龍江省大豆種植情況略顯不合理。

圖2 東北三省大豆種植適宜性等級

圖3 東北三省大豆種植區(qū)

圖4 東北三省大豆適宜性等級區(qū)域與實際種植區(qū)域?qū)Ρ?/p>

4 討論

(1)在文章中討論的適宜性等級是相對概念，黑龍江省很多地區(qū)被定性為不適宜區(qū)并不是指黑龍江省不適宜種植大豆，而是相對于吉林省和遼寧省，其適宜性等級稍低;更具體地說，相對于松遼平原，黑龍江三江平原山前臺地地區(qū)種植適宜性略差。適宜性評價與評價空間尺度及評價指標(biāo)體系是否統(tǒng)一息息相關(guān)，如果在大空間尺度上，如全國尺度，黑龍江省的大豆適宜性未必為非適宜。東北三省作為我國最重要的糧食主產(chǎn)區(qū)之一及商品糧基地，其作物種植結(jié)構(gòu)的合理性直接關(guān)切到該區(qū)域的糧食生產(chǎn)效率及我國糧食安全。因此，在研究區(qū)域僅限于東北三省的前提下，其又是比較典型的完整的研究區(qū)域，研究大豆種植空間結(jié)構(gòu)相對合理性有其必然意義。

(2)從全球空間尺度到農(nóng)地空間尺度，無法應(yīng)用統(tǒng)一的適宜性評價體系，原因大致可以歸為兩點:一是數(shù)據(jù)無法獲取，二是評價體系指標(biāo)面臨尺度推繹過程中的一致性問題。不同的指標(biāo)在不同的空間尺度上對于適宜性的影響有所不同。因此，在評價作物種植適宜性時，不但要用到自然領(lǐng)域通用的因子如土壤、地形因子，而且要根據(jù)實際情況提出相應(yīng)的社會經(jīng)濟影響因子。該文提出應(yīng)用大豆反應(yīng)函數(shù)計算氣候三要適宜度并加權(quán)獲取氣候適宜度，可以精確評價氣候因子的影響;在此基礎(chǔ)上，綜合考慮農(nóng)民意愿(農(nóng)民耕作半徑)、自然災(zāi)害率及其他社會經(jīng)濟指標(biāo)，可以為大豆種植適宜性精細(xì)評價提供有益參考。在研究中一些以行政區(qū)域邊界空間化的指標(biāo)無疑將會對評價結(jié)果產(chǎn)生一定的影響，為提高空間數(shù)據(jù)準(zhǔn)確度，未來社會經(jīng)濟指標(biāo)空間化問題依然是地理科學(xué)研究的重點問題和熱點。

(3)平原地區(qū)通常是作物種植適宜區(qū)，例如松嫩平原，而山區(qū)或山前臺地很多為不適宜種植區(qū)域，例如吉林省山區(qū)為不適宜區(qū)域。但是根據(jù)圖2，發(fā)現(xiàn)非常有趣的現(xiàn)象是:三江平原是大豆種植非適宜區(qū)的集中地帶，這主要是因為大豆是光感作物，三江平原地區(qū)在大豆第一階段種子萌發(fā)和出苗期、第二階段幼苗生長期和第五階段鼓粒成熟期日照不充足有關(guān)，這對于大豆種植的適宜性產(chǎn)生了不良的影響，另外灌溉條件、農(nóng)民耕作意愿及表土含鹽量也是大豆種植適宜性的負(fù)面影響因素。

(4)從圖4可以看出，松嫩平原很多大豆種植適宜區(qū)并沒有種植大豆，而黑龍江很多相對不適宜種植大豆的地區(qū)種植了大豆，大豆種植適宜度的高低并不能決定大豆空間分布和結(jié)構(gòu)。大豆種植從20世紀(jì)80年代到目前，面積波動較大，主要由農(nóng)民利益最大化的思維驅(qū)動，而農(nóng)民利益最大化的思維通常受到大豆價格與耕作成本的左右，而這一因素又受到大豆種植政策、市場、原油價格甚至國內(nèi)外游資對大豆種植的影響。在不同年份比較效益發(fā)生波動及基本農(nóng)田保護的嚴(yán)厲政策執(zhí)行的前提下，農(nóng)民根據(jù)盈利的需要播種不同的作物，致使作物種植空間競爭較為激烈。該文分析大豆種植適宜性目的在于為未來合理的作物結(jié)構(gòu)調(diào)整提供思路和參考，但并不主張在未考慮農(nóng)民種植意愿前提下改種或不種大豆，只有在農(nóng)民意愿得到充分滿足的提前下，作物結(jié)構(gòu)調(diào)整才能得以順利實施。

5 結(jié)論

該文建立了一套涵蓋自然和經(jīng)濟影響因子的大豆適宜性評價指標(biāo)體系，并以東北三省為例進行了大豆種植適宜性的研究;指標(biāo)體系中不但考慮了溫度、降水和日照氣候三個要素在不同生長期對于大豆種植適宜性的影響，而且還引進了自然災(zāi)害率作為自然評價指標(biāo);此外，在社會經(jīng)濟因子類別中還應(yīng)用了基礎(chǔ)設(shè)施、科技進步及農(nóng)民意愿等因子，進行大豆適宜性評價。在獲取東北三省大豆適宜性等級圖基礎(chǔ)上，將其與東北三省實際大豆種植區(qū)域進行對比，結(jié)果發(fā)現(xiàn)松嫩平原是最適合于大豆生長的區(qū)域，黑龍江大豆種植略欠合理性，而三江平原由于日照等因素并非大豆種植的絕佳地帶。在研究結(jié)果基礎(chǔ)上，建議農(nóng)業(yè)及相關(guān)部門在保護農(nóng)民切身利益的情況下適當(dāng)調(diào)整大豆種植的空間結(jié)構(gòu);該文也為其他相關(guān)研究提供參考。

［1］凌云川.土地適宜性評價理論與方法研究.現(xiàn)代農(nóng)業(yè)科技，2007(18):191～194

［2］倪紹祥.土地類型與土地評價.北京:高等教育出版社，1999

［3］倪紹樣，陳傳康.我國土地評價研究近今進展.地理學(xué)報，1993，48(1):63～75

［4］鄭宇，胡業(yè)翠，劉彥隨，等.山東省土地適宜性空間分析及其優(yōu)化配置研究.農(nóng)業(yè)工程學(xué)報，2005，21(1):60～65

［5］邢世和，毛艷鈴，周碧青，等.福建省農(nóng)用地資源評價及其利用保護對策.福建農(nóng)業(yè)大學(xué)學(xué)報，2000，29(3):356～362

［6］謝樹春，趙玲.基于GIS的湘中紫色土丘陵地區(qū)土宜評價——以衡南縣譚子山鎮(zhèn)紫色土綜合治理實驗區(qū)為例.經(jīng)濟地理，2005，25(1):101～105

［7］白淑英，張樹文，寶音，等.遙感和GIS在土宜評價研究中的應(yīng)用——以呼和浩特武川縣為例.水土保持學(xué)報，2003，17(6):18～22

［8］佟寶全，陳才，寶音，等.北方農(nóng)林牧交錯帶土地利用空間結(jié)構(gòu)調(diào)整研究——以內(nèi)蒙古阿榮旗為例.干旱區(qū)資源與環(huán)境，2006，20(1):113～118

［9］倪宏伶，劉甫增，滕建珍，等.喀斯特地區(qū)土宜評價研究.貴州科學(xué)，1996，14(1):63～67

［10］張建明，郁光磊，璩向?qū)?寧夏南部黃土高原地區(qū)土宜評價研究.寧夏大學(xué)學(xué)報 (自然科學(xué)版)，1999，20(2):154～156

［11］張友眾，周澤福，程金花.黃土丘陵溝壑區(qū)上地適宜性評價研究——以山西省中陽縣屹針耳流域為例.水土保持學(xué)報，2003，17(1):93～97

［12］彭補拙，李春華，濮勵杰.中亞熱帶北緣青梅土宜評價方法探討.自然資源學(xué)報，1994(2):14～21

［13］華熙成.浙閩山區(qū)茶葉生產(chǎn)基地評價的探討.地理研究，1994，13(3):98～103

［14］張紅旗.GIS支持下的縣級區(qū)域柑桔土地適宜性綜合評價.資源科學(xué)，1998，20(1):62～70

［15］王桂芝.基于GIS的三亞市熱作土宜評價模型的建立.測繪信息與工程，1997(2):23～25

［16］邱炳文，池天河，王欽敏.基于GIS和多目標(biāo)評價方法的果樹適宜性評價 (英文).農(nóng)業(yè)工程學(xué)報，2005，21(6):96～100

［17］唐嘉平，劉釗.基于GIS的特色經(jīng)濟作物種植適宜性評價系統(tǒng).農(nóng)業(yè)系統(tǒng)科學(xué)與綜合研究，2002，18(1):9～12

［18］Allen，R.B.，Hewitt，A.E.，Partridge，T.R.Predicting land use suitability from vegetation and landform in depleted semiarid grassland，New Zealand.Landscape and Urban Planning，1995，32:31～42

［19］Steiner，F(xiàn).，McSherry，L.，Cohen，J.Land suitability analysis for the upper Gila River watershed.Landscape and Urban Planning，2000，50:199～214

［20］Becker，A.，Bugmann，H.(Eds.).Global Change and Mountain Regions:The Mountain Research Initiative.IGBP Report No.49，Stockholm，Sweden，2001，86pp

［21］Kalogirou，S.Expert systems and GIS:An application of land suitability evaluation computers.Environment and Urban Systems，2002，26:89～112

［22］Wu，L.X.，Sun，B.，Zhou，S.L.，et al.A new fusion technique of remote sensing images for land use/cover.Pedosphere，2004，14:187～194

［23］Liu，Y.S.，Hu，Y.C.，Peng，L.Y.Accurate quantification of grassland covers density in an alpine meadow soil based on remote sensing and GPS.Pedosphere，2005，15(6):778～783

［24］Bhatia，V.S.，Piara，S.，Wani，S.P.，et al.Analysis of potential yields and yield gaps of rainfed soybean in India using CROPGRO-Soybean model.Agriculture and Forest Meteorology，2008，148:1252 ～1265

［25］Collins，M.G.，Steiner，F(xiàn).R.，Rushman，M.J.Land-use suitability analysis in the United States:historical development and promising technological achievements.Environmental Management，2001，28:611 ～621

［26］Cutforth，H.W.，Shaykewich，C.F.1990.A temperature response functions for corn development.Agriculture and Forest Meteorology，1990，50:159～171

［27］Xu，X.X.，Gao，P.，Jiang，D.S.Fuzzy analysis on the suitability of precipitation to the agro2crop growing in Yanan City.Research of Soil and Water Conservation，2000，7:73～76

［28］Pan，T.F.，Zhang，D.R.，Zhang，W.G.，et al.Study on climatic biology of soybean in northeast China.Soybean Science，1982，1:17～28

［29］趙峰，千懷遂，焦士興.農(nóng)作物氣候適宜度模型研究——以河南省冬小麥為例.資源科學(xué)，2003，25(6):77～82

［30］He，Y.B.，Yao，Y.M.，Chen，Y.Q.，et al.Using an integrated response-function method to explore agro-climatic suitability for spring soybean growth in north China.Journal of Applied Meteorology and Climatology，2011，DOI:10.1175/2010JAMC2577.1(in press)

［31］FAO.Global Agro-ecological Assessment for Agriculture in the 21st century:Methodology and Results.Edited by Günther，F(xiàn).，Harrij，V.V.，Mahendra，S.，F(xiàn)reddy，N.2002，F(xiàn)AO，Rome，Italy

［32］楊力，張民，萬連步.大豆優(yōu)質(zhì)高效栽培.山東科學(xué)技術(shù)出版社，2006

［33］陳佑啟，姚艷敏，何英彬，等.中國區(qū)域性耕地資源變化影響評價與糧食安全預(yù)警研究.北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社，2009

［34］韓天富.大豆優(yōu)質(zhì)高產(chǎn)栽培技術(shù)指南.北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社，2005

［35］何英彬，唐華俊，楊鵬，等.不同政策情景下我國荒漠化地區(qū)土地耕作適宜性評價.農(nóng)業(yè)工程學(xué)報，2010，26(10):319～324

［36］何英彬，姚艷敏，唐華俊，等.基于尺度推繹方法的生態(tài)脆弱區(qū)耕作適宜性研究.中國農(nóng)學(xué)通報，2011，27(01):263～269

［37］石淑芹，陳佑啟，姚艷敏，等.東北三省耕地變化對糧食生產(chǎn)能力的影響評價.地理學(xué)報，2008，63(6):574～586

［38］石淑芹，陳佑啟，姚艷敏，等.東北三省耕地自然質(zhì)量和利用質(zhì)量評價.資源科學(xué)，2008，30(3):378-384

［39］周新安，張曉娟.大豆優(yōu)質(zhì)高產(chǎn)栽培技術(shù).北京:中國農(nóng)業(yè)科學(xué)技術(shù)出版社，2001

［40］黃青，唐華俊，周清波，等.東北三省主要作物種植結(jié)構(gòu)遙感提取及長勢監(jiān)測.農(nóng)業(yè)工程學(xué)報，2010，26(9):218～223