基于陸海同步調(diào)查的萊州灣營養(yǎng)鹽時空分布及限制因子分析?

梁生康， 李姍姍， 馬浩陽， 楊燕群， 呂浩然， 許澤昊， 段曉萌

(中國海洋大學(xué) 1.化學(xué)化工學(xué)院； 2.海洋化學(xué)理論與工程技術(shù)教育部重點實驗室， 山東 青島 266100)

營養(yǎng)鹽是浮游植物生命活動的物質(zhì)基礎(chǔ)，浮游植物按一定比例(Redfield比)攝取營養(yǎng)物質(zhì)維持自身物質(zhì)循環(huán)和能量代謝，營養(yǎng)鹽水平會影響海域的浮游植物生物量及初級生產(chǎn)力水平；而營養(yǎng)鹽結(jié)構(gòu)會影響浮游植物的組成[1-3]。適宜的營養(yǎng)鹽水平及結(jié)構(gòu)能促進浮游植物的生長和繁殖，反之則會影響海洋生態(tài)系統(tǒng)結(jié)構(gòu)，甚至引發(fā)富營養(yǎng)化、赤潮及水體缺氧等諸多生態(tài)環(huán)境問題，并影響食物網(wǎng)及漁業(yè)資源[4-6]。

萊州灣是渤海的三大海灣之一，其生態(tài)環(huán)境質(zhì)量對環(huán)萊州灣社會經(jīng)濟可持續(xù)發(fā)展具有重要的支撐作用。萊州灣受人類活動影響顯著，通過小清河、濰河、白浪河及黃河等十余條河流承接了山東省境內(nèi)的濟南、淄博、濰坊等14個地級市排放的氮、磷營養(yǎng)物質(zhì)及其他污染物。自20世紀(jì)70年代末以來，隨著環(huán)萊州地區(qū)社會經(jīng)濟的高速發(fā)展，萊州灣環(huán)境質(zhì)量持續(xù)惡化，突出表現(xiàn)在水質(zhì)污染嚴重[2]，富營養(yǎng)化加劇，營養(yǎng)鹽結(jié)構(gòu)失衡，浮游生態(tài)系統(tǒng)物種多樣性降低等方面[7]，對環(huán)萊州灣地區(qū)的社會經(jīng)濟高質(zhì)量發(fā)展形成剛性制約。自2000年代以來，圍繞萊州灣環(huán)境質(zhì)量改善，相關(guān)部門相繼開始實施入海河流整治、海岸線控制等綜合整治措施[8-12]，特別是于2011年開始實行陸源氨氮(NH4-N)入海排放總量控制措施，按照國家行政區(qū)等比例減排考核指標(biāo)，年均減排率約3%[13]。然而，據(jù)《2018年中國海洋生態(tài)環(huán)境質(zhì)量狀況公報》顯示，近年來萊州灣水質(zhì)呈現(xiàn)出惡化勢頭趨緩，固穩(wěn)向好趨勢卻并不明顯，溶解無機氮(DIN)超標(biāo)嚴重，富營養(yǎng)化水平居高不下；這說明萊州灣生態(tài)環(huán)境問題依然嚴重且難以破解，其中重要原因是陸源營養(yǎng)物質(zhì)的排放與海域響應(yīng)關(guān)系不明[14]。

已有眾多學(xué)者就萊州灣營養(yǎng)鹽入海通量、時空分布及對海域生態(tài)環(huán)境的影響進行了較為深入研究[15-17]。調(diào)查監(jiān)測數(shù)據(jù)顯示，近40年來萊州灣陸源營養(yǎng)物質(zhì)的入海排放通量發(fā)生顯著變化，突出表現(xiàn)在總?cè)芙鈶B(tài)氮(TDN)入海排放通量逐年增加而總?cè)芙鈶B(tài)磷(TDP)入海排放通量則呈現(xiàn)由1980—1995年間增大然后自1995年后逐步減小的跡象[2]；與此同時，萊州灣海域營養(yǎng)鹽濃度和結(jié)構(gòu)也發(fā)生顯著變化，突出表現(xiàn)在DIN濃度持續(xù)增大而活性磷酸鹽(PO4-P)濃度有所降低[2,18-19]，DIN/PO4-P比例由1980年代的23∶1左右持續(xù)增大到2010年代中期的110∶1左右[2]，遠超Redfield比值，結(jié)果可能對萊州灣浮游植物生物量和優(yōu)勢種產(chǎn)生顯著影響。由于陸源入海通量和近海輸運-遷移-轉(zhuǎn)化過程都處在動態(tài)變化，只有采用陸海同步調(diào)查的精細化數(shù)據(jù)，才有可能全面系統(tǒng)剖析陸源營養(yǎng)物質(zhì)輸運與近海營養(yǎng)物質(zhì)濃度的關(guān)系。然而，由于陸海調(diào)查監(jiān)測時間不同步和調(diào)查監(jiān)測站位布設(shè)不銜接等原因，結(jié)果導(dǎo)致陸源營養(yǎng)物質(zhì)入海通量監(jiān)測數(shù)據(jù)與近海營養(yǎng)鹽時空分布的海上調(diào)查數(shù)據(jù)之間缺乏匹配性，難以進行量化分析。

本文通過2018—2021年度不同季節(jié)在萊州灣開展的陸海時間同步空間對接的調(diào)查，查明萊州灣陸源主要入海河流及直排海口的營養(yǎng)鹽入海通量和海域營養(yǎng)鹽的分布和組成的季節(jié)變化特征，解析萊州灣浮游植物生長的營養(yǎng)鹽限制性因素，為萊州灣富營養(yǎng)化的陸海統(tǒng)籌調(diào)控管理提供直接調(diào)查監(jiān)測依據(jù)。

1 材料和方法

1.1 萊州灣海域及其高聯(lián)陸域簡況

萊州灣是典型的半封閉海灣，面積約7 000 km2，約占渤海海域總面積的10%，水深較淺，大部分海域水深小于15 m。萊州灣灣內(nèi)水動力較弱，灣內(nèi)常年存在一個逆時針環(huán)流，水交換時間長，灣底平均交換時間為40 d以上[20]。入灣河流包括黃河、小清河、濰河等20余條，在山東境內(nèi)集水面積約為46 056 km2，涵蓋山東境內(nèi)的菏澤、聊城、濟寧、德州、泰安、濟南、淄博、濰坊、濱州、東營、煙臺、青島、日照13個地級市。

1.2 調(diào)查站位和時間

通過陸海時間同步、空間對接調(diào)查，獲得環(huán)萊州灣河流營養(yǎng)鹽入海通量和海域營養(yǎng)鹽分布的高質(zhì)量數(shù)據(jù)。在調(diào)查站位布設(shè)方面，入海河口及直排海口與海上水質(zhì)響應(yīng)區(qū)在空間上有機對接[21]。其中，陸上監(jiān)測站位布設(shè)于黃河、廣利河、小清河、彌河、白浪河、虞河、濰河、膠萊河、界河、泳汶河這10條河流入?？诩皟蓚€直排?？?萊州灣污水處理廠和龍口市第二污水處理廠排放口)，這些河流和直排海口的徑流量和污染物入海通量之和達到了萊州灣陸源排放總量的95%以上。在萊州灣海域共布設(shè)48個站位，不僅覆蓋整個萊州灣海域，而且考慮到海上水質(zhì)響應(yīng)區(qū)生態(tài)環(huán)境要素變化梯度大，從而進行加密設(shè)置。其中R站位位于水質(zhì)響應(yīng)區(qū)邊界點，M站位位于混合區(qū)邊界點，H站位位于污染熱點區(qū)域邊界，S為加密站位(見圖1)。另外，在調(diào)查時間方面，陸源河流及直排?？诘恼{(diào)查與海上調(diào)查的時間準(zhǔn)同步。在調(diào)查要素方面，河流及直排?？跔I養(yǎng)鹽濃度與流量同步，海上水文、化學(xué)和生物生態(tài)要素同步。

(①Offshore monitoring stations;②River monitoring stations;③Main hydrologic response units;④Yellow River hydrologic response unit;⑤Xiaoqing River hydrologic response unit;⑥Wei River hydrologic response unit;⑦Mi River hydrologic response unit;⑧Jiaolai River hydrologic response unit;⑨Bailang River hydrologic response unit;⑩Yu River hydrologic response unit;Jie River hydrologic response unit;Laizhou waste water treatment facility hydrologic response unit;Longkou waste water treatment facility hydrologic response unit;Main river;Municipal administrative.R站位位于水質(zhì)響應(yīng)區(qū)邊界點，M站位位于混合區(qū)邊界點，H站位位于污染熱點區(qū)域邊界，S為加密站位。Station R were located at the boundary point of the water quality response areas, station M were located at the boundary point of the mixing areas, station H were located at the boundary of the pollution hot spot areas, and station S were encrypted stations.)

陸海同步調(diào)查分別于2018年10月24—25日、2019年5月10—12日、2019年8月18—21日和2021年3月22—27日開展。其中2019年8月份采樣時間為強降雨后7天，受臺風(fēng)“利奇馬”影響，自2019年8月11日開始，環(huán)萊州灣地區(qū)普降大到暴雨，平均降雨量達158 mm，濟南、濰坊、濱州等萊州灣的集水區(qū)出現(xiàn)極端強降雨天氣，最大降雨量達到346 mm，破山東過程降雨量歷史記錄[22]。

1.3 樣品采集與前處理方法

按照《海洋調(diào)查規(guī)范》(GB/T 12763.3—2007)采集水樣，水樣經(jīng)過孔徑為0.45 μm醋酸纖維濾膜過濾后，裝入經(jīng)酸洗的聚乙烯瓶中，-20 ℃下保存，用于營養(yǎng)鹽測定。玻璃濾器在使用前先經(jīng)過鹽酸浸泡24 h，然后用超純水清洗3～6次。按照《地表水和污水監(jiān)測技術(shù)規(guī)范標(biāo)準(zhǔn)(HJ-T92002)》采集環(huán)萊州灣河流和直排?？诘臉悠?，并同步監(jiān)測流量，水樣盛于2 L高密度聚乙烯塑料桶中4 ℃下冷藏保存，6 h內(nèi)運回實驗室中進行過濾和分析。

1.4 樣品測定和分析

河流和排污口流量采用直讀轉(zhuǎn)子流速儀(ZSX-5，中國)測量。海水溫度和鹽度用多參數(shù)儀(Manta 3.0,Eureka，美國)測定。營養(yǎng)鹽通過營養(yǎng)鹽自動分析儀(QuAAtro Applications，德國)測定；NO2-N、NO3-N、NH4-N、PO4-P、SiO3-Si檢測限分別為0.01、0.02、0.03、0.01、0.14 μmol·L-1；DIN為NO3-N、NO2-N和NH4-N之和。葉綠素a(Chla)通過丙酮萃取后應(yīng)用熒光法(F-4500型熒光分光光度計，日立)測定[23]，檢測限為0.01 μg·L-1。

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

應(yīng)用公式(1)計算各項營養(yǎng)鹽的入海通量：

Fj=Cj·Qj·M·10-9·24·31·3 600。

(1)

其中：Fj是在本文中指調(diào)查營養(yǎng)鹽的入海月通量，單位為t·month-1(噸·月-1)；Cj是營養(yǎng)鹽的濃度,單位μmol·L-1；Qj是河流月徑流量或直排?？谖鬯铝髁?，單位為m3·s-1；M為元素的相對原子質(zhì)量。

應(yīng)用SPSS 25.0進行數(shù)據(jù)相關(guān)性分析、主成分分析；相關(guān)性分析采用Pearson相關(guān)系數(shù)(two-tailed,α=0.05)判別，P<0.05時說明變量之間顯著相關(guān)；統(tǒng)計變量之間的差異性檢驗使用Mann-Whitney U-test(two-tailed,α=0.05)，P<0.05說明變量之間存在顯著性差異。

2 結(jié)果與分析

2.1 2018—2021年萊州灣營養(yǎng)鹽入海通量及比值

萊州灣各入海河流和直排?？诹髁咳鐖D2所示，夏季流量高于其它三個季節(jié)，其中，黃河和小清河三個季節(jié)徑流量遠高于萊州灣其他入海河流，春季、夏季、秋季和冬季兩者徑流量合計占比達環(huán)萊州灣河流徑流量的99%、90%、97%和96%。

圖2 2018—2021年度萊州灣入海河流及直排?？诹髁俊DN、TDP和SiO3-Si通量

2018—2021年春季、夏季、秋季和冬季萊州灣入海河流及直排海口TDN的入海通量分別為1.3×104、5.7×104、1.5×104和1.4×104t·month-1，其中DIN的平均占比分別為46%、23%、37%和56%；TDP的入海通量分別3.1×102、7.2×102、1.9×102和1.6×102t·month-1，其中PO4-P的平均占比分別為81%、21%、30%和44%；SiO3-Si的入海通量則分別為1.7×104、4.1×104、1.9×104和0.7×104t·month-1；春、夏、秋和冬季TDN/TDP通量比分別為111、212、199和197，均值為180，TDN/SiO3-Si通量比則分別3、12、5和28，均值為12。季節(jié)分布上，TDP和SiO3-Si入海通量受徑流量影響在夏季最高，冬季最低，TDN通量在夏季最高，但在冬季與春季相差不大。這表明，萊州灣營養(yǎng)鹽入海通量受徑流量的控制更為明顯。各河流及排污口其分擔(dān)率有顯著差異，基本按黃河、小清河、濰河、虞河、廣利河、白浪河、膠萊河、泳汶河、界河的次序降低。

2.2 萊州灣水文分布特征

受陸源淡水輸入影響，萊州灣鹽度分布基本與海岸平行，呈現(xiàn)由灣頂向灣口增大的特征，其中夏季由于河流徑流量增大使得萊州灣內(nèi)鹽度最低，為28.14；春季最高，為29.68。萊州灣水溫以8月份最高，為25.64 ℃；3月最低，為5.93 ℃。秋季、春季、夏季和冬季水體表層鹽度均值較底層分別低2.56、1.29、3.48和1.56，溫度分別相差-0.68、2.44、0.69、1.32 ℃，春季和夏季會存在一定程度的層化現(xiàn)象(見表1)。

2.3 萊州灣營養(yǎng)鹽濃度

2018—2021年春、夏、秋和冬季萊州灣水體DIN分布基本上呈現(xiàn)由灣頂向灣中央和灣口逐步降低的變化趨勢，在小清河和黃河等河口區(qū)出現(xiàn)明顯高值區(qū)(見圖3)。其中，春季和夏季DIN濃度均值相近，分別為(33.73±20.07)和(34.00±42.39) μmol·L-1，均高于國家Ⅲ類但滿足Ⅳ類海水水質(zhì)標(biāo)準(zhǔn)，秋季和冬季較低，為(20.76±53.28)和(16.65±14.70) μmol·L-1，均高于國家Ⅰ類但滿足Ⅱ類海水水質(zhì)標(biāo)準(zhǔn)。NO3-N為DIN主要的賦存形態(tài)，春季、夏季、秋季和冬季平均占比分別為(69±17)%、(86±7)%、(59±16)%和(67±12)%(見表2)。

圖3 2018—2021年不同季節(jié)萊州灣營養(yǎng)鹽濃度分布圖

萊州灣海域不同季節(jié)PO4-P濃度變幅較大，春、夏、秋和冬季分別為0.02～1.81、0.02～4.95、0.01～1.66和0.01～0.51 μmol·L-1，均值分別為(0.13±0.31)、(0.28±0.75)、(0.15±0.28)和(0.07±0.08) μmol·L-1；季節(jié)上，夏季最高，秋季和春季次之，冬季最低，均滿足國家Ⅰ類海水水質(zhì)標(biāo)準(zhǔn)(見表2)；春、夏、秋和冬季PO4-P濃度低于浮游植物生長的絕對閾值(0.1 μmol·L-1)的海域面積占比分別為85.1%、49.5%、40.7%和58.7%，均值為58.5%。

萊州灣水體SiO3-Si濃度基本呈現(xiàn)灣頂河口區(qū)高、灣口低的分布特征。受黃河及灣頂部入海河流的影響[24]，河口區(qū)出現(xiàn)高值。萊州灣SiO3-Si的濃度均值具有明顯的季節(jié)差異，呈現(xiàn)出夏季(20.21±30.29) μmol·L-1最高，秋季(12.14±10) μmol·L-1其次，春季(6.75±7.29) μmol·L-1再次，冬季(3.46±12.30) μmol·L-1最低的變化特征。

2.4 萊州灣營養(yǎng)鹽結(jié)構(gòu)

DIN/PO4-P比值具有較大空間差異，基本呈現(xiàn)近岸河口海域逐漸向灣口降低的趨勢(見圖4)。萊州灣整個海域DIN/PO4-P比值均值呈現(xiàn)春季(744.74)最高，夏季(404.64)其次，冬季(330.06)再次，秋季(235.11)最低的變化趨勢，均遠高于Redfield比。不同季節(jié)萊州灣DIN/SiO3-Si均值差異明顯，其中，冬季最高，達到13.08，春季其次，為7.41；夏季再次，為4.39；秋季最小，為1.60。

圖4 2018—2021年不同季節(jié)萊州灣DIN/PO4-P、DIN/SiO3-Si分布圖

2.5 萊州灣葉綠素

萊州灣表層和底層Chla整體呈現(xiàn)由灣頂向灣口逐步降低的變化趨勢，表層Chla濃度高于底層(見圖5)。Chla的濃度均值具有明顯的季節(jié)差異，呈現(xiàn)出夏季(5.70±5.37) μg·L-1最高，冬季(5.56±5.19) μg·L-1其次，秋季(4.08±3.54) μg·L-1再次，春季(1.35±0.87) μg·L-1最低的變化特征(見表2)。其中，夏季灣南部和東南部及灣中央部分海域Chla含量高于10 μg·L-1，高于赤潮發(fā)生標(biāo)準(zhǔn)，其面積約占海域總面積的60%。

圖5 2018—2021年不同季節(jié)萊州灣海域Chl a分布圖

3 討論

3.1 萊州灣營養(yǎng)鹽時空分布特征及影響因素

萊州灣河流各季節(jié)的TDN、TDP和SiO3-Si通量呈現(xiàn)出一致的規(guī)律，即黃河及萊州灣西南岸的小清河的TDN、TDP和SiO3-Si營養(yǎng)鹽通量占比均高達85%以上，這直接影響了萊州灣海域DIN、PO4-P和SiO3-Si分布規(guī)律，使得濃度高值區(qū)集中在黃河口及小清河河口附近海域(見圖3)，與之前的研究一致[25-28]。值得注意的是夏季小清河TDP通量高于黃河的TDP通量(見圖2)，這是由于2019年8月11日“利奇馬”臺風(fēng)登陸后，濟南濰坊等地出現(xiàn)極端強降雨天氣，導(dǎo)致河流徑流量顯著升高，吸附于泥沙顆粒表面的顆粒磷,經(jīng)過雨水沖刷作用使得河流中總磷含量增加，并能夠向溶解態(tài)磷的形式轉(zhuǎn)變[29]，后經(jīng)河流輸運進入萊州灣內(nèi)，導(dǎo)致萊州灣西南岸小清河口附近PO4-P出現(xiàn)高值區(qū)(見圖3)。夏季萊州灣海域內(nèi)營養(yǎng)鹽的高值區(qū)與鹽度低值區(qū)基本吻合，這是由于監(jiān)測時間在臺風(fēng)“利奇馬”所帶來的強降雨后1周，河流徑流短時期內(nèi)顯著增大，不僅使得萊州灣表層水體鹽度降低，而且徑流攜帶大量陸源營養(yǎng)物質(zhì)輸送到萊州灣內(nèi)，導(dǎo)致海域營養(yǎng)鹽濃度升高[30]。此外，小清河流域除農(nóng)業(yè)耕地外，城市和工業(yè)建設(shè)用地占比也較高(11.4%)[31]，河流中TDN和TDP含量較高，同樣也使得萊州灣西南部近岸海域DIN和PO4-P濃度較高。黃河SiO3-Si通量在各個季節(jié)中占比均為最高，黃河上游巖石類型以硅酸鹽為主，下游巖石類型主要為黏土類硅酸鹽[32]，所以在黃河口附近海域存在一個SiO3-Si含量的高值區(qū)。萊州灣DIN、PO4-P和SiO3-Si濃度與鹽度大都呈顯著負相關(guān)關(guān)系(見圖6)，表明陸源輸入是該海域營養(yǎng)鹽分布的決定性因素之一[33-34]。各入海河流輸入的營養(yǎng)鹽比例不同是導(dǎo)致灣內(nèi)DIN/PO4-P、DIN/SiO3-Si時空異質(zhì)性明顯的直接原因。萊州灣入海河流TDN/TDP和TDN/SiO3-Si通量的比例均值分別180和12，大大偏離Redfield比值，其中白浪河TDN/TDP通量比均值最高，為307.94；其次為萊州污水處理廠，為277.28，這與對應(yīng)的萊州灣海域的DIN/PO4-P高值分布區(qū)基本一致(見圖4)，膠萊河TDN/SiO3-Si通量比均值最高，為35.08；其次為廣利河，為35.11，也與萊州灣海域的DIN/SiO3-Si的高值分布區(qū)基本對應(yīng)(見圖4)，陸源輸入營養(yǎng)鹽的不均衡導(dǎo)致了萊州灣營養(yǎng)鹽結(jié)構(gòu)呈現(xiàn)“氮高、硅少、磷寡”的特征。

圖6 不同季節(jié)萊州灣DIN、PO4-P和SiO3-Si與鹽度相關(guān)性分析圖

水動力條件也是影響萊州灣營養(yǎng)鹽空間分布格局的重要因素之一。萊州灣東部海域刁龍嘴以西存在一個逆時針環(huán)流區(qū)[35]；秋季Chla、DIN，夏季DIN、PO4-P和SiO3-Si，春季DIN濃度較高，與逆時針環(huán)流區(qū)域相吻合，營養(yǎng)鹽在此區(qū)域停留時間較長，凈遷移距離短，難以擴散[35]；此外該區(qū)域圍填?；顒蛹觿∈沟玫撞克w交換不暢，進一步降低了萊州灣的流速及水體交換速率，導(dǎo)致營養(yǎng)鹽在此積累[36-38]。夏季黃河入海徑流增大，科氏力加強，在萊州灣灣底形成窄的浮力沿岸流，部分余流進入萊州灣內(nèi)，從黃河口開始沿著海岸順時針向東岸運移后再出海灣，加之其他河流和排污口輸入量的增大，河口余流導(dǎo)致河口海域底部營養(yǎng)鹽濃度顯著高于其他海域[39]。

萊州灣各季節(jié)DIN與Chla呈顯著正相關(guān)，PO4-P濃度在春季和冬季也與Chla呈正相關(guān)關(guān)系，SiO3-Si除春季外其他三個季節(jié)也與Chla呈現(xiàn)正相關(guān)關(guān)系(見圖7)，夏季表層DIN和SiO3-Si在灣西南部和南部近岸海域的高值區(qū)與夏季表層Chla高值區(qū)吻合(見圖3，5)，表明除陸源輸入這一主要因素外，當(dāng)外界達到適宜的條件(溫度、光照等)，浮游植物初級生產(chǎn)也會對海域營養(yǎng)鹽時空分布具有調(diào)控作用[40]。

圖7 不同季節(jié)萊州灣DIN、PO4-P和SiO3-Si及Chl a相關(guān)性分析圖

3.2 萊州灣浮游植物生長限制性因子

萊州灣PO4-P濃度低于其閾值的海域分布于灣中央和灣東南部，四個季節(jié)其面積占比均高達40%以上，并且DIN/PO4-P和SiO3-Si/PO4-P比值則遠高于22(見表2)，這表明磷為萊州灣浮游植物生長的限制因子，這與姜會超等[36]研究結(jié)果一致。與膠州灣[47]、大亞灣[48]、欽州灣和廉州灣[49]等海域相比，萊州灣磷限制更為嚴重。春季萊州灣內(nèi)磷限制最顯著，這是由于該季節(jié)陸源磷的入海通量顯著低于夏季和秋季，且浮游植物生長旺盛，即按照16∶1的比例吸收DIN和PO4-P，也會加劇磷受限的程度[33,50]；秋季垂直混合加強，營養(yǎng)鹽通過海底沉積物(間隙水)向上覆水釋放，表層營養(yǎng)鹽得到補充[51]，緩解了萊州灣磷限制的狀況；夏季強降雨期間，萊州灣營養(yǎng)物質(zhì)入海通量顯著增大，浮游植物生長的磷限制狀況有所改善，適宜的光照、溫度及營養(yǎng)條件使得該季節(jié)浮游植物快速繁殖增長甚至暴發(fā)赤潮。

自改革開放以來隨著沿海經(jīng)濟高速發(fā)展和人口快速增加，人們生產(chǎn)生活所產(chǎn)生的氮、磷等營養(yǎng)物質(zhì)數(shù)量顯著增大，這些營養(yǎng)物質(zhì)最終入海結(jié)果導(dǎo)致近海富營養(yǎng)化程度加劇。1990—2000年代中期，由于氮肥使用量居高不下且削減總量有限，而生產(chǎn)生活中禁止含磷物資的使用大大減少了磷的排放，特別由于黃河經(jīng)常性斷流，導(dǎo)致泥沙攜帶的磷大大減少，結(jié)果導(dǎo)致2000年后萊州灣入海河流的氮磷通量比值不斷增大[2,54]，萊州灣的營養(yǎng)狀況從20世紀(jì)90年代前的氮限制變?yōu)?1世紀(jì)初的磷限制，2010年達到最高值367.57(見圖8)；隨著2011年國家開始實施以NH4-N為約束性指標(biāo)的陸源污染物總量減排計劃，環(huán)萊州灣地區(qū)氮減排效果顯著，萊州灣入海河流的氮磷通量比值開始明顯下降[19](見圖8)。黃河筑壩工程對各項營養(yǎng)鹽的截留作用有顯著差異，滯留效率硅>磷>氮[55-56]。研究表明，近40年來，黃河氮輸送通量除1990年代后期黃河斷流外，呈現(xiàn)不斷升高的趨勢，而總磷入海通量先下降然后進入2000年代中期后緩慢增長，這加劇了萊州灣氮多磷寡的營養(yǎng)鹽失衡狀況。

圖8 萊州灣1959—2019年DIN/PO4-P變化圖[52-53]

營養(yǎng)鹽結(jié)構(gòu)的不對稱變化會導(dǎo)致浮游植物群落發(fā)生演替，低水平的PO4-P有利于非硅藻的生長[57-58]，渤海近50年間氮磷比不斷升高，導(dǎo)致浮游植物群落結(jié)構(gòu)由硅藻為主逐漸演變?yōu)楣柙搴图自骞餐刂芠59]。Xin等[2]指出渤海1990年代以前浮游植物主要以硅藻為主，Chla含量相對較低，浮游動物種群豐富，底棲動物生物量較高，赤潮罕見；而1990s后硅藻優(yōu)勢度降低；進入21世紀(jì)后，Chla含量增加且浮游植物中非硅藻比例增加，浮游動物數(shù)量減少，大型底棲動物小型化加劇，赤潮災(zāi)害頻發(fā)，赤潮致災(zāi)種呈多樣化和微型化變化趨勢；這些變化進而改變了魚類群落的結(jié)構(gòu)和多樣性，導(dǎo)致漁業(yè)資源下降。

4 結(jié)論

(1)2018—2021年各季節(jié)萊州灣入海河流及直排海口的營養(yǎng)鹽入海通量和營養(yǎng)鹽濃度時空分布具有良好的響應(yīng)關(guān)系。入萊州灣TDN、TDP、SiO3-Si排放通量，整體呈現(xiàn)夏季>秋季>春季>冬季的變化特征，海域DIN、PO4-P、SiO3-Si濃度均值也基本呈現(xiàn)夏季>秋季>春季>冬季的變化特征。從空間分布看，不僅DIN、PO4-P、SiO3-Si濃度高值區(qū)與鹽度低值區(qū)基本吻合，而且營養(yǎng)鹽分布呈現(xiàn)灣西南部近岸海域高于東部海域的特征，這與灣西部和西南部的黃河、小清河入海的營養(yǎng)物質(zhì)排放通量顯著高于灣東部和東南部入海通量直接關(guān)聯(lián)。除陸源輸入外，海域水動力和浮游物植物生長對營養(yǎng)鹽分布的調(diào)控作用也不容忽視。

(2)萊州灣陸源氮磷營養(yǎng)鹽入海通量比例對海域營養(yǎng)鹽結(jié)構(gòu)具有決定性影響。萊州灣各季節(jié)陸源TDN/TDP入海通量比值均值高達180；相應(yīng)的，萊州灣海域各季節(jié)DIN/PO4-P比值均值高達428，都遠高于Redfield比值，表明萊州灣營養(yǎng)鹽結(jié)構(gòu)嚴重失衡；各季節(jié)PO4-P濃度低于浮游植物生長閾值的海域面積占比均值高達58.5%，海域浮游植物生長磷限制嚴重，可能導(dǎo)致萊州灣浮游植物優(yōu)勢種由硅藻為主逐漸演變?yōu)楣杓自骞餐純?yōu)的趨勢,對萊州灣生態(tài)環(huán)境的影響不容忽視。