雄安新區富硒土地資源分布特征及開發利用評價

周亞龍,郭志娟*,王喬林,王成文,宋云濤,劉 飛

雄安新區富硒土地資源分布特征及開發利用評價

周亞龍1,2,3,郭志娟1,2,3*,王喬林1,2,3,王成文1,2,3,宋云濤1,2,3,劉 飛1,2,3

(1.中國地質科學院地球物理地球化學勘查研究所,河北 廊坊 065000；2.中國地質調查局土地質量地球化學調查評價研究中心,河北 廊坊 065000；3.中國地質科學院地球表層碳－汞地球化學循環重點實驗室,河北 廊坊 065000)

采集雄安新區表層土壤樣品(0～20cm)、大宗農作物(小麥籽實)及根系土樣品,利用GIS空間分析、相關分析等方法對土壤硒含量分布、富集特征及影響因素進行分析,基于生態位適宜性評價理論模型開展富硒土地資源開發利用適宜性評價.結果表明:新區表層富硒土地16993hm2,土壤中硒賦存形態以強有機結合態為主,其次為殘渣態和腐殖酸結合態,三者占土壤硒總量的87%以上,生物可利用態硒含量所占比例相對較小(水溶態硒占比4.44%、離子交換態硒占比3.43%).小麥籽實樣品富硒率28.45%,籽實中硒含量與根系土壤中硒含量具有顯著相關性(=0.77,<0.01).土壤硒富集主要受第四紀沉積地貌環境的影響,此外灌溉水、化肥和大氣干濕沉降(燃煤降塵)等外源輸入也影響貧化富集.土壤硒含量與沉積地貌環境“標志性”土壤指標Al2O3、Fe2O3、MgO和K2O含量,土壤有機碳(Corg.)、陽離子交換量(CEC)、黏粒和粉粒含量等具有顯著正相關關系,與SiO2、Na2O含量和土壤砂粒含量等具有顯著負相關關系(<0.01).基于生態位評價理論模型,新區富硒土地資源開發利用適宜性評價結果表明:高適宜區(Ⅰ級)和中適宜區(Ⅱ級)主要位于新區西北容城縣南張鎮和小里鎮、安新縣寨里鄉和老河頭鎮等地,土地面積為14397.35hm2,占新區規劃總面積的8.13%.

硒；土壤；根系土；農作物；適宜性評價；雄安新區

硒是生命必須微量元素之一,與人體健康息息相關[1-2].人體硒的營養水平主要來源于日常飲食中硒的攝入量,土壤中的硒是植物硒的主要來源[3-4].因此,近年來富硒土地資源的開發利用研究受到廣泛關注[5-7].目前,我國已完成的101.7萬km2耕地質量1:25萬地球化學調查顯示,處于低硒含量水平(Se< 0.2mg/kg)的耕地比例高達44.9%.富硒耕地(Se> 0.4mg/kg)比例僅10%[8].基于我國耕地普遍缺硒的現狀,開展天然富硒土地資源評價,促進天然富硒土地資源高效利用,對于提升我國人群硒攝入量、構建健康中國具有重要意義.

如何科學高效地開發利用富硒土地資源,是富硒產業發展的關鍵前提[9].本文利用雄安新區土地質量地球化學調查的土壤與農作物中Se等元素含量數據,對新區富硒土地資源分布特征、土壤硒來源與影響因素、土壤硒形態與生物有效性的分析,基于生態位適宜性評價理論模型開展富硒土地資源開發利用適宜性評價,以期為新區富硒土地資源科學高效利用提供依據,有效支撐“綠色生態宜居”的雄安新區建設,為新區國土空間規劃提供重要參考.

1 材料與方法

1.1 研究區概況

研究區地處華北平原,其地理坐標38°43′～ 39°10′N,115°38′～116°20′E,屬暖溫帶季風型半濕潤半干旱氣候.研究區地勢平坦,地面高程5～26m,坡降0.2‰～0.7‰[10].以低海拔平原和洼地地貌為主,區內第四系地層覆蓋,以沖洪積、沖湖積和沖積為主,主要由砂、黏土和礫石組成,層厚348～437m[11].研究范圍覆蓋新區1770km2,其土地利用類型主要包括耕地、園地、林地、草地、交通運輸用地、水域及水利設施用地和城鎮村住宅用地及工礦用地8種類型(圖1),其中耕地面積占比53.33%[12].土壤類型以潮土為主,容城北部零散分布少量褐土,白洋淀周邊地區發育少量沼澤土.

1.2 樣品采集與測試

1.2.1 樣品采集 在系統分析研究區土地利用類型數據的基礎上,根據研究區耕地、園地、林地、草地和建設用地等土地利用類型圖斑面積情況,采用差異化布點的方式(表1)分類布設樣點[13],系統采集表層土壤樣品37818件.采集大宗農作物(小麥籽實)及對應根系土樣品136件,其中按照4km×4km網格采樣密度在全區采集116件,在富硒土壤區采集20件.采集灌溉水樣品101件.根據研究區內不同沉積地貌單元布置土壤垂直剖面8條,其中沖洪積平原亞區布置5條(PM01、PM02、PM03、PM04和PM07),沖積平原亞區布置3條(PM05、PM06和PM08).土壤垂直剖面樣品采集以1個/10cm的密度等間距采樣,采樣深度0～200cm.野外樣品采集加工方法執行DZ/T 0295-2016規范[14].

圖1 研究區地理位置及其土地利用類型

表1 差異化分類布設樣點控制數

1.2.2 樣品測試 土壤樣品和灌溉水樣Se元素分析測試由河南省巖石礦物測試中心(BAF-2000原子熒光光度計)和中國地質科學院地球物理地球化學勘查研究所中心實驗室(XGY-2020A原子熒光光譜儀)采用氫化物-原子熒光光譜法(AFS)測定,土壤樣品分析測試檢出限為0.01mg/kg,灌溉水樣分析測試檢出限為0.0001mg/L,農作物樣品分析測試檢出限為0.005mg/kg.根系土壤Se形態分析采用逐步浸提方法依次提取土壤硒水溶態、離子交換態、碳酸鹽結合態、腐殖酸結合態、鐵錳氧化物結合態、強有機結合態和殘渣態等.土壤理化性質等相關指標: 采用電感耦合等離子體發射光譜法(ICP-OES)測定MgO、K2O和Na2O,X射線熒光光譜法(XRF)測定Al2O3、Fe2O3和SiO2,容重法(VOL)測定土壤有機碳(Corg.),離子選擇性電極法(ISE)測定pH值.依據LY/T 1243-1999標準[15],采用容量法(VOL)測定CEC,依據NY/T 1121.3-2006標準[16],采用比重計法測定質地.樣品分析測試采用國家一級標準土樣監控分析測試的準確度,采用重復樣監控分析測試的精密度.所有分析測試指標一級標準物質合格率為100%.土壤重復性樣品檢驗Se相對雙差合格率99.65%,異常點樣品檢驗Se相對雙差合格率96.89%,樣品分析測試質量均滿足相關規范要求,分析數據可靠.

1.3 評價方法與標準

1.3.1 富硒評價標準 根據DZ/T 0380-2021標準[17],中酸性土壤(pH￡7.5)中土壤富硒標準閾值為30.40mg/kg,堿性土壤(pH>7.5)中土壤富硒標準閾值為30.30mg/kg.目前沒有小麥籽實富硒國家,本次小麥籽實富硒標準參照已有的相關大宗農作物中富硒水稻籽實硒含量的規定[18],富硒小麥硒含量標準為0.04～0.30mg/kg.

硒輸入量=++(1)

其中:為農田土壤硒的輸入通量,g/(hm2·a);為大氣干濕沉降對農田土壤硒的輸入通量,g/(hm2·a);為灌溉水對農田土壤硒的輸入通量,g/(hm2·a);為化肥對農田土壤硒的輸入通量,g/(hm2·a).

灌溉水中硒輸入通量=?.其中,為灌溉水樣中硒的含量(mg/L),為每年單位面積農田灌溉用水量(m3/(hm2·a)).大氣干濕沉降和化肥輸入通量引用相關研究結果.

2 結果與討論

2.1 土壤硒含量與空間分布

如圖2所示,研究區土壤Se元素含量區間為0.01～5.70mg/kg,平均值為0.22mg/kg,中位數為0.20mg/kg,標準差為0.20mg/kg,元素變異系數為93.94%,屬高度變異[19],其區域分布具有明顯差異.雄安新區表層土壤Se 含量背景值(剔除3倍離差以后的數據平均值0.20mg/kg)略低于全國土壤背景值(0.29mg/kg)[20],略高于河北省海河平原北部土壤背景值(0.18mg/kg)[21].

圖2 雄安新區表層土壤硒(Se)含量統計直方圖

研究區圈定表層富硒土地面積為16993hm2,占總面積的11.05%,其中無重金屬污染的富硒耕地9380hm2.主要分布在容城縣西部的南張鎮、小里鎮,安新縣寨里鄉、老河頭鎮等地(圖3).此外,研究區還存在5.54%的土壤硒缺乏地塊.

圖3 土壤硒(Se)含量等級

2.2 土壤硒賦存形態特征

土壤中硒的賦存形態特征是研究土壤-作物系統中硒含量關系的關鍵因素之一,其生物有效性是評價土壤硒有效性的重要指標之一[22].研究區小麥籽實根系土壤中Se賦存形態均以強有機結合態為主(圖4),其次為殘渣態和腐殖酸結合態,所占比例均值分別為45.60%、30.11%和11.29%,其他4種形態的硒含量均值所占比例相對較小,水溶態硒占比4.44%,離子交換態硒占比3.43%,碳酸鹽態硒占比2.96%,鐵錳結合態態硒占比2.17%.

2.3 農作物硒含量與硒富集評價

雄安新區116件小麥籽實樣品中Se含量介于0.001～1.152mg/kg,平均值為0.074mg/kg,達到富硒小麥標準的樣品33件,研究區小麥籽實富硒率達28.45%,此外5件小麥籽實樣品含量超過富硒標準上限(0.30mg/kg).新區富硒小麥籽實主要位于研究區容城西部、安新西北部、安新老河頭鎮等地區.

富硒土壤區20件小麥籽實樣品中Se含量區間0.076～1.858mg/kg,平均值為0.555mg/kg,顯著高于全區116件小麥籽實樣品Se含量均值.富硒土壤區小麥籽實樣品富硒率100%,此外11件小麥籽實樣品含量超過富硒標準上限.

圖4 根系土壤中硒的賦存形態分布特征

2.4 土壤富硒成因與影響因素探討

土壤中硒元素含量主要受地質背景(成土母質、地貌、土壤類型等)自然因素[2,23-24]和人類活動(燃料燃燒、金屬冶煉、化肥等)[7,25]的影響.

王慧卿介紹，這項工作首先是根據醫院實際情況選擇合適的崗位評價方法——要素計點法，接著初步擬訂崗位模型的要素框架，并篩選崗位評價模型的評價要素和評價維度。項目組最終選定基于七因素的崗位評價模型（圖1），每個因素又設有2～3個維度，每個要素具有不同評分權重。以業務支撐要素為例，在影響層次和貢獻度兩個維度交叉聯系的基礎上，項目組確定了從68分至400分的不同的評分標準。

2.4.1 地質因素控制富集因子 研究區位于華北平原海河流域,地表出露地層為第四系沖洪積松散地層,以沖積、洪積和湖積為主,地層厚348～437m,無基巖出露.研究區沉積地貌特征反映其土壤中硒的富集不是來源于基底巖石原地風化成因.研究區土壤垂直剖面硒的分布模式也證實了該特征.不同地貌單元垂直剖面中土壤硒的含量隨深度變化如圖5所示,各垂直剖面中土壤硒元素含量的分布模式總體上均表現為線性表聚型,即從表層到深部硒含量逐漸降低,一定程度上反映研究區土壤硒含量僅在表層富集,該分布模式可能反映其來源可能未受成土母質的影響,有可能受到外源輸入的影響.

質地是土壤一種重要的自然屬性,可反映母質來源及成土過程中的地球化學特征[26].研究區土壤垂直剖面中硒含量與土壤質地中粉粒和黏粒含量均呈現顯著正相關(= 0.48和0.60,<0.01),而與土壤質地中的砂粒含量呈現顯著負相關(=-0.59,< 0.01),反映土壤黏粒對Se具有較強的吸附富集作用,土壤中黏粒含量越多,硒含量越高.

除地質背景影響外,土壤硒含量還與成土過程中硒的次生富集作用有關.土壤物理、化學和生物作用對硒的吸附固定、活化遷移有重要影響,決定了土壤硒的富集貧化及生物有效性[7].土壤中SiO2、Al2O3、Fe2O3、MgO、CaO、Na2O和K2O等指標含量組合特征對地貌環境具有“指示”或“標志性”意義[27-28].相關分析顯示(表2),研究區表層土壤硒含量與地貌環境“標志性”指標Al2O3、Fe2O3、MgO和K2O含量具有顯著正相關關系,與SiO2和Na2O具有顯著負相關關系(<0.01).表層土壤硒含量與土壤理化性質中的有機碳(Corg.)、陽離子交換量(CEC)和土壤黏粒和粉粒含量等具有顯著正相關關系,與pH值和土壤砂粒含量具有顯著負相關關系(<0.01).

研究區土壤硒含量與地貌環境“標志性”指標、土壤理化性質指標的顯著相關性反映土壤硒的富集貧化與土壤黏土礦物、有機質等吸附載體密切相關.在滯水環境形成的湖相或沖湖積沉積物上發育的土壤中,由于土壤偏黏性,富含黏土礦物、有機質等吸附載體,Al2O3和Fe2O3等易發生富集[28];在沖洪積平原故河道高地、扇上平地或緩斜地等沉積環境中,長石類礦物發育,土壤以砂土、粉砂土為主,Na2O和SiO2呈現高值特征.

圖5 土壤垂直剖面中硒含量分布與土壤理化性質關系

表2 土壤硒含量與礦物含量和土壤理化性質的相關系數

注:*在0.01水平(雙側)上顯著相關.

2.4.2 人類活動影響 河北平原區土壤硒異常成因研究證實,灌溉、化肥和大氣干濕沉降是其主要輸入途徑,其中燃煤、有色金屬冶煉等產生的大氣降塵是形成土壤硒異常的主要因素[29-31].

雄安新區101件灌溉水樣中Se含量均值為0.0029mg/L.根據2019年度河北省水資源公報,河北省農田灌溉用水總量為102.91億m3,國家統計局數據顯示,2019年河北省全省有效灌溉面積448.22萬hm2,河北省農田灌溉年需水量約為2295.99m3/ (hm2·a),新區農田區灌溉水Se輸入量約為6.66g/ (hm2·a),略高于保定市農田灌溉水Se輸入量(5.17g/ (hm2·a))[31]和河北省平原區農田灌溉水Se輸入量(3.17g/(hm2·a))[30].新區4件大氣干濕沉降樣品中Se輸入通量均值為3.18mg/kg[32],略低于保定市大氣干濕沉降Se輸入量(5.15g/(hm2·a))[31]和河北省平原區大氣干濕沉降Se輸入通量均值(5.88g/(hm2·a))[30],大氣干濕沉降硒輸入通量的降低可能與近年來新區冬季燃煤減少,部分冶煉企業關停等人類活動有關.根據保定市化肥樣品中Se含量數據和單位面積年施肥量,其化肥Se元素輸入量約為1.4g/(hm2·a)[32],高于河北省平原區Se化肥輸入通量均值(0.81g/ (hm2·a))[30].研究區土壤Se外源輸入主要為灌溉水和大氣干濕沉降輸入,其合計占比達87.54%,化肥輸入僅占比12.46%.

2.4.3 農作物富硒成因與影響因素 土壤-農作物系統中Se的吸收、運移與土壤Se含量水平、存在形態、生物有效性及土壤理化性質(pH值、氧化還原電位、有機質)等眾多因素有關[33].雄安新區小麥籽實Se生物富集系數(農作物Se含量/根系土Se含量)位于0.005～2.373之間,平均值為0.216.從富硒小麥籽實空間分布位置看,大部分富硒小麥都分布在調查圈定的富硒耕地或者Se含量相對較高耕地中.小麥籽實中Se含量與對應根系土樣品中Se含量相關性統計顯示(圖6a),在置信度為0.01時,兩者顯著相關,Spearman相關系數為0.79.土壤中硒的水溶態和可交換態含量之和作為植物可利用態硒即有效態硒含量[33].研究區根系土壤中硒有效態含量與小麥籽實中硒含量具有顯著相關性,Spearman相關系數為0.75(<0.01)(圖6b).從相關性和空間位置的耦合關系看,研究區小麥籽實中Se主要來源于根系吸收、轉運土壤中生物可利用態硒.

圖6 農作物根系土與作物籽實硒含量相關性(N=116)

2.5 富硒土地資源開發利用適宜性評價

生態位適宜性評價是定量表述生物對其生境資源條件的適宜程度,其核心是用N維空間區域描述生態系統中每種生物生存所必需的生境最小閾值.生態位理論廣泛應用在土地整治[34]、建設用地適宜性評價[35]、高標準農田建設適宜性評價[36]、土地開發適宜性評價[37]等土地資源相關的領域.

2.5.1 生態位適宜性評價方法 本次生態位富硒土地資源適宜性評價的基本單元以新區全國第二次土地調查形成的圖斑為研究單元.根據研究區土地利用現狀,從富硒土地開發利用的主導因素、區位因素和限制因素3個維度構建評價指標因子體系[9].指標因子體系及指標效應見表3.

綜合考慮新區富硒土地資源開發利用主要在縣域范圍以內的中小尺度上進行,由于新區范圍內光照、溫度和氣候等無明顯差別,地勢平坦,因此從富硒土壤開發利用的自身條件、土地利用和生態環境3方面綜合考慮,選擇土壤硒含量、土壤有機質、土壤質地、土壤酸堿度、土地利用類型和土壤綜合污染指數等6個關鍵因子,構建富硒土地資源開發利用生態位適宜性評價指標體系.綜合污染指數是土壤中As、Cr、Cd、Cu、Hg、Ni、Pb和Zn 8種重金屬綜合污染等級,評價標準采用GB15618-2018標準中的篩選值[38].綜合污染指數采用內梅羅指數法計算,具體計算方法詳見文獻[12].

各指標生態位適宜度值是通過評判評價單元實現生態位與最適值的貼近程度確定,設定生態位適宜性指數是[0,1]之間.正向效應指標,其指標值越大則生態位值越大,其富硒土地資源開發利用適宜性越強,當評價因子為負向效應因子時,其指標值越小則生態位值越大,其富硒土地資源開發利用適宜性越強.各指標因子計分規則參照GB/T 28407-2012標準[39].

表3 富硒土地資源開發利用適宜性評價指標因子及指標效應

各指標因子權重采用層次分析法和熵權法確定.研究區6項因子指標利用層次分析法確定的指標權重W1,C=0.048,CR=0.038<0.1,分析結果滿足一致性檢驗;利用熵權法確定為W2(表4).采用平均法確定綜合指標權重W=(W1+W2)/2.

表4 富硒土地資源生態位適宜性評價指標權重

2.5.2 富硒土地資源適宜性評價 綜合考慮適宜度指數值,結合自然斷裂點分級結果,將=0.5作為是否適宜的閾值,研究區富硒土地資源開發利用適宜性劃分為4個等級.>0.5的區域劃為適宜區,由高到低依次劃分I級區域(>0.85),Ⅱ級區域(0.75<￡0.85)和Ⅲ級區域(0.5<￡0.75);將適宜性￡0.5的區域劃為Ⅳ級區域(即不適宜區).研究區各評價單元富硒土地資源開發利用適宜性等級評價劃分結果見圖7.

圖7 富硒土壤資源開發利用適宜性評價

新區富硒土地資源開發利用高、中適宜區土地面積為14397.35hm2,占新區規劃總面積的8.13%.高適宜區(Ⅰ級)土地面積為5314.45hm2,占新區規劃總面積的3.00%,其中耕地面積為5068.06hm2,主要位于新區西北容城縣南張鎮和小里鎮、安新縣寨里鄉等地.該高適宜區域主導因子較好(土壤硒含量富集),土壤肥力、土壤酸堿度、土壤質地等區位因素適宜,土壤環境質量等限制因素對該區域制約作用較小.中適宜區(Ⅱ級)土地面積為9082.90hm2,占新區規劃總面積的5.13%,其中耕地面積為5547.29hm2,主要位于新區西南部安新縣老河頭鎮、蘆莊鄉等地,該適宜區域土壤環境質量等限制因素對其產生了一定的制約作用.低適宜區(Ⅲ級)土地面積96392.31hm2,占新區規劃總面積的54.43%,主要是受主導因子(土壤硒含量)等影響.不適宜區(Ⅳ級)土地面積66292.93hm2,占新區規劃總面積的37.44%,不適宜區主要受主導因子(土壤硒含量)和區位因素(土地利用類型)等因素的共同控制.

3 結論

3.1 雄安新區土壤硒含量均值為0.22mg/kg,圈定表層富硒土地面積16993hm2,其中無重金屬污染的富硒耕地9380hm2.土壤硒賦存形態以強有機結合態為主,其次為殘渣態和腐殖酸結合態的硒,三者占土壤硒總量的87%以上;而生物可利用態硒含量均值所占比例相對較小(水溶態硒占比4.44%,離子交換態硒占比3.43%).全區小麥籽實樣品富硒率28.45%,小麥籽實中硒含量與根系土壤中硒含量具有顯著相關性(=0.77,<0.01).

3.2 新區土壤硒富集主要受第四紀沉積地貌環境的影響.土壤硒含量與沉積地貌環境“標志性”指標Al2O3、Fe2O3、MgO和K2O含量顯著正相關,與SiO2和Na2O含量顯著負相關(<0.01).土壤硒含量與有機碳(Corg.)、陽離子交換量(CEC)、黏粒和粉粒含量顯著正相關,與土壤砂粒含量顯著負相關(<0.01).此外,新區灌溉水、化肥和大氣干濕沉降(燃煤降塵)等外源輸入影響其土壤硒異常分布.

3.3 新區富硒土地資源開發利用高適宜區(Ⅰ級)和中適宜區(Ⅱ級)主要位于新區西北容城縣南張鎮和小里鎮、安新縣寨里鄉和老河頭鎮等地,土地面積為14397.35hm2,占新區規劃總面積的8.13%.

[1] Rayman M P. The importance of selenium to human health [J]. Lancet,2000,356(9225):233-241.

[2] 王學求,柳青青,劉漢糧,等.關鍵元素與生命健康:中國耕地缺硒嗎? [J]. 地學前緣,2021,28(3):412-423.

Wang X Q,Liu Q Q,Liu H L,et al. Key elements and human health: Is China’s arable land selenium-deficient? [J]. Earth Science Frontiers,2021,28(3):412-423.

[3] Rayman M P. Food-chain selenium and human health: emphasis on intake [J]. British Journal of Nutrition,2008,100(2):254-268.

[4] 趙中秋,鄭海雷,張春光,等.土壤硒及其與植物硒營養的關系 [J]. 生態學雜志,2003,22(1):22-25.

Zhao Z Q,Zheng H L,Zhang C G,et al. Advances in the studies on selenium in soil and selenium biological effect [J]. Chinese Journal of Ecology,2003,22(1):22-25.

[5] 尹 炳,夏沙沙,霍天滿,等.旱作區土壤硒富集特征及綜合利用分區 [J]. 干旱區資源與環境,2021,35(1):84-91.

Yin B,Xia S S,Huo T M,et al. Characteristics of Se enrichment and its comprehensive utilization in rainfed farm land [J]. Journal of Arid Land Resources and Environment,2021,35(1):84-91.

[6] 喬新星,晁 旭,任 蕊,等.陜西關中富硒土壤研究及開發利用——以三原—閻良地區為例 [J]. 物探與化探,2021,45(1):230-238.

Qiao X X,Chao X,Ren R,et al. Research,development and utilization of selenium-rich soil of Shaanxi: A case study of Sanyuan- Yanliang area [J]. Geophysical and Geochemical Exploration,2021,45(1):230-238.

[7] 周國華.富硒土地資源研究進展與評價方法 [J]. 巖礦測試,2020,39(3):319-336.

Zhou G H. Research progress of Selenium-enriched land resources and evaluation methods [J]. Rock and Mineral Analysis,2020,39(3): 319-336.

[8] 王惠艷,曾道明,郭志娟,等.天然富硒土地劃定的富硒閾值 [J]. 環境科學,2021,42(1):333-342.

Wang H Y,Zeng D M,Guo Z J,et al. Selenium threshold for the delimitation of natural selenium-enriched land [J]. Environmental Science,2021,42(1):333-342.

[9] 蔡海生,陳 藝,張學玲.基于生態位理論的富硒土壤資源開發利用適宜性評價及分區方法 [J]. 生態學報,2020,40(24):9208-9219.

Cai H S,Chen Y,Zhang X L. Suitability evaluation and zoning method research on development and utilization of selenium-rich soil resources based on niche theory [J]. Acta Ecologica Sinica,2020,40(24):9208-9219.

[10] 馬 震,夏雨波,李海濤,等.雄安新區自然資源與環境-生態地質條件分析 [J]. 中國地質,2021,48(3):677-696.

Ma Z,Xia Y B,Li H T,et al. Analysis of natural resources and environment eco-geological conditions in the Xiong'an New Area [J]. Geology in China,2021,48(3):677-696.

[11] 何登發,單帥強,張煜穎,等.雄安新區的三維地質結構:來自反射地震資料的約束 [J]. 中國科學.地球科學,2018,48(9):1207-1222.

He D F,Shan S Q,Zhang Y Y,et al. 3-D geologic architecture of Xiong’an New Area: constraints from seismic reflection data [J]. Science China Earth Sciences,2018,48(9):1207-1222.

[12] 郭志娟,周亞龍,王喬林,等.雄安新區土壤重金屬污染特征及健康風險 [J]. 中國環境科學,2021,41(1):431-441.

Guo Z J,Zhou Y L,Wang Q L,et al. Characteristics of soil heavy metal pollution and health risk in Xiong'an New District [J]. China Environmental Science,2021,41(1):431-441.

[13] 李 括,彭 敏,趙傳冬,等.全國土地質量地球化學調查二十年 [J]. 地學前緣,2019,26(6):128-158.

Li K,Peng M,Zhao C D,et al. Vicennial implementation of geochemical survey of land quality in China [J]. Earth Science Frontiers,2019,26(6):128-158.

[14] DZ/T 0295-2016 土地質量地球化學評價規范 [S].

DZ/T 0295-2016 Specification of land quality geochemical assessment [S].

[15] LY /T 1243-1999 森林土壤陽離子交換量的測定[S].

LY /T 1243-1999 Determination of cation exchange capacity in forest soil [S].

[16] NY /T 1121.3-2006 土壤檢測第3部分:土壤機械組成的測定[S].

NY /T 1121.3-2006 Soil teating Part 3: Method for determination of mechanical composition [S].

[17] DZ /T 0380 -2021 天然富硒土地劃定與標識[S].

DZ /T 0380 –2021 Delimitation and logo of natural selenium- enriched land [S].

[18] GB /T 22499-2008 富硒稻谷[S].

GB /T 22499-2008 Rich selenium paddy [S].

[19] 張仁鐸.空間變異理論及應用 [M]. 北京:科學出版社,2005:23-39

Zhang R D. Theory and application of spatial variation [M]. Beijing: Science Press,2005:23-29.

[20] 中國環境監測總站.中國土壤元素背景值 [M]. 北京:中國環境科學出版社,1990:2-257.

China National Environmental Monitoring Center. The Background Concentrations of Soil Elements in China [M]. Beijing: China Environmental Science Press,1990.2-257.

[21] 曹 峰,李瑞敏,王 軼,等.海河平原北部地區土壤地球化學基準值與環境背景值 [J]. 地質通報,2010,29(8):1215-1219.

Cao F,Li R M,Wang Y,et al. Soil geochemical baseline and environmental background values in northern Haihe plain,China [J]. Geological Bulletin of China,2010,29(8):1215-1219.

[22] Supriatin S,Weng L,Comans R N J. Selenium speciation and extractability in Dutch agricultural soils [J]. Science of the Total Environment,2015,532(1):368-382.

[23] Wang Z J,Gao Y X. Biogeochemical cycling of Selenium in Chinese environments [J]. Applied Geochemistry,2001,16(11/12):1345-1351.

[24] Dinh Q T,Cui Z W,Huang J,et al. Selenium distribution in the Chinese environment and its relationship with human health: A review [J]. Environment International,2018,112:294-309.

[25] Hartikainen H. Biogeochemistry of Selenium and its impact on food chain quality and human health [J]. Journal of Trace Elements in Medicine and Biology,2005,18(4):309-318.

[26] 張秀芝,趙相雷,李 波,等.基于區域土壤元素地球化學的河北平原土壤質地類型劃分 [J]. 第四紀研究,2017,37(1):25-35.

Zhang X Z,Zhao X L,Li B,et al. The classifying of soil texture types based on the regional soil geochemical elements in Hebei plain [J]. Quaternary Sciences,2017,37(1):25-35.

[27] 朱立新,馬生明,周國華,等.沖積平原區土壤元素組成特征及其示蹤作用 [J]. 地質與勘探,2002,38(4):56-59.

Zhu L X,Ma S M,Zhou G H,et al. Characters of soil elements and its tracing application in alluvial plain [J]. Geology and Exploration,2002,38(4):56-59.

[28] 張秀芝,馬忠社,魏 靜,等.利用元素地球化學特征研究河北平原第四紀沉積環境的探討與實踐 [J]. 地學前緣,2012,19(4):194-205.

Zhang X Z,Ma Z S,Wei J,et al. The use of geochemical characteristics of elements in the study of Quaternary sedimentary environment in the Hebei Plain [J]. Earth Science Frontiers,2012,19(4):194-205.

[29] 張秀芝,馬忠社,王蔭楠,等.河北平原土壤硒異常成因及其生態效應 [J]. 地球與環境,2012,40(4):541-547.

Zhang X Z,Ma Z S,Wang Y N,et al. The origin and ecological effects of Selenium abnormity in soil in Hebei Plain [J]. Earth and Environment,2012,40(4):541-547.

[30] 宋澤峰,李振寧,王豐翔,等.基于GIS的河北平原土壤Se年增長量估算 [J]. 福建師范大學學報(自然科學版),2019,35(4):77-85.

Song Z F,Li Z N,Wang F X,et al. Se annual increment estimated data in soils of the Hebei Plain basd on GIS technique [J]. Journal of Fujian Normal University (Natural Science Edition),2019,35(4):77-85.

[31] 李振寧.河北省平原區土壤中硒異常源追蹤及生態效應評價 [D]. 石家莊:石家莊經濟學院,2010.

Li Z N. The study on source tracking of Se anomaly and ecological appraisal in Plain Terrain of Hebei [D]. Shijiazhuang: Shijiazhuang University of Economics,2010.

[32] 張素榮.京津冀魯耕地區土地質量地球化學調查成果報告 [R]. 天津:中國地質調查局天津地質調查中心,2019.

Zhang S R. Report on the results of land quality geochemical survey in Beijing-Tianjin-Hebei-Shandong cultivated land area [R]. Tianjin: Tianjin Center,China Geological Survey,2019.

[33] Li Z,Liang D L,Peng Q,et al. Interaction between selenium and soil organic matter and its impact on soil selenium bioavailability: A review [J]. Geoderma,2017,295:69-79.

[34] 郭笑東,陳利根,畢如田,等.基于生態位理論的黃土丘陵區耕地整治優先度及模式研究 [J]. 水土保持通報,2019,39(1):184-190.

Guo X F,Chen L G,Bi R T,et al. priority and mode of cultivated land consolidation in loess hilly region based on niche-fitness model [J]. Bulletin of Soil and Water Conservation,2019,39(1):184-190.

[35] 張海龍,郭彥龍,高 蓓.基于生態位理論的多模型秦嶺山區建設用地適宜性評價——以商洛市商州區為例 [J]. 地理與地理信息科學,2016,32(4):83-89.

Zhang H L,Guo Y L,Gao B. Multi-model suitability assessment of construction land in Tsinling Mountains based on the niche theory: A case study of Shangzhou,Shangluo [J]. Geography and Geo- Information Science,2016,32(4):83-89.

[36] 趙素霞,牛海鵬,張捍衛,等.基于生態位模型的高標準基本農田建設適宜性評價 [J]. 農業工程學報,2016,32(12):220-228.

Zhao S X,Niu H P,Zhang B W,et al. Suitability evaluation on high quality capital farmland consolidation based on niche-fitness model [J]. Transactions of the Chinese Society of Agricultural Engineering,2016,32(12):220-228.

[37] 張榮群,王大海,艾 東,等.基于生態位和“反規劃”思想的城市土地開發適宜性評價 [J]. 農業工程學報,2018,34(3):258-264.

Zhang R Q,Wang D H,Ai D,et al. Suitability evaluation of urban land development based on niche and anti-planning [J]. Transactions of the Chinese Society of Agricultural Engineering,2018,34(3):258-264.

[38] GB 15618-2018 土壤環境質量農用地土壤污染風險管控標準(試行) [S].

GB 15618-2018 Soil environmental quality risk control standard for soil contamination of agricultural land [S].

[39] GB /T 28407-2012 農用地質量分等規程[S].

GB /T 22499-2008 Regulation for gradation on agriculture land quality [S].

Distribution characteristics and utilization evaluation of selenium-rich land resources in Xiong’an New District.

ZHOU Ya-long1,2,3,GUO Zhi-juan1,2,3*,WANG Qiao-lin1,2,3,WANG Cheng-wen1,2,3,SONG Yun-tao1,2,3,LIU Fei1,2,3

(1.Institute of Geophysical & Geochemical Exploration,Chinese Academy of Geological Sciences,Langfang 065000,China；2.Research Center of Geochemical Survey and Assessment on Land Quality,China Geological Survey,Langfang 065000,China；3.Key Laboratory of Geochemical Cycling of Carbon and Mercury in the Earth’s Critical Zone,Chinese Academy of Geological Sciences,Langfang 065000,China),2022,42(8)：3913～3921

Surface soil samples (0～20cm),bulk crops (wheat seeds) and root soil samples were collected from Xiong’an New District for characterizing the distribution and enrichment of soil selenium through GIS spatial analysis,correlation analysis,etc. Then,the suitability assessment of development and utilization of selenium-rich land resources was carried out with the theory model of niche suitability evaluation. The results show that the surface selenium-rich land in the new area is about 16993hm2; and the soil selenium is mainly in the form of strong organic binding,then in the residue and humic acid binding state,together accounting for more than 87% of the total pool. The proportion of average bioavailable selenium content is relatively small (water-soluble selenium accounts for 4.44%,ion-exchange selenium accounts for 3.43%). The selenium-enriched rate in wheat seed samples is 28.45% and significantly correlated with that in root soil (=0.77,<0.01). The content of soil selenium has a significant positive correlation with “signature” soil indicators of sedimentary geomorphological environment,such as the content of Al2O3,Fe2O3,MgO,K2O,soil organic carbon (Corg.),cation exchange capacity (CEC),clay and silt,and negatively correlated with the content of SiO2、Na2O and soil sand (<0.01). Obviously,the enrichment of soil selenium is mainly affected by Quaternary sedimentary geomorphic environment even though its anomaly could be related to external sources such as irrigation water,chemical fertilizers,atmospheric dry and wet deposition (coal burning and dust reduction). And the niche theory modeling results show that the high-suitability areas (level I) and the moderately suitable areas (level II) are mainly present in Nanzhang and Xiaoli Town of Rongcheng County,Zhaili and Laohetou Town of Anxin County in the northwest of the new area,covering an area of 14397.35hm2,or accounting for 8.13% of the area.

selenium；soil；root soil；crop；suitability evaluation；Xiongan New District

X53

A

1000-6923(2022)08-3913-09

2022-01-04

中國地質調查局地調項目(DD20189123,DD20190518)

* 責任作者,高級工程師,gzhijuan@mail.cgs.gov.cn

周亞龍(1984-),男,湖北武漢人,高級工程師,碩士,主要從事環境地球化學研究.發表論文40余篇.