凍融期土壤溫度對有機污染物遷移行為的影響

崔 嵩, 付 強, 李天霄, 李一凡

(1.東北農(nóng)業(yè)大學 水利與建筑學院, 哈爾濱 150030; 2.東北農(nóng)業(yè)大學 國際持久性有毒物質(zhì)聯(lián)合研究中心,哈爾濱 150030; 3.哈爾濱工業(yè)大學 國際持久性有毒物質(zhì)聯(lián)合研究中心, 哈爾濱 150090)

?

凍融期土壤溫度對有機污染物遷移行為的影響

崔 嵩1,2, 付 強1,2, 李天霄1,2, 李一凡3

(1.東北農(nóng)業(yè)大學 水利與建筑學院, 哈爾濱 150030; 2.東北農(nóng)業(yè)大學 國際持久性有毒物質(zhì)聯(lián)合研究中心,哈爾濱 150030; 3.哈爾濱工業(yè)大學 國際持久性有毒物質(zhì)聯(lián)合研究中心, 哈爾濱 150090)

為了解季節(jié)性凍融期土壤垂直分層溫度變化對有機污染物遷移擴散行為的影響。通過冬季大田試驗分別測定了裸地，5，10，15 cm厚度秸稈覆蓋條件下3，10，20，40，60，100 cm深度土壤溫度，并應用土—氣交換層及相鄰土壤層遷移擴散通量系數(shù)，深入研究了土壤溫度變化對有機污染物遷移行為的影響，并系統(tǒng)地分析了不同溫度數(shù)據(jù)對模擬結果產(chǎn)生的差異性。研究結果表明：不同覆蓋條件下深層土壤有機氯農(nóng)藥α-HCH遷移擴散通量系數(shù)數(shù)值的變化與土壤層厚度有關；土壤垂直分層溫度的變化與α-HCH的遷移擴散通量系數(shù)呈顯著正相關，且隨著土壤深度的增加相關性逐漸減弱；日平均氣溫變化與土—氣交換層α-HCH的遷移擴散通量系數(shù)的變化具有明顯的一致性(R=0.999，p=0.000)；日平均氣溫與土壤垂直分層溫度觀測數(shù)據(jù)建立的回歸方程能夠較好地與0—40 cm土壤深度實測溫度相吻合：y=-0.788+0.818x(3 cm，R=0.964，p<0.01)，y=-1.214+0.705x(10 cm，R=0.942，p<0.01)，y=0.912+0.474x(20 cm，R=0.836，p<0.01)，y=1.004+0.361x(40 cm，R=0.714，p<0.01)；40 cm深度土壤預測溫度計算α-HCH遷移擴散通量系數(shù)產(chǎn)生的相對誤差小于使用日平均氣溫作為替代數(shù)據(jù)的計算結果。本研究結果將可能為有機污染物在土壤中的遷移擴散行為及相關數(shù)值模擬研究提供參考價值。

凍融期; 土壤溫度; 有機污染物; 遷移行為

土壤因其具有相對穩(wěn)定的性質(zhì)而成為有機污染物在陸地生態(tài)環(huán)境系統(tǒng)重要的貯存場所，同時由于土壤所擁有較大的貯存能力，以及有機污染物在土壤中較長的半衰期，從而使得土壤成為有機污染物在全球范圍內(nèi)進行重新分配的重要主導力量[1-3]。有機污染物在土壤中的環(huán)境監(jiān)測研究雖能較好地代表監(jiān)測范圍內(nèi)土壤的污染狀況，但是限于大尺度范圍內(nèi)土壤樣品采集的不現(xiàn)實性，運用多介質(zhì)環(huán)境數(shù)值模型或擴散模型，進行數(shù)值模擬研究則能夠很好地解決這一問題。Harner等[4]對農(nóng)藥土—氣交換及土壤中賦存狀態(tài)的研究表明，細化土壤垂直分層能夠更好地描述有機污染物的環(huán)境行為。另有研究表明，環(huán)境溫度是土壤水熱狀況及有機污染物在土壤中殘留濃度最重要的影響因素之一[5-7]。同時有關土壤溫度變化的研究多見于凍融期不同覆蓋條件下[8-10]，而對有機污染物在凍融期的遷移擴散行為的研究卻不多見。然而，現(xiàn)有的對土壤垂直分層較為全面且細化的多介質(zhì)環(huán)境模型或擴散模型[11-13]，由于缺乏土壤垂直分層不同深度土壤溫度的觀測數(shù)據(jù)，多以日平均氣溫作為不同深度土壤溫度的替代數(shù)據(jù)，這可能會給有機污染物土壤殘留狀況的數(shù)值模擬研究結果帶來較大的不確定性。

基于此，本研究將以我國北方高寒地區(qū)典型城市——哈爾濱市為研究對象，測定凍融期不同覆蓋條件下土壤垂直分層溫度，通過土—氣交換層和相鄰土壤層的遷移擴散通量系數(shù)，了解土壤垂直分層溫度變化對凍融期有機污染物遷移擴散行為的影響。

1　材料與方法

1.1數(shù)據(jù)來源

研究區(qū)為我國北方高寒地區(qū)典型城市——哈爾濱市，試驗于2013年11月—2014年4月在東北農(nóng)業(yè)大學節(jié)水灌溉試驗場進行為期6個月的現(xiàn)場觀測。試驗區(qū)位于北緯45°44′24″，東經(jīng)126°43′07″，屬溫帶大陸性季風氣候，冬長夏短，全年平均氣溫約2～4℃，冬季1月氣溫最低，平均氣溫為-19℃；夏季7月氣溫最高，平均氣溫為23℃。研究區(qū)域面積為1.0萬m2，主要種植旱田作物，地表類型分別設置為裸地，5，10，15 cm秸稈覆蓋。土壤垂直層深度劃分為3，10，20，40，60，100 cm，且溫度采用JL—04型6路地表溫度記錄儀進行自動采集。2013年11月—2014年10月日平均氣溫數(shù)據(jù)采用TRM-ZS1型氣象生態(tài)環(huán)境監(jiān)測系統(tǒng)進行自動記錄。試驗區(qū)土壤凍結期為11月9日—翌年3月8日，融化期為3月9日—4月28日。有機氯農(nóng)藥α-HCH在我國雖已停止使用多年，但其較為活躍的物理化學性質(zhì)使其易于在土壤及大氣環(huán)境介質(zhì)中遷移擴散，故本研究選取持久性有機污染物的典型代表α-HCH作為目標研究對象，其物理化學參數(shù)來自于文獻[14—16]。

1.2模型概述

有機污染物在土壤中的環(huán)境行為包括擴散、揮發(fā)、降解及滲濾等過程，而在土壤中的遷移(擴散)行為過程則由土壤—大氣界面的交換及相鄰土壤層之間的擴散過程構成。本文以土壤模型[1]中涉及的有機污染物在土壤層間及土壤—大氣層交換的物理過程為基礎，研究有機污染物在土壤中的遷移行為。通常有機污染物在相鄰土壤層之間的遷移擴散通量系數(shù)是由這兩層土壤的空氣和水的擴散過程組成，可由下式表示[4]：

(1)

式中：n，n+1——相鄰的兩個土壤層；Dv，Dv(n,n+1)——有機污染物的遷移擴散D值。而有機污染物在土壤交換層與大氣間的遷移行為過程模塊由遷移(擴散)通量系數(shù)Dv[mol/(Pa·h)]表示[4]：

(2)

Da=ABAEZa/YD

(3)

Dw=ABWEZw/YD

(4)

De=AKVZa

(5)

式中：De——空氣邊界層的揮發(fā)速率系數(shù)[mol/(Pa·h)]；Da，Dw——土壤中空氣和水的擴散速率系數(shù)[mol/(Pa·h)]；A——研究區(qū)域面積(m2)；KV——邊界層質(zhì)量遷移系數(shù)(m/h)，KV=BA/0.00475；BAE——土壤中空氣的有效擴散系數(shù)(m2/h)；BWE——土壤中水的有效擴散系數(shù)(m2/h)；YD——土壤的對數(shù)平均擴散距離(m)，可由(dl-du)/ln(dl/du)計算得出，dl和du——下層和上層土壤深度；Za，Zw——有機污染物在大氣和土壤中的逸度容量[mol/(m3·Pa)]。

進而利用Millington-Quirk方程計算出土壤中空氣和水的有效擴散系數(shù)BAE和BWE：

(6)

(7)

式中：fsa，fsw——土壤中空氣和水的體積分數(shù)。

空氣分子擴散系數(shù)BA(m2/h)和水相分子擴散系數(shù)BW(m2/h)可以用下式計算[18-20]：

(8)

(9)

式中：T——絕對溫度(K)；Ma——空氣的摩爾質(zhì)量(g/mol)；M——模擬化合物的摩爾質(zhì)量(g/mol)；P——大氣壓(1 atm)；Va——空氣的摩爾體積(cm3/mol)；Vm——模擬化合物的摩爾體積(cm3/mol)；μ——水的粘度系數(shù)。

1.3模擬試驗方案設計

有機污染物在大氣和土壤中的環(huán)境過程包括大氣向土壤的干濕沉降、雨水淋洗和擴散過程，土壤向大氣的揮發(fā)(擴散)過程，土壤層間擴散、滲濾和徑流過程，以及有機污染物在大氣和土壤中的降解行為，可以說有機污染物在土壤中的環(huán)境行為是比較復雜的，所有這些環(huán)境過程構成了有機污染物在大氣—土壤及不同土壤層間的遷移轉化行為。已有研究表明不同土壤類型(旱田、水田、林地、草地、荒地和城市土壤)的含水和含氣體積分數(shù)及土壤有機質(zhì)質(zhì)量分數(shù)均會對有機污染物的遷移擴散行為產(chǎn)生影響[13,21]。因研究期間處于溫度相對較低的冬季，而α-HCH作為有機氯農(nóng)藥類化學品，在土壤中主要通過微生物的降解動力學過程來進行轉化的行為相對較弱，因此在不考慮α-HCH大氣和土壤中化學過程，同時假定土壤有機質(zhì)、土壤無機質(zhì)、土壤水和氣的體積分數(shù)一定的情況下，單純從其物理過程即遷移(擴散)行為的角度，通過擴散D值來分析土壤垂直分層溫度的變化對旱田土壤α-HCH遷移擴散行為的影響。

2　結果與分析

2.1不同覆蓋條件下溫度對α-HCH遷移行為的影響分析

2.1.1不同覆蓋條件下土壤溫度對土壤垂直分層α-HCH遷移行為的影響由圖1可知，在不考慮地表覆蓋類型的情況下，α-HCH的遷移擴散D值隨著土壤深度的增加呈現(xiàn)遞減的趨勢，且3 cm和100 cm深度相差近兩個數(shù)量級，而這種差異則隨著土壤深度的增加逐漸減小，這可能受土壤對數(shù)平均擴散距離(YD)的影響，YD與土壤層厚度有關，隨著土壤層厚度的增加而增大，從而導致土壤中空氣和水的擴散速率系數(shù)減小，進而會使α-HCH遷移擴散D值變小。然而，將考慮的6種土壤垂直分層的YD值與α-HCH的遷移擴散D值作相關性分析發(fā)現(xiàn)，二者之間并沒有統(tǒng)計學意義上的相關性；而將20，40，60，100 cm深度的YD值與遷移擴散D值作相關性分析則表明，二者之間存在顯著負相關，R=-0.962，p<0.05，這說明土壤對數(shù)平均擴散距離對20—100 cm深度土壤α-HCH的遷移擴散行為具有影響，同時從遷移擴散D值的數(shù)值變化來看，這種差異性會隨著土壤深度的增加而減小。

秸稈覆蓋相對于裸地來說，對土壤溫度的變化具有雙向作用，即試驗初期當溫度下降時秸稈覆蓋能夠阻止土壤熱量的散失，而試驗末期當溫度上升時也會影響太陽輻射和氣溫對土壤溫度的傳遞。由圖1可以看出，裸地土壤垂直分層α-HCH遷移擴散D值隨溫度變化的波動最大且隨著土壤深度的增加波動振幅逐漸減弱；對于10 cm和15 cm秸稈覆蓋，不同深度土壤α-HCH的遷移擴散D值變化基本一致且振幅較??；裸地不同深度土壤α-HCH遷移擴散D值的變化則較大；同時與裸地相比，所有覆蓋類型α-HCH的遷移擴散D值的變化由淺層土壤(3，10，20 cm)至深層土壤(40，60，100 cm)隨相應土壤深度溫度的變化均逐漸趨于平穩(wěn)，且這種變化存在明顯的滯后效應，即當試驗期間日平均氣溫由下降到上升期間，在秸稈覆蓋條件下，α-HCH的遷移擴散D值隨著土壤深度的增加，受溫度變化影響而減小或增大的趨勢逐漸變緩，即由于秸稈覆蓋的存在，凍結期α-HCH遷移擴散D值高于裸地，而融化期(除100 cm深度外)則呈現(xiàn)相反狀態(tài)。為進一步探尋土壤垂直分層溫度與遷移擴散D值的關系，將裸地3，10，20，40，60，100 cm深度土壤溫度分別與相應層的遷移擴散D值作相關性分析得到，其相關系數(shù)分別為：R=0.851；R=0.882；R=0.943；R=0.977；R=0.987；R=0.956，且均在0.01水平上顯著相關。由此可知，不同覆蓋條件下，深層土壤α-HCH遷移擴散D值數(shù)值的變化與土壤層厚度有關，而遷移擴散D值曲線變化的波動振幅則與溫度的變化有關，同時土壤垂直分層溫度將會影響遷移擴散D值的變化，與裸地相比不同秸稈覆蓋厚度會使遷移擴散D值曲線的變化相對平穩(wěn)且在凍結期高于裸地的遷移擴散D值。

圖1　不同覆蓋條件下土壤垂直分層擴散D值隨時間的變化特征

2.1.2日平均溫度變化對土—氣交換層α-HCH環(huán)境行為的影響分析為了多介質(zhì)環(huán)境數(shù)值模型模擬有機污染物環(huán)境行為研究的需要，通常將土壤垂直分層劃分為交換層(深度0.1 cm)、緩沖層(深度1 cm)、耕作層(深度20 cm)和貯蓄層(深度30 cm)，而土壤交換層則主要是用于描述模擬有機污染物在界面交換時的過渡過程。

因秸稈覆蓋對土壤熱量的傳遞和散失有阻滯作用，且由于覆蓋層的存在也會阻礙污染物在大氣與土壤交換層之間的遷移擴散行為，故選取裸地來研究污染物的遷移擴散行為。由圖2和圖1比較可知，土壤交換層的遷移擴散D值高出3～100 cm土壤深度約1～3個數(shù)量級，這表明土壤交換層對有機污染物的土—氣交換及土壤垂直分層中有機污染物遷移擴散行為均具有重要的影響，即為了保持平衡狀態(tài)，當污染物的大氣濃度高于土壤濃度時，其環(huán)境行為會表現(xiàn)為通過土壤交換層向下層土壤進行垂直遷移，而當大氣濃度低于土壤濃度時深層土壤中的污染物會通過擴散行為對土壤交換層進行補給。通過2013年11月-2014年10月裸地土壤交換層遷移擴散D值與日平均氣溫的相關性分析得到，二者之間具有顯著相關性，R=0.999，p=0.000，隨著溫度升高會增大α-HCH的遷移擴散D值，每1℃溫度的變化可以引起α-HCH遷移擴散D值的變化量為0.55 mol/(Pa·h)。進一步探討日溫度變化對α-HCH遷移擴散D值的影響，由圖3可以看出，2013年11月1日0—24時，二者之間的變化具有明顯的一致性(R=1.000，p=0.000)，溫度在-1.4～10.9℃范圍可以引起α-HCH遷移擴散D值在188.07～194.93 mol/(Pa·h)之間變化，即每1℃溫度的變化可以引起α-HCH遷移擴散D值的變化量為0.56 mol/(Pa·h)。由此可知，日平均氣溫的變化將會強烈影響土—氣交換層α-HCH的環(huán)境行為。

圖2　日平均氣溫與土-氣交換層擴散D值的關系

2.2日平均氣溫變化對土壤垂直分層溫度的影響

土壤垂直分層由于溫度的差異性會對α-HCH的

遷移擴散D值的變化產(chǎn)生影響。與裸地相比，3種秸稈覆蓋條件下的α-HCH遷移擴散D值的波動均較為平穩(wěn)，裸地土壤因與大氣環(huán)境直接接觸而受環(huán)境溫度變化的影響較為劇烈，從而會引起α-HCH的遷移擴散D值的強烈波動。為了進一步探尋日平均氣溫對土壤垂直分層溫度變化的影響，將試驗期日平均氣溫與觀測獲取的不同土壤深度的溫度以及相鄰土壤層間的溫度數(shù)據(jù)進行相關性分析，試圖獲得二者之間是否存在必然的聯(lián)系。由表1可以看出，除100 cm深度的土壤溫度外，其他深度土壤與日平均氣溫之間均存在顯著正相關關系，且相關系數(shù)隨著土壤深度的增加逐漸降低。然而，在相鄰土壤層間土壤溫度則具有較高的相關程度，這表明相鄰土壤層間存在較強烈的溫度傳遞關系，而大氣環(huán)境溫度的傳遞作用則隨著土壤深度的增加明顯減弱。

圖3　2013年11月9日環(huán)境溫度變化與擴散D值

日平均氣溫土壤分層溫度3cm10cm20cm40cm60cm100cm日平均氣溫1.0003cm0.964**1.00010cm0.942**0.994**1.000土壤分20cm0.836**0.929**0.955**1.000層溫度40cm0.714**0.832**0.872**0.962**1.00060cm0.578**0.709**0.756**0.893**0.967**1.000100cm0.020 0.135*0.202**0.415**0.608**0.765**1.000

注：**表示p<0.01,*表示p<0.05。

在研究有機污染物在土壤中遷移擴散行為時，由于實際觀測數(shù)據(jù)的缺失，往往忽略土壤垂直分層溫度的差異，而以日平均氣溫作為不同深度土壤溫度的替代數(shù)據(jù)，這將可能會增加有機污染物土壤殘留濃度的模擬預測結果的不確定性。

為此，本研究試圖通過日平均氣溫與模型研究需要的土壤垂直分層3，10，20，40 cm深度溫度的相關關系建立預測模型，將日平均氣溫(變化范圍-23.11～20.85℃)作為自變量，分別與其作回歸分析得到預測方程：

y=-0.788+0.818x(3 cm，R=0.964，p<0.01)

(10)

y=-1.214+0.705x(10 cm，R=0.942，p<0.01)

(11)

y=0.912+0.474x(20 cm，R=0.836，p<0.01)

(12)

y=1.004+0.361x(40 cm，R=0.714，p<0.01)

(13)

根據(jù)預測方程計算不同深度土壤溫度的預測值，與實際觀測值進行對比分析(圖4)可以看出，由日平均氣溫得到的預測值與實際觀測值基本吻合，并能反映出不同深度土壤溫度隨環(huán)境溫度的變化趨勢，而40 cm土壤深度的預測值與實際觀測值略有偏差，但也處于實際觀測最大值與最小值區(qū)間范圍之內(nèi)，這表明本研究根據(jù)日平均氣溫建立的回歸方程能夠較好地反映不同深度土壤溫度的變化情況，可以為考慮北方高寒地區(qū)土壤垂直分層溫度變化條件下的多介質(zhì)環(huán)境數(shù)值模擬模型的構建及優(yōu)化提供一定參考價值。

圖4　土壤實測溫度與預測溫度對比

2.3不同溫度數(shù)據(jù)對40 cm土壤深度α-HCH遷移行為的差異性分析

由圖4可知，40 cm土壤深度溫度預測值與實測值的偏差略大，為進一步驗證這種偏差對α-HCH遷移擴散行為帶來的影響，本研究將分別用40 cm深度土壤溫度的實測值、預測值與日平均氣溫來計算該深度α-HCH的遷移擴散D值，結果如圖5所示。應用實測與預測溫度值計算40 cm土壤深度α-HCH的遷移擴散D值基本一致，而應用日平均氣溫作為替代數(shù)據(jù)的計算值則存在強烈的波動，這表明根據(jù)日平均氣溫預測的40 cm土壤深度的溫度數(shù)值能夠較好地用于α-HCH遷移擴散D值的計算。為了能夠更為直觀地體現(xiàn)出日平均氣溫、40 cm深度預測土壤溫度、40 cm深度實測土壤溫度對α-HCH的遷移擴散行為的差異性，分別計算預測溫度與日平均氣溫對真實結果的相對誤差。

由圖6可以看出，應用預測溫度計算的α-HCH遷移擴散D值的誤差要遠小于應用日平均氣溫計算的誤差。通常模型模擬的運行周期一般長達幾年甚至幾十年，這種誤差也將會逐步放大，進而會對污染物土壤殘留狀況數(shù)值模擬結果的準確性產(chǎn)生較大的影響。由此說明，根據(jù)日平均氣溫預測不同深度土壤溫度可以進一步應用到北方高寒區(qū)有機污染物遷移行為的數(shù)值模擬研究中，且具有一定的可靠性。

圖6　日平均氣溫與預測溫度對實際擴散D值的相對誤差

3　結 論

(1) 不同覆蓋條件下，土壤垂直分層深層土壤α-HCH遷移擴散D值數(shù)值的變化與土壤層厚度有關，土壤垂直分層溫度的變化與α-HCH的遷移擴散通量系數(shù)呈顯著正相關，且隨著土壤深度的增加相關性逐漸減弱。在凍結期秸稈覆蓋會使α-HCH遷移擴散D值大于裸地，且變化相對平穩(wěn)，而融化期(除100 cm土壤深度外)則呈相反趨勢。

(2) 日平均氣溫的變化將會強烈影響土—氣交換層α-HCH的環(huán)境行為，溫度每1℃的變化可以引起α-HCH的遷移擴散D值的變化量約為0.56 mol/(Pa·h)。

(3) 日平均氣溫會影響土壤垂直分層溫度的變化，以日平均氣溫作為自變量分別與模擬研究需要的土壤垂直分層溫度觀測數(shù)據(jù)作回歸分析得到預測方程，能夠較好地與不同土壤層的實測溫度吻合。本研究雖具有一定的局限性，但通過研究可以證實土壤垂直分層溫度的變化將會影響有機污染物的遷移擴散行為。

(4) 本研究根據(jù)日平均氣溫建立的回歸方程可以為北方高寒地區(qū)有機污染物在土壤中的遷移擴散行為研究提供一定的參考價值。通過對旱田土壤垂直分層溫度變化對有機污染物遷移擴散行為的研究也可為進一步對不同土壤類型(水田、草地、林地、城市土壤及荒地等)垂直分層土壤溫度的觀測試驗，以及相應土壤類型有機污染物的遷移擴散行為的研究提供借鑒。

[1]Mackay D. Multimedia environmental models: the fugacity approach[M]. Boca Raton: CRC Press,2001.

[2]Wang X, Sheng J, Gong P, et al. Persistent organic pollutants in the Tibetan surface soil: spatial distribution, air-soil exchange and implications for global cycling[J]. Environmental Pollution,2012,170(8):145-151.

[3]Liu X, Ming L L, Nizzetto L, et al. Critical evaluation of a new passive exchange-meter for assessing multimedia fate of persistent organic pollutants at the air-soil interface[J]. Environmental Pollution, 2013,181(6):144-150.

[4]Harner T, Bidleman T F, Jantunen L M M, et al. Soil-air exchange model of persistent pesticides in the United States cotton belt[J]. Environmental Toxicology and Chemistry,2001,20(7):1612-1621.

[5]Tripathi V, Abhilash P C, Singh H B, et al. Effect of temperature variation on lindane dissipation and microbial activity in soil[J]. Ecological Engineering,2015,79:54-59.

[6]何匯虹,劉文兆.黃土塬區(qū)農(nóng)田土壤溫度變化特征研究[J].水土保持研究,2013,20(5):124-128.

[7]Zhang Q Q, Ying G G, Chen Z F, et al. Multimedia fate modeling and risk assessment of a commonly used azole fungicide climbazole at the river basin scale in China[J]. Science of the Total Environment,2015,520:39-48.

[8]陳軍鋒,鄭秀清,秦作棟,等.凍融期秸稈覆蓋量對土壤剖面水熱時空變化的影響[J].農(nóng)業(yè)工程學報,2013,29(20):102-110.

[9]胡銘,劉志輝,陳凱,等.雪蓋影響下季節(jié)性凍土消融期的土壤溫度特征分析[J].水土保持研究,2013,20(3):39-43.

[10]邢述彥,劉虎,鄭秀清,等.秸稈覆蓋厚度對凍融期土壤溫度的影響[J].太原理工大學學報,2013,43(6):741-744.

[11]Ma J, Daggupaty S, Harner T, et al. Impacts of lindane usage in the Canadian prairies on the Great Lakes ecosystem. 1. Coupled atmospheric transport model and modeled concentrations in air and soil[J]. Environmental Science & Technology,2003,37(17):3774-3781.

[12]Tian C, Liu L, Ma J, et al. Modeling redistribution ofα-HCH in Chinese soil induced by environment factors[J]. Environmental Pollution, 2011,159(10):2961-2967.

[13]崔嵩,楊萌,李一凡.不同土壤類型多氯聯(lián)苯土壤殘留特征變化分析[J].吉林大學學報:地球科學版,2012,42(6):1880-1888.

[14]Zhang N, Yang Y, Liu Y, et al. Determination of octanol-air partition coefficients and supercooled liquid vapor pressures of organochlorine pesticides[J]. Journal of Environmental Science and Health Part B,2009,44(7):649-656.

[15]Mackay D, Shiu W, Ma K. Physical-chemical properties and environmental fate handbook[CD]. Boca Raton: CRC Press,1999.

[16]Shoeib M, Harner T. Using measured octanol-air partition coefficients to explain environmental partitioning of organochlorine pesticides[J]. Environmental Toxicology and Chemistry,2002,21(5):984-990.

[17]Jury W A, Spencer W F, Farmer W J. Behavior assessment model for trace organics in soil: I. Model description[J]. Journal of Environmental Quality,1983,12(4):558-564.

[18]Schwarzenbach R P, Gschwend P M, Imboden D M. Environmental Organic Chemistry[M]. New York: John Wiley & Sons Inc.,2003.

[19]Hayduk W, Laudie H. Prediction of diffusion coefficients for nonelectrolytes in dilute aqueous solutions[J]. AIChE Journal, 1974,20(3):611-615.

[20]Severinsen M, Jager T. Modelling the influence of terrestrial vegetation on the environmental fate of xenobiotics[J]. Chemosphere,1998,37(1):41-62.

[21]崔嵩,劉麗艷,李一凡.土壤有機質(zhì)影響多氯聯(lián)苯土壤殘留變化特征[J].中南大學學報:自然科學版,2013,44(11):4759-4767.

Effect of Soil Temperature on Migration Behavior of Organic Pollutants During Freeze-Thaw Period

CUI Song1,2, FU Qiang1,2, LI Tianxiao1,2, LI Yifan3

(1.SchoolofWaterConservancy&CivilEngineering,NortheastAgriculturalUniversity,Harbin150030,China;2.InternationalJointResearchCenterforPersistentToxicSubstances,NortheastAgriculturalUniversity,Harbin150030,China;3.InternationalJointResearchCenterforPersistentToxicSubstances,HarbinInstituteofTechnology,Harbin150090,China)

In order to understand soil effects of temperatures in different vertical layers on migration and diffusion behavior of organic pollutants during freeze-thaw period, the coefficient of migration and diffusion flux was used to assess the migration behavior of the soil-air exchange layer and adjacent soil layers. We deeply discussed the soil temperature effect on migration of organic pollutants and systematically analyzed differences on simulation results applying different temperature data. The results show that the migration and diffusionDvalue ofα-HCH in deep soil was related with the thickness of soil layer; the soil temperatures of different vertical layers affected the migration and diffusion behavior ofα-HCH; the change of diurnal average air temperature strongly affected the environmental behavior ofα-HCH at soil-air exchange layer (R=0.999,p=0.000); the established regression equation between diurnal average air temperature and soil vertical layer temperature could well coincide with measured temperatures of different soil depths:y=-0.788+0.818x(3 cm,R=0.964,p<0.01)；y=-1.214+0.705x(10 cm,R=0.942,p<0.01);y=0.912+0.474x(20 cm,R=0.836,p<0.01)；y=1.004+0.361x(40 cm,R=0.714,p<0.01)；The relative error of migration and diffusion D value used to predict temperature was less than using diurnal average air temperature as surrogate data. The established regression equation applying diurnal average air temperature in this study has great significance for examining the soil temperature effects of different vertical layers on migration and diffusion behavior of organic pollutants, and can further improve and enhance the accuracy of simulation results applying multimedia environmental modeling and diffusion modeling.

freeze-thaw period; soil temperature; organic pollutants; migration behavior

2015-05-17

2015-6-16

國家自然科學基金(51279031,41401550);黑龍江省自然科學基金(E201241);黑龍江省博士后資助項目(LBH-Z13029);教育部新世紀優(yōu)秀人才支持計劃(NCET-11-0952);黑龍江省杰出青年基金(JC201402)；東北農(nóng)業(yè)大學“青年才俊”基金(14QC49)

崔嵩(1981—),男,黑龍江寶清人,博士,副教授,主要從事農(nóng)業(yè)水土資源環(huán)境效應及持久性有毒物質(zhì)數(shù)值模擬研究。E-mail:cuisong-bq@163.com

付強(1973—),男,黑龍江哈爾濱人,博士,教授,博士生導師,主要從事農(nóng)業(yè)水土資源系統(tǒng)分析和節(jié)水灌溉研究。E-mail:fuqiang0629@126.com

X592; S152.8

A

1005-3409(2016)03-0314-07