湖泊富营养化问题在全球范围内日益严重, 已经引发了水华暴发和水环境质量恶化等一系列问题.氮是控制初级生产力的限制性营养物质, 也是造成富营养化的主要原因之一.已有研究表明, 氮在沉积物中的转化包括氨化, 硝化, 厌氧氨氧化, 反硝化和异化NO3--N还原(DNRA)等过程, 这些过程影响了水生生态系统中氮的转化和去除[1, 2].其中沉积物中的反硝化作用可去除82%左右的硝态氮[3].水华的暴发使沉积物中氮的转化过程更为复杂, 并影响着水体中氮的去除[4].随着水华的衰亡, 水华残体沉降到沉积物中可以改变水体与沉积物中营养物质的比例, 同时藻类在降解过程中会消耗大量的溶解氧, 为反硝化过程提供缺氧条件, 从而影响氮转化过程[5, 6].但目前关于水华对湖泊中氮转化和去除影响的认识仍然不足, 这限制了人们对湖泊中氮转化过程的理解以及对湖泊氮含量长期变化及水华暴发的预测能力.
太湖中存在着以水华和水生植物为主的不同的生态类型, 具有高度的异质性, 是研究水华对氮转化和去除影响的良好天然实验场.本研究通过对水华频发和大型水生植物为主的两个湖区的沉积物样品的分析, 分析水华对沉积物中氮含量的影响, 探讨水华对氮转化和去除的影响.为分析水华对沉积物中氮转化的直接和间接影响途径, 采用结构方程模型对水华影响沉积物氮转化和去除的关键因素进行了分析, 以期了解水华对沉积物中氮转化的影响, 并为合理选择氮削减策略提供依据.
1 材料与方法 1.1 采样点和样本采集本研究采样点分别位于西太湖的N1点(31°18′14.0″N, 119°58′17.0″E)和贡湖湾的N3点(31°27′52.5″N, 120°10′36.9″E).自2011年5~12月, 每月采集沉积物样品, 所有样品在4℃条件下立即运送到实验室进行分析.
1.2 沉积物和水样的化学分析100 g沉积物4 000 r·min-1离心30 min, 所得上清为孔隙水, 过0.22 μm滤膜后, 分析孔隙水中总溶解氮(TDN)、铵态氮(NH4+-N)、硝态氮(NO3--N)和亚硝态氮(NO2--N)[7].使用TOC分析仪(Elementar, 德国)测定总有机碳(TOC)含量[8].
1.3 DNA提取和氮转化相关功能基因丰度分析每个样品取0.5 g沉积物用MoBio PowerSoilTM DNA Isolation Kit(San Diego, CA)提取沉淀物中的总DNA.总DNA经电泳和NanoDrop ND 2000(NanoDrop Technologies, 美国)检测后, 通过定量PCR进行细菌16S rRNA(338F/518R)、厌氧氨氧化16S rRNA(AMX809F/1066R[9])、氨氧化细菌(amoA1F/amoA2R -TC)、氨氧化古菌(Arch-amoAF/Arch-amoAR)、反硝化菌(nirS, cd3af/r3cd)、反硝化菌(nirK, nirK-1F/nirK-5R)基因丰度的测定4[4].
1.4 统计分析不同采样位点和时间点之间的差异采用单因素方差分析(ANOVA)和HSD事后检验分析.氮循环基因的丰度与环境变量之间的关系通过SPSS进行Spearman相关分析.水华对沉积物中氮转化和去除之间的直接和间接影响通过结构方程模型(SEM)分析[10].
2 结果与讨论 2.1 沉积物孔隙水中不同形态氮含量的变化两个采样点沉积物中TDN和NH4+-N的浓度均在夏季达到最高, 这与夏季水华的暴发相一致, 表明水华的暴发可能引发了大量的氮输入至沉积物中.与N3点相比, N1点的沉积物中总溶解性氮(TDN)、NH4+-N和NO3--N浓度更高, 而C/N相对较低, 表明与N3相比, N1点的污染更为严重[11]. N1点沉积物中氮浓度较高, 可能是由藻华碎屑沉积造成的.藻华碎屑蛋白质含量高、C/N低, 很容易被沉积物中的微生物分解[12], 可能是沉积物中TDN和NH4+-N在夏季浓度较高的原因之一.同时, 沉积物中的低氧环境可能减缓藻华碎屑的分解速率, 使得N1点沉积物中有机物和氮的含量更高[13].与TDN和NH4+-N相比, NO3--N的浓度呈波动性变化.两个地点的NO2--N浓度在5~7月较高, 而在上覆水中的NO2--N浓度较低(图 1), 这表明沉积物中的NO2--N虽有可能来自上覆水的扩散, 但也可能是由于沉积物中的硝化作用大于反硝化作用[14]所造成的.
![]() |
图 1 沉积物孔隙水与上覆水中TDN、NO3--N、NH4+-N、NO2--N、TOC的浓度和C/N的变化 Fig. 1 Monthly variation in TDN, NO3--N, NH4+-N, NO2--N, and TOC concentration and C/N in the pore water of sediments and overlaying water at the two sampling sites |
图 2结果表明, 不同形态氮转化率在N1与N3点的变化趋势大致相同. 7月和8月的N1点NH4+-N分别增加了396.20%和189.90%, 表明有机氮的矿化在氮转化过程中占主导地位.有机氮矿化作用中释放的NH4+-N易于结合到大生物和微生物的生物质中[15], 在氮的再循环中起着核心作用, 释放的NH4+-N可能会影响包括沉积物硝化和反硝化的整个氮转化过程. TDN在9月去除率最高为71%, 此时C/N为1.24[图 2(a)和图 1(f)], 表明沉积物中在这种情况下, C/N比对反硝化没有显著影响. 7~9月, NH4+-N的转化率呈下降趋势, NO3--N呈上升趋势, 与N1和N3点的潜在硝化速率(PNR)在夏季较高相符(图 3).在整个研究过程中NO2--N积累量较低, 这表明沉积物中的反硝化活性能够去除两个地点的NO3--N和NO2--N.
![]() |
图 2 沉积物中TDN, NH4+-N, NO3--N和NO2--N的氮转化率 Fig. 2 Nitrogen transformation ratio of TDN, NH4+-N, NO3--N, and NO2--N in sediments at the two sampling sites |
![]() |
图 3 两个采样点沉积物PNR的变化 Fig. 3 Monthly variation of PNR in sediments at the two sampling sites |
图 4(a)结果显示, 在采样期间氨氧化细菌amoA(AOB)丰度相对稳定, 氨氧化古菌amoA(AOA)丰度存在显著的季节性变化(P < 0.05), 并在冬季最高.这可能是因为冬季太湖沉积物中氨氮含量相对较低, AOA适合在氨氮含量较低时生长[16]. AOA和AOB丰度在N1点整体上均高于N3(P < 0.05). 5~8月, N1和N3的amoA基因分别增加7.88倍和5.89倍, 氨氧化古菌占据优势[图 4(a)].同时, 在N1和N3点中均检测到了厌氧氨氧化菌, 但两个位点之间的丰度不存在显著性差异, 且其丰度均低于氨氧化古菌.
![]() |
图 4 amoA、厌氧氨氧化16S rRNA和nirK/nirS基因丰度的变化 Fig. 4 Abundance of the amoA gene from ammonia-oxidizing archaea and ammonia-oxidizing bacteria andthat of the 16S rRNA gene of ANAMMOX and nirK gene, and nirS gene at the two sampling sites |
图 4(b)显示nirK和nirS基因的丰度在5~8月增加, 这可能导致沉积物反硝化能力的增加. 9月nirK和nirS的基因丰度降低.与前人的研究结果类似, nirS基因丰度在N1和N3点均高于nirK基因(P < 0.01)[17, 18], 这也表明nirS型的反硝化菌是太湖沉积物中的主要生态类型.
2.4 水华对氮转化和去除的影响结构方程模型分析结果表明水华对沉积物氮循环的直接和间接影响对总溶解性氮去除的解释度为65.1%(图 5).水华对沉积物中氮去除没有直接影响, 但水华对沉积物中的TOC(λ=0.34; P < 0.001)、TDN(λ=0.73; P < 0.001)、厌氧氨氧菌16S rDNA丰度(λ=0.53; P < 0.001)和nirS/nirK基因丰度(λ=0.49; P < 0.001)呈显著的正相关.
![]() |
图中仅显示了显著的效果;实线箭头表示正效应, 虚线箭头表示负效应;箭头的宽度表示因果关系的强度;箭头上的数字是标准化路径系数;内生变量的百分比(R2)表示模型解释的方差 图 5 水华对沉积物中氮转化和去除影响的结构方程模型分析 Fig. 5 Structural equation model for the effects of water bloom on nitrogen transformation and removal in sediment |
沉积物中TDN的含量对氮的矿化作用表现出强烈的正效应(λ=0.48; P < 0.05).厌氧氨氧化菌(λ=0.31; P < 0.001)和amoA(λ=-0.28; P < 0.001)均对氨氮浓度有很大影响, 氨氮浓度则对nirS/nirK基因丰度有很大影响(λ=0.75; P < 0.001). nirS/nirK基因丰度, NH4+-N去除率和上覆水中NO3--N浓度对脱氮效率有直接影响, 分别为0.423、-0.674和-0.447(图 5).水华可以直接提高TDN和TOC的浓度, 厌氧氨氧化菌和nirS/nirK的丰度, 并间接提高沉积物氨氮和硝态氮浓度, 从而影响氮的去除.首先, 水华生物量可以通过直接促进细胞的存活和维持, 同时减少细胞生长的滞后时间来加速分解过程[19, 20].其次, 在水华暴发和沉积期间, 与有机物质降解有关的功能基因, 微生物酶活性和微生物都受到刺激[21, 22], 导致有机质矿化能力增强, 氮利用率增高[23].
amoA基因丰度对厌氧氨氧化菌的丰度具有显著正效应(λ=0.37; P < 0.001), 表明NH4+-N氧化所产生的NO2--N可能是激发沉积物中厌氧氨氧化潜能的重要因素(图 5), 这表明了在沉积物中硝化与厌氧氨氧化作用耦合脱氮的可能性.这两个途径的耦合已经在污水处理厂[24]、人工湿地[25]以及海洋[26]和淡水生态系统[27]中的脱氮过程中都发挥着重要作用, 但需要更多地实证研究来证实太湖沉积物中这种耦合作用对沉积物中氮的转化合脱除的作用.
此外, 结构方程模型分析显示nirS/nirK丰度是唯一与TDN去除正相关的因素, nirK/nirS丰度可以解释沉积物中总溶解性氮去除的42.3%, 表明反硝化是太湖沉积物脱氮的主要过程(图 5).反硝化和厌氧氨氧化都消耗了通过amoA硝化或通过narG/napA反硝化产生的NO2--N[28].因此, 厌氧氨氧化菌和反硝化细菌之间存在资源利用竞争.然而, 在本研究中, nirS/nirK和厌氧氨氧化菌基因之间并没有显著的负相关关系, 表明在太湖沉积物中厌氧氨氧化和反硝化之间的竞争是可以忽略的.这可能是由于nirK/nirS的丰度(7.91×106~2.43×108copies·g-1)远高于厌氧氨氧化菌的基因丰度(4.44×104~4.99×105copies·g-1)造成的(图 4).水华与厌氧氨氧化基因丰度(λ=0.53, P < 0.001)和nirK /nirS基因丰度(λ=0.49, P < 0.001)的显著正相关表明, 水华可通过增强厌氧氨氧化和反硝化过程加速沉积物中氮的去除.
3 结论N1和N3点的沉积物中TDN和NH4+-N的浓度在夏季达到最高, 并且N1点沉积物中氮浓度高于N3点, 表明水华的暴发与沉积造成沉积物中氮的输入.结构方程分析结果表明水华可直接对沉积物中的TDN和TOC浓度产生正向影响, 从而间接影响沉积物中NH4+-N浓度.在太湖沉积物中硝化作用通过与厌氧氨氧化耦合促进了氮去除过程.然而, 与厌氧氨氧化相比, 通过nirK/nirS进行的反硝化作用是太湖沉积物中氮去除的主要途径.
[1] | Rios-Del Toro E E, Valenzuela E I, Ramirez J E, et al. Anaerobic ammonium oxidation linked to microbial reduction of natural organic matter in marine sediments[J]. Environmental Science & Technology Letters, 2018, 5(9): 571-577. |
[2] | Tan E H, Hsu T C, Huang X F, et al. Nitrogen transformations and removal efficiency enhancement of a constructed wetland in subtropical Taiwan[J]. Science of the Total Environment, 2017, 601-602: 1378-1388. DOI:10.1016/j.scitotenv.2017.05.282 |
[3] | Kreiling R M, Richardson W B, Cavanaugh J C, et al. Summer nitrate uptake and denitrification in an upper Mississippi River backwater lake:the role of rooted aquatic vegetation[J]. Biogeochemistry, 2011, 104(1-3): 309-324. DOI:10.1007/s10533-010-9503-9 |
[4] | Peng Y K, Liu L, Jiang L J, et al. The roles of cyanobacterial bloom in nitrogen removal[J]. Science of the Total Environment, 2017, 609: 297-303. DOI:10.1016/j.scitotenv.2017.03.149 |
[5] | Ma J R, Qin B Q, Paerl H W, et al. The persistence of cyanobacterial (Microcystis spp.) blooms throughout winter in Lake Taihu, China[J]. Limnology and Oceanography, 2016, 61(2): 711-722. DOI:10.1002/lno.10246 |
[6] | Ding S M, Chen M S, Gong M D, et al. Internal phosphorus loading from sediments causes seasonal nitrogen limitation for harmful algal blooms[J]. Science of the Total Environment, 2018, 625: 872-884. DOI:10.1016/j.scitotenv.2017.12.348 |
[7] | 国家环境保护总局. 水和废水监测分析方法[M]. 第四版. 北京: 中国环境科学出版社, 2002. |
[8] | Li H L, Chen J J, Wu W, et al. Distribution of polycyclic aromatic hydrocarbons in different size fractions of soil from a coke oven plant and its relationship to organic carbon content[J]. Journal of Hazardous Materials, 2010, 176(1-3): 729-734. DOI:10.1016/j.jhazmat.2009.11.095 |
[9] | Tsushima I, Kindaichi T, Okabe S. Quantification of anaerobic ammonium-oxidizing bacteria in enrichment cultures by real-time PCR[J]. Water Research, 2007, 41(4): 785-794. DOI:10.1016/j.watres.2006.11.024 |
[10] | Xiong Z Q, Guo L D, Zhang Q F, et al. Edaphic conditions regulate denitrification directly and indirectly by altering denitrifier abundance in wetlands along the Han River, China[J]. Environmental Science & Technology, 2017, 51(10): 5483-5491. |
[11] | Wu Z, Liu Y, Liang Z Y, et al. Internal cycling, not external loading, decides the nutrient limitation in eutrophic lake:a dynamic model with temporal Bayesian hierarchical inference[J]. Water Research, 2017, 116: 231-240. DOI:10.1016/j.watres.2017.03.039 |
[12] | Kitidis V, Tait K, Nunes J, et al. Seasonal benthic nitrogen cycling in a temperate shelf sea:the Celtic Sea[J]. Biogeochemistry, 2017, 135(1-2): 103-119. DOI:10.1007/s10533-017-0311-3 |
[13] | Kinsman-Costello L E, Sheik C S, Sheldon N D, et al. Groundwater shapes sediment biogeochemistry and microbial diversity in a submerged Great Lake sinkhole[J]. Geobiology, 2017, 15(2): 225-239. DOI:10.1111/gbi.12215 |
[14] | Yang J H, Zhang X W, Xie Y W, et al. Ecogenomics of zooplankton community reveals ecological threshold of ammonia nitrogen[J]. Environmental Science & Technology, 2017, 51(5): 3057-3064. |
[15] | Gold A C, Thompson S P, Piehler M F. Nitrogen cycling processes within stormwater control measures:a review and call for research[J]. Water Research, 2019, 149: 578-587. DOI:10.1016/j.watres.2018.10.036 |
[16] | Verhamme D T, Prosser J I, Nicol G W. Ammonia concentration determines differential growth of ammonia-oxidising archaea and bacteria in soil microcosms[J]. The ISME Journal, 2011, 5(6): 1067-1071. DOI:10.1038/ismej.2010.191 |
[17] | Yin C, Fan F L, Song A L, et al. Denitrification potential under different fertilization regimes is closely coupled with changes in the denitrifying community in a black soil[J]. Applied Microbiology and Biotechnology, 2015, 99(13): 5719-5729. DOI:10.1007/s00253-015-6461-0 |
[18] | Han H Y, Li Z K. Effects of macrophyte-associated nitrogen cycling bacteria on ANAMMOX and denitrification in river sediments in the Taihu Lake region of China[J]. Ecological Engineering, 2016, 93: 82-90. DOI:10.1016/j.ecoleng.2016.05.015 |
[19] | Sun R, Sun P F, Zhang J H, et al. Microorganisms-based methods for harmful algal blooms control:a review[J]. Bioresource Technology, 2018, 248: 12-20. DOI:10.1016/j.biortech.2017.07.175 |
[20] | Hutchins D A, Fu F X. Microorganisms and ocean global change[J]. Nature Microbiology, 2017, 2(6): 17058. DOI:10.1038/nmicrobiol.2017.58 |
[21] | Bennett J A, Maherali H, Reinhart K O, et al. Plant-soil feedbacks and mycorrhizal type influence temperate forest population dynamics[J]. Science, 2017, 355(6321): 181-184. DOI:10.1126/science.aai8212 |
[22] | Kyaschenko J, Clemmensen K E, Hagenbo A, et al. Shift in fungal communities and associated enzyme activities along an age gradient of managed Pinus sylvestris stands[J]. The ISME Journal, 2017, 11(4): 863-874. DOI:10.1038/ismej.2016.184 |
[23] | Alvarez G, Shahzad T, Andanson L, et al. Catalytic power of enzymes decreases with temperature:new insights for understanding soil C cycling and microbial ecology under warming[J]. Global Change Biology, 2018, 24(9): 4238-4250. DOI:10.1111/gcb.14281 |
[24] | Du R, Cao S B, Li B K, et al. Performance and microbial community analysis of a novel DEAMOX based on partial-denitrification and anammox treating ammonia and nitrate wastewaters[J]. Water Research, 2017, 108: 46-56. DOI:10.1016/j.watres.2016.10.051 |
[25] | Li H Y, Tao W D. Efficient ammonia removal in recirculating vertical flow constructed wetlands:complementary roles of anammox and denitrification in simultaneous nitritation, anammox and denitrification process[J]. Chemical Engineering Journal, 2017, 317: 972-979. DOI:10.1016/j.cej.2017.02.143 |
[26] | Hellemann D, Tallberg P, Bartl I, et al. Denitrification in an oligotrophic estuary:a delayed sink for riverine nitrate[J]. Marine Ecology Progress Series, 2017, 583: 63-80. DOI:10.3354/meps12359 |
[27] | Xiao K, Wu J P, Li H L, et al. Nitrogen fate in a subtropical mangrove swamp:Potential association with seawater-groundwater exchange[J]. Science of the Total Environment, 2018, 635: 586-597. DOI:10.1016/j.scitotenv.2018.04.143 |
[28] | Ma B, Qian W T, Yuan C S, et al. Achieving mainstream nitrogen removal through coupling anammox with denitratation[J]. Environmental Science & Technology, 2017, 51(15): 8405-8413. |