祝亞飛，劉 俊，王 陽(yáng)，齊曉波,汪建飛

(1.安徽理工大學(xué)地球與環(huán)境學(xué)院，安徽 淮南 232001；2.中國(guó)科學(xué)院地理科學(xué)與資源研究所，北京 100101；3.安徽科技學(xué)院資源與環(huán)境學(xué)院，安徽 滁州 233100)

硒(Selenium,Se)作為人體正常生長(zhǎng)發(fā)育所需要的基本微量元素，可以有效防止細(xì)胞因氧化而引起的衰老，并具有有效預(yù)防癌癥、提高人體免疫力、降低重金屬毒性等功能[1-2]。硒過(guò)量和缺乏都會(huì)導(dǎo)致人體產(chǎn)生一些疾?。何^(guò)量可導(dǎo)致智力低下、皮膚褪色等； 硒缺乏可引發(fā)克山病、 白肌病、 大骨節(jié)病等[3-5]。 由于硒無(wú)法靠人體自身合成，而是通過(guò)從外部攝取來(lái)獲取，土壤便是人體和動(dòng)植物體從外部獲取硒的最根本來(lái)源[6]， 區(qū)域內(nèi)的土壤硒資源對(duì)人體健康有重要意義。 通過(guò)研究土壤硒的地球化學(xué)特征及其影響因素，明確區(qū)域土壤硒資源的空間分布特征，為發(fā)展區(qū)域特色富硒農(nóng)業(yè)提供理論依據(jù)[7]。

近些年來(lái)，各個(gè)領(lǐng)域都在重視有關(guān)土壤硒資源開(kāi)發(fā)利用方面的研究。 國(guó)內(nèi)的研究方向主要集中于硒的地球化學(xué)特征、 硒和土壤其他理化性質(zhì)以及重金屬元素之間的拮抗協(xié)同機(jī)制等方面[8-9]；國(guó)外的研究方向則更側(cè)重于探究土壤-作物系統(tǒng)中的不同形態(tài)硒的遷移轉(zhuǎn)化機(jī)制和規(guī)律、硒元素在區(qū)域生態(tài)系統(tǒng)中的各種循環(huán)機(jī)理以及硒的生物有效性等[10]。目前的研究大多針對(duì)富硒土壤，而對(duì)區(qū)域硒資源在足硒層次以及非富硒區(qū)域如何發(fā)展富硒農(nóng)業(yè)的研究較少[11-13]。 因此， 本文基于生態(tài)位理論，通過(guò)研究硒和其他理化性質(zhì)在土壤中的富集特征，分別進(jìn)行生態(tài)位適宜性評(píng)價(jià)，最后利用地統(tǒng)計(jì)學(xué)方法，結(jié)合研究區(qū)生態(tài)敏感性分布結(jié)果，將全區(qū)進(jìn)行硒資源綜合利用分區(qū)[14-16]。最后結(jié)合綜合利用分區(qū)結(jié)果，進(jìn)一步探究如何科學(xué)地開(kāi)發(fā)和利用區(qū)域土壤硒資源，為高效利用硒的經(jīng)濟(jì)價(jià)值和發(fā)展特色富硒農(nóng)業(yè)提出科學(xué)建議[17]。

1 材料與方法

1.1 研究區(qū)概況及樣點(diǎn)布設(shè)

研究區(qū)(39°26'～39°50'N，116°13'～116°43'E)位于北京市南郊， 地理位置為華北平原的東北部， 整體區(qū)域西北較高、 東南較低， 總體地勢(shì)較為平坦， 海拔高度在15～45m之間，自然坡度為0.08%至0.01%，基礎(chǔ)農(nóng)業(yè)條件較好，適合發(fā)展各類(lèi)特色農(nóng)業(yè)。研究區(qū)土壤母質(zhì)主要以沖積物為主，其次是風(fēng)砂沖積物，洪沖積物最少；質(zhì)地構(gòu)型以海綿型為主，夾層型最少。

圖1 研究區(qū)位置和樣點(diǎn)分布

1.2 數(shù)據(jù)獲取

按1km×1km網(wǎng)格采樣，取表層土(0～20cm)，共采集167份土樣。本文利用SPSS 22.0、Excel 2019、ArcGIS 10.3和ENVI 5.3進(jìn)行數(shù)據(jù)分析和相關(guān)圖形繪制。

1.3 研究方法

1)單因子生態(tài)位適宜性分析

a. 主導(dǎo)因子

(1)

式中：Ni為評(píng)價(jià)因子i的生態(tài)位適宜度(i∈[1，n])；Xi為現(xiàn)實(shí)生態(tài)位；Diopt為最適生態(tài)位；Dimin為最小生態(tài)位。

b. 區(qū)位因子

(2)

式中：Ni為生態(tài)位適宜度(i∈[1，n])；Xi為現(xiàn)實(shí)生態(tài)位；Diopt為最適生態(tài)位；Dimin為最小生態(tài)位；Dimax為最大生態(tài)位。

c. 限制因子

(3)

式中：Ni為生態(tài)位適宜度(i∈[1，n])；Xi為現(xiàn)實(shí)生態(tài)位；Diopt為最適生態(tài)位；Dimin為最小生態(tài)位；Dimax為最大生態(tài)位。

2)生態(tài)位適宜性綜合分析

a. 生態(tài)位適宜性綜合指數(shù)分析

(4)

式中：N為生態(tài)位適宜性綜合指數(shù)；Ni為生態(tài)位適宜度(i∈[1，n])；n為評(píng)價(jià)因子個(gè)數(shù)。

b. 生態(tài)位適宜性綜合分析

(5)

式中：S為生態(tài)位適宜性加權(quán)求和綜合指數(shù)；Ni為生態(tài)位適宜性綜合指數(shù)；Wi為生態(tài)位適宜性權(quán)重；n為生態(tài)位維數(shù)。

結(jié)合生態(tài)位適宜性分析結(jié)果，通過(guò)GIS軟件將適宜性綜合評(píng)價(jià)結(jié)果分為：高適宜區(qū)、中適宜區(qū)、低適宜區(qū)和不適宜區(qū)[18]。

2 結(jié)果與討論

2.1 土壤硒分布特征分析

1)土壤硒描述性統(tǒng)計(jì)分析 對(duì)研究區(qū)土壤硒含量、有機(jī)質(zhì)和pH進(jìn)行描述性統(tǒng)計(jì)，研究區(qū)表層土壤硒含量介于0.09～0.41mg·kg-1之間，其算術(shù)平均值為0.19mg·kg-1，相較于我國(guó)平均土壤硒含量0.29mg·kg-1偏低；土壤有機(jī)質(zhì)含量介于2.84～39.60g·kg-1之間，其算術(shù)平均值為15.36g·kg-1；土壤pH值介于7.65～8.89之間，其算術(shù)平均值為8.41，整體偏堿性。研究區(qū)土壤硒和有機(jī)質(zhì)均為中度變異(37.44%、43.73%)土壤pH值呈弱變異性(3.33%)，這可能與研究區(qū)內(nèi)不同區(qū)域的種植模式、土地利用類(lèi)型不同及環(huán)境因素差異有關(guān)。

2)土壤硒空間分布特征分析 空間自相關(guān)性是探究相關(guān)元素之間是否存在顯著關(guān)系的基礎(chǔ)， 是進(jìn)行地統(tǒng)計(jì)數(shù)學(xué)處理的重要依據(jù)。圖2為全局莫蘭指數(shù)散點(diǎn)圖和局部莫蘭指數(shù)聚類(lèi)圖，研究區(qū)土壤硒全局莫蘭指數(shù)均值為0.848；局部莫蘭指數(shù)聚類(lèi)圖表表明，全區(qū)高-高聚類(lèi)主要分布在中北部，低-低聚類(lèi)主要分布在西部和東部地區(qū)。利用Geoda軟件檢驗(yàn)全局莫蘭指數(shù)特征：Z=3.29>1.96；P<0.01。符合空間自相關(guān)分析的要求。

(a) 全局莫蘭指數(shù)散點(diǎn)圖 (b) 局部莫蘭指數(shù)聚類(lèi)圖

研究區(qū)土壤硒與土壤其他理化性質(zhì)的相關(guān)性分析如圖3所示。結(jié)果表明，土壤硒與有機(jī)質(zhì)、pH和重金屬Hg等呈不同程度的相關(guān)。硒與pH呈正相關(guān)(P<0.05)。研究區(qū)表層土壤pH平均值為8.41， 為中、 堿性土壤。 硒與有機(jī)質(zhì)呈顯著正相關(guān)(P<0.05)，表明有機(jī)質(zhì)對(duì)硒具有一定的吸附和固定作用，這與文獻(xiàn)[19]研究認(rèn)為硒是典型親生物元素，在富含有機(jī)質(zhì)的土壤環(huán)境中易富集的結(jié)論一致，但有關(guān)研究表明，土壤中大量N、P的存在不利于土壤硒向植物遷移轉(zhuǎn)化。硒與Hg呈正相關(guān)(P<0.05)，這與它們同屬于親硫元素，具有相似的地球化學(xué)性質(zhì)有關(guān)。硒與As、Ni 呈正相關(guān)(P<0.05)，但相關(guān)系數(shù)較小，而與Cd無(wú)相關(guān)關(guān)系，表明土壤中硒的自然來(lái)源與土壤母質(zhì)有關(guān)。

圖3 硒與土壤理化性質(zhì)的關(guān)系

為進(jìn)一步揭示土壤硒與理化性質(zhì)及環(huán)境要素之間的關(guān)系，利用方差分解探究各指標(biāo)對(duì)硒的影響，結(jié)果顯示有機(jī)質(zhì)、pH和Hg對(duì)土壤硒積累的貢獻(xiàn)作用較大，與相關(guān)性分析結(jié)果相同。而根據(jù)以往研究，在酸性和堿性環(huán)境下表層土壤硒主要與pH呈不同程度負(fù)相關(guān)關(guān)系[20]；在中性環(huán)境下，呈較弱正相關(guān)關(guān)系[21]。本文結(jié)果可能與樣點(diǎn)分布以農(nóng)田為主，受土地利用方式及種植管理模式影響有關(guān)[22]；Cr、Ni和Cd對(duì)土壤硒積累無(wú)貢獻(xiàn)作用，說(shuō)明土壤硒含量受Cr、Ni和Cd的影響不大。

2.2 生適宜性評(píng)價(jià)

1)土壤硒生態(tài)位適宜性評(píng)價(jià)體系構(gòu)建 土壤硒資源開(kāi)發(fā)利用受多方面因子的影響，如天然的成土母質(zhì)會(huì)影響土壤硒含量的高低，區(qū)域的基礎(chǔ)農(nóng)業(yè)設(shè)施和條件會(huì)影響富硒農(nóng)業(yè)的發(fā)展，還要考慮發(fā)展富硒產(chǎn)業(yè)所投入的經(jīng)濟(jì)占比及經(jīng)濟(jì)收益。因此本文從主導(dǎo)因子、區(qū)位因子和限制因子3個(gè)方面建立指標(biāo)體系，通過(guò)熵值法確定各評(píng)價(jià)因子權(quán)重，如表1所示。本文富硒分級(jí)標(biāo)準(zhǔn)參考文獻(xiàn)[23-24]劃分。

2)主導(dǎo)因子適宜性分析 本文的主導(dǎo)因子區(qū)域土壤硒的含量，也是發(fā)展當(dāng)?shù)靥厣晃r(nóng)業(yè)的基礎(chǔ)。根據(jù)評(píng)價(jià)體系標(biāo)準(zhǔn)對(duì)因子進(jìn)行重分類(lèi)賦值，得出主導(dǎo)因子適宜性評(píng)價(jià)結(jié)果如圖4所示。按照評(píng)價(jià)因子現(xiàn)實(shí)生態(tài)位等級(jí)將研究區(qū)分為高適宜(Ⅰ級(jí))、中適宜(Ⅱ級(jí))、低適宜(Ⅲ級(jí))和不適宜(Ⅳ級(jí))4種類(lèi)型區(qū)?？梢钥闯?，研究區(qū)表層土壤硒評(píng)價(jià)以低適宜(足硒)為主，主要集中在研究區(qū)西北地區(qū)和南部地區(qū)；高適宜(富硒)區(qū)主要分布在研究區(qū)中部，表層土壤硒含量較高，適合發(fā)展富硒產(chǎn)業(yè)；中適宜區(qū)(潛在富硒)較少，主要分布在中北部地區(qū)；不適宜區(qū)(缺硒)主要分布在研究區(qū)西部和東部地區(qū)，由于土壤硒含量偏低，不適合發(fā)展富硒產(chǎn)業(yè)。

表2 土壤硒資源開(kāi)發(fā)利用適宜性評(píng)價(jià)指標(biāo)體系及權(quán)重

圖4 主導(dǎo)因子適宜性評(píng)價(jià)

3)區(qū)位因子適宜性分析 本文選擇區(qū)域最典型區(qū)位因子地形坡度、土壤肥力和土地利用類(lèi)型進(jìn)行評(píng)價(jià)，根據(jù)評(píng)價(jià)體系標(biāo)準(zhǔn)對(duì)因子進(jìn)行重分類(lèi)賦值，得出相應(yīng)適宜性評(píng)價(jià)結(jié)果(見(jiàn)圖5)。

圖5 有機(jī)質(zhì)、坡度、土地利用因子適宜性評(píng)價(jià)

土壤有機(jī)質(zhì)適宜性以中適宜(Ⅱ級(jí))為主，主要集中在中部和南部地區(qū)；高適宜(Ⅰ級(jí))區(qū)主要集中在研究區(qū)西部和東南部，這部分地區(qū)以菜地和糧地為主，在長(zhǎng)期耕作與施肥的影響下土壤有機(jī)質(zhì)含量較高，評(píng)價(jià)結(jié)果與實(shí)際相符，適合農(nóng)業(yè)生產(chǎn)；低適宜(Ⅲ級(jí))和不適宜區(qū)(Ⅳ級(jí))主要集中在研究區(qū)北部和東北部，這部分地區(qū)城鎮(zhèn)化程度較高，農(nóng)業(yè)發(fā)展程度低，土壤有機(jī)質(zhì)含量偏低，不適合農(nóng)業(yè)生產(chǎn)。

地形坡度適宜性評(píng)價(jià)以高適宜(Ⅰ級(jí))為主，由于研究區(qū)處于華北平原，地形坡度變化較小，研究區(qū)南部和中部及部分北部地區(qū)均為高適宜區(qū)，基礎(chǔ)農(nóng)業(yè)設(shè)施建設(shè)較為完善、交通方便，適合發(fā)展農(nóng)業(yè)；低適宜(Ⅲ級(jí))和不適宜區(qū)(Ⅳ級(jí))主要在西北部和中西部地區(qū)，這類(lèi)地區(qū)地勢(shì)較高，農(nóng)業(yè)設(shè)施建設(shè)較為薄弱、交通不便，不適合農(nóng)業(yè)生產(chǎn)。

土地利用類(lèi)型適宜性評(píng)價(jià)以高適宜為主(Ⅰ級(jí))，主要集中在中部和南部，研究區(qū)耕地較多，且大部分地區(qū)為平原，較適合農(nóng)業(yè)生產(chǎn)；低適宜(Ⅲ級(jí))和不適宜區(qū)(Ⅳ級(jí))主要集中在北部，這部分地區(qū)主要為城鄉(xiāng)、工礦和居民用地，結(jié)合城市規(guī)劃發(fā)展考慮，不適合發(fā)展農(nóng)業(yè)。

4)限制因子適宜性分析 本文采用生態(tài)敏感性作為生態(tài)保護(hù)建設(shè)的限制因子，生態(tài)環(huán)境敏感性是指生態(tài)系統(tǒng)對(duì)區(qū)域內(nèi)自然和人類(lèi)活動(dòng)干擾的敏感程度。

生態(tài)敏感性評(píng)價(jià)結(jié)果如圖6所示，利用GIS技術(shù)對(duì)DEM數(shù)據(jù)進(jìn)行高程、坡度、坡向和水系數(shù)據(jù)提取，疊加分析將研究區(qū)生態(tài)敏感性分為Ⅰ級(jí)、Ⅱ級(jí)、Ⅲ級(jí)和Ⅳ級(jí)。Ⅰ級(jí)為高敏感區(qū)，主要為水源地和生態(tài)保護(hù)建設(shè)區(qū)，不適合農(nóng)業(yè)生產(chǎn)；Ⅱ級(jí)為敏感區(qū)，距離水源地和生態(tài)保護(hù)建設(shè)區(qū)較近，存在人為活動(dòng)影響到生態(tài)環(huán)境的可能；Ⅲ級(jí)為較敏感區(qū)，正常的人為活動(dòng)和農(nóng)業(yè)生產(chǎn)不會(huì)影響到生態(tài)環(huán)境；Ⅳ級(jí)為低敏感區(qū)，人為活動(dòng)對(duì)生態(tài)環(huán)境影響較小，適合農(nóng)業(yè)生產(chǎn)。研究區(qū)生態(tài)敏感性以低敏感性為主，大部分地區(qū)適合進(jìn)行人為活動(dòng)和農(nóng)業(yè)生產(chǎn)。

圖6 生態(tài)敏感性分級(jí)

5)生態(tài)位適宜性綜合分析 綜合疊加計(jì)算得出生態(tài)位適宜性評(píng)價(jià)結(jié)果如圖7所示，研究區(qū)分為高適宜、中適宜、低適宜和不適宜4個(gè)類(lèi)別。高適宜區(qū)主要分布在中部和西南部，適宜發(fā)展富硒產(chǎn)業(yè)；中適宜區(qū)主要在南部和東部地區(qū)，為潛在富硒產(chǎn)業(yè)發(fā)展區(qū)，可作為土壤硒資源利用后備區(qū)；低適宜區(qū)主要集中在水源地和生態(tài)保護(hù)建設(shè)地周?chē)鶕?jù)生態(tài)保護(hù)原則，此類(lèi)地區(qū)不適合發(fā)展富硒產(chǎn)業(yè)；不適宜區(qū)主要在西北部，這部分地區(qū)地形坡度偏高、交通不便、基礎(chǔ)農(nóng)田設(shè)施建設(shè)滯后，不適合發(fā)展農(nóng)業(yè)以及富硒產(chǎn)業(yè)。

圖7 生態(tài)適宜性評(píng)價(jià)

2.3 土壤硒資源綜合利用分區(qū)

基于生態(tài)位適宜性評(píng)價(jià)結(jié)果，結(jié)合農(nóng)用地等級(jí)，得出區(qū)域土壤硒資源綜合利用分區(qū)結(jié)果。為科學(xué)地為區(qū)域特色富硒產(chǎn)業(yè)發(fā)展提出建議，將研究區(qū)劃分為富硒產(chǎn)業(yè)發(fā)展區(qū)、重點(diǎn)農(nóng)業(yè)生產(chǎn)區(qū)和綜合產(chǎn)業(yè)發(fā)展區(qū)(見(jiàn)圖8)。

圖8 土壤硒資源綜合利用分區(qū)

富硒產(chǎn)業(yè)發(fā)展區(qū)主要在研究區(qū)中部，這部分區(qū)域的農(nóng)業(yè)基礎(chǔ)設(shè)施較好，耕地質(zhì)量等級(jí)較高，土壤硒資源較為豐富，特色富硒農(nóng)業(yè)發(fā)展前景廣闊，可充分發(fā)揮農(nóng)業(yè)基礎(chǔ)優(yōu)勢(shì)，以重點(diǎn)開(kāi)發(fā)綠色富硒農(nóng)產(chǎn)品為核心，構(gòu)建富硒小麥、西甜瓜、蔬菜和玉米等特色產(chǎn)品農(nóng)業(yè)產(chǎn)業(yè)體系，走特色富硒生態(tài)農(nóng)業(yè)發(fā)展道路。具體開(kāi)發(fā)過(guò)程中，應(yīng)調(diào)查清楚各類(lèi)農(nóng)產(chǎn)品含硒值，掌握具體地塊土壤硒分布情況，堅(jiān)持以當(dāng)?shù)靥厣晃r(nóng)產(chǎn)品為基礎(chǔ)，有針對(duì)性地打造富硒農(nóng)產(chǎn)品。

重點(diǎn)農(nóng)業(yè)生產(chǎn)區(qū)主要集中在研究區(qū)中部和南部部分區(qū)域，這部分區(qū)域產(chǎn)業(yè)基礎(chǔ)好、交通便利，基礎(chǔ)農(nóng)業(yè)發(fā)展較好，土壤硒資源為潛在富硒和足硒。應(yīng)堅(jiān)持基礎(chǔ)農(nóng)業(yè)發(fā)展為主導(dǎo)，并加強(qiáng)人工富硒相關(guān)技術(shù)的研究，通過(guò)外源富硒技術(shù)加強(qiáng)和改良農(nóng)產(chǎn)品籽粒果實(shí)硒含量，綜合發(fā)展現(xiàn)代基礎(chǔ)農(nóng)業(yè)和特色富硒產(chǎn)業(yè)[25]。

綜合產(chǎn)業(yè)發(fā)展區(qū)主要在研究區(qū)中西部，該部分區(qū)域耕地質(zhì)量等級(jí)相對(duì)較低，應(yīng)有針對(duì)性地發(fā)展特色產(chǎn)業(yè)，加強(qiáng)農(nóng)田監(jiān)管和保護(hù)，改良農(nóng)田土壤障礙因子，優(yōu)化土地生產(chǎn)要素投入結(jié)構(gòu)，提高耕地質(zhì)量等級(jí)，提升生產(chǎn)能力和生態(tài)效益。結(jié)合當(dāng)?shù)赝临|(zhì)狀況，優(yōu)化種植模式和農(nóng)業(yè)種植結(jié)構(gòu)，通過(guò)合理施用復(fù)合肥、改進(jìn)耕作措施、集中安置管理糧田和農(nóng)產(chǎn)品品種改良等措施挖掘農(nóng)田發(fā)展?jié)摿26]。

3 結(jié)論

本文基于生態(tài)位理論對(duì)土壤硒資源開(kāi)發(fā)利用進(jìn)行適宜性評(píng)價(jià)，通過(guò)分區(qū)結(jié)果得出以下結(jié)論：

(1)研究區(qū)土壤硒含量在0.09～0.41mg·kg-1之間，均值為為0.19mg·kg-1，總體屬于足硒，全區(qū)土壤硒全局莫蘭指數(shù)均值為0.848，存在顯著的空間正相關(guān)性；

(2)全區(qū)表層土壤硒低適宜(足硒)主要集中在研究區(qū)西北地區(qū)和南部地區(qū)，高適宜(富硒)區(qū)主要分布在研究區(qū)中部，表層土壤硒含量較高；土壤有機(jī)質(zhì)以中適宜(Ⅱ級(jí))為主，主要集中在中部和南部地區(qū)，地形坡度和土地利用類(lèi)型主要以高適宜(Ⅰ級(jí))為主，生態(tài)敏感性以低敏感性為主，地理?xiàng)l件與農(nóng)業(yè)生產(chǎn)基礎(chǔ)條件較好，適合構(gòu)建以農(nóng)業(yè)產(chǎn)業(yè)為主導(dǎo)的生產(chǎn)體系；

(3)本文將生態(tài)位適宜性評(píng)價(jià)結(jié)果和農(nóng)用地土地利用等別相結(jié)合，將研究區(qū)劃分為富硒產(chǎn)業(yè)發(fā)展區(qū)、重點(diǎn)農(nóng)業(yè)生產(chǎn)區(qū)和綜合產(chǎn)業(yè)發(fā)展區(qū)。結(jié)合不同類(lèi)型的發(fā)展區(qū)實(shí)際狀況給出相應(yīng)的發(fā)展建議，富硒產(chǎn)業(yè)發(fā)展區(qū)應(yīng)充分發(fā)揮農(nóng)業(yè)基礎(chǔ)優(yōu)勢(shì)，重點(diǎn)開(kāi)發(fā)綠色富硒農(nóng)產(chǎn)品；重點(diǎn)農(nóng)業(yè)生產(chǎn)區(qū)應(yīng)堅(jiān)持基礎(chǔ)農(nóng)業(yè)發(fā)展為主導(dǎo)，并加強(qiáng)人工富硒相關(guān)技術(shù)的研究；綜合產(chǎn)業(yè)發(fā)展區(qū)應(yīng)加強(qiáng)農(nóng)田監(jiān)管和保護(hù)，提高耕地質(zhì)量等級(jí)。