黃玉潔，張銀龍，李海東，孫雅輝，張 波，楊均科

（1．南京林業(yè)大學(xué) 森林資源與環(huán)境學(xué)院，南京210037；2．環(huán)境保護(hù)部 南京環(huán)境科學(xué)研究所，南京210042）

水體富營(yíng)養(yǎng)化是一種水體衰老現(xiàn)象，系指氮、磷等植物性營(yíng)養(yǎng)物大量進(jìn)入湖泊、海灣等相對(duì)封閉、水流緩慢的水體，引起藻類和其它水生植物大量繁殖，水體透明度和溶解氧下降，水質(zhì)惡化，水生生物大量死亡，水體生態(tài)系統(tǒng)和水功能受到阻礙和破壞的現(xiàn)象［1－3］。國(guó)內(nèi)外學(xué)者對(duì)高等水生植物凈化富營(yíng)養(yǎng)化水體進(jìn)行了大量研究［4－9］，結(jié)果顯示，高等水生植物能夠大量吸收水體中的氮、磷營(yíng)養(yǎng)物質(zhì)，能有效地凈化富營(yíng)養(yǎng)化水體，有助于水體生態(tài)系統(tǒng)恢復(fù)正常的功能。

太湖流域隨著工農(nóng)業(yè)生產(chǎn)的迅速發(fā)展以及城市化進(jìn)程的加快，大量含有氮、磷等營(yíng)養(yǎng)元素工業(yè)廢水和生活污水直接排入到河流、湖泊中去，使得藻類大量繁殖，造成水體富營(yíng)養(yǎng)化。湖泊濕地在物質(zhì)循環(huán)、凈化水質(zhì)、降解環(huán)境污染、調(diào)節(jié)氣候等方面，具有其它生態(tài)系統(tǒng)不可替代的作用［10－11］。實(shí)踐證明，人工濕地植物在凈化受污染水體，恢復(fù)和重建濕地生態(tài)系統(tǒng)已經(jīng)獲得了廣泛共識(shí)。已有文獻(xiàn)對(duì)沉積物中的氮、磷含量進(jìn)行了報(bào)道、研究［12－16］，結(jié)果表明：濕地植物的莖葉可減小水流速度，降低水的輸送能力，進(jìn)而將氮、磷元素通過(guò)促淤效應(yīng)沉積于底泥中，其經(jīng)化學(xué)和生物學(xué)過(guò)程轉(zhuǎn)化，最終被植物吸收、儲(chǔ)存起來(lái)；沉積物中氮、磷含量的高低往往與濕地植物類型、密度和生物量密切相關(guān)。因此，研究不同濕地植物群落的建植對(duì)氮、磷的攔截、沉積作用效果，對(duì)太湖流域濕地植被恢復(fù)工程建設(shè)和水體富營(yíng)養(yǎng)化防治具有重要意義。

1 研究區(qū)域與研究方法

1．1 研究區(qū)概況

太湖為我國(guó)第三大淡水湖，位于長(zhǎng)江下游地區(qū)（33°55′42″－31°33′50″N，119°53′45″－120°36′15″E），全流域面積36 900 k m2，屬亞熱帶季風(fēng)氣候區(qū)，四季分明，雨水豐沛，熱量充裕，多年平均氣溫15～17℃，降水量1 177 mm。

研究區(qū)位于太湖百瀆港，是太湖西北角上唯一的河港，其周邊城區(qū)和鄉(xiāng)鎮(zhèn)工廠遍布，大量工業(yè)廢水、生活污水排入其中，經(jīng)監(jiān)測(cè)上游太滆運(yùn)河屬于劣Ⅴ類水質(zhì)［17］，是研究區(qū)外源污染物輸入的主要通道。2005年當(dāng)?shù)卣槍?duì)水域污染問(wèn)題實(shí)施了濕地植被恢復(fù)工程，工程先進(jìn)行了一定程度的清淤，后以蘆葦（Phr ag mites austr alis）、茭 草 （Zizania caducif lor a）、香蒲（Typha l atif olia）、荷花（Nel umbo nucifer a）和沉水植物作為建群種，植物群落從西南向東北大體分布依次為：沉水植物群落、香蒲群落、茭草群落、蘆葦群落和荷花群落，并在植物群落周圍設(shè)置了木樁消浪帶，使區(qū)域內(nèi)除河口位置外水流比較平緩，對(duì)植物的生長(zhǎng)起到很好的保護(hù)作用。經(jīng)過(guò)幾年后現(xiàn)已基本形成穩(wěn)定的濕地植物群落，蘆葦群落平均高度為2．8 m，平均生物量（鮮重）4．70 kg／m2，茭草群落平均高度為1．6 m，平均生物量3．26 kg／m2，香蒲群落平均高度為2．3 m，平均生物量3．68 kg／m2，荷花群落平均高度為1．3 m，平均生物量1．25 kg／m2，沉水植物包括馬來(lái)眼子菜（Pota mogeton mal aianus）、菹草（Potamogeton crispus）、苦草（Vallisneria spiralis）、金魚藻（Cer atophyll u m demersu m）、輪葉黑藻（Hydrill a verticill ata）等，平均生物量1．13 kg／m2。

1．2 研究方法

于2009年7月，在濕地植被建植區(qū)（由于植物比較茂密，一些區(qū)域難以進(jìn)入，只能粗略勾繪出采樣點(diǎn)的植被分布情況），采用GPS定位共設(shè)58個(gè)采樣點(diǎn)（圖1），樣點(diǎn)的布設(shè)充分考慮到濕地植物的類型、群落的分布情況，使各個(gè)群落內(nèi)的采樣點(diǎn)基本上能代表該區(qū)域沉積物的狀況。其中，蘆葦群落共有26個(gè)采樣點(diǎn)（樣點(diǎn)號(hào)為1－5，16－18，31－32，35－40，44－47，52－56，58），茭草群落共有10個(gè)采樣點(diǎn)（樣點(diǎn)號(hào)為6－9，26－28，49－50，57）；香蒲群落共有5個(gè)采樣點(diǎn)（樣點(diǎn)號(hào)為10－13，48），荷花群落共有7個(gè)采樣點(diǎn)（樣點(diǎn)號(hào)為14－15，33－34，41－43），沉水植物群落共有7個(gè)采樣點(diǎn)（樣點(diǎn)號(hào)為19－22，29－30，51）；并在分布少量植物的開(kāi)闊水面設(shè)3個(gè)采樣點(diǎn)（樣點(diǎn)號(hào)為23－25）作為對(duì)照區(qū)。

圖1 太湖百瀆港采樣點(diǎn)分布

用直徑10 c m的PVC塑料管采集濕地表層10 c m的沉積物，樣品帶回實(shí)驗(yàn)室，所有樣品經(jīng)風(fēng)干、剔除雜質(zhì)、研磨和過(guò)篩后，用于測(cè)定沉積物的全氮、全磷和有效磷。其中，全氮采用半微量凱氏法；全磷采用硫酸－高氯酸混合酸做消化劑，鉬銻抗比色法測(cè)定全磷含量；有效磷采用碳酸氫鈉浸提，鉬銻抗比色法測(cè)定［18］。

數(shù)據(jù)處理主要采用地統(tǒng)計(jì)學(xué)方法的克里格（Kriging）插值法，Kriging插值方法是利用原始數(shù)據(jù)和半方差函數(shù)的結(jié)構(gòu)性，對(duì)未采樣點(diǎn)的區(qū)域化變量進(jìn)行最優(yōu)無(wú)偏估值［19－20］。

2 結(jié)果與分析

2．1 沉積物中氮、磷經(jīng)典統(tǒng)計(jì)特征

研究區(qū)沉積物中全氮、全磷和有效磷含量的經(jīng)典描述性統(tǒng)計(jì)見(jiàn)表1，可以看出，全氮、全磷和有效磷含量變化范圍分別為0．37～5．42 g／kg，1．01～9．06 g／kg，150～760 mg／kg，最大值分別是最小值的14．65倍、8．97倍和5．07倍，差異較大。全氮平均值和中值相差較大，全磷和有效磷相差較小，有效磷均值和中值幾乎相等。從變異系數(shù)來(lái)看，全氮變異系數(shù)為79．5%，遠(yuǎn)大于全磷（33．8%）和有效磷（28．5%），按照變異系數(shù)的劃分等級(jí)：弱變異性，CV＜10%；中等變異性，CV＝10%～100%；強(qiáng)變異性，CV＞100%［20］；全氮、全磷和有效磷均表現(xiàn)為中等程度的變異性。

從統(tǒng)計(jì)平均值來(lái)看，研究區(qū)全氮的平均值為1．25 g／kg，低于加拿大安大略省環(huán)境和能源部發(fā)布的具有嚴(yán)重級(jí)別生態(tài)風(fēng)險(xiǎn)效應(yīng)沉積物中全氮濃度（4．80 g／kg）［21］，最大值為5．42 g／kg也低于富營(yíng)養(yǎng)化嚴(yán)重的滇池沉積物中全氮的最高值（9．10 g／kg）［22］，可見(jiàn)沉積物中全氮含量不是很高，同時(shí)亦表明濕地水體氮素的富營(yíng)養(yǎng)化程度也不是特別嚴(yán)重。全磷和有效磷的平均值分別為5．32 g／kg和420 mg／kg，根據(jù)加拿大安大略省環(huán)境和能源部發(fā)布的沉積物中能引起最低級(jí)別生態(tài)毒性效應(yīng)的全磷濃度及具有嚴(yán)重級(jí)別生態(tài)毒性效應(yīng)的全磷濃度標(biāo)準(zhǔn)（分別為0．6 g／kg和2．0 g／kg）［23］，本研究中全磷的平均值遠(yuǎn)超過(guò)嚴(yán)重級(jí)別生態(tài)毒性效應(yīng)的濃度，全磷的最小值也高于最低級(jí)別生態(tài)毒性效應(yīng)的濃度，全磷變化范圍亦比其他非富營(yíng)養(yǎng)化湖泊沉積物中全磷的變化范圍高得多（洱海0．74～1．60 g／kg［24］、瓦埠湖0．29～0．87 g／kg［25］、梁子湖0．635～0．824 g／kg［26］），可見(jiàn)沉積物中全磷含量很高，同時(shí)亦表明濕地水體磷素的富營(yíng)養(yǎng)化程度特別嚴(yán)重。從全氮和全磷在沉積物中的含量來(lái)看，其主要來(lái)源于外源輸入，且對(duì)水體的富營(yíng)養(yǎng)化狀況全磷的貢獻(xiàn)要比全氮大。

表1 沉積物中全氮、全磷和有效磷的描述性統(tǒng)計(jì)

2．2 正態(tài)分布性檢驗(yàn)及半方差函數(shù)分析

數(shù)據(jù)的正態(tài)分布性是使用地統(tǒng)計(jì)學(xué)方法進(jìn)行空間分析的前提，為了檢驗(yàn)這些實(shí)驗(yàn)數(shù)據(jù)是否符合地統(tǒng)計(jì)學(xué)插值條件，采用偏斜度法和非參數(shù)Kol mogor ov－Smironov（K－S）正態(tài)分布檢驗(yàn)概率對(duì)全氮、全磷和有效磷進(jìn)行正態(tài)分布檢驗(yàn)（見(jiàn)表1），對(duì)不符合正態(tài)分布的變量數(shù)據(jù)進(jìn)行對(duì)數(shù)轉(zhuǎn)換后再進(jìn)行檢驗(yàn)。全氮的偏度系數(shù)絕對(duì)值大于1，具有偏態(tài)趨勢(shì)。經(jīng)K－S檢驗(yàn)，全氮不能滿足正態(tài)分布的要求，全磷和有效磷均為正態(tài)分布形式。對(duì)全氮含量進(jìn)行對(duì)數(shù)轉(zhuǎn)換后繼續(xù)檢驗(yàn)，達(dá)到正態(tài)分布的要求。

半方差函數(shù)模型的選擇是對(duì)沉積物全氮、全磷和有效磷空間變異結(jié)構(gòu)分析的關(guān)鍵，采用GS＋軟件進(jìn)行半方差理論函數(shù)模型選擇，經(jīng)比對(duì)不同模型產(chǎn)生的R2、RSS等參數(shù)，分別選擇球狀模型、指數(shù)模型為全氮、全磷和有效磷的半方差函數(shù)理論模型，見(jiàn)表2。按照區(qū)域化變量空間相關(guān)性程度的分級(jí)標(biāo)準(zhǔn)，當(dāng)C0／C0＋C）＜25%，表明變量的空間變異以結(jié)構(gòu)性變異為主，變量具有強(qiáng)烈的空間相關(guān)性；在25%～75%之間，說(shuō)明變量具有中等程度的空間相關(guān)性；當(dāng)比值＞75%時(shí)，表明變量以隨機(jī)性變異為主，變量的空間相關(guān)性很弱；如果比值接近于l，則說(shuō)明變量在整個(gè)尺度上具有恒定的變異［27］。由表2可以看出，全氮的C0／（C0＋C）值為20．10%，小于25%，說(shuō)明結(jié)構(gòu)性因素對(duì)其空間變異分布影響較大。全磷和有效磷的C0／（C0＋C）分別為42．31%和37．97%，在25%～75%之間，表明全磷和有效磷含量的空間分布具有中等程度空間相關(guān)性，隨機(jī)效應(yīng)與結(jié)構(gòu)效應(yīng)相接近。決定系數(shù)R2較大、殘差RSS較小表明半方差函數(shù)理論模型擬合效果較好，適合于采用Kriging局部?jī)?nèi)插法進(jìn)行空間預(yù)測(cè)。

表2 全氮、全磷和有效磷的半方差函數(shù)理論模型及其相關(guān)參數(shù)

2．3 沉積物中氮、磷空間分布特征

選擇上述半方差模型及相關(guān)參數(shù)，根據(jù)有限的空間采樣點(diǎn)，運(yùn)用Kriging插值方法對(duì)未抽取樣點(diǎn)的氮、磷含量進(jìn)行最優(yōu)無(wú)偏估計(jì)，繪制出研究區(qū)全氮、全磷和有效磷空間插值分布圖（圖2）。

研究區(qū)沉積物中全氮含量自南向北逐漸增加，低值出現(xiàn)在西南部，高值出現(xiàn)在西北部，并形成了幾個(gè)比較明顯的島狀高值區(qū)域（見(jiàn)圖2a），以45，46，58號(hào)采樣點(diǎn)和31，52號(hào)采樣點(diǎn)周邊區(qū)域最高，其次為17，52，53號(hào)和35，37，56，55號(hào)采樣點(diǎn)周邊區(qū)域，最低區(qū)域主要分布在23，24，25號(hào)采樣點(diǎn)周邊。太滆運(yùn)河河口位于1號(hào)采樣點(diǎn)附近，該點(diǎn)附近全氮含量并不是最高，因此處水流較快，不利于元素的沉積，而在離河口稍遠(yuǎn)，水流變平緩的45，46，58號(hào)和31，52號(hào)采樣點(diǎn)周邊的蘆葦群落分布區(qū)域含量最高；同時(shí)分布著蘆葦群落的35，55，56號(hào)采樣點(diǎn)周邊區(qū)域，離河口位置較遠(yuǎn)，全氮含量卻明顯高于離河口較近的其他群落分布區(qū)，但低于位于河口附近的蘆葦群落分布區(qū)；由此可以得出影響研究區(qū)全氮含量的主要因素有兩個(gè)，一是濕地植物群落與河口（污染源）的位置關(guān)系，靠近河口的植物能先攔截、沉積外源污染物，對(duì)比圖1中茭草群落的兩個(gè)分布區(qū)域也能驗(yàn)證了此結(jié)論，離河口較近的茭草群落區(qū)大于離河口遠(yuǎn)的茭草群落區(qū)；二是濕地植物的類型，蘆葦相比其他幾種植物，其莖葉茂密、根系發(fā)達(dá)、生物量大，這就使蘆葦群落對(duì)全氮的攔截、沉積能力最強(qiáng)，進(jìn)而造成研究區(qū)蘆葦群落沉積物中全氮含量整體最高。荷花群落分布區(qū)與離河口位置較近的茭草群落分布區(qū)域（6－9號(hào)采樣點(diǎn)周邊區(qū)域）含量大體相當(dāng)，高于香蒲群落、沉水植物群落分布區(qū)和離河口位置較遠(yuǎn)的茭草群落分布區(qū)域（26－28號(hào)、49－50號(hào)采樣點(diǎn)周邊區(qū)域），后三者之間區(qū)分的不是很明顯，這都是植物密度、分布位置及自身生物學(xué)等因素綜合作用的結(jié)果。對(duì)照區(qū)植物較少，攔截、沉積效果差，且離河口最遠(yuǎn)，其含量自然最低。

全磷總體分布與全氮具有一定的相似性，從西南向東北先增加后降低，在蘆葦群落分布區(qū)達(dá)到最高，之后在荷花群落分布區(qū)又有所降低；最高值也形成兩個(gè)島狀區(qū)域，分別是4，5，18，54號(hào)采樣點(diǎn)落的在西北部島狀區(qū)域和35，36，56采樣號(hào)點(diǎn)落在的東南部島狀區(qū)域；最低值主要在24，25號(hào)兩采樣點(diǎn)周邊區(qū)域（見(jiàn)圖2b）。由此可見(jiàn)，濕地植物對(duì)全磷的攔截、沉積效果最好的仍是蘆葦群落，其次茭草群落。全氮和全磷分布情況的基本一致性，反映了濕地植物對(duì)沉積物中氮、磷含量高低作用原理的共性。

圖2 研究區(qū)沉積物全氮、全磷和有效磷插值分布

有效磷在空間分布上與全氮和全磷相差較大，基本趨于相反的方向，從南向北成金字塔形依次減少，高值出現(xiàn)在23，24，25號(hào)采樣點(diǎn)周邊區(qū)域（見(jiàn)圖2c）。有效磷與全磷、全氮相關(guān)分析結(jié)果表明，全磷和有效磷在α＝0．01的水平上具有極顯著負(fù)相關(guān)，相關(guān)系數(shù)為－0．316，說(shuō)明有效磷隨著全磷的增加而逐漸減少。結(jié)合濕地植物分布情況，從茭草群落、沉水植物群落、香蒲群落、蘆葦群落到荷花群落有效磷含量大體依次減少，對(duì)照區(qū)最高。分析原因，有效磷是可被植物吸收的磷組分，在全磷中的比例很小，對(duì)照區(qū)植物分布較少，對(duì)有效磷的吸收能力很弱，有效磷占全磷的比例較大，致使該區(qū)含量最高；濕地植物群落分布區(qū)全磷含量較高，并且植物根際微生物很活躍，對(duì)磷的活化作用也很強(qiáng)，但采樣正值植物生長(zhǎng)最旺盛的7月份，處于對(duì)有效磷吸收最強(qiáng)的時(shí)期，因此致使?jié)竦刂参锶郝浞植紖^(qū)低于對(duì)照區(qū)。特別是蘆葦群落分布區(qū)，全磷含量最高而有效磷較低，說(shuō)明蘆葦群落對(duì)有效磷的吸收能力最強(qiáng)。

3 討 論

通過(guò)上述分析，人工濕地植物群落的建植明顯增大沉積物中氮、磷的含量，濕地植物通過(guò)攔截、沉積作用，將氮、磷聚集促積在沉積物中，再經(jīng)生物和化學(xué)轉(zhuǎn)化過(guò)程，主要是微生物的降解作用，轉(zhuǎn)化成能夠被植物吸收的組分，最后儲(chǔ)存在植物體中，起到凈化水環(huán)境的功能。但氮、磷被攔截、沉積并被吸收、儲(chǔ)存到植物體后，并不能保證其就可以從水體中去除，首先，在植物生長(zhǎng)周期過(guò)程中，凋落、腐爛的代謝體能將被吸收氮、磷重新釋放到水體中，最終會(huì)被再次沉積到沉積物中，通過(guò)調(diào)查發(fā)現(xiàn)，自濕地植被恢復(fù)工程以來(lái)，一直未對(duì)植物進(jìn)行過(guò)收獲，這就造成了研究區(qū)氮、磷在沉積物中一直處于積累狀態(tài)；其次，當(dāng)?shù)⒘缀砍^(guò)臨界飽和值時(shí)，沉積物反而作為內(nèi)源污染源，會(huì)將它們釋放到水體中，引起水體的富營(yíng)養(yǎng)化。所以，及時(shí)更新、收獲濕地植物，源源不斷地把氮、磷元素轉(zhuǎn)移出濕地，并采取疏浚清理沉積物的方式是防止水體再次被污染的重要措施，這樣才能達(dá)到最終凈化濕地水體，恢復(fù)濕地生態(tài)功能的目的。

4 結(jié) 論

（1）全氮、全磷和有效磷含量的變化范圍分別為0．37～5．42 g／kg，1．01～9．06 g／kg和150～760 mg／kg，差異較大，在空間上均呈中等程度的變異性。沉積物中氮素含量相對(duì)不是特別高，全氮平均值為1．25 g／kg。磷素含量相對(duì)很高，全磷平均值為1．25 g／kg，有效磷平均值為420 mg／kg。

（2）全氮自南向北逐漸增加，低值出現(xiàn)在西南部，高值出現(xiàn)在西北部；全磷從西南向東北先增加后降低，總體分布與全氮基本一致，說(shuō)明了濕地植物對(duì)氮、磷含量高低作用原理的共性；有效磷分布與全氮和全磷相差較大，基本趨于相反的方向，從南向北成金字塔形依次減少。

（3）人工濕地植物群落的建植明顯增大了沉積物中氮、磷含量，蘆葦群落對(duì)氮、磷的攔截、沉積、吸收作用最強(qiáng)，茭草群落、香蒲群落、荷花群落和沉水植物群落由于受分布密度、位置以及自身生物學(xué)等因素的影響，表現(xiàn)出一定的差異性，但都優(yōu)于對(duì)照區(qū)。

［1］ 高建平，王珊玲．水體富營(yíng)養(yǎng)化評(píng)價(jià)和防治的一些進(jìn)展［J］．農(nóng)村生態(tài)環(huán)境，1989（3）：55－60．

［2］ 古濱河．美國(guó)Apopka湖的富營(yíng)養(yǎng)化及其生態(tài)修復(fù)［J］．湖泊科學(xué)，2005，17（1）：1－8．

［3］ 史丹．我國(guó)湖泊富營(yíng)養(yǎng)化問(wèn)題及防治對(duì)策［J］．資源開(kāi)發(fā)與市場(chǎng)，2005，21（1）：17．

［4］ Holland D．Down the riverside［J］．Landscape Design，2002，11：65－67．

［5］ 謝建華，楊華．不同植物對(duì)富營(yíng)養(yǎng)化水體凈化的靜態(tài)試驗(yàn)研究［J］．工業(yè)安全與環(huán)保，2006，32（6）：23－25．

［6］ Robyn A O，David L P．The uptake of uraniu m by Eleocharis dulcis （Chinese water chestnut）in the Ranger Uranium Mine constructed wetland filter［J］．Environ mental Pollution，2004，132：307－320．

［7］ Groudeva V I，Groudev S N，Doycheva A S．Bioremediation of waters contaminated with cr ude oil and toxic heavy metals［J］．Inter national Jour nal of Mineral Processing，2001，62：293－299．

［8］ Todd L I，Lawrence A R．Nitrate removal in wetland microcos ms［J］．Water Research，1998，32（3）：677－684．

［9］ 孔繁翔．大型淺水富營(yíng)養(yǎng)化湖泊中藍(lán)藻水華形成機(jī)理的思考［J］．生態(tài)學(xué)報(bào)，2005，25（3）：589－595．

［10］ 萬(wàn)曉紅，周懷東，劉玲花，等．白洋淀湖泊濕地中氮素分布的初步研究［J］．水土保持學(xué)報(bào)，2008，22（2）：166－169．

［11］ van Hees P A W，Jones D L，Nyberg L，et al．Modelling low molecular weight organic acid dynamics in forest soils［J］．Soil Biology ＆ Biochemistr y，2005，37：517－531．

［12］ 朱元榮，張潤(rùn)宇，吳豐昌．滇池沉積物生物有效性氮和磷的分布及相互關(guān)系［J］．環(huán)境科學(xué)研究，2010，23（1）：993－998．

［13］ 李文朝．東太湖水生植物的促淤效應(yīng)與磷的沉積［J］．環(huán)境科學(xué)，1997，18（3）：11－13．

［14］ Luijn F V，Boers P C M，Lijklema L，et al．Nitrogen fluxes and pr ocesses in sandy and muddy sedi ments fr om a shalow lake［J］．Water Research，1999，33（1）：33－42．

［15］ Sundby B，Gobeil C，Silberberg N．The phosphorcycle in coastal marine sedi ment［J］． Li mnology and Ocanogrphy，1992，37（6）：1129－1145．

［16］ 熊漢鋒，譚啟玲，王運(yùn)華．梁子湖沉積物中氮磷分布特征研究［J］．華中農(nóng)業(yè)大學(xué)學(xué)報(bào)，2008，27（2）：235－238．

［17］ 孫衛(wèi)紅，張利民，劉偉京，等．漕橋河區(qū)域水環(huán)境質(zhì)量現(xiàn)狀分析［J］．環(huán)境科學(xué)與管理，2009，34（1）：48－51．

［18］ 劉焱．中國(guó)林業(yè)標(biāo)準(zhǔn)匯編（營(yíng)造林卷）［M］．北京：中國(guó)標(biāo)準(zhǔn)出版社，1998．

［19］ 雷志棟，楊詩(shī)秀，許志榮，等．土壤特性空間變異性初步研究［J］．水利學(xué)報(bào)，1985（9）：10－21．

［20］ 李海東，張波，沈渭?jí)?，等．蘇南山丘區(qū)小流域土壤養(yǎng)分特性空間分布［J］．長(zhǎng)江流域資源與環(huán)境，2009，18（9）：831－835．

［21］ 陳永川，湯利，張德剛，等．滇池沉積物總氮的時(shí)空變化特征研究［J］．土壤，2007，39（6）：879－883．

［22］ 岳維忠，黃小平，孫翠慈．珠江口表層沉積物中氮、磷的形態(tài)分布特征及污染評(píng)價(jià)［J］．海洋與湖沼，2007，38（2）：111－116．

［23］ 孫惠民，何江，高興東，等．烏梁素海沉積物中全磷的分布特征［J］．沉積學(xué)報(bào)，2006，24（14）：579－584．

［24］ 王雨春，此里能布，馬根連，等．洱海沉積物磷的化學(xué)賦存形態(tài)研究［J］．中國(guó)水利水電科學(xué)研究院學(xué)報(bào)，2005，3（2）：150－154．

［25］ 潘成榮，李凌，葉琳琳，等．瓦埠湖沉積物中氮與磷賦存形態(tài)分析［J］．水資源保護(hù)，2007，23（4）：10－14．

［26］ 熊漢鋒，陳治平，黃世寬，等．梁子湖水體P的季節(jié)變化與沉積物P釋放初步研究［J］．濕地科學(xué)，2006，4（3）：174－178．

［27］ 鄭海龍，陳杰，鄧文靖，等．城市邊緣帶土壤重金屬空間變異及其污染評(píng)價(jià)［J］．土壤學(xué)報(bào)，2006，43（1）：40－44．