北京市掛甲峪地區土壤全硒含量分布特征①

崔 萌，孫向陽*，姜新福，李素艷，楊少斌，李 松，龔小強，張 樂

北京市掛甲峪地區土壤全硒含量分布特征①

崔 萌1，孫向陽1*，姜新福2，李素艷1，楊少斌1，李 松1，龔小強1，張 樂1

(1 北京林業大學林學院，北京 100083；2 北京市通州區林業種子苗木服務中心，北京 101100)

為調查北京市掛甲峪地區土壤全硒含量分布特征，探討不同土地利用方式和海拔高度對土壤全硒含量的影響，采集掛甲峪地區表層土壤(0 ~ 20 cm)樣品122個，采用GIS方法對土壤全硒含量進行了空間分布特征研究。結果表明：掛甲峪地區土壤全硒含量變幅為0. 06 ~ 0.76 mg/kg，平均值為0. 30 mg/kg。86.88% 的土壤全硒含量處于足硒及以上水平(> 0.175 mg/kg)，其中有20.49% 的土壤屬于富硒含量(0.40 ~ 0.30 mg/kg)。不同土地利用類型間以草地土壤全硒含量平均值最高(0.40 mg/kg)，其次是林地(0.37 mg/kg)，最低為園地(0.25 mg/kg)。表層土壤全硒含量隨海拔高度的增加呈上升趨勢。土壤硒元素的塊金值為0.001 4，塊金系數為88.10%，空間相關性弱。掛甲峪地區整體上呈北西向帶狀分布，存在富硒土壤。

掛甲峪；空間相關性；硒

硒(Se)是生物體必需的微量元素之一，是環境中重要的生命元素[1-2]。土壤中的硒元素被植物吸收后通過食物鏈進入人體，發揮增強免疫、抗氧化和重金屬解毒等多種生物學作用[3]，但其過量或缺乏都會導致人體產生多種病癥[4-5]。研究表明，貫穿我國從東北到西南的低硒帶與我國克山病和大骨節病帶的分布完全吻合[6-7]。自然土壤中的硒元素主要來源于土壤母質和巖石，不同地區土壤中硒的含量分布極不均勻，從世界范圍來看，全球地殼中硒的平均豐度為0.13 mg/kg，普遍缺乏[8]。我國是世界上公認的嚴重缺硒國家之一，約有72% 的地區處于缺硒和低硒狀態，其中約有30% 以上為嚴重缺硒[9]。

在以往研究中，北京地區土壤硒含量分布規律差異較大。戴昭華和黃衍初[10]調查發現，北京地區表層土壤硒含量在0.05 ~ 0.48 mg/kg，平均值為0.275 mg/kg。郭莉等[11]調查發現，北京平原地區有高硒土壤，沿平原區西部山前斷續出現，并在市區零星分布，硒在土壤表層富集。鄧厚培和陳慶中[12]對北京東郊順義調查發現，土壤硒含量在0.041 ~ 0.133 mg/kg，含量較低。黃淇等[13]對北京市房山平原區調查發現，該區域存在富硒土壤，土壤硒含量受山區成土母質的影響較大。李娜等[14]對北京周口店地區調查發現土壤環境呈低硒狀態，且當地居民血清硒水平普遍處于缺乏狀態。可見，當前針對北京土壤硒的含量及分布已有較多相關研究，北京地區土壤硒含量整體上屬于缺乏水平，局部地區硒含量差異較大。但對硒元素存在的地區空間分布特征及系統性深入調查還不夠充分。

平谷是北京地區桃的主產區，平谷掛甲峪具有“平谷大桃村”之稱，果品年總產量可達150萬kg，是北京市重要的水果產出地區之一，“掛甲峪名果”等名牌享譽國內外[15]。本研究擬以北京掛甲峪地區為例，采集該區域表層土壤(0 ~ 20 cm)，測定土壤全硒含量，分析硒元素的空間分布特征，為地方發展富硒農產品以及合理利用硒資源等提供理論依據。

1 材料與方法

1.1 研究區概況

掛甲峪地區(117°05′27.90 ~ 117°06′32.50″E，40°15′08.15″ ~ 40°16′15.11″N)位于北京市平谷區北部，屬于暖溫帶半濕潤大陸性季風氣候，海拔123.18 ~ 371.62 m，年均氣溫10 ~ 20 ℃，年均降水量644 mm，面積2.411 km2[15]。該區域北、東、南三面環山，中間為一狹小的丘陵盆地。土壤巖性以砂巖、砂質礫巖、砂質及粉砂質頁巖、碳酸鹽巖、黏土質巖石等沉積巖為主。土壤母質為洪積沖積物、殘積物和坡積物，土壤類型為褐土。土地利用方式主要為林地、園地和草地，林地大部分位于中低山區，園地主要分布在淺山區，草地分布以山前區、川谷及地勢較平緩的崗臺地帶為主。區內林果業發達，林木覆蓋率達到70% 以上。

1.2 樣品采集與測定

2016年8月對調查區0 ~ 20 cm土壤進行樣品采集，采用100 m×100 m的網格，避開房屋及道路等設施，共采集樣品122個。采樣點分布如圖1所示。采用GPS定位記錄取樣點空間坐標。土壤樣品經風干、磨碎后過100目篩，土壤全硒含量測定采用NY/ T1104—2006氫化物發生-原子熒光光譜法(HG-AFS)。

圖1 研究區采樣點分布

1.3 數據處理

數據采用SPSS 21.0和Excel 2007進行整理與統計分析，采用GS+7.0對硒元素進行空間結構性分析。采用單樣本非參數檢驗法(One-Sample Kolmogorov- Smirnov Test)對樣本數據進行正態分布檢驗(＜0.05，2-tailed＝[16]，硒元素服從K-S正態分布。以殘差(RSS)最小和決定系數(2)最大為原則，選擇最佳的半變異函數擬合模型[17]。其中，塊金值(C0)通常由測量誤差和取樣引起的隨機變異組成；塊金系數(C0/(C0+C))表示隨機部分引起的空間異質性占系統總變異的比例[18]。利用ArcGIS 10.2.2軟件對土壤硒元素進行克里金空間插值。

2 結果與分析

2.1 掛甲峪地區不同土地利用方式土壤全硒含量比較

不同土地利用方式的土壤全硒含量見表1，林地土壤全硒含量變化范圍為0.10 ~ 0.76 mg/kg，均值為0.37 mg/kg；園地變化范圍為 0.06 ~ 0.56 mg/kg，均值為0.25 mg/kg；草地變化范圍為0.11 ~ 0.57 mg/kg，均值為0.40 mg/kg。不同土地利用方式中草地全硒含量平均值最高，其次是林地，最低為園地。全硒含量最高值出現在林地，最低值出現在園地。全硒含量變異系數為草地(0.40)>林地(0.35)>園地(0.33)，分布較均勻。

表1 掛甲峪地區不同土地利用方式土壤全硒含量比較

注：本研究參照土地利用現狀分類標準GB/T21010—2007[19]對土壤樣品進行分類。

2.2 掛甲峪地區不同海拔土壤全硒含量比較

不同海拔土壤全硒含量統計分析結果見表2，海拔在100 ~ 200 m的土壤全硒含量變化范圍為0.11 ~ 0.39 mg/kg，均值為0.23 mg/kg；海拔在200 ~ 300 m的土壤全硒含量變化范圍為0.06 ~ 0.76 mg/kg，均值為0.30 mg/kg；海拔在300 ~ 400 m的土壤全硒含量變化范圍為0.10 ~ 0.57 mg/kg，均值為0.37 mg/kg。土壤全硒平均含量表現為海拔300 ~ 400 m > 海拔200 ~300 m > 海拔100 ~ 200 m。

對掛甲峪地區土壤全硒含量與海拔做相關性分析(圖 2)，結果表明，土壤全硒含量與海拔高度呈顯著正相關關系(<0.01)，即土壤全硒含量隨著海拔增高有增加趨勢。

2.3 掛甲峪地區土壤全硒含量總體特征

由掛甲峪地區與其他地區表層土壤硒含量比較(表3)可知，該地區土壤全硒含量平均值為 0.30 mg/kg，變幅為0.06 ~ 0.76 mg/kg，平均值高于北京市平原區[11]、內蒙古武川縣[20]、黑龍江[21]、廣東[22]和三峽區(重慶段)[23]等地區，同時高于我國土壤背景值[24]和世界土壤平均硒含量[28]，但低于我國江西省豐城市[24]及香港等富硒地區[25]。

表2 掛甲峪地區不同海拔土壤全硒含量比較

圖2 表層土壤全硒含量與海拔的相關關系

表3 掛甲峪地區與其他地區表層土壤硒含量比較(mg/kg)

根據譚見安[27]的土壤硒含量分級方法，對采集的土壤樣品全硒含量劃分為5個等級(表4)，調查樣點中有20.49% 的土壤樣品全硒含量處于富硒水平，有66.39% 的土壤樣品處于足硒水平，10.66% 的土壤樣品處于潛在硒不足，2.46% 的土壤樣品處于硒含量缺乏水平，所采集樣品中未發現硒含量過剩土壤。可見，掛甲峪地區存在大面積足硒土壤和富硒土壤。

2.4 掛甲峪地區土壤硒元素的空間變異特征

利用GS+7.0對硒元素進行空間結構性分析，不同模型擬合參數見表5。不同模型中殘差(RSS)最小為0.000 008 15，決定系數(2)最大為0.912，即指數模型為硒元素的最佳半變異函數擬合模型。該模型中，塊金值為0.001 4，塊金值較小，表明采樣方案可以較好地反映土壤全硒元素的變異程度；塊金系數大于75%，表明硒元素空間相關性弱，其原有的空間分布格局受到人為因素的干擾。

采用ArcGis10.2.2軟件克里金插值方法繪制土壤全硒含量分布圖(圖3)。由圖3可知，掛甲峪地區土壤硒含量分布不均勻，整體上呈北西向帶狀分布。高硒土壤主要分布在區域西南角(林全亭的西北方向)、北側山坡上和東側山坡上，空間分布比較集中。

表4 掛甲峪地區土壤硒含量分級

表5 硒元素不同模型擬合參數結果

圖3 掛甲峪地區土壤全硒空間分布圖

3 討論

掛甲峪地區土壤全硒含量范圍為0.06 ~ 0.76 mg/kg，硒含量高的區域巖性以灰巖、泥巖、泥灰巖互層和砂巖等沉積巖為主，這些沉積巖中含有豐富的黏土礦物，有利于硒元素在土壤中吸附固定。土壤母質影響土壤硒元素的含量與分布，這與前人調查一致。王美珠和章明奎[28]及陳俊堅等[22]調查發現，土壤母質是造成自然土壤中硒含量高低的主要因素，土壤硒元素的含量和分布在很大程度上取決于成土母質的組成和性質。遲鳳琴等[21]對黑龍江省土壤硒的調查發現，土壤中豐富的黏土礦物和有機質有利于硒酸鹽等的累積。

掛甲峪地區林地土壤全硒含量平均值為0.37 mg/kg，是園地平均值的1.48倍。草地土壤硒全量平均值為0.40 mg/kg，是園地平均含量的1.60倍。說明不同土地利用方式和人類活動強弱會影響硒元素分布的不均一性和多變性[29]。林地和草地分別位于海拔較高的中低山區和地勢較平緩的地區，硒含量受人類干擾較少，生態系統中的硒通過山地侵蝕、風化、腐殖化、有機物礦化等作用向表層土壤富集，其中少量硒元素被植物吸收利用[16]，這可能是引起林地和草地土壤富集硒元素的原因。園地土壤中有機肥和化肥的定期施用會增加土壤中硒元素向水溶態轉化的速度，使得土壤中的硒元素易于向植物體遷移轉化，這可能是引起該類型土壤硒元素較低的原因。

掛甲峪地區土壤硒含量隨海拔的升高有增加趨勢，這與王曉杰等[20]、章海波等[25]、張曉平和張玉霞[30]、武少興等[31]調查結果一致。原因可能是：該地區氣候屬大陸性季風氣候，降水較少，淋溶作用較弱，硒元素的遷移作用減弱，且北、東、南三面環山，隨著海拔的增加溫度降低，有機質分解減慢，有機復合態的硒與有機質結合減少，導致植物吸收的水溶態硒和土壤中被淋溶的硒含量減少，使得硒元素在高海拔地區更易于富集。

4 結論

掛甲峪地區表層土壤硒含量高于我國和世界土壤硒含量平均值，調查點有20.49% 的土壤處于富硒水平(0.40 ~ 3.00 mg/kg)，有66.39% 的土壤處于足硒狀態(0.175 ~ 0.400 mg/kg)，僅2.46% 的土壤處于硒含量缺乏水平(<0.125 mg/kg)，未發現硒含量過剩土壤，即該地區存在大面積足硒土壤和富硒土壤。不同的土地利用方式表層土壤硒含量存在差異，草地土壤硒水平最高，林地次之，園地最低。表層土壤全硒含量隨海拔的增加而增加。掛甲峪地區硒元素整體上呈北西向帶狀分布，高硒土壤的空間分布較集中。土壤母質對全硒含量及分布的影響最大，硒元素空間變異受人為因素、土地利用方式、地形條件和人為管理方式的影響。

[1] Mao J, Pop V J, Bath S C, et al. Effect of low-dose selenium on thyroid autoimmunity and thyroid function in UK pregnant women with mild-to-moderate iodine deficiency[J]. European Journal of Nutrition, 2014, 55(1): 1–7

[2] Lenz M, Lens P N. The essential toxin: The changing perception of selenium in environmental sciences[J]. Sci. Total Environ., 2009, 407(12): 3620–3633

[3] 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[J]. Environmental Earth Sciences, 2011, 62(7): 1551–1560

[4] Wedepohl K H. The Composition of the continental crust[J]. Geochimica et Cosomochimica Acta, 1995,59(7): 1217–1232

[5] 張寶軍, 鐘松臻, 龔如雨, 等. 贛南低丘紅壤水稻土硒及其生物有效形態的組成與分布[J]. 土壤, 2017, 49(1): 150–154

[6] 廖金鳳. 土壤環境中的硒對人和動物健康的影響[J]. 廣東微量元素科學, 2002, 9(3): 20–23

[7] 周鑫斌, 賴凡, 張城銘, 等. 不同形態硒向水稻籽粒轉運途徑及品種差異[J]. 土壤學報, 2017, 54(5): 1251–1258

[8] 譚見安. 生命元素硒的地域分異與健康[J]. 中國地方病學雜志, 1996(2): 67

[9] 王海軍. 北京市大骨節病輕病區32年監測報告[J]. 中國民康醫學, 2011, 23(12): 1481–1482

[10] 戴昭華, 黃衍初. 北京、南京地區土壤中硒和鈹的分布[J]. 環境科學, 1984, 5(3): 26–28

[11] 郭莉, 楊忠芳, 阮起和, 等. 北京市平原區土壤中硒的含量和分布[J]. 現代地質, 2012, 26(5): 859–864

[12] 鄧厚培, 陳慶中. 北京郊區環境中的硒、鈣、氟與碘缺乏病關系初探[J]. 地方病通報, 1995, 10(1): 17–18

[13] 黃淇, 成杭新, 陳出新, 等. 北京市房山區富硒土壤調查與評價[J]. 物探與化探, 2013, 37(5): 889–894

[14] 李娜, 高志東, 羅德剛, 等. 北京周口店地區環境硒含量及居民血清硒水平分析[J]. 首都公共衛生, 2007, 1(1): 13–17

[15] 牛然. 不同資源類型的鄉村旅游規劃研究[D]. 北京: 首都師范大學, 2008

[16] 鄭袁明, 陳煌, 陳同斌, 等. 北京市土壤中Cr, Ni含量的空間結構與分布特征[J]. 第四紀研究, 2003, 23(4): 436–445

[17] 王幼奇, 白一茹, 王建宇. 基于GIS的銀川市不同功能區土壤重金屬污染評價及分布特征[J]. 環境科學, 2016, 37(2): 710–716

[18] 郭旭東, 傅伯杰. 河北省遵化平原土壤養分的時空變異特征:變異函數與Kriging插值分析[J]. 地理學報, 2000, 1(5): 555–566

[19] 國土資源部. 土地利用現狀分類(GB/T 21010—2007)[S]. 北京: 中國標準出版社, 2007

[20] 王曉杰, 孟凡喬, 吳文良. 內蒙古武川縣土壤硒分布特性研究[J]. 土壤通報, 2016, 47(3): 624–629

[21] 遲鳳琴, 徐強, 匡恩俊, 等. 黑龍江省土壤硒分布及其影響因素研究[J]. 土壤學報, 2016, 53(5): 1262–1274

[22] 陳俊堅, 張會化, 余煒敏, 等. 廣東省土壤硒空間分布及潛在環境風險分析[J]. 生態環境學報, 2012, 21(6): 1115– 1120

[23] 羅友進, 韓國輝, 孫協平, 等. 三峽庫區(重慶段)土壤硒分布特征及影響因素[J]. 土壤, 2018, 50(1): 131–138

[24] 魏復盛, 陳靜生. 中國土壤環境背景值研究[J]. 環境科學, 1991, 12(4): 12–19

[25] 章海波, 駱永明, 吳龍華, 等. 香港土壤研究Ⅱ.土壤硒的含量、分布及其影響因素[J]. 土壤學報, 2005, 42(3): 404–410

[26] Lisk D J. Trace metals in soils, plants, and animals[J]. Advances in Agronomy,1972, 24: 267–325

[27] 譚見安. 環境生命元素與克山病[M]. 北京: 中國醫藥科技出版社, 1996

[28] 王美珠, 章明奎. 我國部分高硒低硒土壤的成因初探[J]. 浙江農業大學學報, 1996, 22(1): 89–93

[29] 戴慧敏, 宮傳東, 董北, 等. 東北平原土壤硒分布特征及影響因素[J]. 土壤學報, 2015, 52(6): 1356–1364

[30] 張曉平, 張玉霞. 西藏土壤中硒的含量及分布[J]. 土壤學報, 2000, 37(4): 558–562

[31] 武少興, 龔子同, 黃標. 我國土壤中的溶態硒含量及其與土壤理化性質的關系[J]. 中國環境科學, 1997, 17(6): 43–46

Distribution Pattern of Soil Total Selenium in Guajiayu Area of Beijing

CUI Meng1, SUN Xiangyang1*, JIANG Xinfu2, LI Suyan1, YANG Shaobin1, LI Song1, GONG Xiaoqiang1, ZHANG Le1

(1 College of Forestry, Beijing Forestry University, Beijing 100083, China; 2 Forest Seed and Seedling Service Center of Tongzhou District in Beijing, Beijing 101100, China)

To illustrate the distribution characteristics of soil total selenium (Se), 122 samples of topsoils (0–20 cm) were collected from Guajiayu area in Beijing to analyze the spatial distribution characteristics of soil total Se content with GIS and its relationship with altitude and land use type. The results showed that total Se content ranged from 0.06 to 0.76 mg/kg with the average value of 0.30 mg/kg. 86.88% of soil samples belonged to Se-sufficient and above levels (>0.175 mg/kg), among of which 20.49% of soil samples belonged to Se-enriched level (0.40 – 0.30 mg/kg). Grassland had the highest total Se content (0.40 mg/kg), followed by woodland (0.37 mg/kg), and the lowest was orchard (0.25 mg/kg). Soil total Se content had a significant positive correlation with altitude. Furthermore, spatial variability indicated that spatial correlation of Se element was weak, with the nugget (0) and the coefficient of nugget (0/(0+)) were 0.001 4 and 88.10%, respectively. The distribution of soil total Se generally showed a northwest zonal distribution trend, and Se-rich soils are obviously existed in the area.

Guajiayu; Spatial correlation; Selenium(Se)

中國工程院咨詢研究項目(2014-XZ-11-02)和北京市科技計劃項目( Z161100001116061 )資助。

(sunxy@bjfu.edu.cn)

崔萌(1993—)，女，河北萬全人，碩士研究生，研究方向為土壤修復與健康。E-mail：974538322@qq.com

S159.2；X53

A

10.13758/j.cnki.tr.2019.04.026