基于通徑分析的土壤性質與硒形態的關系①
——以黑龍江省主要類型土壤為例

徐 強，遲鳳琴*，匡恩俊，張久明，宿慶瑞，韓錦澤，王 偉，魏 丹

（1 黑龍江省農業科學院土壤肥料與環境資源研究所，黑龍江省土壤環境與植物營養重點實驗室，哈爾濱 150086；2 東北農業大學資源與環境學院，哈爾濱 150086）

基于通徑分析的土壤性質與硒形態的關系①
——以黑龍江省主要類型土壤為例

徐 強1，2，遲鳳琴1，2*，匡恩俊1，張久明1，宿慶瑞1，韓錦澤1，2，王 偉1，魏 丹1

（1 黑龍江省農業科學院土壤肥料與環境資源研究所，黑龍江省土壤環境與植物營養重點實驗室，哈爾濱 150086；2 東北農業大學資源與環境學院，哈爾濱 150086）

黑龍江省是缺硒（Se）比較嚴重的省份之一，位于全國低Se帶的始端。本文選擇黑龍江省不同類型土壤，采用連續浸提法測定了土壤中Se的形態，并運用通徑分析法研究了土壤理化性質對土壤Se形態、全Se的影響。結果表明：總Se含量在不同類型土壤中差異較大，整體上處于中等偏低水平；各結合態Se中以有機態和殘渣態為主，分別占8.16% ～ 50.5% 和26.32% ～ 70.90%，酸溶態占0.80% ～ 33.97%，而水溶態與交換態僅占0.70% ～ 7.18% 和0.75% ～9.37%。通徑分析顯示，土壤理化性質通過直接和間接作用共同影響著土壤各賦存形態Se的含量和分解轉化方向，但其各自的作用機理和影響強度不同。在土壤各性質中，土壤有效鐵、錳和黏粒以及它們間的共同作用決定了土壤Se形態分布，其中，土壤有效鐵對于除殘渣態Se之外的4種形態Se的生成轉化起到了主導作用，而土壤有機碳、pH等其他性質主要通過正或負的間接作用影響Se形態。土壤總Se方面，土壤有效鐵、錳及黏粒含量對總Se有較強的富集作用，但土壤有機碳和pH等因素間的相互作用也不可忽略。

土壤硒；形態；理化性質；通徑分析；黑龍江省

硒（Se）是一種稀有分散元素，介于氧族非金屬元素硫（S）和金屬元素碲（Te）之間。世界上共40多個國家和地區缺Se，我國也有3/4以上的國土面積缺Se，Se攝取不足的人口達1億之多［1］。土壤中的Se主要來源于火山噴發及巖石風化，在地殼中的分布極不均勻［2］。在我國，既有黑龍江克山與河北張家口等低Se病區，也有陜西紫陽和湖北恩施等Se毒地區［3］。基于此，研究Se的賦存形態及其與土壤理化性質的關系對于Se的富集、遷移和轉化過程及人體健康安全等問題具有重要意義。目前，國內外多地已開展Se形態的相關研究，而在中國東北此方面的報道卻十分有限。

各結合態Se與土壤各個組分如有機質、鐵錳氧化物及黃銅礦（CuFeS2）、黃鐵礦（FeS2）等之間不斷發生著吸附-解吸、沉淀-溶解等地球化學過程［4］，土壤理化性質更是因為對Se在土壤中的賦存形態具有重要作用而受到研究者的廣泛關注［5］。目前對于土壤Se總量及各賦存形態的研究主要針對成土母質、土壤類型、土地利用方式和土壤理化性質等環境因素，所采用的方法多為簡單相關、多元回歸分析等［6-7］。然而，簡單相關分析不能全面考察變量間的相互關系，其結果往往帶有一定的片面性；多元回歸分析雖然在一定程度上能夠消除變量之間的多重共線性，能真實地反映出各個自變量與因變量的關系，但卻無法進一步深入分析相關變量之間對結果的共同影響。因此兩種方法都帶有一定的局限性。通徑分析在多元回歸的基礎上將相關系數分解為直接通徑系數（某一自變量對因變量的直接作用）和間接通徑系數（該自變量通過其他自變量對因變量的間接作用），分析結果更符合實際，是研究多個相關變量之間關系的更為有效的方法［8］。但迄今為止，該方法在土壤性質對土壤Se影響的研究中仍未見報道。為此，本研究選擇了黑龍江省具有代表性的土壤樣品，研究Se賦存形態在不同土壤類型中的分配規律，并運用通徑分析方法對土壤Se與土壤理化性質之間的關系進行了探討。研究結果能夠進一步闡釋環境因素對土壤Se賦存形態的作用，為揭示Se在土壤中的地球化學行為提供科學依據。

1 材料與方法

1.1 研究區概況

黑龍江省位于中國最東北部（43°26′～53°33′ N，121°11′～135°05′ E），土地總面積47.3萬km2，占全國陸地總面積的4.9%，人口3 833萬，全省耕地面積11.78萬km2，約占全國耕地面積的9%，居全國第1位。該省處于溫帶大陸性季風氣候區，是國內氣溫最低的省份，冬長夏短，年日照時數為2 621.4 h，年均降雨量在370 ～ 670 mm。全省地形地貌復雜，西有大興安嶺，北有小興安嶺，東南部為山地，中西部為松嫩平原，東北為三江平原，主要的土壤類型有：黑土、白漿土、暗棕壤、黑鈣土、草甸土、火山灰土、風沙土、鹽堿土及泥炭土等。

1.2 樣品采集與處理

2013年7月，對上述5個自然地理區域具有代表性的土壤進行了采樣。對表層土壤以梅花形多點采樣法，垂直采集地表至20 cm 深處的土柱，保證上下均勻采集，去除樣品中草根、礫石磚塊等雜物，混勻并用四分法取約2 kg；土壤剖面采用人工挖井，井深一般為2 m，按照土壤剖面層次和機械分層結合的方式采集樣品。原則上，耕作層獨立采樣，厚達10 cm者獨立采樣，＜10 cm者并層采樣，樣品質量2 kg左右。全部帶回實驗室，風干，磨碎，分別過不同孔徑尼龍篩后裝瓶備用。

1.3 分析項目與方法

土壤Se形態分級：參照瞿建國等［9］的逐級連續化學浸提技術，將土壤中Se分為水溶態 （0.25 mol/L KCl 浸提）、交換態（0.7 mol/L pH 5.0 KH2PO4浸提）、酸溶態（2.5 mol/L HCl 浸提）、有機態（5% K2S2O8+ 1∶1 HNO3浸提） 和殘渣態（HCl-HNO3混酸浸提）。

土壤總Se：土壤總Se的測定采用中華人民共和國國家環境保護標準HJ680-2103［10］推薦的HCl-HNO3法消解。稱取風干、過篩（0.149 mm）的土壤樣品0.1 ～ 0.5 g（精確至0.000 1 g）置于聚四氟乙烯消解內罐中，于通風櫥中依次順序加入6 ml HCl和2 ml HNO3，若有劇烈反應可加彎口漏斗靜置一宿，次日將消解內灌嵌入消解外灌和保護支架，一并置于微波消解儀（上海新拓XT-9916型密閉式智能微波消解/萃取儀）爐腔中。Se測定采用氫化物發生-原子熒光光譜法，儀器為北京海光AFS-230E雙道原子熒光光度計。用國家標準物質（土壤GSS-10）檢驗結果的可靠性。

土壤理化性質：土壤粒徑分布（particle size distribution，PSD）用Malvern 2000激光粒度儀測定［11］；土壤總有機碳（soil organic carbon，SOC）用HT1300干法測定［12］，儀器為德國耶拿（Analytikjena）公司生產的multi N/C 2100型號總有機碳/總氮分析儀。土壤pH用無CO2水浸提，pH計測定；土壤電導率（electric conductivity，EC）用電導儀測定；土壤有效態鐵、錳、銅、鋅（分別用aFe、aMn、aCu、aZn表示）采用DTPACaCl2-TEA體系浸提，原子吸收分光光度計測定；土壤有效硫（用aS表示）采用Ca（H2PO4）2·H2O-CH3COOH體系浸提，BaCl2比濁法測定［13］；土壤理化性質見表1。

1.4 數據分析

數據處理和統計分析運用Microsoft Office Excel 2013和SPSS 19.0軟件完成，采用逐步回歸分析方法建立環境因子對土壤總Se和各形態Se含量的多元回歸方程，并對其顯著性影響因子進行識別，然后根據回歸分析結果對識別出的顯著性影響因子進行通徑分析。作圖采用Origin 8.0 軟件完成。

表1 供試土壤的物理化學性質Table1 Physico-chemical properties of tested soils

2 結果與分析

2.1 土壤中總硒及各形態硒分布

由表2可知：黑龍江省主要類型土壤總Se含量變化較大，按照Tan等［14］對我國表層土壤總Se含量的分級標準：缺Se土壤（＜0.125 mg/kg）、潛在缺Se土壤（0.125 ～ 0.175 mg/kg）、中等含Se土壤（0.175 ～ 0.450 mg/kg）、富Se土壤（0.450 ～ 2.000 mg/kg）、高Se土壤（2.000 ～ 3.000 mg/kg）、過量Se土壤（＞3.000 mg/kg），黑龍江省土壤Se含量整體上處于中等偏低水平。在石灰巖、頁巖等沉積巖上發育的土壤，因受到母質的影響，Se含量較高，不同土類中以小興安嶺山地富含的暗棕壤總Se含量最高，而松嫩平原的風沙土及鹽堿土總Se含量偏低。

表2 供試土壤各形態Se含量Table2 Contents of various Se fractions in tested soils

研究區表土層水溶態Se含量極低且變異性較大，最大值 （暗棕壤，19.77 μg/kg） 與最小值 （火山灰土，1.21 μg/kg） 之間相差16倍之多，但不同類型土壤水溶態Se含量占總Se的百分比差異并不大，為0.70% ～ 7.18%。交換態Se與水溶態Se規律相似，其中最大值 （28.50 μg/kg） 與最小值 （1.31 μg/kg）相差近22倍，其含量占總Se的百分比為0.75% ～9.37%。酸溶態Se含量在不同類型土壤中占總Se的百分比差異較大，為0.80% ～ 33.97%；不同類型土壤中的Se以有機態與殘渣態為主，分別占總硒的8.16%～ 50.50%和26.32% ～ 70.90%。全部樣品中有63% 以殘渣態Se為主，占總Se的34.17% ～ 70.90%；31% 的土壤以有機態Se為主，占總Se的40.16% ～ 50.50%，而鹽堿土主要以酸溶態Se存在，占總Se的33.97%（圖1）。

2.2 土壤各形態硒與土壤理化性質的通徑分析

2.2.1 相關分析 結果表明（表3），土壤有機碳及陽離子交換量與土壤水溶態Se、交換態Se及酸溶態Se含量呈負相關關系，而與有機態Se和殘渣態Se均呈正相關關系，但相關性均未達到顯著水平；土壤有效態鐵、錳、銅、鋅含量與土壤各形態Se含量關系密切，其中土壤有效鐵對土壤水溶態Se、交換態Se、有機態Se及總Se的正相關關系達到極顯著水平；土壤各形態Se含量與土壤有效硫、pH、土壤電導率及砂粒含量存在不同程度的負相關關系，與黏粒含量呈正相關關系，其中有機態Se及殘渣態Se與土壤有效鐵、黏粒和砂粒的相關性分別達到顯著或極顯著水平。土壤總Se方面，土壤有效硫、pH、電導率及砂粒含量與土壤總Se呈負相關關系，其中與砂粒達到顯著相關水平；其他土壤理化性質與總Se均呈現正相關關系，其中土壤有效鐵、黏粒含量與之均達到極顯著相關水平。

圖1 供試土壤各形態Se含量百分比Fig. 1 Percentage of various Se fraction in tested soils

表3 土壤Se與土壤理化性質間的相關系數Table3 Correlation coefficients between various Se fractions and soil properties

2.2.2 逐步多元回歸分析 分別以土壤各形態Se、總Se為因變量，土壤理化性質為自變量，進行逐步多元回歸分析，得到回歸方程如表4所示。

2.2.3 通徑分析 表5表明，土壤有效鐵含量對水溶態Se含量表現出明顯的直接正效應（0.969 6）和較小的間接負效應（-0.207）。土壤有效錳量對水溶態Se含量有一定的直接負效應影響，為 -0.321 1，但其通過土壤有效鐵對水溶態Se含量也產生了具有正效應的間接作用（0.456 9）。土壤有機碳含量和有效鐵及有效錳含量相比，無論直接作用還是間接作用都表現最低。此外，所選環境因子對水溶態Se含量的決定系數R2為0.728，剩余因子的通徑系數e為0.521 5，該值較大，說明所選環境因子未能充分地解釋土壤水溶態Se含量的變異。這可能是由于農業土壤中常年種植的作物不斷從土壤中提取水溶態Se，加上外源Se肥甚至是人類農業活動帶來的有機、無機Se污染物的介入［15］，使得水溶態Se含量總是處于動態變化之中所致。

表4 土壤Se與土壤理化性質的回歸分析Table4 Regression analysis between Se contents and soil properties

表5 土壤性質對水溶態Se含量影響的通徑系數Table5 Path analysis coefficients between soil properties and SOL-Se content

表6表明，與水溶態Se相似，土壤有效鐵含量對交換態Se含量也表現出明顯的直接正效應（0.728 2）和較小的間接負效應（-0.049 8）。土壤有機碳含量對交換態Se含量影響的直接負效應（-0.303 6）和間接正效應（0.119 5）有所抵消，表明土壤有機碳不是影響土壤交換態Se含量的主要因素。另外，所選環境因子對交換態Se含量的決定系數R2為0.550，剩余因子的通徑系數e為0.670 8，該值較大，說明所選環境因子未能較為充分地解釋土壤水溶態Se含量的變異。交換態Se是植物有效Se，Wang等［16］對于克山病區土壤的研究發現，土壤有效Se和氣候因素 （海拔和降水） 關系密切，共同影響著該地區居民Se的攝入量。而Na-和Mg-硫酸鹽、氧化還原電位和含水量等因素共同控制著局部土壤Se的沉淀和溶解，進而影響土壤Se對植物的有效性［17］，而涉及上述因素的影響機制需要進一步研究。

表6 土壤性質對交換態Se含量影響的通徑系數Table6 Path analysis coefficients between soil properties and EX-Se content

表7表明，影響土壤酸溶態Se的因素較多，其中，土壤電導率對酸溶態Se含量具有明顯的直接負效應（-1.708 1），但同時其通過有效硫含量等因子對酸溶態Se含量也產生了強烈的具有正效應的間接作用（1.555 9）。土壤有效硫含量對酸溶態Se含量具有明顯的直接正效應（1.358 1），但同時其通過電導率等因子對酸溶態Se含量也產生了強烈的具有負效應的間接作用（-1.482 3），使得土壤有效硫含量與酸溶態Se含量在表觀上顯示出一種負相關關系。因此，若由相關系數 -0.124 （表4） 就簡單地認為有效硫對土壤酸溶態Se僅具有負效應是不恰當的。土壤有效銅、鐵含量對酸溶態Se含量的直接正效應較大 （0.807 2和0.522 5），而間接負效應較小（-0.231 7和-0.187 6）。土壤有效錳與陽離子交換量對酸溶態Se含量均有一定的直接效應影響（-0.436 4和-0.243 5），但幾乎被其間接效應（0.380 7和0.223 9）所抵消，表明兩種因子不是影響酸溶態Se含量的主要因素。

表7 土壤性質對酸溶態Se含量影響的通徑系數Table7 Path analysis coefficients between soil properties and FMO-Se content

表8表明，土壤有效鐵與黏粒含量對有機態Se量的直接通徑系數分別為0.728 5和0.524 0，這說明有效鐵對有機態Se的直接正效應最大，黏粒次之。土壤有效鋅對有機態Se的直接負效應（-0.492 2）幾乎被間接正效應（0.455 7）所抵消，表明土壤有效鋅不是影響有機態Se含量的主要因素。

表8 土壤性質對有機態Se含量影響的通徑系數Table8 Path analysis coefficients between soil properties and OM-Se content

表9 土壤性質對殘渣態Se含量影響的通徑系數Table9 Path analysis coefficients between soil properties and RES-Se content

表9表明，兩種環境因子無論直接效應（-0.454 9和0.331 3）還是間接效應（-0.122 7和0.168 5）對殘渣態Se的貢獻都不大，這可能是由于殘渣態Se是研究區土壤硒的主要賦存形態（表2和圖1），殘渣態Se的多少取決于礦物的天然組成，Se常與硫化物礦等共生，牢固地結合在晶格中，不因土壤理化性質的改變而變異，在自然環境條件下極難轉化成植物吸收利用的Se形態［18］。此外，所選環境因子對殘渣態Se的決定系數R2=0.428 0，剩余因子的通徑系數e為0.756 3，該值較大，說明對殘渣態Se含量影響較大的一些因素還沒有考慮到，如成土母質、巖石風化過程及土地利用方式等因素，有待進一步研究。

表10表明，土壤有效鐵對土壤總Se含量表現出明顯的直接正效應（0.728 2）和較小的間接負效應（-0.049 8）。土壤有效錳的間接通徑系數大于其直接通徑系數，表明其對土壤總Se含量的主要貢獻表現為通過影響土壤有效鐵及砂粒含量等環境因子而產生的間接作用效應。土壤質地對總硒含量也有一定的影響，具體表現為土壤黏粒含量通過土壤有效鐵等因子對總Se產生間接正效應，而通過土壤有效錳對總Se含量產生間接負效應。土壤砂粒含量通過土壤有效鐵等因子對總Se含量產生間接負效應，而通過土壤有效錳對總Se含量產生間接正效應。由此可知，若僅從相關系數0.672和-0.542 （表4） 就簡單地認為土壤質地與總Se含量之間僅具有顯著的直接影響是不恰當的。

表10 土壤性質對總Se含量影響的通徑系數Table10 Path analysis coefficients between soil properties and total Se content

3 討論

有學者指出，在全國范圍內存在一條從中國東北地區向西南方向，經過黃土高原再向西南延伸到西藏高原的低Se帶，其Se含量均值僅為0.1 mg/kg。黑龍江省位于全國低Se帶的始端，早在1935年黑龍江省克山縣就有關于人類克山病的報道，而家畜的白肌病 （畜禽硒-維生素E 缺乏綜合癥） 又相繼在1936年發生［19］。由于黑龍江省地形地貌復雜，地球化學條件具有明顯差別，這導致本實驗所采集的土壤樣品總Se含量變幅較大。其中，泥頁巖黑土由于受到了母質的影響總Se含量較高（0.291 4 mg/kg），而在石灰巖等沉積巖上發育的暗棕壤也較高（0.462 8 mg/kg）。風沙土及鹽堿土由于pH較高（＞8.5）而有機質含量較低等特點，使Se的遷移能力增強而具有較低的含Se量。其他農田土壤Se含量偏低，屬于中等及缺乏水平。

土壤Se形態方面，研究結果表明，殘渣態Se是研究區土壤中Se的主要存在形態，63% 的土樣以殘渣態Se為主，占總Se的34.17% ～ 70.9%。這部分Se是以較穩定的化合物和晶格的形態存在，難溶解于水和一般的酸堿性溶液，植物難以吸收利用，但這部分Se是土壤Se的重要儲備庫源；其次是有機態Se，有31% 的樣品以有機態Se為主，占土壤總Se的40.16% ～ 50.50%，有機態Se可以礦化成硒酸鹽或亞硒酸鹽被植物吸收，屬于土壤潛在有效Se。只有鹽堿土主要以酸溶態Se存在，占總Se的33.97%，一般認為這部分Se極難轉化為植物可吸收利用的Se形態，而被視為無效Se。植物可吸收利用的水溶態Se與交換態Se含量較低，在不同類型土壤間的差異較大，但占土壤總Se含量的百分比差異不大，分別為0.70% ～ 7.18% 和0.75% ～ 9.37%，這也是導致黑龍江省土壤Se生物有效性低的主要原因。

目前，國內外學者在土壤Se形態及其影響因素之間關系的研究中多以簡單相關分析和回歸分析加以探討，而結果不盡相同，主要原因在于影響Se形態間轉化的因素太多、太復雜，很多時候單純的相關或回歸分析會因沒有考慮到土壤性質之間內在的相互作用而得出錯誤的結論，而實際上變量之間的關聯關系往往是一個復雜的傳遞過程。通徑分析將自變量與因變量之間的相關系數分解為直接通徑系數和間接通徑系數，這樣不僅可以分析自變量對因變量的直接影響，還可以分析其通過其他相關的自變量對因變量的間接影響，從而能夠透過簡單相關的表面深入研究自變量與因變量之間的因果關系，揭示各個因素對結果的相對重要性。

通徑分析結果表明，土壤理化性質通過直接和間接作用共同影響著土壤各賦存形態Se的含量和分解轉化方向，但其各自的作用機理和影響強度不同。土壤有效鐵對于除殘渣態Se之外的4種形態Se均具有較強的直接正效應，在土壤各性質中，土壤有效鐵、錳和黏粒以及它們間的共同作用決定了土壤Se在各個形態中的分配，是影響其變化的主導因素，而土壤有機碳、pH等其他性質主要通過正或負的間接作用影響Se形態。土壤有效鐵主要指來自游離鐵中的活性鐵，其中無定形氧化鐵和絡合態鐵是其主要供給形態。這可以從Wang等［20］的研究中得到驗證，認為土壤鐵、錳等氧化物對于土壤Se的吸附和固持具有重要貢獻，土壤無定形鐵在對于除水溶態Se之外的4種Se形態均為正影響，且其影響作用大于土壤有機質。研究表明，土壤黏粒對于Se有較強的富集作用，黏粒含量越高，土壤的保肥性越好，能有效地減少Se的流失。土壤顆粒對Se的吸附量主要與其邊面結構有關，故其粒徑大小和穩定性直接影響土壤中Se的含量。Wang和Chen［21］研究指出，與土壤細粉粒及黏粒相結合的土壤金屬氧化物及有機物質富集Se的效果顯著。Kausch和Pallud［22-23］通過建立二維動態模型研究發現，對于直徑大于1 μm的黏粒，隨著黏粒粒徑的增大，其對Se的固定能力越強，并據此提出了通過提高土壤團聚情況來增強保Se能力的措施。徐文等［24-25］發現，粒徑 ＞1 mm 的顆粒對Se基本沒有吸附固定能力，而粒徑 ＜0.025 mm 的顆粒對亞硒酸鈉有較強的吸附固定作用。這些都與本研究通過通徑分析得到的黏粒對Se含量影響的結果相吻合。土壤總Se方面，土壤金屬氧化物及黏土礦物依然是制約土壤總Se變異的主要因素。另外，前人研究表明土壤中的碳、硫、磷和氮元素與Se元素均有較好的相關性，而在表淺環境下涉及成土過程和生物代謝過程的元素，其地球化學行為對Se的分布與遷移具有深刻的影響［26］。本研究盡管通徑分析未把土壤有機碳及pH納入模型之中，但這兩種因子對于土壤總Se的作用不可忽略。黑龍江省地勢地貌復雜，加之氣候上的干濕交替和凍融作用，使得不同地區土壤酸堿度及有機質的積累具有明顯的差異，黑土、暗棕壤及草甸土往往具有較高的土壤有機質含量而pH較低，風沙土及鹽堿土則相反。研究表明富含有機質的土壤當土壤pH較低（＜7.5）時，土壤中鐵錳氧化物和有機質更趨向富集Se，使Se的遷移能力減弱，本研究的測定結果與之相吻合（表2）。土壤性質間對Se的作用及影響機理往往相輔相成，更多時候不能僅僅考慮其中某一個因素的影響［27-30］。今后可在相對控制某一變量的前提下，采用定量模型技術、同位素示蹤等手段分析Se的有效性、形態的真實影響因素和定量規律［31］。

4 結論

對黑龍江省典型土壤Se 含量和分布的研究發現，供試土壤Se 含量總體上處于中等偏低水平，不同土類中以小興安嶺山地富含的暗棕壤總Se含量最高，而松嫩平原的風沙土及鹽堿土總Se含量偏低。利用連續浸提技術對土壤Se形態的分析表明，殘渣態Se及有機態Se是研究區低Se土壤中Se 的主要存在形態， 其次是酸溶態Se，植物可吸收利用的水溶態Se和交換態Se含量很低，這是導致黑龍江省土壤Se生物有效性低的主要原因。通徑分析結果表明，在土壤各性質中，土壤有效鐵、錳和黏粒以及它們間的共同作用決定了土壤Se形態分布，其中，土壤有效鐵對于除殘渣態Se之外的4種形態Se的生成轉化起到了主導作用，而土壤有機碳、pH等其他性質主要通過正或負的間接作用影響Se形態。土壤總Se方面，土壤有效鐵、錳及黏粒含量對總Se有較強直接影響作用，是其變異的主要因素，但土壤有機碳和pH等因素的作用也不可忽略。

［1］ 中國環境監測總站. 中國土壤元素背景值［S］. 北京： 中國環境科學出版社， 1990： 134-137， 370-373

［2］ Geerke H F， Gabriela R R. Selenium in volcanic environments： A review［J］. Applied Geochemistry， 2012，27（3）： 517-531

［3］ Yang G Q， Wang Z Z， Zhou R H， et al. Endemic selenium intoxication of humans in China［J］. The American Journal of Clinical Nutrition， 1983（37）： 872-881

［4］ 朱建明， 梁小兵， 凌宏文， 等. 環境中硒存在形式的研究現狀［J］. 礦物巖石地球化學通報， 2003， 22（1）： 75-81

［5］ Antana A， Jadvyg L， Sarunas A， et al. Selenium concentration dependence on soil properties［J］. Journal of Food， Agriculture & Environment， 2008， 1（6）： 163-167

［6］ 田應兵， 陳芬， 熊明彪， 等. 若爾蓋高原濕地土壤硒的數量、形態與分布［J］. 水土保持學報， 2004， 18（3）： 66-70

［7］ 張艷玲， 潘根興， 胡秋輝， 等. 江蘇省幾種低硒土壤中硒的形態分布及生物有效性［J］. 植物營養與肥料學報. 2002， 8（3）， 355-359

［8］ 蓋鈞鎰. 試驗統計方法［M］. 4版. 北京： 中國農業出版社，2013： 1-150

［9］ 瞿建國， 徐伯興， 龔書椿. 連續浸提技術測定土壤和沉積物中硒的形態［J］. 環境化學， 1997， 16（3）： 277-283

［10］ 環境保護部科技標準司. 土壤和沉積物 汞、砷、硒、鉍、銻的測定微波消解/原子熒光法（HJ680-2013）［S］. 北京：中國標準出版社， 2014

［11］ Malvern Instruments Ltd. MASTERSIZER 2000 MANUAL［S］. United Kingdom： 2007

［12］ Analytik Jena AG. MULTI N/C 2100 /2100S 中文說明書［S］. Germany： 2000

［13］ 魯如坤. 土壤農業化學分析方法［M］. 北京： 中國農業科技出版社， 2000

［14］ Tan J A， Zhu W Y， Wang W Y， et al. Selenium in soil and endemic diseases in China［J］. Science of the Total Environment， 2002， 284（1）： 227-235

［15］ 王松山， 梁東麗， 魏 威， 等. 基于路徑分析的土壤性質與硒形態的關系［J］. 土壤學報， 2011， 48（4）： 823-830

［16］ Wang J， Li H.， Li Y， et al. Speciation， distribution， and bioavailability of soil selenium in the Tibetan Plateau Kashin-Beck disease area—A case study in Songpan County， Sichuan Province， China. Biological Trace Element Research， 2013， 156（1/2/3）： 367-375

［17］ Gerla P J， Sharif M U， Korom S F. Geochemical processes controlling the spatial distribution of selenium in soil and water， west central South Dakota， USA. Environmental Earth Sciences， 2011， 62（7）： 1 551-1 560

［18］ 王子健， 孫喜平， 孫景芳. 土壤樣品中Se 的結合態分析［J］. 中國環境科學， 1988， 8（6）： 51- 54

［19］ 郭慶雨， 張佳誼. 黑龍江省及內蒙呼倫貝爾盟地區畜禽常用植物性飼料及土壤含硒量的調查［J］. 黑龍江八一農墾大學學報， 1982（2）： 31-51

［20］ Wang S S， Liang D L， Wang D. Selenium fractionation and speciation in agriculture soils and accumulation in corn（Zea mays L.） under field conditions in Shaanxi Province，China［J］. Science of the Total Environment， 2012， 427/428：159-164

［21］ Wang M C， Chen H M. Forms and distribution of selenium at different depths and among particle size fractions of three Taiwan soils［J］. Chemosphere， 2003，52（3）： 585-593

［22］ Kausch M F， Pallud C E. Modeling the impact of soil aggregate size on selenium immobilization. Biogeosciences，2013， 10（3）： 1 323-1 336

［23］ Tolu J， Hecho L I， Bueno M， et al. Selenium speciation analysis at trace level in soils. Analytica Chimica Acta， 2011， 684（1/2）： 126-133

［24］ 徐文， 唐文浩， 鄺春蘭， 等. 海南省土壤中硒含量及影響因素分析［J］. 安徽農業科學， 2010（6）： 3 026-3 027

［25］ 周鑫斌， 于淑慧， 謝德體. pH和三種陰離子對紫色土亞硒酸鹽吸附-解吸的影響［J］. 土壤學報， 2015， 52（5）： 1 069-1 077

［26］ Li Y H， Wang W Y， Luo K L. Environmental behaviors of selenium in soil of typical selenosis area， China［J］. Journal of Environmental Sciences， 2008， 20： 859-864

［27］ Elrashidi M A， Adriano D C， Workman S M， et al. Chemical equilibria of selenium in soils： A theoretical development. Soil Science， 1987， 144（2）： 141-152

［28］ 戴慧敏， 宮傳東， 董北， 等. 東北平原土壤硒分布特征及影響因素［J］. 土壤學報， 2015， 52（6）： 154-161

［29］ 鄭甲成， 劉婷. 不同濃度硒肥對秈稻硒含量和產量的影響［J］. 土壤， 2014， 46（1）： 88-93

［30］ 張 木， 胡承孝， 孫學成， 等. 鉬硒配施對小白菜鉬硒形態及硒價態的影響［J］. 2014， 46（5）： 894-900

［31］ Mitchell K， Raoul M C， Thomas M J， et al. Selenium sorption and isotope fractionation： Iron（III） oxides versus iron（II） sulfides［J］. Chemical Geology， 2013， 342： 21-28

Relationship Between Soil Physico-chemical Properties and Selenium Species Based on Path Analysis

XU Qiang1，2， CHI Fengqin1，2*， KUANG Enjun1， ZHANG Jiuming1， SU Qingrui1，HAN Jinze1，2， WANG Wei1， WEI Dan1
（1 Institute of Soil Fertilizer and Environment Resources， Heilongjiang Academy of Agriculture Sciences， Key Laboratory of Soil Environment and Plant Nutrition of Heilongjiang Province， Harbin 150086， China； 2 College of Resources and Environment，Northeast Agricultural University， Harbin 150086， China）

Samples of different kinds of soils were collected from Heilongjiang Province for analysis of Selenium （Se）species using a sequential extraction method. Se in soils may be partitioned into： 1） water soluble Se （SOL-Se）， 2） exchangeable Se and Se bound to carbonate （EX-Se）， 3） Acid soluble Se （FMO-Se）， 4） Se bound to organic matters （OM-Se）， 5） residual Se（RES-Se）. Relationships between soil physico-chemical properties and Se species were investigated with path analysis. Results showed that total Se content in the soils from the province was low， but varied with soil types. In terms of Se species， Se in the soils existed mainly in the forms of OM-Se and RES-Se， which accounted for 8.16%-50.5% and 26.32%-70.90% of total Se，respectively. FMO-Se accounted for 0.80%-33.97%， but SOL-Se and EX-Se only accounted for 0.70%-7.18% and 0.75%-9.37% of total Se in the soils， respectively. The results of path analysis indicated that soil properties affected transformation direction and contents of Se species directly or indirectly， but their influencing mechanisms were different. Soil available iron （aFe） showed stronger positive effects on MO-Se， FMO-Se， SOL-Se and EX-Se. The contents of aFe， available Mn （aMn） and clay as well as the interaction among them determined largely the distribution of Se species in the soils. Soil organic carbon （SOC） and pH had indirect impacts on Se forms positively or negatively. The contents of aFe， aMn and clay also showed stronger direct impacts on soil total Se content. The interaction between soil SOC and pH could have certain influence on soil total Se.

Soil selenium； Species； Physico-chemical properties； Path analysis； Heilongjiang Province

S151.9

10.13758/j.cnki.tr.2016.05.022

公益性行業（農業）科研專項（201303106）資助。

*通訊作者（fqchi2013@163.com）

徐強（1990—），男， 黑龍江哈爾濱人，碩士研究生，主要從事土壤硒和養分空間變異方面的研究。E-mail： qiangxu2014@163.com