朱 偉,尤星宇,段志鵬①,呂 藝,侯 豪
(1.河海大學(xué)環(huán)境學(xué)院,江蘇 南京 210098;2.河海大學(xué)土木與交通學(xué)院,江蘇 南京 210098)
太湖是我國第三大淡水湖,屬超大型淺水湖泊,其水質(zhì)環(huán)境在空間上存在顯著差異。例如,太湖東南部區(qū)域的水質(zhì)要明顯優(yōu)于西北部和湖心區(qū)域[1],總體呈“東清西濁”特征[2]。以胥口灣為代表的太湖東部是草型清水區(qū),可提供優(yōu)質(zhì)的飲用水,對維系區(qū)域水資源安全具有重要意義[3]。監(jiān)測顯示,近年來太湖北部水質(zhì)環(huán)境發(fā)生了顯著變化[1],可能影響總體水質(zhì)格局。然而,胥口灣缺乏水質(zhì)跟蹤監(jiān)測,導(dǎo)致難以全面掌握太湖水質(zhì)總體變化趨勢,不利于水資源管理。
經(jīng)典湖泊穩(wěn)態(tài)理論表明,處于穩(wěn)態(tài)的淺水湖泊生態(tài)系統(tǒng)能夠強(qiáng)化該區(qū)域的水質(zhì)特征[4]。例如,草型清水區(qū)的水生植物既能吸收營養(yǎng)物質(zhì)[5-6],又能有效削減風(fēng)浪對底泥的擾動以減少底泥再懸浮,降低水體營養(yǎng)鹽濃度[3]。此外,水生植物還能分泌化感物質(zhì)直接抑制藻類生長[7-8],共同維系草型清水狀態(tài)。然而,極端天氣如持續(xù)高溫或由其引起的流域洪水可能導(dǎo)致草型系統(tǒng)發(fā)生災(zāi)變。研究表明,持續(xù)高溫可能對水生植物造成高溫脅迫,抑制水生植物生長和繁殖[9]。同時(shí),高溫還能增加底泥微生物活性,驅(qū)動底泥有機(jī)質(zhì)加速分解和氮磷等營養(yǎng)物質(zhì)快速釋放[10]。流域洪水還可能帶來上游含高濃度營養(yǎng)鹽物質(zhì)的水團(tuán),直接導(dǎo)致草型湖區(qū)的水質(zhì)指標(biāo)惡化[11]。流域洪水導(dǎo)致的水位快速上漲也會嚴(yán)重削弱水下光場,可顯著抑制水生植物的光合作用,進(jìn)而導(dǎo)致水生植物衰亡,誘導(dǎo)草型淺水湖泊生態(tài)系統(tǒng)的災(zāi)變發(fā)生[12]。因此,明確太湖典型草型清水湖區(qū)水質(zhì)環(huán)境變化趨勢及對水溫、水位的響應(yīng),對區(qū)域水資源管理具有重要意義。
近年來,太湖流域發(fā)生多次大洪水和持續(xù)高溫事件。2016年夏季發(fā)生特大洪水[13-14],隨后2020年夏季又發(fā)生大洪水[15-16],2022年夏季太湖流域乃至整個(gè)長江流域發(fā)生罕見的持續(xù)高溫和干旱事件[17],尤其是8月整個(gè)太湖流域出現(xiàn)大范圍、長時(shí)間的高溫干旱天氣,使得太湖夏季水位相比歷史同期下降約45%[18]。極端天氣及水溫、水位變化可能顯著影響太湖水質(zhì)格局[19],威脅區(qū)域生態(tài)環(huán)境安全。
研究選取太湖典型草型清水湖區(qū)——胥口灣,分析自2016年以來的水質(zhì)變化規(guī)律,補(bǔ)充胥口灣近年來較長時(shí)間序列的水質(zhì)監(jiān)測數(shù)據(jù),有利于深化對太湖水環(huán)境總體格局的認(rèn)識。同時(shí),探究胥口灣水質(zhì)環(huán)境與水溫、水位的變化關(guān)系,力圖探明水質(zhì)環(huán)境對流域內(nèi)極端天氣及其引起的流域洪水的響應(yīng)特征,以期為草型淺水湖泊生態(tài)系統(tǒng)對水溫、水位的響應(yīng)研究提供理論參考,為水資源管理提供借鑒。
研究在太湖胥口灣設(shè)有2個(gè)環(huán)境監(jiān)測站點(diǎn)(圖1)。其中,樣點(diǎn)1處于胥口灣與湖心區(qū)和北部灣區(qū)的交界區(qū)域,屬于草型區(qū)向藻型區(qū)過渡區(qū)域,水質(zhì)特征受上游來水影響可能更加顯著[2]。樣點(diǎn)2處于胥口灣內(nèi)部中心區(qū),是典型的草型區(qū),水動力和生態(tài)環(huán)境均較為穩(wěn)定,屬于飲用水源地保護(hù)區(qū)。
圖1 太湖胥口灣水質(zhì)環(huán)境監(jiān)測站點(diǎn)
監(jiān)測時(shí)間為2016年1月1日至2022年12月31日,每天采用HACH-GLI在線pH/ORP分析儀(sc200P53P33PRO-P3,美國)實(shí)時(shí)監(jiān)測水溫和pH值,并通過邁德施在線葉綠素檢測儀(MDS-B500CPL,中國)估算葉綠素a(Chl-a)濃度。同時(shí),在表層水下0.5 m處采集水樣,收集于500 mL塑料樣品瓶中,并在4 ℃下避光快速轉(zhuǎn)運(yùn)至實(shí)驗(yàn)室,24 h內(nèi)測定氮磷濃度。需要注意的是,由于儀器等原因,樣點(diǎn)1缺乏2022年8月以后的水質(zhì)指標(biāo)數(shù)據(jù),樣點(diǎn)2缺乏2016年水溫的日均數(shù)據(jù)(僅存月均數(shù)據(jù))、部分營養(yǎng)鹽和Chl-a濃度數(shù)據(jù)。
水樣氮磷濃度測定包含總氮(TN)、氨氮(NH3-N)和總磷(TP)指標(biāo),同時(shí)測定水樣的高錳酸鹽指數(shù)(CODMn)以表征水體中的有機(jī)物質(zhì)含量。上述水質(zhì)指標(biāo)的測定方法均參照文獻(xiàn)[20]中的標(biāo)準(zhǔn)方法進(jìn)行。每個(gè)水質(zhì)指標(biāo)的月均數(shù)據(jù)是將該指標(biāo)的每日數(shù)據(jù)累計(jì)后取平均值。
太湖水位數(shù)據(jù)來源于水利部太湖流域管理局(http:∥www.tba.gov.cn/),包含2016年1月1日至2022年12月31日太湖水位的2 551個(gè)日均數(shù)據(jù)點(diǎn)。月均水位數(shù)據(jù)是將該月的每日水位數(shù)據(jù)累計(jì)后取平均值。
將水溫觸及30 ℃且后續(xù)水溫連續(xù)3 d高于該溫度時(shí)定義為高水溫時(shí)段或高溫氣候事件,這是因?yàn)樗疁爻^30 ℃通常會顯著改變藻類和水生植物的生理生化活性及物種組成[9]。太湖水位受流域降水和人工調(diào)蓄共同作用,流域極端降雨導(dǎo)致的洪水對湖泊生態(tài)系統(tǒng)可能具有顯著影響。然而,流域降雨對湖泊生態(tài)系統(tǒng)的影響通常具有滯后性,水位能更直接地表征流域極端降雨對湖泊系統(tǒng)的影響。基于此,將太湖的設(shè)計(jì)洪水警戒水位3.8 m設(shè)為高水位,當(dāng)太湖水位觸及3.8 m且后續(xù)水位連續(xù)3 d高于該水位時(shí)定義為高水位時(shí)段或洪峰時(shí)段。
所有統(tǒng)計(jì)分析均在R語言(R version 4.1.0)中進(jìn)行。通過R語言自帶的aov函數(shù)進(jìn)行單因素方差分析,對比分析高水溫、高水位時(shí)段胥口灣水質(zhì)環(huán)境與非高水溫、非高水位水質(zhì)環(huán)境指標(biāo)日均數(shù)據(jù)的差異顯著性,顯著性水平設(shè)為P<0.05。此外,采用R語言自帶的prcomp函數(shù)對胥口灣水質(zhì)環(huán)境指標(biāo)月均數(shù)據(jù)的季節(jié)性變化和空間差異性進(jìn)行主成分分析(PCA),并采用ggplot2程序包中的autoplot函數(shù)繪制主成分分析數(shù)據(jù)圖。同時(shí)在R語言中運(yùn)行vegan程序包中的rda函數(shù),通過冗余分析(RDA)和分步冗余分析(pRDA)定量分析月均水溫和水位對胥口灣水質(zhì)環(huán)境指標(biāo)月均數(shù)據(jù)波動的貢獻(xiàn)率,并采用permutest函數(shù)對自變量貢獻(xiàn)值的顯著性進(jìn)行置換檢驗(yàn)(999次置換),顯著性水平設(shè)為P<0.05。
2016—2022年太湖胥口灣2個(gè)樣點(diǎn)處的TN濃度、CODMn和pH值水平和變化趨勢均較為接近(圖2)。
圖2 2016—2022年胥口灣TN、TP、NH3-N、Chl-a、CODMn和pH值的變化趨勢
其中,2016—2018年TN質(zhì)量濃度在1.5 mg·L-1左右,TN質(zhì)量濃度于2018年春季達(dá)到峰值(1.8~2.3 mg·L-1)后快速下降,2019—2022年處于0.6~1.5 mg·L-1之間(均值約1.0 mg·L-1)。而CODMn在2016—2018年有快速下降趨勢,隨后變化趨于平緩,在3.0 mg·L-1左右小幅波動。近7 a來胥口灣的pH值穩(wěn)定在8.2左右,但樣點(diǎn)2處的pH值浮動更為劇烈。
相比之下,TP、NH3-N和Chl-a濃度在樣點(diǎn)1和樣點(diǎn)2之間有顯著差異。其中,樣點(diǎn)1處的TP和Chl-a濃度要顯著高于樣點(diǎn)2,而NH3-N濃度在2個(gè)樣點(diǎn)之間的差異相反(圖2)。樣點(diǎn)1處TP質(zhì)量濃度在2016—2022年總體上表現(xiàn)出下行趨勢,但在2016年夏季和2022年春季分別出現(xiàn)低谷值(0.04 mg·L-1)和高峰值(0.09 mg·L-1);而樣點(diǎn)2處的TP質(zhì)量濃度在0.03 mg·L-1左右浮動。雖然樣點(diǎn)1處Chl-a濃度維持在較低水平(4.3 μg·L-1),但常在冬春季出現(xiàn)峰值;而樣點(diǎn)2處的Chl-a質(zhì)量濃度常年低于3 μg·L-1。樣點(diǎn)1處的NH3-N質(zhì)量濃度常年低于0.1 mg·L-1,且波動幅度小;而樣點(diǎn)2處的NH3-N濃度要明顯高于樣點(diǎn)1,且在2017—2019年出現(xiàn)極大波動,最大值達(dá)0.6 mg·L-1。
近年來太湖胥口灣水溫呈快速上升趨勢,2016—2022年水溫年均值上升近1 ℃(圖3)。夏季是太湖流域溫度最高的季節(jié),就月均值而言,2016—2022年中只有2021年全年水溫月均值均低于30 ℃。2017和2022年7、8月水溫月均值均高于30 ℃,分別為31.6、30.7 ℃以及31.2、31.6 ℃。就日均值而言,2017和2022年高水溫時(shí)段天數(shù)分別為44和54 d,遠(yuǎn)高出2016—2022年高水溫時(shí)段天數(shù)均值(35.1 d)。
圖3 2016—2022年胥口灣水溫及水位的月均值、年均值變化
近年來太湖水位波動幅度也較大,2016年是特大洪水年,水位年均值接近3.6 m。2020年也是洪水年,水位年均值接近3.4 m。月均值上,2016年6、7、10月水位均超過警戒水位(3.8 m)。日均值上,2016年高水位(>3.8 m)時(shí)段天數(shù)為89 d,2020年高水位時(shí)段天數(shù)為47 d,其余年份的高水位時(shí)段天數(shù)均小于5 d。
通過PCA分析發(fā)現(xiàn),太湖胥口灣水質(zhì)環(huán)境存在一定的季節(jié)差異和顯著的空間差異(圖4)。主成分第一軸(PC1)對太湖胥口灣水質(zhì)數(shù)據(jù)變異性的解釋率為32.9%,第二軸(PC2)的解釋率為19.2%。其中,造成胥口灣水質(zhì)環(huán)境季節(jié)差異的主要因素為TN、水溫和水位,其他因素?zé)o顯著貢獻(xiàn)??臻g差異性分析顯示,樣點(diǎn)1和樣點(diǎn)2存在顯著差異,且樣點(diǎn)1和樣點(diǎn)2在第一軸上差異化分布的主要貢獻(xiàn)水質(zhì)指標(biāo)為NH3-N、Chl-a、TP濃度和CODMn。相反,TN濃度、水位和水溫的空間差異性較低。
圖4 太湖胥口灣水質(zhì)環(huán)境季節(jié)差異和空間差異的主成分分析
通過RDA分析發(fā)現(xiàn),水位和水溫變化對胥口灣水質(zhì)環(huán)境的影響非常有限,總解釋率不足5%(圖5)。其中,水溫顯著貢獻(xiàn)RDA1(R2=0.11,P=0.001),水位顯著貢獻(xiàn)RDA2(R2=0.08,P=0.003)。分步冗余分析(pRDA)分析顯示,水溫對水質(zhì)環(huán)境有顯著影響(P=0.026),但水位的影響不顯著(P=0.397)。
圖5 水溫和水位對胥口灣水質(zhì)指標(biāo)影響的冗余分析
通過方差分析高水溫期間(>30 ℃)與非高水溫期間水質(zhì)指標(biāo)之間的差異,發(fā)現(xiàn)高水溫對胥口灣水質(zhì)環(huán)境有顯著影響(圖6)。具體而言,高水溫期的TN、NH3-N濃度和CODMn分別比非高水溫期降低8.8%(P<0.001)、12.9%(P<0.05)和4.6%(P<0.001),而pH值上升0.65%(P=0.003)。然而,高水溫對胥口灣TP和Chl-a濃度均未產(chǎn)生顯著影響(P>0.05)。
圖6 高水溫(>30 ℃)對胥口灣水質(zhì)指標(biāo)的影響
方差分析顯示,流域洪水或高水位(>3.8 m)同樣顯著影響胥口灣水質(zhì)環(huán)境(圖7)。洪水期間,NH3-N濃度下降35.6%(P<0.001),TN濃度降低9.8%(P=0.004),TP濃度下降6.1%(P=0.02),pH值下降1.3%(P<0.001),而Chl-a濃度上升65.5%(P<0.001)。相較而言,流域大洪水對胥口灣CODMn未產(chǎn)生顯著影響(P=0.08)。
圖7 流域高水位(>3.8 m)對胥口灣水質(zhì)指標(biāo)的影響
監(jiān)測顯示,2016—2022年胥口灣水質(zhì)指標(biāo)總體穩(wěn)定,表現(xiàn)出改善趨勢(圖2)。這既展現(xiàn)出草型湖區(qū)生態(tài)系統(tǒng)的穩(wěn)態(tài)性[4],也受益于近年來環(huán)太湖大力度的水環(huán)境治理措施。通過流域綜合治理,環(huán)太湖入湖河道氮磷等污染物輸入量已經(jīng)快速降低,水質(zhì)指標(biāo)持續(xù)改善[21]。其中, 2020—2021年TN和TP平均濃度相較于2007年分別降低47.5%和18.0%[22]。然而,胥口灣水質(zhì)分布存在顯著的空間異質(zhì)性(圖4)。相比于樣點(diǎn)2(灣區(qū)內(nèi)部),樣點(diǎn)1(西部外沿區(qū)域)的TP和Chl-a濃度較高且波動較大。其中,TP濃度月均值在2022年2月達(dá)近年來的峰值(約0.1 mg·L-1),而Chl-a濃度也常在冬春季出現(xiàn)峰值。盡管如此,TP與Chl-a濃度之間沒有顯著相關(guān)性(R2=0.004,P>0.05),表明樣點(diǎn)1的冬春季藻類峰值并非TP濃度波動引起,亦或是TP濃度波動與藻類繁殖關(guān)系不大。由于樣點(diǎn)1臨近湖心區(qū)和北部灣區(qū)(圖1),因而受上游來水的影響可能較為顯著。樣點(diǎn)1位于北部區(qū)域水團(tuán)進(jìn)入胥口灣的主要通道[2],上游水團(tuán)進(jìn)入胥口灣后有較長的停留時(shí)間[2],這可能解釋了樣點(diǎn)1的TP和Chl-a濃度較高、變幅較大。NH3-N濃度表現(xiàn)出相反趨勢,樣點(diǎn)2處的NH3-N濃度是樣點(diǎn)1處的2~3倍,這可能與區(qū)域水生植物生長過程相關(guān)。樣點(diǎn)2處的水草茂盛[23],但水生植物衰敗過程中可能會釋放NH3-N[24]。這表明適度的水草管理措施(分區(qū)收割等)可能有利于維持胥口灣NH3-N濃度穩(wěn)定。
此外,胥口灣水質(zhì)指標(biāo)的季節(jié)性變化不明顯。除水溫和TN濃度以外,其余水質(zhì)指標(biāo)均不能有效分辨出季節(jié)性變化規(guī)律。胥口灣TN濃度通常冬春季較高,而夏秋季出現(xiàn)顯著下降,這種變化趨勢可能主要由微生物的反硝化作用引起。研究表明,反硝化是太湖水體脫氮的主要過程,且隨著水溫的上升,水體微生物反硝化作用逐漸加強(qiáng)[25]。這可能是TN濃度冬春高、夏秋低的原因之一。在草型湖區(qū),水生植物的衰敗過程也可能顯著影響NH3-N濃度。但筆者研究表明,胥口灣內(nèi)部區(qū)域多數(shù)情況下NH3-N濃度占TN濃度的比值低于20%,這表明NH3-N濃度波動并沒有顯著影響TN濃度的變化趨勢。
水溫是湖泊生態(tài)系統(tǒng)中的重要物理參數(shù),直接影響水生動植物的新陳代謝和微生物對有機(jī)質(zhì)的分解等湖泊物質(zhì)能量交換速率[26-27]。因而,水溫變化可能會對湖泊生態(tài)系統(tǒng)造成復(fù)雜的影響[28]。筆者研究顯示,2016—2022年太湖胥口灣水溫年均值出現(xiàn)快速上升趨勢(圖3,R2=0.89,P<0.01)。需要指明的是,由于缺乏日均數(shù)據(jù),該研究中2016年水溫年均值通過月均數(shù)據(jù)計(jì)算獲得。月均水溫?cái)?shù)據(jù)可能存在采樣日期代表性差的問題,難以準(zhǔn)確計(jì)算水溫年均值,這可能是水溫年均值在2016年出現(xiàn)最低值(18.48 ℃)的原因。然而,即便剔除2016年的數(shù)據(jù),2017—2022年的水溫年均值同樣表現(xiàn)出顯著上升趨勢(R2=0.87,P<0.01),上升幅度達(dá)0.65 ℃,遠(yuǎn)高于全球大部分湖泊表面水溫的上升速率〔0.34 ℃·(10 a)-1〕[28]。彭宗琪[29]發(fā)現(xiàn),2000—2018年我國主要大型湖泊(巢湖、滇池、太湖、洪澤湖、鄱陽湖和洞庭湖)的全年平均水溫上升速率約為0.3 ℃·(10 a)-1,且夏季上升更加顯著。這表明太湖胥口灣水溫年均值可能進(jìn)入了快速上升期。
湖泊升溫可能會加劇湖泊富營養(yǎng)化進(jìn)程,改變水質(zhì)環(huán)境健康狀態(tài)[30-31]。RDA分析發(fā)現(xiàn),雖然水溫變化對胥口灣水質(zhì)指標(biāo)波動的貢獻(xiàn)較低,但水溫顯著影響了胥口灣的水質(zhì)環(huán)境(圖5)。pRDA分析顯示,高水溫(>30 ℃)期間TN、NH3-N濃度和CODMn顯著下降(圖6)。水生植物的適宜生長溫度通常低于25 ℃。研究表明,高溫會顯著抑制水生植物的生長[9],影響其對氮磷等元素的吸收。胥口灣的優(yōu)勢水生植物種類為苦草(Vallisneriaspiralis)、金魚藻(Ceratophyllumdemersum)和狐尾藻(Myriophyllumspicatum)[32]。其中,苦草和金魚藻最適生長溫度均可達(dá)30 ℃,且能適應(yīng)更廣的溫度范圍[33-34]。高溫期間胥口灣水生植物對氮的吸收利用可能會增加,導(dǎo)致TN和NH3-N濃度降低。此外,高水溫會增加微生物活性,加快有機(jī)物質(zhì)分解及反硝化,提高脫氮速率[35],有利于降低水體中氮和有機(jī)質(zhì)含量。
水位也是影響湖泊生態(tài)系統(tǒng)的重要物理參數(shù)。研究顯示,草型淺水湖泊生態(tài)系統(tǒng)對水位變化極為敏感[36]。因?yàn)楦咚粫@著削弱水下光強(qiáng),抑制水生植物生長或誘導(dǎo)水生植物大面積衰亡,最終可能引發(fā)系統(tǒng)性災(zāi)變[37]。然而,筆者研究發(fā)現(xiàn),水位不能有效解釋胥口灣水質(zhì)環(huán)境的年際變化趨勢(圖5)。雖然每年太湖水位會在夏秋季出現(xiàn)峰值,但水位波動并沒有引起胥口灣水質(zhì)指標(biāo)的規(guī)律性波動(圖3)。這可能是由于胥口灣多數(shù)年份的水位均較低,且波動幅度小(圖3),未引起水生植物的顯著變化。盡管如此,流域洪水或高水位顯著影響了胥口灣水質(zhì)環(huán)境(圖7)。其中, NH3-N和Chl-a濃度受高水位的影響最為顯著,分別下降35%和上升65%。同時(shí),上游來水還可能攜帶太湖北部灣區(qū)的藍(lán)藻顆粒進(jìn)入胥口灣,這或許是Chl-a濃度顯著上升的原因之一。
此外,高水位期間胥口灣TP濃度出現(xiàn)顯著下降趨勢(圖7)。尤其是2016年7月太湖流域特大洪水期間,胥口灣西部外沿區(qū)域的TP濃度出現(xiàn)極低值(0.004 mg·L-1)。前人研究表明,2016年的特大洪水導(dǎo)致太湖TP入湖通量顯著高于平水年[38],但洪水期間高磷濃度水團(tuán)由太湖西北部向胥口灣等東南部遷移過程中,水體磷濃度出現(xiàn)了顯著降低趨勢[11]。這是因?yàn)榻?jīng)由河道等陸源輸入太湖的磷以顆粒態(tài)磷為主,外源輸入的顆粒態(tài)磷進(jìn)入湖泊后會快速沉降進(jìn)入底泥[2],這使得洪水期間高磷濃度水體不能到達(dá)胥口灣等太湖東南部區(qū)域[11],相反可能對胥口灣內(nèi)部的TP進(jìn)行稀釋,從而導(dǎo)致水體TP濃度出現(xiàn)下降趨勢。綜上所述,雖然近7 a水位變化對胥口灣水質(zhì)環(huán)境年際變化的總體趨勢無顯著影響,但流域洪水或高水位造成了NH3-N、Chl-a和TP濃度等關(guān)鍵水質(zhì)指標(biāo)的顯著波動。
胥口灣是太湖水質(zhì)環(huán)境最優(yōu)質(zhì)的區(qū)域之一,雖然2016—2022年間太湖流域發(fā)生了持續(xù)高溫的極端天氣或特大洪水事件,但胥口灣的水質(zhì)環(huán)境總體穩(wěn)定,氮磷濃度波動向好的趨勢不變。然而,胥口灣水質(zhì)存在顯著的空間異質(zhì)性。其中,胥口灣西部外沿區(qū)域的水質(zhì)要明顯劣于水生植被更加繁茂的內(nèi)部區(qū)域,但內(nèi)部區(qū)域的NH3-N濃度要顯著高于西部邊沿區(qū)。因此,NH3-N濃度應(yīng)成為胥口灣內(nèi)部的重點(diǎn)監(jiān)測指標(biāo)。此外,雖然水溫和水位對胥口灣水質(zhì)環(huán)境波動的貢獻(xiàn)較低,但高水位和高水溫對多數(shù)水質(zhì)指標(biāo)有顯著影響,表明流域洪水和持續(xù)高溫等極端天氣可能成為破壞胥口灣水環(huán)境穩(wěn)定的重要不確定性因素。該研究闡明了近年來胥口灣水質(zhì)特征,識別出對流域洪水和持續(xù)高溫高度敏感的關(guān)鍵水質(zhì)指標(biāo),可為胥口灣水資源風(fēng)險(xiǎn)評估和高效管理提供理論參考和借鑒。