魏西會(huì)，房俊旭

（1.青島市生態(tài)環(huán)境局城陽分局，山東 青島266109；2.青島市生態(tài)環(huán)境局即墨分局，山東 青島，266200）

0 引言

墨水河是青島市北部區(qū)域的一條入海河流，發(fā)源于三標(biāo)山，流經(jīng)城陽區(qū)和即墨區(qū)，最終匯入膠州灣，全長(zhǎng)41.5 km，流域總面積317.2 km2。墨水河枯水期和豐水期水量變化顯著，是一條季節(jié)性河流[1]。2021年墨水河5—9月為豐水期、1月、2月和12月為枯水期，3月、4月及10月為平水期[2]。

河流不僅可以為流域提供豐富的水源，還可以收納流域內(nèi)陸源入河的各種污染物。長(zhǎng)期以來，國(guó)內(nèi)外學(xué)者對(duì)河流中氮磷變化特征及影響因素進(jìn)行了廣泛研究，不過研究大多集中在南方河流，針對(duì)北方河流的研究較少。目前公開發(fā)表的主要有對(duì)密云水庫(kù)入庫(kù)河流[3-4]、奎河[5]、等的研究，奎河氮磷濃度明顯高于長(zhǎng)江、黃河和珠江等南方河流。部分學(xué)者對(duì)墨水河水質(zhì)亦有研究[1]，但是主要集中在河流水質(zhì)評(píng)價(jià)方面，缺少對(duì)氮磷營(yíng)養(yǎng)鹽的變化特征及影響因素的研究。文章以墨水河為研究對(duì)象，分析了該河流氮磷元素含量的分布特征，探討了影響其變化主要因素，為墨水河污染治理及水環(huán)境安全提供科學(xué)依據(jù)。

1 實(shí)驗(yàn)與方法

1.1 水樣的采集與分析

監(jiān)測(cè)點(diǎn)設(shè)在墨水河入?？跀嗝妫∟120.3704° E36.3123°），每月開展1次監(jiān)測(cè)，采集表層（水面下0.5 m）樣品，水溫、pH值、溶解氧及電導(dǎo)率現(xiàn)場(chǎng)測(cè)定，其余項(xiàng)目帶回實(shí)驗(yàn)室分析。樣品采集、運(yùn)輸、交接及質(zhì)量控制均按照《HJ/T 91-2002 地表水和污水監(jiān)測(cè)技術(shù)規(guī)范》要求進(jìn)行?？偭撞捎肎B 11893-1989 鉬酸氨分光光度法、總氮采用 HJ 636-2012 堿性過硫酸鉀消解紫外分光光度法，氨氮采用HJ 535-2009 納氏試劑分光光度法、高錳酸鹽指數(shù)采用GB/T 11892-1989 高錳酸鹽指數(shù)的測(cè)定法、亞硝酸鹽采用GB/T 7493-1987 分光光度法、硝酸鹽、亞硝酸鹽及氟離子均采用HJ 84-2016 離子色譜法進(jìn)行測(cè)定。

1.2 數(shù)據(jù)分析

2021年共獲得12組監(jiān)測(cè)數(shù)據(jù)，每組數(shù)據(jù)包含9個(gè)參數(shù)。運(yùn)用WPS Office對(duì)氮磷進(jìn)行分布特征分析；利用IBM SPSS 26軟件，通過Pearson分析不同環(huán)境因子的相關(guān)性。

2 結(jié)果與分析

2.1 墨水河氮磷含量及時(shí)間分布特征

2021年監(jiān)測(cè)結(jié)果顯示，墨水河水質(zhì)中總氮和總磷含量變化幅度較大。從圖1可以看出，總氮有明顯的時(shí)間變化規(guī)律，冬季＞春秋季＞夏季，最大值出現(xiàn)在2月，最小值出現(xiàn)在6月?？偭讜r(shí)間變化亦有規(guī)律，夏季＞春秋季＞冬季，最大值出現(xiàn)在7月，最小值出現(xiàn)在2月。從圖2可以看出，不同水期總氮含量也有差別，枯水期高于豐水期；總磷變化規(guī)律與總氮相反，豐水期高于枯水期。

圖1 墨水河總氮、總磷含量月度變化圖

圖2 墨水河總氮、總磷含量水期變化圖

圖3 墨水河無機(jī)氮含量水期變化圖

2.2 墨水河不同形態(tài)氮含量時(shí)間分布特征

對(duì)墨水河水體中不同水期、不同形態(tài)的無機(jī)氮進(jìn)行分析，發(fā)現(xiàn)硝酸鹽氮、氨氮、亞硝酸鹽氮雖然含量相差較大，但是有同樣的時(shí)間變化特征：冬季＞春秋季＞夏季，枯水期高于豐水期，這和總氮的變化特征一致。進(jìn)一步分析發(fā)現(xiàn)，硝酸鹽氮在總氮的占比最大，在31.2%～91.6%之間，年均值為62.8%；其次為氨氮，年均占比為8.0%；亞硝酸鹽氮最小，年均占比為2.0%。

2.3 墨水河氮磷比值時(shí)間分布特征

由于水中氮磷含量比例會(huì)影響水體藻類生長(zhǎng)，因此水體中氮磷之比常被用來判別浮游藻類的營(yíng)養(yǎng)限制類型。Redfield提出氮磷比為16∶1時(shí)，可以滿足藻類快速生長(zhǎng)對(duì)氮磷營(yíng)養(yǎng)的的需求比例[6]。由圖4可以看出，墨水河氮磷比值有春冬高、夏秋低的分布特征，比值范圍為16～80，年均值為39。這說明2021全年磷都是墨水河藻類生長(zhǎng)潛在的限制性營(yíng)養(yǎng)鹽。

圖4 墨水河氮磷含量比值含量月度變化圖

3 討論

3.1 墨水河氮磷時(shí)間變化影響因素

墨水河水體中氮磷含量的時(shí)間變化規(guī)律明顯，總氮冬季高于夏季，總磷夏季高于冬季。墨水河水體中氮磷含量時(shí)間變化特征與密云水庫(kù)入庫(kù)河流[3-4]、奎河[5]、流溪河[7]、贛江[8]、龍江河[9]一致。這說明造成墨水河氮磷含量變化特征的原因可能是冬季水溫低，微生物活動(dòng)受限，分解和吸收氮的能力降低，使得水體氮濃度升高；到了夏季水溫升高，浮游植物進(jìn)入生長(zhǎng)活躍期，氮元素被吸收，加之夏季降雨頻繁，流量增加對(duì)水體中的氮含量有一定的稀釋作用。因?yàn)槟铀w長(zhǎng)期處于磷限值狀態(tài)，所以磷含量與浮游植物活躍性關(guān)系不大，主要受降雨沖刷和河水稀釋效應(yīng)雙重影響。夏季強(qiáng)降雨形成的徑流沖刷地表并攜帶泥沙中的顆粒態(tài)磷以及溶解態(tài)磷進(jìn)入墨水河，其效果遠(yuǎn)大于稀釋作用，所以夏季水體中總磷含量升高。

研究發(fā)現(xiàn)，氮磷濃度與流量的關(guān)系能反映它們的主要來源[7，10]?？傮w上，墨水河水體中氮元素含量枯水期高于豐水期，而磷元素則是豐水期高于枯水期。進(jìn)一步分析發(fā)現(xiàn)，墨水河水體中總氮含量隨流量增加而降低而后趨于穩(wěn)定，磷濃度則是先隨流量增加而增加，而后同樣趨于穩(wěn)定。墨水河總氮、總磷含量變化特征與Kalkhoff等人對(duì)面源污染河流氮磷含量變化特征的研究結(jié)果基本一致[10]。這說明，墨水河氮磷元素可能主要來自面源污染。

表1 墨水河環(huán)境指標(biāo)與氮磷元素相關(guān)性分析

墨水河氮污染嚴(yán)重，全年總氮濃度均超過國(guó)家地表水V類標(biāo)準(zhǔn)。墨水河水體中不同形態(tài)無機(jī)氮含量相差較大，硝酸鹽氮是氮元素的主要存在形式（占比年均值為62.8%），這一結(jié)果和密云水庫(kù)入庫(kù)河流[3-4]、奎河[5]、龍江河[9]一致。這進(jìn)一步表明，墨水河水質(zhì)受到農(nóng)業(yè)面源污染，此外可能還受到河兩岸點(diǎn)源污染排放、河底底泥釋放、水生生物作業(yè)等影響。

3.2 墨水河氮磷及環(huán)境因子相關(guān)性分析

氮磷是浮游植物生長(zhǎng)所必須的營(yíng)養(yǎng)元素，氨氮和高錳酸鹽指數(shù)是水體有機(jī)污染物的重要指標(biāo)[11]。選取總氮、總磷、氨氮、高錳酸鹽指數(shù)以及環(huán)境指標(biāo)進(jìn)行相關(guān)性分析，結(jié)果顯示高錳酸鹽指數(shù)與水溫呈現(xiàn)顯著正相關(guān)（P＜0.05），與總氮呈現(xiàn)顯著負(fù)相關(guān)（P＜0.01）；電導(dǎo)率與總磷呈現(xiàn)顯著負(fù)相關(guān)（P＜0.05）；水溫與氨氮、總氮呈現(xiàn)顯著負(fù)相關(guān)（P＜0.01），與總磷呈現(xiàn)顯著正相關(guān) （P＜0.01），這表明水溫對(duì)于自然界中氮磷的循環(huán)產(chǎn)生顯著影響?？偟桶钡嚓P(guān)性不明顯，說明總氮和氨氮可能來源不同。

4 結(jié)論

（1）墨水河水體中氮磷含量的時(shí)間變化規(guī)律明顯。總氮規(guī)律為冬季＞春秋季＞夏季，總磷則相反：夏季＞春秋季＞冬季。不同水期總氮含量也有差別，枯水期高于豐水期，總磷為豐水期高于枯水期。墨水河氮磷元素可能主要來自面源污染。

（2）墨水河水體中硝酸鹽氮、氨氮、亞硝酸鹽氮時(shí)間變化特征與總氮相似：冬季＞春秋季＞夏季，枯水期高于豐水期，但含量相差較大。硝酸鹽氮在總氮的占比最大，年均值為62.8%；其次為氨氮，亞硝酸鹽氮最小。

（3）水溫與氨氮、總氮呈現(xiàn)顯著負(fù)相關(guān)，與總磷呈現(xiàn)顯著正相關(guān)，這表明水溫對(duì)于自然界中氮磷的循環(huán)產(chǎn)生顯著影響?？偟桶钡嚓P(guān)性不明顯，說明總氮和氨氮可能來源不同。