蘭州某鋼廠周邊土壤中磁化率積累速率研究

呂玉甲，戴　霜，2*，張有賢，2，栗澤苑，丁　雷（.蘭州大學(xué)資源環(huán)境學(xué)院，甘肅 蘭州 730000；2.蘭州大學(xué)西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，甘肅 蘭州 730000）

蘭州某鋼廠周邊土壤中磁化率積累速率研究

呂玉甲1，戴霜1，2*，張有賢1，2，栗澤苑1，丁雷1（1.蘭州大學(xué)資源環(huán)境學(xué)院，甘肅 蘭州 730000；2.蘭州大學(xué)西部環(huán)境教育部重點(diǎn)實(shí)驗(yàn)室，甘肅 蘭州 730000）

為了尋找一種簡(jiǎn)便快捷的積累速率的研究方法，通過(guò)對(duì)蘭州某鋼廠近地表土壤進(jìn)行高密度分層采樣和磁化率測(cè)量，發(fā)現(xiàn)各土層磁化率值隨著采樣深度的逐漸增加而迅速減小，并且這種關(guān)系可以用可靠度很高的方程式來(lái)定量描述.在此基礎(chǔ)上，通過(guò)對(duì)該方程式積分處理，計(jì)算了土層中磁化率的新增積累量，然后結(jié)合鋼廠的排污年限計(jì)算了磁化率在各土層中的積累速率.結(jié)果表明磁化率值在1cm以上的土層中富集程度較高，約占整個(gè)土層新增磁化率總積累量的70.24%～97.56%；1cm深度以下土層中磁化率值的變化幅度較小；各采樣點(diǎn)新增磁化率總積累速率分別為6.5472×10-8m3/（kg·a），40.4178×10-8m3/（kg·a），49.7683×10-8m3/（kg·a），31.8679×10-8m3/（kg·a），且地表1cm以上土層的積累速率最大，向下直至5cm迅速變小.同時(shí)發(fā)現(xiàn)下風(fēng)向的積累速率約為上風(fēng)向的4.87～7.60倍.因此，利用土層磁化率值縱向變化規(guī)律來(lái)計(jì)算積累速率的方法是可行、可信的，可以反映土壤中重金屬的動(dòng)態(tài)積累過(guò)程.

磁化率；重金屬；積累速率；土壤污染；鋼廠

鋼鐵廠在生產(chǎn)過(guò)程中排放的污染物會(huì)對(duì)周邊的環(huán)境造成較大的危害，特別是重金屬污染物會(huì)通過(guò)大氣沉降等方式富集在土壤表層，并在雨水下滲等作用下逐漸遷移至土層較深處，從而導(dǎo)致土壤質(zhì)量嚴(yán)重受損，并對(duì)當(dāng)?shù)氐纳鷳B(tài)環(huán)境造成較大破壞.研究重金屬在土壤中的分布和動(dòng)態(tài)積累過(guò)程，對(duì)于認(rèn)識(shí)土壤重金屬污染和環(huán)境修復(fù)具有重要意義.目前，國(guó)內(nèi)外已有許多針對(duì)鋼廠周邊土壤重金屬污染監(jiān)測(cè)及治理方面研究成果.但是這些研究主要集中在土壤重金屬污染的空間分布［1］、遷移規(guī)律［2］及污染評(píng)估［3］等方面.而重金屬在土壤中積累速率的相關(guān)報(bào)道較少.王國(guó)梁等［4］曾采用空間序列代替時(shí)間序列的方法研究了土壤中不同重金屬隨著作物種植年限增加而變化的過(guò)程和發(fā)展趨勢(shì).Rohit等［5］以磚塊為載體研究了自1910年以后印度阿格拉地區(qū)土壤中銅元素的積累變化過(guò)程.這些方法需要大量的監(jiān)測(cè)數(shù)據(jù)，且監(jiān)測(cè)周期較長(zhǎng).而我國(guó)環(huán)境監(jiān)測(cè)相對(duì)國(guó)外起步較晚，數(shù)據(jù)積累較少，特別是大中型鋼廠對(duì)周邊土壤環(huán)境影響的數(shù)據(jù)十分罕見(jiàn).另外，依賴于化學(xué)分析的傳統(tǒng)方法在大面積布點(diǎn)、多層次采樣、快捷動(dòng)態(tài)地監(jiān)測(cè)土壤重金屬污染狀況方面存在不足.因此，尋求簡(jiǎn)便快捷的污染監(jiān)測(cè)手段，適合的研究對(duì)象，并結(jié)合已有的研究成果對(duì)重金屬積累速率的研究方法及思路進(jìn)行優(yōu)化，對(duì)構(gòu)建鋼廠對(duì)周邊環(huán)境污染的過(guò)程具有十分重要的意義.

環(huán)境磁學(xué)在重金屬污染研究方面具有簡(jiǎn)便、快捷、經(jīng)濟(jì)、靈敏度高和非破壞性等特點(diǎn)［6-7］.磁化率是土壤中磁性礦物類型、含量、大小的實(shí)驗(yàn)室表征參數(shù)［8-9］，是對(duì)土壤受重金屬污染狀況的總體反映.目前，國(guó)內(nèi)外有許多關(guān)于鋼廠排放物環(huán)境磁學(xué)效應(yīng)的研究以及利用該方法代替?zhèn)鹘y(tǒng)化學(xué)分析方法來(lái)反映土壤中重金屬污染程度可行性的報(bào)道.倪劉建等［10］研究發(fā)現(xiàn)鋼廠工業(yè)燃煤、冶煉等過(guò)程釋放的氣態(tài)污染物、工業(yè)粉塵中磁性礦物和重金屬元素組合基本一致，這反映了土壤中Cu，F(xiàn)e，Pb，Mn和Zn等元素與磁性礦物具有類似的來(lái)源或地球化學(xué)行為.段雪梅等［11］通過(guò)研究鋼廠周邊受污染土壤的縱剖面之后，認(rèn)為磁指標(biāo)（χ，ARM和SIRM）與重金屬（Ni，Cu，F(xiàn)e，Pb，V和Zn）呈顯著相關(guān)關(guān)系.Lu等［12］在研究了杭州市“城郊污染梯度”之后，認(rèn)為環(huán)境磁學(xué)可以作為一種替代方法監(jiān)測(cè)土壤中重金屬含量.閆倩等［13］之前在本研究的鋼廠周邊采樣，對(duì)比分析后發(fā)現(xiàn)，在 α=0.01的顯著性水平下，各重金屬含量和磁化率值之間呈十分顯著相關(guān)關(guān)系，相關(guān)系數(shù)分別為0.837（Pb），0.857（Zn），0.783（Cu），0.951（Fe），0.882（Mn），其據(jù)此提出用磁化率濃集因子（S）代替重金屬含量來(lái)評(píng)價(jià)土壤受污染程度.這表明磁化率指標(biāo)可以在一定程度上代替重金屬含量，來(lái)反映土壤受污染程度.總之，近年來(lái)有很多利用磁化率來(lái)監(jiān)測(cè)土壤重金屬污染狀況可行性的研究，以及環(huán)境磁學(xué)應(yīng)用于金屬元素在不同環(huán)境介質(zhì)中遷移擴(kuò)散的研究成果［14-21］，都為本研究利用磁化率積累速率來(lái)反映重金屬動(dòng)態(tài)積累過(guò)程提供了理論基礎(chǔ).

本文所研鋼廠位于蘭州市東北郊一個(gè)小型盆地內(nèi).前人對(duì)該區(qū)域的研究表明土壤的磁化率值隨深度向下呈明顯的減小趨勢(shì)［13］，但沒(méi)有對(duì)這種變化趨勢(shì)進(jìn)行進(jìn)一步的分析.本文通過(guò)進(jìn)一步的采樣分析，旨在建立磁化率隨深度變化的定量表達(dá)式，探討利用這種定量關(guān)系來(lái)描述污染物在土壤中分布的情況，并為積累速率的研究提供一種新的研究思路.

1　樣品采集及測(cè)定

鋼廠所在區(qū)域位于蘭州市東北方向約 25km處，年降水量為330mm左右.自2003年12月正式投產(chǎn)以來(lái)年產(chǎn)鋼鐵190余萬(wàn)t、型材95萬(wàn)t.廠區(qū)周邊均為荒地和農(nóng)田；無(wú)其他較大重金屬排放源；表層土為馬蘭黃土，土層結(jié)構(gòu)和成分較為單一.

樣品全部來(lái)源于鋼廠周邊的山坡及殘存黃土梁峁上.考慮到該地區(qū)常年盛行東風(fēng)，夏秋季以東南風(fēng)為主，降水主要集中在夏秋兩季，在上風(fēng)向設(shè)置采樣點(diǎn)A1，在下風(fēng)向設(shè)置采樣點(diǎn)A2、A3、A4，具體見(jiàn)圖1.

每個(gè)采樣點(diǎn)的樣品采集深度為 15cm，采樣方法為刻槽采樣，0～1cm為不等間距采樣，分別將0～0.2cm、0.2～0.5cm及0.5～1cm的土層作為一個(gè)樣品；1～2cm以0.5cm為間距分別采樣；2～10cm以1cm為間距分別采樣；最后以 1cm為間距采集14～15cm.每個(gè)采樣點(diǎn)每次共采集14個(gè)樣品.土樣顏色自上而下逐漸由灰黑色變?yōu)橥咙S色.于2013 年4月份分別對(duì)A1、A2、A3及A4點(diǎn)進(jìn)行了樣品采集，采得的樣品裝入密封袋帶回實(shí)驗(yàn)室.樣品在實(shí)驗(yàn)室中用烘箱烘干，然后用瑪瑙研缽研細(xì)后裝于2cm×2cm×2cm樣品盒（聚乙烯材質(zhì)）中以備檢測(cè)，稱重后利用英國(guó) Bartington儀器公司生產(chǎn)的MS2磁化率儀測(cè)得樣品的低頻磁化率，最后通過(guò)計(jì)算可得樣品的質(zhì)量磁化率.

圖1　采樣布點(diǎn)（地形據(jù)文獻(xiàn)［13］修改）Fig.1　Sketch map showing the location of the steel plant and samples （landforms following reference［13］）

2　數(shù)據(jù)處理

2.1最小可靠度確定

線性、非線性回歸分析是科學(xué)研究中十分常用的數(shù)據(jù)統(tǒng)計(jì)方法，即利用最小二乘法的原理定量地描述一組數(shù)據(jù)中不同變量之間相互作用的變化關(guān)系.據(jù)此，若自變量的數(shù)值發(fā)生改變，則可給出因變量的預(yù)測(cè)值.但是僅僅給出預(yù)測(cè)值，在統(tǒng)計(jì)學(xué)上是沒(méi)有意義的，還必須給出某一顯著性水平α下的可靠度.當(dāng)給定顯著性水平α?xí)r，其對(duì)應(yīng)某一確定的Fα值.α一般取0.05或0.01，1-α表示數(shù)據(jù)的可靠度.若 F＜F0.05，則稱兩變量之間沒(méi)有明顯的關(guān)系，回歸方程不可信；若F0.01＞F＞F0.05，則稱變量之間有顯著關(guān)系；若 F＞F0.01，則稱兩變量之間有十分顯著的關(guān)系.

但是重金屬在土壤中的積累速率較為緩慢，影響因素較多，并且考慮到采樣、檢測(cè)等環(huán)節(jié)不可避免地存在誤差，所以顯著性水平α取0.01就顯得誤差較大.為進(jìn)一步減小預(yù)測(cè)值與真實(shí)值之間的誤差，參考飲用水中金屬元素的檢測(cè)方法來(lái)確定最小可信度.

在利用分光光度法［22］檢測(cè)水中金屬元素含量時(shí)，若以5個(gè)標(biāo)準(zhǔn)液擬合所得的標(biāo)準(zhǔn)曲線的相關(guān)系數(shù)r≥0.999則認(rèn)為預(yù)測(cè)值與真實(shí)值之間的誤差可以忽略不計(jì).經(jīng)計(jì)算［23］，當(dāng)r=0.999其顯著性水平α=0.000037955，可靠程度為99.9962%.雖然實(shí)驗(yàn)室相較于野外環(huán)境具有更高的可控性，但是考慮到本文數(shù)據(jù)處理的實(shí)際需要，最小可靠程度最終確定為99.9999%.

2.2樣品擬合分析

圖2　各樣點(diǎn)樣品數(shù)據(jù)擬合曲線Fig.2　Sample data fitting curve

利用軟件 Orign8.0，對(duì)數(shù)據(jù)進(jìn)行非線性回歸分析.分別用方程y=a-b· cx（1）及方程y=a（1-e-bx）c（2）對(duì)原始數(shù)據(jù)進(jìn)行擬合處理（其中a、b、c為擬合參數(shù)），所得結(jié)果見(jiàn)圖2.通過(guò)相關(guān)系數(shù)來(lái)判斷，式（2）對(duì)各樣點(diǎn)樣品數(shù)據(jù)的擬合效果都要略優(yōu)于式（1），并且在兩種擬合方式中，下風(fēng)向采樣點(diǎn)的擬合效果都要優(yōu)于上風(fēng)向采樣點(diǎn)，各樣點(diǎn)的磁化率數(shù)據(jù)幾乎都分布在擬合曲線上，相關(guān)系數(shù)都大于0.995，A2點(diǎn)的相關(guān)系數(shù)甚至都大于0.999，而相比之下A1點(diǎn)則都小于0.99.這可能是由于A1點(diǎn)位于鋼廠的上風(fēng)向位置不利于鋼廠排放的污染物沉降，且該地區(qū)降水主要集中在夏秋季節(jié)，強(qiáng)烈的下滲作用容易產(chǎn)生斷層等不均勻遷移現(xiàn)象，從而使其擬合效果較差.

利用 SPSS軟件，對(duì)回歸分析所得方程進(jìn)行顯著性檢驗(yàn)，主要信息見(jiàn)表1.其中，顯著性水平α 取 0.000001，其所對(duì)應(yīng)的 F0.000001（3，11）為 50.63，其遠(yuǎn)小于兩種擬合方程中A1點(diǎn)的F因子，當(dāng)然更小于其他各點(diǎn).因此各樣點(diǎn)擬合方程的可靠度都大于99.9999%.

綜上所述，式（1）及式（2）的可靠度都滿足統(tǒng)計(jì)學(xué)的要求，其曲線上的預(yù)測(cè)值可以代表土壤磁化率值的真實(shí)含量，這在一定程度上也反映了重金屬在土壤中縱向空間的分布情況.

表1　回歸分析參數(shù)及顯著性檢驗(yàn)Table 1　Regression analysis parameters and significance test

2.3磁化率在土層中的累積

采樣區(qū)域土層中土樣的磁化率值隨深度呈明顯的指數(shù)分布，且可以用可靠度很高的方程式來(lái)定量描述，所以通過(guò)對(duì)該方程式進(jìn)行積分處理來(lái)計(jì)算土層中磁化率的累積量（即橫、縱坐標(biāo)及擬合曲線所圍的面積）是可行、可信的.

土層中新增磁化率累積量可用式（3）計(jì)算：

式中：Q為新增磁化率積累總量；h1、h2表示土層的不同深度；F（x）表示磁化率值隨土層深度變化的定量關(guān)系式，m為背景值，本文以各采樣點(diǎn)土層中14～15cm深度處的磁化率值為準(zhǔn).

雖然式（2）較式（1）具有更高的可靠度，但是由于其形式較為復(fù)雜，在利用 MATLAB處理時(shí)不易得到其原函數(shù)，也不利于后續(xù)的進(jìn)一步處理，所以相比較而言式（1）具有更高的實(shí)用價(jià)值.所以本文以式（1）來(lái)計(jì)算新增磁化率積累量.推導(dǎo)過(guò)程如下：

表2　不同采樣點(diǎn)新增磁化率積累量（×10-8m3/kg）隨深度變化Table 2　Depth-dependent accumulation of magnetic susceptibility on the sampled sections （×10-8m3/kg）

利用式（4），以1cm為間隔分別求出0～5cm土層及整個(gè)土層（0～5cm）中新增磁化率累積量，計(jì)算結(jié)果詳見(jiàn)表2.另外由圖2和表2可知，A4點(diǎn)較A2點(diǎn)更接近鋼廠但其磁化率值卻小于A2點(diǎn)，這主要是由于采樣點(diǎn)的布置受到地形的限制，鋼廠位于河谷盆地內(nèi)，采樣點(diǎn)應(yīng)盡量設(shè)置在當(dāng)?shù)鼐用窕顒?dòng)較少的山坡或者高地上，以避免汽車尾氣及礦渣等對(duì)于采樣點(diǎn)土壤的影響.A4位于一個(gè)比較狹窄的谷地，常年有風(fēng)活動(dòng)，可能對(duì)粉塵沉降有影響，因而A4的磁化率值比更遠(yuǎn)處的A2低.

2.4磁化率值在土層中的累積速率

盡管在某些環(huán)境介質(zhì)（冰芯、湖泊和淺海沉積物）中重金屬元素在100a以上尺度的變化趨勢(shì)呈指數(shù)關(guān)系［24］.但是鋼廠排放的顆粒污染物經(jīng)過(guò)大氣沉降到達(dá)土壤表層，然后經(jīng)過(guò)降水下滲，土壤阻滯等作用逐漸縱向遷移至較深土層.以土壤作為介質(zhì)的縱向遷移是較為緩慢的，特別是當(dāng)積累尺度較短時(shí)，可以將重金屬元素在土層的積累近似地視為一種線性積累的過(guò)程［25］.所研鋼廠正式投產(chǎn)時(shí)間較短，其所排放的污染物在周邊土壤中的積累也應(yīng)滿足這種線性關(guān)系.即 Q=k·t，其中 Q為土層中新增磁化率積累總量；k為積累速率；t為積累時(shí)間.計(jì)算結(jié)果詳見(jiàn)表3.

表3　磁化率在土壤中的積累速率［×10-8m3/（kg·a）］Table 3　The accumulation rate of magnetic mineral inevery soil layer［×10-8m3/（kg·a）］

3　結(jié)果與討論

3.1磁化率隨土壤深度的變化

由圖3可知，各樣點(diǎn)樣品的磁化率值變化趨于一致.在5～15cm深度，磁化率值幾乎保持不變，平均值為 34.3×10-8～57.1×10-8m3/kg，代表了本地磁化率背景值.在1～5cm土層中，1cm處磁化率值是5cm處的1.13～1.71倍，說(shuō)明在深度1～5cm土層中，隨著深度增大，磁化率值逐漸降低.而0～1cm土層中，表層土的磁化率值達(dá)到最大，為107×10-8～716×10-8m3/kg，但隨著深度的增加，磁化率值驟然減小，表層土約為5cm以下土層磁化率值的3.12～12.54倍；且該段土壤為灰黑色，與其下（深度5cm以下）棕黃色-土黃色土壤層顏色迥然不同，這種磁化率值和顏色的顯著差異說(shuō)明磁性礦物主要富集在表層1cm的土層中.

圖3　各樣點(diǎn)樣品磁化率值隨深度變化趨勢(shì)Fig.3　Every point the sample magnetic susceptibility values with depth change trend

由表2可知，各采樣點(diǎn)0～1cm土層中新增磁化率 總 累積 量 達(dá)41.3911×10-8～420.7614× 10-8m3/kg，約占新增磁化率總累積量的 70.24%～97.56%，當(dāng)土層深度大于1cm時(shí)，各土層（以1cm為間隔）中新增磁化率積累量變化較小，這也說(shuō)明磁性礦物主要集中在1cm以上的土層中，并且該現(xiàn)象在下風(fēng)向樣點(diǎn)中更為明顯.下風(fēng)向樣點(diǎn)0～1cm土層中新增磁化率積累量占新增磁化率總累積量的比重平均為上風(fēng)向樣點(diǎn)的1.36倍.這主要是因?yàn)?A1點(diǎn)位于上風(fēng)向位置，不利于鋼廠排放的顆粒物沉降積累，在相同降雨強(qiáng)度的作用下，有較多磁性礦物遷移至土層較深處.

綜上所述，通過(guò)定性和定量的方法來(lái)研究磁化率值在土壤中的縱向變化情況，發(fā)現(xiàn)磁性礦物主要集中在地表1cm以上的土層中，這一研究結(jié)果與閆倩等［13］之前對(duì)該區(qū)域的研究結(jié)論有所區(qū)別，這主要是由于其是通過(guò)定性描述的方法來(lái)研究磁化率值在土壤中的變化情況，誤差較大；1～5cm深度的土層受鋼廠的干擾較?。?～15cm深度的土層基本不受干擾，代表了本地磁化率背景值.這種上高下低的分布模式與淋溶模式中污染物上低下高的分布明顯不同，說(shuō)明土壤表層磁化率值較高不是自然因素的結(jié)果，而是由于鋼廠在生產(chǎn)過(guò)程中持續(xù)排放的污染物所造成的.

3.2土層中磁化率累積速率

由表3可見(jiàn)，各采樣點(diǎn)整個(gè)土層（0～5cm）中新增磁化率積累速率依次為6.5472×10-8m3/ （kg·a），40.4178×10-8m3/（kg·a），49.7683×10-8m3/（kg· a）及31.8679×10-8m3/（kg·a）；其中2cm以下的土層中，各土層（以1cm為間隔）磁化率的年累積速率幾乎保持不變，平均為 0.2458×10-8～0.7513× 10-8（kg·a），1～2cm 土層中磁化率累積速率為0.2496×10-8～0.7633×10-8（kg·a），0～1cm 深度土層的累積速率最大，可達(dá) 4.5990×10-8～46.7513× 10-8m3/（kg·a）.

積累速率隨土層深度變化的趨勢(shì)與磁化率縱向空間分布狀況都呈明顯的指數(shù)模式，這可能是由于鋼廠排放的顆粒物沉降到地表后，其在降水的作用下有了向下遷移的動(dòng)力，但是在土壤本身機(jī)械截留的作用下，較大顆粒物質(zhì)迅速沉積在0～1cm土層中，而較小顆粒及可溶性物質(zhì)隨著降水滲透到土層較深處，并在有機(jī)質(zhì)固定及膠體吸附等［26］作用下被逐漸固定在土層中.處于下風(fēng)向采樣點(diǎn)的整個(gè)土層（0～5cm）中新增磁化率積累速率為上風(fēng)向采樣點(diǎn)積累速率的4.87～7.60倍，表層（0～1cm）積累速率可達(dá)6.58～10.17倍.這表明風(fēng)向不僅會(huì)直接影響磁化率的橫向空間分布，也會(huì)間接影響其在土壤中的縱向空間分布.

綜上所述，A3點(diǎn)的積累速率最大，A1點(diǎn)的積累速率最小，在同一采樣點(diǎn)中表層（0～1cm）土壤的積累速率最大；風(fēng)向不僅影響磁性礦物的積累速率，還影響其空間分布.

4　結(jié)論

4.1蘭州某鋼廠近地表土層磁化率較背景值有明顯的積累，且磁化率主要集中在 1cm以上的土層中，其約占整個(gè)土層磁化率積累量的70.24%～ 97.56%，在1～5cm土層中，磁化率值逐漸降低，在5～15cm土層中，磁化率值基本保持不變，代表了本地磁化率背景值；各采樣點(diǎn)的總積累 速 率 分 別6.5472×10-8m3/（kg·a），40.4178× 10-8m3/（kg·a），49.7683×10-8m3/（kg·a），31.8679× 10-8m3/（kg·a）.其中以 A3點(diǎn)的積累速率最大，A1點(diǎn)的積累速率最小，且縱向變化過(guò)程中以地表1cm土層的積累速率最大.下風(fēng)向的積累速率約為上風(fēng)向的4.87～7.60倍.

4.2該鋼廠近地表土層中磁化率值隨著采樣深度的逐漸增加而迅速減少，并且這種關(guān)系可以用可靠度高于99.9999%的方程式來(lái)定量描述.因此，利用土壤磁化率值縱向變化規(guī)律來(lái)計(jì)算積累速率的方法是可行、可信的，可以反映土壤中重金屬的動(dòng)態(tài)積累過(guò)程.

［1］ 陳曉晨，崔巖山.城市表層土壤中重金屬的小尺度空間分布——以首鋼廠區(qū)附近小區(qū)域?yàn)槔?［J］. 中國(guó)科學(xué)院研究生院學(xué)報(bào)，2010，27（2）：176-183.

［2］ 薛粟尹，李萍，王勝利，等.工礦型綠洲農(nóng)田土壤中氟和鎘的剖面行為研究 ［J］. 中國(guó)環(huán)境科學(xué)，2014，34（1）：195-200.

［3］ 肖雪，趙南京，袁靜，等.淮河流域某鎮(zhèn)農(nóng)業(yè)土壤重金屬含量特征及污染評(píng)價(jià) ［J］. 光譜學(xué)與光譜分析，2014，34（7）：1785-1788.

［4］ 王國(guó)梁，周生路，趙其國(guó)，等.菜地土壤剖面上重金屬元素含量隨時(shí)間的變化規(guī)律研究 ［J］. 農(nóng)業(yè)工程學(xué)報(bào)，2006，22（1）：79-84.

［5］ Rohit S，Sanjay K M，Shobhit S，et al. Bricks as histor-ical record of metals fallout： study on copper accumulation in Agra soils since 1910. Environmental Monitoring and Assessment，1996，40：271-278.

［6］ Schmidt A，Yarnold R，Hill M，et al. Magnetic suscepti-bility as proxy for heavy metal pollution： a site study ［J］. Jo-urnal of Geochemical Exploration，2005，85：109-117.

［7］ Boyko T，Scholger R，Stanjek H，et al. Topsoil magnetic susceptibility mapping as a tool for pollution monitoring：r-epeatability of in situ measurements ［J］. Journal of Applied Geophysics，2004，55：249-259.

［8］ 王君蘭，李暉，鄧偉，等.內(nèi)蒙古嘎順諾爾湖泊沉積物磁化率與粒度的古環(huán)境意義 ［J］. 中國(guó)沙漠，2012，25（3）：661-668.

［9］ 王建，劉澤純，姜文英.磁化率與粒度、礦物的關(guān)系及其古環(huán)境意義 ［J］. 地理學(xué)報(bào)，1996，51（2）：155-163.

［10］ 倪劉建，張甘霖，楊金玲，等.鋼鐵工業(yè)區(qū)降塵對(duì)周邊土壤的影響［J］. 土壤學(xué)報(bào)，2007，44（4）：637-642.

［11］ 段雪梅，胡守云，閆海濤，等.南京某鋼鐵公司周邊耕作土壤的磁學(xué)性質(zhì)與重金屬污染的相關(guān)性研究 ［J］. 地球科學(xué)，2009，39（9）：1304-1312.

［12］ Lu Shenggao，Wang Hongyan，Bai Shiqiang. Heavy me-tal contents and magnetic susceptibility of soils along urbanr-ural gradient in rapidly growing city of eastern china ［J］. En-vironmental Monitoring and Assessment，2009，155：91-101.

［13］ 閆倩，戴霜，劉海嬌，等.蘭州某鋼廠附近土壤磁化率特征及其環(huán)境意義 ［J］. 中國(guó)環(huán)境科學(xué)，2011，31（10）：1732-1737.

［14］ Thompson R，Oldfield F. Environmental Magnetism ［M］. London：Allen and Unwin，1986：227.

［15］ Mayer T，Morris W A，Versteeg K J. Feasibility of usi-ng magnetic properties for assessment of particle-associated contaminant transport ［J］. Water Quality Research Journal，1996，31（4）：741-752.

［16］ Maher B A，Thompson R. Quaternary climates，environments and magnetism ［M］. Cambridge： Cambridge University Press，1999：279-322.

［17］ Chester R，Sharples E J，Sanders G，et al. The distribu-tion of nature and non-crustal ferromagnetic minerals in soil-size particulates from the Mediterranean atmosphere ［J］. Water Air Soil Pollution，1984，23：25-35.

［18］ Hunt A，Jones J，Oldfield F. Magnetic measurements and heavy metals in atmospheric particulates of anthropogenic origin ［J］. Science of the Total Environment，1984，33：129-139.

［19］ Hoffmann V，Knab M，Apple E. Magnetic susceptibility napping of roadside pollution ［J］. Journal of Geochemical Exploration，1999，66：313-326.

［20］ Evans M E，Heller F. Environmental magnetism： Princi-ples and applications of enviromagnetics ［M］. Oxford （UK）： Academic press，2003.

［21］ 閆海濤，朱育新，胡守云，等.西笤溪流域重金屬污染磁學(xué)研究及多元分析 ［J］. 中國(guó)環(huán)境科學(xué)，2004，24（4）：385-389.

［22］ 國(guó)家環(huán)?？偩?水和廢水監(jiān)測(cè)分析方法 ［M］. 第四版（增補(bǔ)版），北京，中國(guó)環(huán)境出版社，2002.

［23］ 邱軼兵，張文利，閔凡飛，等.試驗(yàn)設(shè)計(jì)與數(shù)據(jù)處理 ［M］. 合肥：中國(guó)科學(xué)技術(shù)大學(xué)出版社，2008，185-195.

［24］ 李戀卿，潘根興，張平究，等.太湖地區(qū)水稻土表層土壤10年重金屬元素積累的估計(jì) ［J］. 環(huán)境科學(xué)，2002，23（3）：119-123.

［25］ 成杰民，張麗娜.近 20年山東省典型農(nóng)田土壤中重金屬含量變化特征及積累速率估計(jì) ［J］. 土壤，2013，45（1）：99-104.

［26］ 阮心玲，張甘霖，趙玉國(guó)，等.基于高密度采樣的土壤重金屬分布特征及遷移速率 ［J］. 環(huán)境科學(xué)，2006，27（5）：1020-1025.

致謝：本實(shí)驗(yàn)的現(xiàn)場(chǎng)采樣工作由馬素輝、張玉麗、張佳佳、甄天坷等同學(xué)協(xié)助完成，在此表示感謝.

Magnetic susceptibility accumulation rates in soils near a steel smelter in Lanzhou.

Lü Yu-jia1，DAI Shuang1，2*，ZHANG You-xian1，2，LI Ze-yuan1，DING Lei1（1.College of Resources and Environment，Lanzhou University，Lanzhou 730000，China；2.Key Laboratory of Western China's Environmental Systems，Ministry of Education，Lanzhou University，Lanzhou 730000，China）.

China Environmental Science，2015，35（7）：2112～2118

In order to establish an easy and fast method to investigate the accumulation rate of heavy metals in soil，samples were collected from the soils on 4sections adjacent a steel smelter near Lanzhou city and their magnetic susceptibility were measured. The data show that the magnetic susceptibilities fast decrease on the section downward. This variation trend can be formulated by a high reliability equation. Using the equation after integrating processing，the newly increased accumulation mass of magnetic susceptibility of soils for each section were calculated. Then the accumulation rate of magnetic susceptibility for each soil layer of the section was calculated according to the running time of the steel smelter. The result shows that the magnetic mineral mainly concentrated in the topmost 1cm soil layer，which accounts for about 70.24%～97.56% of the newly increased accumulation mass of magnetic susceptibility in soil. The magnetic susceptibilities of the soils below 1cm on the section are slightly various. The newly increased accumulation rate of magnetic susceptibility for 4sections are respectively 6.5472×10-8m3/（kg·a），40.4178×10-8m3/（kg·a），49.7683× 10-8m3/（kg·a） and 31.8679×10-8m3/（kg·a）. The value is maximum in the topsoil above 1cm depth，and it decreases fast down to the 5cm depth on each section. In horizontal，the accumulation rate of magnetic susceptibility of soils in the downwind area is as 4.87 to 7.60 times as that in the upwind area. Therefore，it is feasible and reliable to obtain the newly increased accumulation rate of soil using the variation trend of magnetic susceptibility on the section. This method can reflect the continuous accumulation process of the heavy metals in the soils.

magnetic susceptibility；heavy metal；accumulation rate；soil pollution；steel smelter

X131.3

A

1000-6923（2015）07-2112-07

2014-12-02

國(guó)家基礎(chǔ)科學(xué)人才培養(yǎng)基金（J1210065）；國(guó)家自然科學(xué)基金（41272127）.

* 責(zé)任作者，教授，daisher@lzu.edu.cn

呂玉甲（1992-），男，甘肅靖遠(yuǎn)人，蘭州大學(xué)資源環(huán)境學(xué)院碩士研究生，主要從事土壤污染研究與治理.