張 巖，李暢游，高 寧，史小紅，喬玲敏

（1：煙臺(tái)大學(xué)土木工程學(xué)院，煙臺(tái)264005）（2：內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院，呼和浩特010018）

結(jié)冰對(duì)烏梁素海水體富營養(yǎng)化的影響?

張 巖1，李暢游2，高 寧1，史小紅2，喬玲敏1

（1：煙臺(tái)大學(xué)土木工程學(xué)院，煙臺(tái)264005）
（2：內(nèi)蒙古農(nóng)業(yè)大學(xué)水利與土木建筑工程學(xué)院，呼和浩特010018）

為研究結(jié)冰對(duì)湖泊富營養(yǎng)化狀態(tài)的影響，以內(nèi)蒙古烏梁素海為研究對(duì)象，用“注射器”式采水器和冰鉆分別采集水樣和冰樣，并測定總氮、總磷、CODCr、葉綠素a濃度、冰厚、水深和透明度等指標(biāo)，運(yùn)用分形理論評(píng)價(jià)模型對(duì)湖泊結(jié)冰前、后的富營養(yǎng)化等級(jí)進(jìn)行綜合評(píng)價(jià).結(jié)果表明：烏梁素海湖泊結(jié)冰厚度與水深存在較好的負(fù)相關(guān)性，其相關(guān)系數(shù)為0.54；湖泊結(jié)冰后，水體的富營養(yǎng)化等級(jí)由富營養(yǎng)變?yōu)橹馗粻I養(yǎng)，表明結(jié)冰過程中，氮、磷等營養(yǎng)元素被排斥至冰下水體，即結(jié)冰過程對(duì)湖泊水體中的營養(yǎng)元素具有濃縮效應(yīng)，導(dǎo)致冰下環(huán)境質(zhì)量惡化，并促進(jìn)營養(yǎng)元素向沉積物轉(zhuǎn)移；但也可以利用該濃縮效應(yīng)，從湖泊冰封期富營養(yǎng)化污染的特殊性入手，運(yùn)用底泥疏浚和冰體-水體分離等工程技術(shù)手段，實(shí)現(xiàn)湖泊富營養(yǎng)化污染的治理.

結(jié)冰；富營養(yǎng)化；分形理論；烏梁素海

湖泊是重要的淡水資源之一，我國湖泊數(shù)量眾多、分布廣泛，但由于水體中氮、磷等濃度增加，導(dǎo)致湖泊富營養(yǎng)化不斷加重.據(jù)統(tǒng)計(jì)，我國有85.4%的湖泊超過了富營養(yǎng)化標(biāo)準(zhǔn)，達(dá)到重富營養(yǎng)化程度的湖泊達(dá)40.1%［1］.為此，眾多學(xué)者在湖泊富營養(yǎng)化方面開展了大量的研究并取得了大量的成果，但這些研究大都基于湖泊的非冰封條件下進(jìn)行的，而對(duì)湖泊冰封條件下湖泊的富營養(yǎng)化特征卻少見報(bào)道.國內(nèi)僅有的幾篇文獻(xiàn)主要研究了氮、磷等在冰-水體系的分配關(guān)系和其在冰生長期的遷移過程［2-5］，而未曾耦合結(jié)冰過程對(duì)總氮（TN）、總磷（TP）、CODCr、葉綠素a（Chl.a）和透明度（SD）等的影響，闡明結(jié)冰過程對(duì)湖泊富營養(yǎng)化狀態(tài)的影響.而國外尚未見湖泊冰封期富營養(yǎng)化方面的報(bào)道.但湖泊冰封期，由于冰層的存在幾乎切斷了大氣-水體物質(zhì)交換和削弱了其間的能量交換，改變了整個(gè)體系的透光性和熱傳導(dǎo)性［6］，水中存在的微生物數(shù)量很少，且微生物的生物降解能力極低；湖水中污染物失去遷移和光解作用，并使自然曝氣形成的復(fù)氧過程幾乎停止，并且光合作用微弱，溶解氧濃度處于極低值，氧化還原電位降低，水體的自凈能力大大減弱.因此，結(jié)冰必將對(duì)湖泊的富營養(yǎng)化狀態(tài)產(chǎn)生影響.

1 材料與方法

1.1 研究區(qū)概況

烏梁素海（40°36′～41°03′N，108°43′～108°57′E）位于內(nèi)蒙古自治區(qū)巴彥淖爾市境內(nèi)，水面面積為283 km2，水深1.4～3.1 m，是國家八大淡水湖之一.該湖泊主要承納河套灌區(qū)的農(nóng)田退水及流域內(nèi)部分工業(yè)廢水和生活污水，湖泊富營養(yǎng)化嚴(yán)重.烏梁素海流域?qū)俦睖貛Т箨懶愿珊担敫珊禋夂?，多年平均氣溫?.6～7.4℃，一般從每年11月開始結(jié)冰，翌年3月末開始融化，冰封期約4～5個(gè)月，冰層的最大厚度可達(dá)0.70 m. 2012實(shí)驗(yàn)?zāi)甓葍?nèi)烏梁素海湖泊的結(jié)冰期為11月25日至次年2月18日.冬季冰蓋形成后，湖泊接近靜水體，環(huán)流作用微弱，與外界水體的交換近乎停止.

1.2 樣品采集

根據(jù)我國水環(huán)境及湖泊調(diào)查規(guī)范，考慮烏梁素海的水文及水環(huán)境特征，設(shè)置采樣點(diǎn)（圖1），并于2012年1月在每個(gè)采樣點(diǎn)上分別采集冰樣和水樣：用冰鉆采集6根冰芯樣品，用電子天平稱量每根冰芯樣品的重量并將其放在切割板上，用冰鋸從上至下每10 cm切割一段冰柱，并將其放在塑料瓶中，在室溫下自然融化后進(jìn)行監(jiān)測；用“注射器”式采水器通過冰孔各抽取冰-水界面、冰下30 cm和水-沉積物界面的水各1000 ml，將上述3個(gè)位置的水樣均勻混合后進(jìn)行監(jiān)測.

圖1 烏梁素海位置及采樣點(diǎn)Fig.1 Location and sampling sites of Lake Ulansuhai

1.3 監(jiān)測方法

各指標(biāo)監(jiān)測方法如表1所示.

1.4 富營養(yǎng)化評(píng)價(jià)方法

選用評(píng)價(jià)結(jié)果分辨率較高的分形理論評(píng)價(jià)模型對(duì)烏梁素海結(jié)冰前后的富營養(yǎng)化等級(jí)進(jìn)行評(píng)價(jià)［7］.

1.4.1 單指標(biāo)分形維數(shù)計(jì)算 （1）選用關(guān)聯(lián)維數(shù)計(jì)算分形維數(shù)，設(shè)研究水質(zhì)指標(biāo)Xi的數(shù)據(jù)向量為：Xi＝（xi1、xi2、…、xin），n為樣本數(shù)據(jù)的總數(shù).

（2）建立1～7維相空間，運(yùn)用式（1）和式（2）分別計(jì)算每維相空間兩點(diǎn)間距離rp，q（s）與平均距離Δxs：

表1 各指標(biāo)監(jiān)測方法Tab.1 Themethods of each item

式中，p，q＝1、2、…、n-s+1；s＝1、2、…、w（w為最大相空間維數(shù)）.

（3）運(yùn)用式（3）和式（4）計(jì)算每維相空間兩點(diǎn)間距離小于rsk的概率Ck（S），從而得到關(guān)聯(lián)積分函數(shù)：

式中，rsk為指定的距離上限；H為Heaviside函數(shù)；p，q＝1、2、…、n-s+1；s＝1、2、…、w；k＝1、2、…、14.

（4）根據(jù)每維相空間求出的一組（k＝14）Ck（s）值，若Ck（s）和rsk雙對(duì)數(shù)圖上為直線，則分形存在，其斜率為分維數(shù)：

若分維數(shù)在相空間維數(shù)升高時(shí)趨向極限，此極限即為空間的分維數(shù)；若分維數(shù)未嚴(yán)格地趨向極限，則需要比較不同相空間的分維數(shù)，選擇相鄰空間分維數(shù)之差滿足一定精度或趨于穩(wěn)定的最大者，將其作為該水質(zhì)指標(biāo)的分維數(shù).

1.4.2 水質(zhì)評(píng)價(jià)指標(biāo)的分維數(shù) 運(yùn)用上述方法可求出各指標(biāo)的分維數(shù)：（Chl.a，TP，TN，CODCr，SD）＝（0.33，1.17，1.00，1.59，1.83），分維數(shù)越大的水質(zhì)指標(biāo)越重要.

1.4.3 評(píng)價(jià)標(biāo)準(zhǔn)的建立 用式（6）建立各級(jí)別的評(píng)價(jià)標(biāo)準(zhǔn)：

式中，F(xiàn)為評(píng)價(jià)標(biāo)準(zhǔn)；i代表第i種指標(biāo)，i＝1、2、…、n；D（i）和S（i）分別為第i種指標(biāo)的分形維數(shù)和實(shí)測值.

表2 水質(zhì)評(píng)價(jià)標(biāo)準(zhǔn)［7］Tab.2 The estimation standard of water quality

所建立的評(píng)價(jià)標(biāo)準(zhǔn)為：（貧營養(yǎng)，貧中營養(yǎng)，中營養(yǎng)，中富營養(yǎng)，富營養(yǎng)和重富營養(yǎng)）＝（52.06，66.30，336.4，569.86，2271.48，6793.32）.

2 結(jié)果與分析

2.1 湖泊水深和冰厚的分布

為了解采樣時(shí)湖泊各采樣點(diǎn)的冰厚和水深等水文特性，對(duì)湖面11個(gè)采樣點(diǎn)的水深、冰厚進(jìn)行監(jiān)測分析（圖2），結(jié)果表明烏梁素海采樣期各采樣點(diǎn)的水深介于0.46～2.88 m之間，其平均值為1.17 m，相對(duì)標(biāo)準(zhǔn)偏差為0.54；冰厚介于0.32～0.55m之間，其平均值為0.43m，相對(duì)標(biāo)準(zhǔn)偏差為0.17.湖泊的南部和北部水深較深，中部水深較淺，而湖泊中部的冰厚較南部和北部大，這表明冰厚與水深存在較好的負(fù)相關(guān)性，其相關(guān)系數(shù)為0.54（圖3），這可從熱力學(xué)對(duì)湖冰形成過程的影響進(jìn)行解釋：整個(gè)湖泊處于相同的環(huán)境條件下，任意湖泊區(qū)域大氣-冰界面的冷通量相等，水深越淺，受冷水的體積越小，冰下水體更易達(dá)到過冷狀態(tài)，進(jìn)而使得冰厚較快增加，即與水深呈負(fù)相關(guān)關(guān)系.

圖2 各采樣點(diǎn)冰厚和水深變化Fig.2 The variety of ice thickness and water depth of each sampling site

圖3 冰厚與水深關(guān)系圖Fig.3 The relationship between ice thickness and water depth

2.2 湖泊富營養(yǎng)化等級(jí)評(píng)價(jià)

運(yùn)用分形理論評(píng)價(jià)模型對(duì)烏梁素海結(jié)冰前水體、結(jié)冰后冰體融水和冰下水的富營養(yǎng)化等級(jí)進(jìn)行綜合評(píng)價(jià)（表3）.

結(jié)冰前，烏梁素海湖泊11個(gè)采樣點(diǎn)均達(dá)到了富營養(yǎng)化程度，已經(jīng)沒有貧營養(yǎng)、貧中營養(yǎng)和中營養(yǎng)的區(qū)域出現(xiàn)（圖4）.這是由于近年來，河套灌區(qū)實(shí)施嚴(yán)格的節(jié)水灌溉制度，入湖水量減少，水位下降，而化肥施用量卻逐年增加.據(jù)統(tǒng)計(jì)，2009年河套灌區(qū)的化肥施用量已達(dá)72×104t，是1978年的10.2倍，其中大約70%的化肥經(jīng)總排干流入湖內(nèi)，并且沿途還匯入了大量的工業(yè)廢水和生活污水［8］，致使烏梁素海湖泊水質(zhì)日益惡化，全湖均達(dá)到了富營養(yǎng)化程度.

結(jié)冰后，冰下水樣中，5個(gè)采樣點(diǎn)達(dá)到了重富營養(yǎng)化程度，占樣品數(shù)量的45.4%，6個(gè)采樣點(diǎn)仍然富營養(yǎng)化，占樣品數(shù)量的54.6%，全湖均值的評(píng)價(jià)結(jié)果為3738.29，評(píng)價(jià)等級(jí)為重富營養(yǎng)；冰層融水樣品中，7個(gè)采樣點(diǎn)為中富營養(yǎng)，占樣品數(shù)量的63.6%，4個(gè)采樣點(diǎn)為富營養(yǎng)化，占樣品數(shù)量的36.4%，全湖均值的評(píng)價(jià)等級(jí)為富營養(yǎng)，但其評(píng)價(jià)結(jié)果為877.14，僅為冰前水評(píng)價(jià)結(jié)果的54.6%（表3）.

2.3 結(jié)冰對(duì)富營養(yǎng)化狀態(tài)的影響

從全湖分形評(píng)價(jià)均值來看，結(jié)冰前后，湖泊水體的富營養(yǎng)化等級(jí)由富營養(yǎng)變?yōu)橹馗粻I養(yǎng)，這表明結(jié)冰過程加劇了水體的富營養(yǎng)化程度，結(jié)冰過程中，氮、磷等營養(yǎng)元素被排斥至冰下水體，即結(jié)冰過程對(duì)營養(yǎng)元素具有濃縮效應(yīng).

湖泊結(jié)冰過程的濃縮效應(yīng)可從湖冰的形成過程進(jìn)行解釋：湖冰的形成是一個(gè)從敞開水域到部分水域，最后到全水域冰覆蓋的過程.當(dāng)氣溫低于0℃時(shí)，湖泊表層達(dá)到過冷狀態(tài)，形成自發(fā)長大的冰核.隨著負(fù)溫度的累積，湖面上形成細(xì)小、平整且不含雜質(zhì)的片狀冰晶，這些冰晶多平臥水面，在生長過程逐漸生長為更加利于釋放潛熱的星形或樹枝形，這些冰晶相互連接便可形成柔性冰，柔性冰經(jīng)凍結(jié)變厚，形成堅(jiān)硬冰蓋，并形成冰-水界面.由于冰蓋的形成，使得冷能通過冰蓋穩(wěn)定的傳遞，而由于氮、磷等污染物的凝固點(diǎn)較低，冷能僅能夠使得在冰-水界面處的水分子達(dá)到凝固點(diǎn)而平穩(wěn)析出，且晶核生長的速度大于其形成的速度［9］，隨著晶核的長大，其將附著在冰蓋的下表面上形成潔凈冰體，而氮、磷等營養(yǎng)元素則被排斥至冰下水體.但當(dāng)溫度較低時(shí)，冰晶體的形成速率加快，在冰的快速生長過程中會(huì)有部分雜質(zhì)不能被及時(shí)排斥至冰下水體，而被俘獲在純凈的冰體之間（圖5）.氮、磷等雜質(zhì)被俘獲的過程與海冰形成過程中鹽的俘獲相似，且該俘獲過程與冰生長速率、冰體結(jié)構(gòu)、冰下水環(huán)流和水體中雜質(zhì)的含量等因素相關(guān)［10］.特別是氣溫驟降時(shí)，冰的生長速率驟增，將會(huì)有更多的雜質(zhì)被俘獲在冰體中［11-12］，但仍有部分雜質(zhì)被排斥至冰下水體，使得冰下水體中氮、磷等污染物濃度仍高于冰體.

表3 實(shí)測數(shù)據(jù)和評(píng)價(jià)結(jié)果Tab.3 Original data and estimating results

圖4 結(jié)冰前水體、冰融水和冰下水體的富營養(yǎng)化等級(jí)（1-中富營養(yǎng)；2-富營養(yǎng)；3-重富營養(yǎng)）Fig.4 The eutrophic state ofwater before icing，icemeltwater and under-ice water

從熱力學(xué)角度來看，由于水和氮、磷等污染物物理性質(zhì)的不同，隨著氣溫的降低，整個(gè)體系的溫度也將下降，使水分子對(duì)氮、磷等污染物分子的溶劑化作用減弱，水分子在氫鍵作用下首先締結(jié)析出，形成純凈冰體，而氮、磷等污染物分子被排斥至冰下水體中.但在湖冰形成初期和氣溫驟降時(shí)，由于氮、磷等污染物的溶解度降低，會(huì)有部分污染物會(huì)從飽和或過飽和的溶液中析出，并可能形成臨界尺寸的冰核，促使水溶液異相成核［13-14］，即部分污染物被凍結(jié)在冰體中，但其在冰體中的濃度仍然小于冰下水體.

從固液平衡理論分析（圖6），曲線ABDE是湖水的冰點(diǎn)線.A點(diǎn)是純水的冰點(diǎn)，但在一定范圍內(nèi)湖水的冰點(diǎn)會(huì)隨著污染物濃度的增加而有所降低.隨著氣溫的逐漸降低，湖水溫度將由T0降至T1，湖水達(dá)到過冷狀態(tài)即開始結(jié)冰，氣溫繼續(xù)降低至T2時(shí)，潔凈冰晶形成、生長并從水體中分離出來，此時(shí)湖水污染物濃度由C0增加至C1.E點(diǎn)為共晶點(diǎn)，T3為共晶點(diǎn)溫度，氣溫繼續(xù)降低時(shí)，水分子以潔凈冰晶的形式從冰下水體中不斷析出，而氮、磷等污染物分子被排斥至冰下水體中，隨著冰層厚度的增加，水體中氮、磷等污染物會(huì)達(dá)到飽和，甚至從水體中析出［15］.

圖5 冰結(jié)構(gòu)剖面示意圖Fig.5 Schematic diagram of ice structure

圖6 溫度-濃度平衡曲線Fig.6 Temperature-concentration balance curve

盡管湖泊冰封期冰下水體氮、磷等營養(yǎng)元素已超過富營養(yǎng)化發(fā)生的臨界濃度，但冰下水體的溫度一直穩(wěn)定在0℃附近，抑制了藻類增殖，同時(shí)由于冰層的覆蓋減弱了冰下水體的光照強(qiáng)度，削弱了水體流速，擾動(dòng)等水動(dòng)力條件［16］，一般不會(huì)出現(xiàn)藻類異常增殖的情況［17］.但根據(jù)Langmuir吸附理論可知，冰下水體氮、磷等濃度的升高，會(huì)打破水-沉積物之間存在著的氮、磷等的濃度平衡，使得氮、磷等污染物在物理、化學(xué)和生物的耦合作用下由水體向沉積物遷移［5］，成為新的內(nèi)源污染，也有研究指出，與非冰封期相比，低溫和冰下水的低擾動(dòng)都將加快污染物向沉積物的遷移［18-19］，而該過程使得冰-水體系中氮、磷等的總量減少，這在一定程度上可以緩解因冷凍濃縮而引起的冰下水體的重富營養(yǎng)狀態(tài).但是，冰體融化后，水體中氮、磷等污染物的濃度迅速下降，打破沉積物-水兩相界面的物質(zhì)平衡，在風(fēng)力以及其他外力的擾動(dòng)下，大量聚集在沉積物中的污染物會(huì)重新懸浮而再次進(jìn)入湖泊水體中［20］.另外，對(duì)部分作為飲用水水源的湖泊而言，結(jié)冰過程使得冰下水體中氮、磷等污染物濃度升高，加之給水廠在冬季的處理效率降低，使得生活飲用水的水質(zhì)下降，甚至難以達(dá)到生活飲用水衛(wèi)生標(biāo)準(zhǔn)（GB 5749-2006）的要求.

2.4 淺水湖泊冰封期富營養(yǎng)化污染治理建議

由前述可知，結(jié)冰過程對(duì)湖泊水體中氮、磷等污染物具有濃縮效應(yīng)，使冰下水體中污染物濃度增加，進(jìn)而使得冰下水體中的污染物向沉積物遷移，并且冰厚與水深存在較好的負(fù)相關(guān)性，因此，可從湖泊冰封期富營養(yǎng)化特征的特殊性入手，提出湖泊特別是淺水湖泊冰封期富營養(yǎng)化的治理建議.

底泥疏浚是治理湖泊內(nèi)源污染的有效措施之一，但由于疏浚過程屬于水上作業(yè)，耗資巨大［21］.而在淺水湖泊冰封期，可以進(jìn)行冰上作業(yè)，從而降低工程成本，可使底泥疏浚工程得以大規(guī)模推廣.此外，湖泊冰封期水體中的污染物向沉積物遷移，致使沉積物中污染物濃度提高，因此可以提高底泥疏浚效果.

由于湖泊冰封期冰體較為潔凈，可采用冰體和水體分離的方法，移出并凈化淺水湖泊冰下較小體積的濃縮液，必要時(shí)可將冰下水體多次凍結(jié)，獲得更小體積和更高濃度的冰下濃縮液［22］，由此可以移除大量污染物.待氣溫回升后，潔凈冰體將自然融化形成低濃度的冰融水.但是該方法會(huì)威脅水生生物群落，破壞水生態(tài)系統(tǒng)的物質(zhì)循環(huán)和能量流動(dòng)過程，降低水體的自凈能力，如能取得突破，可為湖泊富營養(yǎng)化污染提供一種經(jīng)濟(jì)與效益并優(yōu)的治理方法.

3 結(jié)論

1）湖泊的結(jié)冰厚度與湖水深度具有良好的負(fù)相關(guān)性，即湖水越深，冰層厚度越小.

2）結(jié)冰過程對(duì)氮、磷等營養(yǎng)元素具有濃縮效應(yīng)，烏梁素海結(jié)冰后，水體的富營養(yǎng)化等級(jí)由富營養(yǎng)變?yōu)橹馗粻I養(yǎng)，即結(jié)冰過程加劇了水體的富營養(yǎng)化程度，并會(huì)對(duì)冰下水環(huán)境質(zhì)量帶來不良影響.

3）利用結(jié)冰過程的濃縮效應(yīng)，理論上可通過底泥疏浚和冰水分離的方式對(duì)富營養(yǎng)湖泊進(jìn)行治理，但在實(shí)際工程中可能存在破壞水體生態(tài)系統(tǒng)等缺點(diǎn)，需為其開展深入研究.

致謝：感謝內(nèi)蒙古農(nóng)業(yè)大學(xué)“河、湖、濕地水環(huán)境保護(hù)與修復(fù)技術(shù)研究團(tuán)隊(duì)”的同學(xué)在野外采樣工作中給予的無私幫助.

［1］ Yang Guishan，Ma Ronghua，Zhang Lu et al.Lake status，major problems and protection strategy in China.J Lake Sci，2010，22（6）：799-810.DOI：10.18307／2010.0601.［楊桂山，馬榮華，張路等.中國湖泊現(xiàn)狀及面臨的重大問題與保護(hù)策略.湖泊科學(xué)，2010，22（6）：799-810.］

［2］ Huang Jiguo，Peng Xiangjie，Yu Shuang et al.Distribution characteristics of nutrient and chlorophyll a during the icebound season.Journal of Jilin University：Science Edition，2008，46（6）：1231-1236.［黃繼國，彭祥捷，俞雙等.水體結(jié)冰期營養(yǎng)鹽和葉綠素a的分布特征.吉林大學(xué)學(xué)報(bào)：理學(xué)版，2008，46（6）：1231-1236.］

［3］ Huang Jiguo，F(xiàn)u Xinting，Wang Xuesong et al.The distribution characteristics of nutrients and phytoplankton during the icebound season in Changchun.Acta Scientiae Circumstantiae，2009，20（8）：1678-1683.［黃繼國，傅鑫廷，王雪松等.湖水冰封期營養(yǎng)鹽及浮游植物的分布特征.環(huán)境科學(xué)學(xué)報(bào)，2009，20（8）：1678-1683.］

［4］ Jiang Huiqin.Experimental study on the spatial distribution of nutrient salts in the ice of lake wuliangsuhai and their release processes duringmelting［Dissertation］.Huhhot：Inner Mongolia Agricultural University，2011.［姜慧琴.烏梁素海營養(yǎng)鹽在冰體中的空間分布及其在凍融過程中釋放規(guī)律的試驗(yàn)研究［學(xué)位論文］.呼和浩特：內(nèi)蒙古農(nóng)業(yè)大學(xué)，2011.］

［5］ Zhang Yan，Li Changyou，Shen Hungtao et al.Total nitrogen migration in Wuliangsuhai Lake during ice growth process. Advances in Water Science，2013，（5）：728-735.［張巖，李暢游，Shen Hungtao等.烏梁素海湖泊冰生長過程中總氮的遷移規(guī)律.水科學(xué)進(jìn)展，2013，（5）：728-735.］

［6］ Huang Wenfeng.Study on lake ice microstructure and its effects on thermal and mechanical parameters［Dissertation］. Dalian：Dalian University of Technology，2013.［黃文峰.湖冰微結(jié)構(gòu)及其及其對(duì)冰熱、力學(xué)參數(shù)影響的研究［學(xué)位論文］.大連：大連理工大學(xué)，2013.］

［7］ Liu Guangping，Du Ping，Wang Kun.Application of fractal theory to evaluation of lake eutrophication.Acta Agriculturae Universitatis Jiangxiensis，2005，27（6）：925-929.［劉光萍，杜萍，王琨.分形理論在湖泊富營養(yǎng)化評(píng)價(jià)中的應(yīng)用.江西農(nóng)業(yè)大學(xué)學(xué)報(bào)，2005，27（6）：925-929.］

［8］ Zhan Yong，Li Changyou，Shi Xiaohong etal.Application of fractal theory to the evaluation ofwater quality in Wuliangsuhai Lake.Journal ofWater Resource and Water Engineering，2012，23（2）：37-43.［詹勇，李暢游，史小紅等.分形理論在烏梁素海湖泊水質(zhì)評(píng)價(jià)中的應(yīng)用.水資源與水工程學(xué)報(bào)，2012，23（2）：37-43.］

［9］ Guillaume G，Olivier L，Aza A et al.Wastewater treatmentby radial freezing with stirring effects.Water Research，2003，37：2520-2524.DOI：10.1016／S0043-1354（03）00020-4.

［10］ WeeksWF，Ackley SF.The growth，structure and properties of sea ice.New York：Plenum Press，1989：9-164.

［11］ Weeks WF，Lee OS.Observation on the physical properties of sea ice at Hopedale，Labrador.Arctic，1958，11（3）：135-155.

［12］ Nakawo M，Sinha NK.Growth rate and salinity profile of first-sea ice in the high Arctic.Journal ofGlaciology，1981，27（96）：315-330.

［13］ GaoW，Smith DW，Sego DC.Ice nucleation in industrial wastewater.Cold Regions Science and Technology，1999，29：121-133.DOI：10.1016／S0165-232X（99）00019-1.

［14］ Song Meifeng，Liu Daoping，Wu Zhimin et al.Study on nucleation and growth mechanism of snow crystal.Journal of Refrigeration，2004，25（3）：46-50.［宋玫峰，劉道平，鄔志敏等.雪晶成核和生長機(jī)理研究.制冷學(xué)報(bào)，2004，25（3）：46-50.］

［15］ Olivier L，Pascal T，Eugenie B etal.Potential of freezing in wastewater treatment：Soluble pollutantapplication.WaterResearch，2001，35（2）：541-547.DOI：10.1016／S0043-1354（00）00287-6.

［16］ Iversen KR，Primicerio R，Larsen A etal.Effectsof small-scale turbulence on lower trophic levels under different nutrient conditions.Journal of Plankton Research，2010，32（2）：197-208.

［17］ Pennak RW.Field and experimentalwinter limnology of three Coloradomountain lakes.Ecology，1968，49（3）：505-520.

［18］ Belzile C，Gibson JAE，VincentWF.Colored dissolved organic matter and dissolved organic carbon exclusion from lake ice：Implication for irradiance transmission and carbon cycling.Limnology and Oceanography，2002，47（5）：1283-1293. DOI：10.4319／lo.2002.47.5.1283.

［19］ He Jiang，Sun Ying，LüChangwei et al.Research on phosphorus release from the surface sediments in the Daihai Lake. Acta Ecologica Sinica，2010，30（2）：389-398.［何江，孫英，呂昌偉等.岱海表層沉積物中內(nèi)源磷的釋放.生態(tài)學(xué)報(bào)，2010，30（2）：389-398.］

［20］ Petterson K.Mechanisms for internal loading of phosphorus in lakes.Hydrobiologia，1998，373／374：21-25.

［21］ Zhong Jicheng，F(xiàn)an Chengxin.Advance in the study on the effectiveness and environmental impactof sediment dredging.J Lake Sci，2007，19（1）：1-10.DOI：10.18307／2007.0101.［鐘繼承，范成新.底泥疏浚效果及環(huán)境效應(yīng)研究進(jìn)展.湖泊科學(xué)，2007，19（1）：1-10.］

［22］ Gao W，Habib M，Smith DW.Removal of organic contaminants and toxicity from industrial effluents using freeing processes.Desalination，2009，245（1）：108-119.

Effect of freezing on eutrophication in Lake Ulansuhai

ZHANG Yan1，LIChangyou2，GAO Ning1，SHIXiaohong2＆QIAO Lingm in1
（1：Civil Engineering College，Yantai University，Yantai264005，P.R.China）
（2：Water Conservancy and Civil Engineering College，Inner Mongolia Agricultural University，Huhhot 010018，P.R.China）

In order to study the effect of freezing on lake eutrophication，Lake Ulansuhai，located in Inner Mongolia，was selected as the study object.The ice and water samplewere collected respectively by syringewater sampler and ice drill.The total nitrogen，total phosphorus，CODCr，chlorophyll-a，ice thickness，water depth and transparency of each sample were measured.The eutrophic state of Lake Ulansuhai before and after freezing were assessed by themeans of fractal theory.The results showed that there was a good negative correlation between water depth and ice thickness in Lake Ulansuhai，and the correlation coefficientwas 0.54.The trophic state ofwater before freezing was eutrophication，and changed to heavy eutrophication after freezing.The results also showed that nutrient elements such as nitrogen and phosphoruswere ejected towater during the period of ice growth，that is to say，the process of freezing had a good concentrating effecton pollutants in water.Thiseffectalso canmake pollutantsmigrate from water to sediment during ice growth process，whichmade the quality of under-ice environment deteriorate.While the concentrating effect can provide theoretical basis and data support for the sediment dredging and ice-water separating during ice-on period.

Freezing；eutrophication；fractal theory；Lake Ulansuhai

DOI 10.18307／2017.0404

?2017 by Journal of Lake Sciences

?山東省優(yōu)秀中青年科學(xué)家科研獎(jiǎng)勵(lì)基金項(xiàng)目（BS2014HZ021）和國家自然科學(xué)基金項(xiàng)目（51609207，51339002）聯(lián)合資助.2016-07-22收稿；2016-10-18收修改稿.張巖（1984～），男，博士，副教授；E-mail：zhangyan-992＠163. com.