环境因子对雁鸣湖沉积物氮磷释放的影响

裴佳瑶, 冯民权. 环境因子对雁鸣湖沉积物氮磷释放的影响[J]. 环境工程学报, 2020, 14(12): 3447-3459. doi: 10.12030/j.cjee.201912021
引用本文: 裴佳瑶, 冯民权. 环境因子对雁鸣湖沉积物氮磷释放的影响[J]. 环境工程学报, 2020, 14(12): 3447-3459. doi: 10.12030/j.cjee.201912021
PEI Jiayao, FENG Minquan. Effects of environmental factors on the release of nitrogen and phosphorus from the sediment of the Yanming Lake, China[J]. Chinese Journal of Environmental Engineering, 2020, 14(12): 3447-3459. doi: 10.12030/j.cjee.201912021
Citation: PEI Jiayao, FENG Minquan. Effects of environmental factors on the release of nitrogen and phosphorus from the sediment of the Yanming Lake, China[J]. Chinese Journal of Environmental Engineering, 2020, 14(12): 3447-3459. doi: 10.12030/j.cjee.201912021

环境因子对雁鸣湖沉积物氮磷释放的影响

    作者简介: 裴佳瑶(1994—),女,硕士研究生。研究方向:底泥营养盐释放。E-mail:447896358@qq.com
    通讯作者: 冯民权(1964—),男,博士,教授。研究方向:环境模拟与污染物治理。E-mail:mqfeng@xaut.edu.cn
  • 基金项目:
    国家自然科学基金资助项目(51679191);陕西省自然科学基础研究计划重点项目(2019JZ-42)
  • 中图分类号: X502

Effects of environmental factors on the release of nitrogen and phosphorus from the sediment of the Yanming Lake, China

    Corresponding author: FENG Minquan, mqfeng@xaut.edu.cn
  • 摘要: 为探明环境因子对雁鸣湖沉积物-上覆水界面间氮磷释放的影响,通过单因子实验和响应面实验分析探究了温度、溶解氧及pH对氮磷释放通量的影响。单因子实验结果表明:雁鸣湖沉积物-上覆水界面间氮磷通量随着温度的升高而增加;随溶解氧浓度的升高而减小;当pH为中性时,氮磷通量最小,氮通量在pH=5时达到最大值,磷通量在pH=9时达到最大值。BBD模型拟合结果显示,各环境因子与氮磷通量间拟合关系均为二次多项式。响应面实验结果表明:温度和pH交互作用对TN通量影响显著;温度和溶解氧交互作用对TP通量影响显著;温度和溶解氧交互作用、温度和pH交互作用对${{\rm{NH}}_4^{+}}$-N通量影响均为显著;温度和溶解氧交互作用对${{\rm{PO}}_4^{3 - }}$-P通量影响极显著;雁鸣湖底泥释放的最不利条件为T=10 ℃、DO=7.87 mg·L−1、pH=7.13,此时的TN交换通量为3.956 mg·(m2·h)−1,TP交换通量为0.471 mg·(m2·h)−1${{\rm{NH}}_4^{+}}$-N交换通量为1.469 mg·(m2·h)−1${{\rm{PO}}_4^{3 - }}$-P交换通量为0.146 mg·(m2·h)−1。以上研究结果可为雁鸣湖富营养化防治提供参考。
  • 加载中
  • 图 1  采样点位置示意图

    Figure 1.  Schematic diagram of sampling point location

    图 2  培养装置示意图

    Figure 2.  Schematic of the culture device

    图 3  温度对营养盐通量的影响

    Figure 3.  Effect of temperature on nutrient flux

    图 4  溶解氧对营养盐通量的影响

    Figure 4.  Effect of dissolved oxygen on nutrient flux

    图 5  pH对营养盐通量的影响

    Figure 5.  Effect of pH on nutrient flux

    图 6  营养盐通量实际值与预测值的对比

    Figure 6.  Comparison of actual and predicted nutrient fluxes

    图 7  X1X3交互作用对TN通量的影响

    Figure 7.  Effect of X1X3 interaction on TN flux

    图 8  X1X2交互作用对TP通量的影响

    Figure 8.  Effect of X1X2 interaction on TP flux

    图 9  X1X2X1X3交互作用对${\rm{NH}}_4^{+}$-N通量的影响

    Figure 9.  Effect of X1X2 and X1X3 interactions on ${\rm{NH}}_4^{+}$-N flux

    图 10  X1X2交互作用对${\rm{PO}}_4^{3-}$-P通量的影响

    Figure 10.  Effect of X1X2 interaction on ${\rm{PO}}_4^{3-}$-P flux

    表 1  响应面实验的因子与水平编码

    Table 1.  Factors and level codes for response surface tests

    编码环境因子水平
    −101
    X1温度/℃517.530
    X2溶解氧/(mg·L−1)369
    X3pH579
    编码环境因子水平
    −101
    X1温度/℃517.530
    X2溶解氧/(mg·L−1)369
    X3pH579
    下载: 导出CSV

    表 2  BBD实验实测值

    Table 2.  Experimental results of BBD

    实验组号TN通量/(mg·(m2·h)−1)TP通量/(mg·(m2·h)−1)${\rm{NH}}_4^{+}$-N通量/(mg·(m2·h)−1)${\rm{PO}}_4^{3-}$-P通量/(mg·(m2·h)−1)
    平行样1平行样2平行样3平行样1平行样2平行样3平行样1平行样2平行样3平行样1平行样2平行样3
    18.4457.7548.5590.7560.7640.6733.3753.1152.9740.2360.2340.217
    25.0125.6895.5410.6230.6450.6512.0451.9972.2890.1860.2150.195
    37.4157.5257.8940.5740.7010.6812.8522.6822.9620.2030.2110.225
    47.5517.0457.5410.6760.6780.6362.9422.9642.6750.1950.2210.215
    54.6045.1974.1410.4680.4510.4311.5891.5761.7860.1370.1510.139
    65.2124.6344.7910.4730.4750.5041.8411.8461.6190.1490.1420.155
    78.1118.6978.1190.6640.7270.6713.0473.2023.1050.2290.2210.197
    85.7895.8656.4360.6120.6270.5892.4122.2762.2480.1970.1890.190
    94.7895.3275.2140.5230.5540.4751.8252.2411.9650.1610.1690.151
    105.5896.5496.9890.6780.6970.6792.7582.3882.5870.2160.1970.214
    114.7074.7844.3140.4670.5040.4411.7721.5901.8100.1420.1450.162
    126.4176.1245.5850.6500.6170.5632.3572.4122.0570.1780.1920.196
    134.3214.0124.3380.4630.5250.5181.6691.5591.4960.1660.1450.154
    148.8077.9818.9130.7110.7890.7653.3213.1973.1240.2240.2450.232
    153.9743.9543.5150.5240.5030.5311.4161.4821.3290.1710.1600.146
    167.2587.5697.8410.6070.5470.6432.7912.7603.0600.1960.1740.200
    177.3057.3947.2110.6780.6970.6112.8042.6492.7170.1960.2210.213
    实验组号TN通量/(mg·(m2·h)−1)TP通量/(mg·(m2·h)−1)${\rm{NH}}_4^{+}$-N通量/(mg·(m2·h)−1)${\rm{PO}}_4^{3-}$-P通量/(mg·(m2·h)−1)
    平行样1平行样2平行样3平行样1平行样2平行样3平行样1平行样2平行样3平行样1平行样2平行样3
    18.4457.7548.5590.7560.7640.6733.3753.1152.9740.2360.2340.217
    25.0125.6895.5410.6230.6450.6512.0451.9972.2890.1860.2150.195
    37.4157.5257.8940.5740.7010.6812.8522.6822.9620.2030.2110.225
    47.5517.0457.5410.6760.6780.6362.9422.9642.6750.1950.2210.215
    54.6045.1974.1410.4680.4510.4311.5891.5761.7860.1370.1510.139
    65.2124.6344.7910.4730.4750.5041.8411.8461.6190.1490.1420.155
    78.1118.6978.1190.6640.7270.6713.0473.2023.1050.2290.2210.197
    85.7895.8656.4360.6120.6270.5892.4122.2762.2480.1970.1890.190
    94.7895.3275.2140.5230.5540.4751.8252.2411.9650.1610.1690.151
    105.5896.5496.9890.6780.6970.6792.7582.3882.5870.2160.1970.214
    114.7074.7844.3140.4670.5040.4411.7721.5901.8100.1420.1450.162
    126.4176.1245.5850.6500.6170.5632.3572.4122.0570.1780.1920.196
    134.3214.0124.3380.4630.5250.5181.6691.5591.4960.1660.1450.154
    148.8077.9818.9130.7110.7890.7653.3213.1973.1240.2240.2450.232
    153.9743.9543.5150.5240.5030.5311.4161.4821.3290.1710.1600.146
    167.2587.5697.8410.6070.5470.6432.7912.7603.0600.1960.1740.200
    177.3057.3947.2110.6780.6970.6112.8042.6492.7170.1960.2210.213
    下载: 导出CSV

    表 3  BBD实验设计及通量平均值

    Table 3.  BBD experiment design and average flux

    实验组号环境因子TN通量/
    (mg·(m2·h)−1)
    TP通量/
    (mg·(m2·h)−1)
    ${\rm{NH}}_4^{+}$-N通量/
    (mg·(m2·h)−1)
    ${\rm{PO}}_4^{3-}$-P通量/
    (mg·(m2·h)−1)
    温度/℃溶解氧/(mg·L−1)pH
    117.5358.2530.7313.1550.229
    217.5955.4140.6402.1100.199
    330377.6110.6912.8320.213
    430657.3790.6632.8600.210
    517.5674.3730.4501.6500.142
    617.5674.8790.4841.7690.149
    730698.3090.6873.1180.216
    817.5996.0300.6092.3120.192
    917.5675.1100.5171.8950.160
    1030976.7690.6852.5780.209
    1117.5674.6020.4711.7240.150
    125696.0420.6102.2750.189
    135974.2240.5021.5750.155
    1417.5398.5670.7553.2140.234
    1517.5673.8140.5191.4090.159
    165657.5560.5992.8700.190
    175377.3030.6622.7230.210
    实验组号环境因子TN通量/
    (mg·(m2·h)−1)
    TP通量/
    (mg·(m2·h)−1)
    ${\rm{NH}}_4^{+}$-N通量/
    (mg·(m2·h)−1)
    ${\rm{PO}}_4^{3-}$-P通量/
    (mg·(m2·h)−1)
    温度/℃溶解氧/(mg·L−1)pH
    117.5358.2530.7313.1550.229
    217.5955.4140.6402.1100.199
    330377.6110.6912.8320.213
    430657.3790.6632.8600.210
    517.5674.3730.4501.6500.142
    617.5674.8790.4841.7690.149
    730698.3090.6873.1180.216
    817.5996.0300.6092.3120.192
    917.5675.1100.5171.8950.160
    1030976.7690.6852.5780.209
    1117.5674.6020.4711.7240.150
    125696.0420.6102.2750.189
    135974.2240.5021.5750.155
    1417.5398.5670.7553.2140.234
    1517.5673.8140.5191.4090.159
    165657.5560.5992.8700.190
    175377.3030.6622.7230.210
    下载: 导出CSV

    表 4  底泥释放通量回归方程方差分析

    Table 4.  Analysis of variance for regression equation of the release flux from sediment

    方差来源TN通量回归分析TP通量回归分析${\rm{NH}}_4^{+}$-N通量回归分析${\rm{PO}}_4^{3-}$-P通量回归分析
    dfFP>F显著性dfFP>F显著性dfFP>F显著性dfFP>F显著性
    模型918.450.000 4**922.260.000 2**920.020.000 3**941.79< 0.000 1**
    X1113.190.008 4**121.680.002 3**114.950.006 2**135.610.000 6**
    X2146.650.000 2**128.420.001 1**144.410.000 3**156.50.000 1**
    X310.0650.806 710.140.717 110.0220.885 210.030.868 2
    X1X215.40.053 118.330.023 4*16.320.040 2*117.130.004 4**
    X1X316.450.038 7*10.0540.822 515.760.047 5*10.320.587 8
    X2X310.0980.762 911.020.345 710.160.699 810.950.362 6
    ${X}_1^{2}$121.530.002 4**115.410.005 7**122.460.002 1**129.290.001**
    ${X}_2^{2}$112.610.009 3**153.20.000 2**114.260.006 9**190.08< 0.000 1**
    ${X}_3^{2}$151.160.000 2**158.570.000 1**161.80.000 1**1120.78< 0.000 1**
    残差7777
    失拟项30.830.54130.550.672 930.940.500 130.230.868 4
    纯误差4444
    总值16161616
    变异系数0.0770.044 40.075 40.032 7
    R20.959 50.966 20.962 60.981 7
    $R_{{\rm{Adj}}}^2$0.907 50.922 80.914 50.958 2
      注:**表示P <0.01下差异显著;*表示P<0.05下差异显著;df表示自由度。
    方差来源TN通量回归分析TP通量回归分析${\rm{NH}}_4^{+}$-N通量回归分析${\rm{PO}}_4^{3-}$-P通量回归分析
    dfFP>F显著性dfFP>F显著性dfFP>F显著性dfFP>F显著性
    模型918.450.000 4**922.260.000 2**920.020.000 3**941.79< 0.000 1**
    X1113.190.008 4**121.680.002 3**114.950.006 2**135.610.000 6**
    X2146.650.000 2**128.420.001 1**144.410.000 3**156.50.000 1**
    X310.0650.806 710.140.717 110.0220.885 210.030.868 2
    X1X215.40.053 118.330.023 4*16.320.040 2*117.130.004 4**
    X1X316.450.038 7*10.0540.822 515.760.047 5*10.320.587 8
    X2X310.0980.762 911.020.345 710.160.699 810.950.362 6
    ${X}_1^{2}$121.530.002 4**115.410.005 7**122.460.002 1**129.290.001**
    ${X}_2^{2}$112.610.009 3**153.20.000 2**114.260.006 9**190.08< 0.000 1**
    ${X}_3^{2}$151.160.000 2**158.570.000 1**161.80.000 1**1120.78< 0.000 1**
    残差7777
    失拟项30.830.54130.550.672 930.940.500 130.230.868 4
    纯误差4444
    总值16161616
    变异系数0.0770.044 40.075 40.032 7
    R20.959 50.966 20.962 60.981 7
    $R_{{\rm{Adj}}}^2$0.907 50.922 80.914 50.958 2
      注:**表示P <0.01下差异显著;*表示P<0.05下差异显著;df表示自由度。
    下载: 导出CSV
  • [1] 中华人民共和国生态环境部. 2018中国生态环境状况公报[EB/OL]. [2019-12-01]. http://www.mee.gov.cn/hjzl/sthjzk/zghjzkgb.
    [2] RUSSEL J M, HOPMANS E C, LOOMIS S E, et al. Distributions of 5- and 6-methyl branched glycerol dialkyl glycerol tetraethers (brGDGTs) in East African lake sediment: Effects of temperature, pH, and new lacustrine paleotemperature calibrations[J]. Organic Geochemistry, 2018, 117: 56-69. doi: 10.1016/j.orggeochem.2017.12.003
    [3] BEUTEL M W, HORNE A J. Nutrient fluxes from profundal sediment of ultra-oligotrophic lake tahoe, california/nevada: Implications for water quality and management in a changing climate[J]. Water Resources Research, 2018, 54(3): 1549-1559. doi: 10.1002/2017WR020907
    [4] 韩宝红, 宋蕾, 李浩, 等. 不同温度条件下稳定剂对沉积物中镉稳定化的影响[J]. 环境科学学报, 2019, 39(8): 2610-2616.
    [5] 李鑫, 耿雪, 王洪伟, 等. 外源输入对底泥疏浚新生表层磷恢复及迁移的影响[J]. 环境科学, 2019, 40(8): 3539-3579.
    [6] OKOGBUE C O, OYESANYA O U, ANYIAM O A, et al. Evaluation of the extent of pollution of discharged oil field brine in the Bonny estuary, Niger Delta, Nigeria[J]. Environmental Earth Sciences, 2018, 77(10): 382-396. doi: 10.1007/s12665-018-7542-z
    [7] MAIERYEMU Y, MAMAT S, NIGELA T, et al. Distribution of heavy metal pollution and assessment of its potential ecological risks in Ugan-Kuqa River Delta of Xinjiang[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(20): 226-233.
    [8] 郭赟, 赵秀红, 黄晓峰, 等. 原位活性覆盖抑制河道底泥营养盐释放的效果研究及工程化应用[J]. 环境工程, 2018, 36(6): 6-11.
    [9] 雷晓玲, 韩亚鑫, 冉兵, 等. 环境因子对三峡库区底泥污染物释放的影响研究[J]. 环境工程, 2016, 34(1): 47-50.
    [10] KAISER D, UNGER D, QIU G, et al. Natural and human influences on nutrient transport through a small subtropical Chinese estuary[J]. Science of the Total Environment, 2013, 450-451: 92-107. doi: 10.1016/j.scitotenv.2013.01.096
    [11] LIANG Z, LIU Z, ZHEN S, et al. Phosphorus speciation and effects of environmental factors on release of phosphorus from sediments obtained from Taihu Lake, Tien Lake, and East Lake[J]. Toxicological & Environmental Chemistry Reviews, 2015, 97(3/4): 335-348.
    [12] LI R F, FENG CH, WANG D X. Role of salinity in the multiphase redistribution of polycyclic aromatic hydrocarbons (PAHs) in sediment suspension[J]. Environmental Earth Sciences, 2016, 75(2): 116-122. doi: 10.1007/s12665-015-5014-2
    [13] 张文斌, 董昭皆, 徐书童, 等. 微生物和藻类分解对荣成天鹅湖沉积物氮磷释放的影响[J]. 海洋环境科学, 2019, 38(4): 561-567. doi: 10.12111/j.mes20190412
    [14] FERENCZ B, TOPOROWSKA M, JAROSLAW D, et al. Hydro-chemical conditions of shaping the water quality of shallow leczna-wlodawa lakes[J]. Clean-Soil Air Water, 2017, 45(5): 356-369.
    [15] DEVORE C L, RODRIGUEZ F L, MEHDI A, et al. Effect of bicarbonate and phosphate on arsenic release from mining-impacted sediments in the Cheyenne River watershed, South Dakota, USA[J]. Environmental Science: Processes and Impacts, 2019, 21(3): 456-468. doi: 10.1039/C8EM00461G
    [16] YAO C, HU X Z, LU S Y, et al. Repression of nitrogen and phosphorus release from lakeshore sediment by five littoral-zone plants[J]. Environmental Science, 2017, 38(2): 589-599.
    [17] 张强, 曹秀芹, 胡明, 等. 扰动对城市河道底泥污染物释放影响[J]. 环境工程, 2019, 37(9): 40-44.
    [18] 张茜. 漳泽水库沉积物和上覆水污染特征及氮磷释放规律研究[D]. 西安: 西安理工大学, 2019.
    [19] 张硕, 方鑫, 黄宏, 等. 基于正交试验的沉积物-水界面营养盐交换通量研究: 以海州湾海洋牧场为例[J]. 中国环境科学, 2017, 37(11): 4266-4276. doi: 10.3969/j.issn.1000-6923.2017.11.032
    [20] 姜伟, 周川, 纪道斌, 等. 三峡库区澎溪河与磨刀溪电导率等水质特征与水华的关系比较[J]. 环境科学, 2017, 38(6): 2326-2335.
    [21] 国家环境保护总局. 水和废水监测分析方法[M]. 4版. 北京: 中国环境科学出版社, 2002.
    [22] 王圣瑞. 湖泊沉积物-水界面过程: 基本理论与常用测定方法[M]. 北京: 科学出版社, 2014.
    [23] JEROEN J M, CISKA C O, CHRISTIAN J J, et al. Effect of temperature on oxygen profiles and denitrification rates in freshwater sediments[J]. Wetlands, 2017, 37(5): 975-983. doi: 10.1007/s13157-017-0933-1
    [24] 张红, 陈敬安, 王敬富, 等. 贵州红枫湖底泥磷释放的模拟实验研究[J]. 地球与环境, 2015, 43(2): 243-251.
    [25] BAREHA Y, GIRAULT R, JIMENEZ J, et al. Characterization and prediction of organic nitrogen biodegradability during anaerobic digestion: A bioaccessibility approach[J]. Bioresource Technology, 2018, 263: 425-436. doi: 10.1016/j.biortech.2018.04.085
    [26] 郭念, 闫金龙, 魏世强, 等. 三峡库区消落带典型土壤厌氧呼吸对铁还原及磷释放的影响[J]. 水土保持学报, 2014, 28(3): 271-276.
    [27] LI Z, WANG S R, WU Z H. Coupling effect of pH and dissolved oxygen in water column on nitrogen release at water-sediment interface of Erhai Lake, China[J]. Estuarine Coastal and Shelf Science, 2014, 149: 178-186. doi: 10.1016/j.ecss.2014.08.009
    [28] 袁和忠, 沈吉, 刘恩峰, 等. 模拟水体pH控制条件下太湖梅梁湾沉积物中磷的释放特征[J]. 湖泊科学, 2015, 21(5): 663-668.
    [29] 王荣杰, 沈本贤, 刘纪昌, 等. 响应面分析法优化不溶性硫黄萃取提浓工艺[J]. 应用化工, 2018, 47(7): 1457-1461. doi: 10.3969/j.issn.1671-3206.2018.07.035
    [30] 闫晓涛, 李杰, 冯淑琪, 等. 响应面法优化混凝处理黄河兰州段低温低浊水[J]. 水资源与水工程学报, 2018, 29(6): 68-74.
    [31] 杜凤龄, 王刚, 徐敏, 等. 新型高分子螯合-絮凝剂制备条件的响应面法优化[J]. 中国环境科学, 2015, 35(4): 158-164.
  • 加载中
图( 10) 表( 4)
计量
  • 文章访问数:  4222
  • HTML全文浏览数:  4222
  • PDF下载数:  65
  • 施引文献:  0
出版历程
  • 收稿日期:  2019-12-03
  • 录用日期:  2020-03-01
  • 刊出日期:  2020-12-10

环境因子对雁鸣湖沉积物氮磷释放的影响

    通讯作者: 冯民权(1964—),男,博士,教授。研究方向:环境模拟与污染物治理。E-mail:mqfeng@xaut.edu.cn
    作者简介: 裴佳瑶(1994—),女,硕士研究生。研究方向:底泥营养盐释放。E-mail:447896358@qq.com
  • 西安理工大学,省部共建西北旱区生态水利国家重点实验室,西安 710048
基金项目:
国家自然科学基金资助项目(51679191);陕西省自然科学基础研究计划重点项目(2019JZ-42)

摘要: 为探明环境因子对雁鸣湖沉积物-上覆水界面间氮磷释放的影响,通过单因子实验和响应面实验分析探究了温度、溶解氧及pH对氮磷释放通量的影响。单因子实验结果表明:雁鸣湖沉积物-上覆水界面间氮磷通量随着温度的升高而增加;随溶解氧浓度的升高而减小;当pH为中性时,氮磷通量最小,氮通量在pH=5时达到最大值,磷通量在pH=9时达到最大值。BBD模型拟合结果显示,各环境因子与氮磷通量间拟合关系均为二次多项式。响应面实验结果表明:温度和pH交互作用对TN通量影响显著;温度和溶解氧交互作用对TP通量影响显著;温度和溶解氧交互作用、温度和pH交互作用对${{\rm{NH}}_4^{+}}$-N通量影响均为显著;温度和溶解氧交互作用对${{\rm{PO}}_4^{3 - }}$-P通量影响极显著;雁鸣湖底泥释放的最不利条件为T=10 ℃、DO=7.87 mg·L−1、pH=7.13,此时的TN交换通量为3.956 mg·(m2·h)−1,TP交换通量为0.471 mg·(m2·h)−1${{\rm{NH}}_4^{+}}$-N交换通量为1.469 mg·(m2·h)−1${{\rm{PO}}_4^{3 - }}$-P交换通量为0.146 mg·(m2·h)−1。以上研究结果可为雁鸣湖富营养化防治提供参考。

English Abstract

  • 湖泊是一种生物与环境、生物与生物之间相互依存与制约的复杂系统。一般情况下,这个系统是处于相对平衡的状态,过于频繁的人类活动会打破原本的平衡状态。2018年中华人民共和国生态环境部发布的关于湖库富营养化现状调查结果[1]显示,在所有被调查的湖库中,中营养占比为61.7%,中营养以上占比为29%,仅有9.3%营养化程度为贫营养。底泥是湖泊水体重要的营养源,在湖泊自净能力范围内,沉积物表现为营养物的“源”,超过其承载力时则表现为“汇”。由于浅水湖泊具有环境比较封闭、流动性差、水深较浅、易污染等特点,所以外界环境极易对湖泊水环境造成干扰,使得沉积物中的污染物释放到水体中形成二次污染。

    氮磷过量输入是导致湖泊富营养化的直接原因,进入湖泊水体的氮磷来源于2个方面:一是外源的直接输入;另一重要来源则是底泥向上覆水的内源释放。近年来,国内外学者在对非洲的东非湖[2]、美国的太浩湖[3]、中国的乌梁素海[4]和太湖梅梁湾[5]等研究中,也证实了底泥污染物质对水质有较大的影响,因此,对因底泥营养盐释放引起水体二次污染的问题进行了极大的关注,并展开了一系列的相关研究。这些研究主要包括沉积物中污染物含量及分布[6]、污染程度评价[7]、底泥的释放机理[8]等。掌握污染物在沉积物-上覆水间的迁移对内源污染防治至关重要,迁移作用主要包括扩散、吸附-解吸、矿化、溶解和分解等。在这些过程中,pH[9]、溶解氧[10]、温度[11]、盐度[12]、微生物或藻类[13]、扰动(包括物理和生物扰动)[14]等因素会影响沉积物-上覆水间物质的交换。目前,关于此方面的研究多是从单因素角度定性地对环境因子和污染物释放间的关系进行分析[15-16]。张强等[17]在研究中发现,高扰动强度增加了TN和COD的释放通量,缩短了TP的吸附-解吸平衡时间,而对氨氮的影响较小。实际上各因子是协同作用的,也有部分学者通过设计正交实验,揭示了多因子共同作用下底泥的释放规律。张茜[18]通过正交实验得出温度、pH、溶解氧对总氮和总磷释放速率有着显著的差异,沉积物中总氮和总磷的释放与环境因子间的显著性排序为pH>温度>溶解氧。张硕等[19]研究发现,温度和溶解氧交互作用对正磷酸盐交换通量影响显著,对氨氮和硝氮交换通量无显著影响。但正交实验也存在一定的局限性,其分析的是离散型数据,所求得的因素最佳组合只能局限在所设计的工况中,但在现实情况中,各因子的变化是连续的,离散型数据会降低所得结果的精度。

    雁鸣湖2#湖(图1)位于西安市,属于浅水人工湖,从浐河引水,属于河道外湖泊。湖西面为黄土源,雨季时雨水夹带着塬上的污染物会流入2#湖,加上受污染的浐河水体会直接进入2#湖,这些因素都会导致湖泊水环境受到污染。在进行本研究前,监测发现,2#湖底泥和水体已经受到一定程度的污染,底泥颜色呈现深褐色至黑色,并散发有难闻气味,水体中TN、TP已经远超过国际上规定的富营养化发生的阈值[20](TN为0.2 mg·L−1,TP为0.025 mg·L−1)。针对雁鸣湖内源污染的研究很少,然而污染现状已不容忽视。本研究以雁鸣湖为研究对象,在单因子实验的基础上,结合响应面法,分别考察了温度、溶解氧、pH对沉积物-上覆水界面氮磷释放规律的影响,此实验方法规避了正交实验工况离散、回归精度较低等不足,分析了单因子作用及交互作用对界面氮磷释放的影响,并通过建立回归模型量化反映环境因子与释放量间的对应关系,研究结果对雁鸣湖富营养化防治具有一定的参考价值。

    • 于2019年秋季进行样品采集,实验所用底泥采自西安市雁鸣湖2#湖中心,具体位置见图1。底泥通过抓斗式采样器获得,上覆水通过有机玻璃采样器获得。将采集的上覆水装入水桶,底泥装进自封袋迅速带回实验室备用。

      在实验室将采回来的泥样平铺于有机玻璃柱底部,将上覆水缓缓加入底泥上部进行实验,培养装置如图2所示。选取的环境因子有温度、溶解氧、pH(X1X2X3)。上覆水温度的调节通过恒温培养箱控制;溶解氧的调节通过向水中充氮气或氧气实现;pH的调节通过向原水中加盐酸或氢氧化钠来调控。

      在单因子实验中,环境因子水平值的设置参考湖泊历史实测值。每个因子设置5个水平值:温度分别为5、10、17.5、25、30 ℃;溶解氧分别为3、4.5、6、7.5、9 mg·L−1;pH分别为 5、6、7、8、9,共15个培养柱。每隔12 h取一次上覆水,检测其中TN、TP、${\rm{NH}}_4^{+}$-N、${\rm{PO}}_4^{3-}$-P的含量,每个指标测定3个平行样,整个实验周期共持续7 d。每次取完水样补充相同体积的上覆水。

      响应面实验根据Box-Behnken模型设计要求设有17个培养柱,为估计实验误差,其中有5组中心点重复实验。所测营养盐种类同单因子实验,每个指标测定3个平行样。因子水平数的设置根据单因子的设置及Box-Behnken模型设计要求决定,最终设置的变量水平及对应编码如表1所示。

    • 上覆水中TN、TP、${\rm{NH}}_4^{+}$-N、${\rm{PO}}_4^{3-}$-P,参照文献中的方法[21]测定,温度、pH、溶解氧用便携式多参数测定仪测定。

      本研究采用实验室培养法进行交换通量的估算,具体计算方法[22]如式(1)所示。

      式中:F为交换通量,mg·(m2·h)−1F>0表示污染物从沉积物向上覆水扩散,否则,从上覆水向沉积物扩散;A为沉积物与水界面营养盐交换面积,m2t为培养时间,h;Mtt时间段内营养盐的质量变化量,mg,计算方法如式(2)所示。

      式中:V为培养柱中上覆水的总体积,L;Ctt时刻测得的上覆水中营养盐浓度,mg·L−1Dt1t−1时刻上覆水中实际的营养盐浓度,mg·L−1,计算方法如式(3)所示。

      式中:V0为每次所取的上覆水体积,L;C0Ct−1分别为原始时刻和t−1时刻上覆水中营养盐的浓度,mg·L−1

    • 图3为氮磷通量随温度变化的趋势。氨氮和总氮、正磷酸盐和总磷表现出一致性的变化规律,通量均随温度增加呈现明显的上升趋势。30 ℃下的氨氮和总氮通量分别为5 ℃时的1.8倍和2.2倍。一方面是由于温度较高时微生物活性有所增强,直接加速有机氮的矿化作用[23];另一方面是由于微生物代谢和运动强度的增加,会消耗环境中氧气,使沉积物-上覆水界面处于厌氧环境,从而增强反硝化作用,进而加速界面间氮扩散。由图3可知,温度为5 ℃时的磷通量最小,正磷酸盐和总磷通量仅为0.166 mg·(m2·h)−1和0.315 mg·(m2·h)−1,随着温度逐渐升高,磷通量也逐渐增大,在30 ℃时正磷酸盐、总磷通量分别增加至0.243 mg·(m2·h)−1和0.566 mg·(m2·h)−1。这是因为升温会增加体系中离子的活性,促进离子交换反应的进行,加快磷酸盐的溶解及扩散作用;同时由于生物活动增强导致界面间耗氧增强,厌氧条件加速了Fe3+还原成Fe2+、Mn4+还原成Mn2+的速率[24],使铁锰结合态磷得到释放。

      在不同溶解氧浓度的条件下,底泥的氮磷释放规律呈现一定的差异性。图4为氮磷通量随溶解氧浓度的变化趋势。氨氮通量为2.102~3.109 mg·(m2·h)−1,总氮通量为3.342~5.942 mg·(m2·h)−1,即溶解氧从最高水平(9 mg·L−1)减小到最低水平(3 mg·L−1)时,氨氮和总氮通量分别增加了1.5倍和1.8倍,在厌氧条件下的氮通量显著有所增加,低溶解氧水平下会促进底泥中有机氮的矿化[25],同时反硝化作用剧烈,无机氮主要以氨氮的形式释放,使得氮通量有明显的增加。对于磷的释放,在高溶解氧(9 mg·L−1)水平下,正磷酸盐通量为实验范围内最小,为0.183 mg·(m2·h)−1;而在厌氧条件下(3 mg·L−1),其通量达到0.333 mg·(m2·h)−1,为好氧条件下的1.8倍,总磷通量也由好氧条件下的0.456 mg·(m2·h)−1增加到0.566 mg·(m2·h)−1。这是因为厌氧条件促使Fe3+还原成Fe2+,使得原本被氢氧化铁胶体吸附的磷释放出来;而在高溶解氧条件下,体系处于氧化状态,Fe2+会被氧化成Fe3+,形成的氢氧化铁胶体会吸附部分磷沉积在底部[26];同时,Fe3+也会直接与部分磷酸盐结合成沉淀吸附在沉积物表层,从而降低磷的交换通量。

      pH对沉积物中氮磷释放的影响如图5所示。可以看出,在偏酸或碱性条件下的氮通量显著大于中性条件下。在pH=7时,氨氮和总氮通量最低;在pH=5时,氨氮和总氮通量达到最大值,均为最小值的1.3倍。在酸性条件下,存在较多H+,会和体系中${\rm{NH}}_4^{+}$竞争吸附在胶体上[27],pH越低,竞争作用越强,氮释放越强;在碱性条件下,上覆水中存在的OH会和${\rm{NH}}_4^{+}$反应生成气态NH3,其可从上覆水中逸出,上覆水中氨氮减少,使得两相间浓度差变大,促进了底泥中氮向上覆水的释放。pH对底泥磷的影响主要通过吸附-解吸和离子交换2种作用得以实现。由图5可知:在中性条件下,正磷酸盐通量最小,在pH=9时达到最大,为中性条件下的1.4倍;随着pH增加,总磷通量也由最小值0.490 mg·(m2·h)−1 (pH=7时)增加到最大值0.536 mg·(m2·h)−1 (pH=9时)。在碱性条件下,磷释放主要以离子交换为主[28],体系中OH会与磷酸根发生离子交换,碱性愈强,交换作用愈强,导致磷通量增加;在酸性条件下,难溶性磷酸盐及吸附了磷的氢氧化物胶体会溶解,使磷脱离底泥进入上覆水,酸性愈强趋势越明显,磷通量越大;而在pH=7时,磷酸盐主要以${\rm{HPO}}_4^{2-}$${{\rm{H}}_2}{\rm{PO}}_4^ - $的形式存在,此时离子易与体系中的金属离子结合被底泥吸附,使得底泥中的磷不容易释放。

    • 1) BBD设计方案及测定结果。根据Box-Behnken模型设计要求,以温度、溶解氧、pH为自变量,以各营养盐交换通量为响应值进行实验设计。设有3组平行样,实测结果如表2所示,实验工况设计及通量平均值如表3所示。

      2)模型的建立和方差分析。利用BBD模型,对表3中的结果进行数据回归拟合,温度(X1)、溶解氧(X2)、pH(X3)与TN通量(FTN)、TP通量(FTP)、${\rm{NH}}_4^{+}$-N通量$({F_{{\rm{NH}}_4^ + - {\rm{N}}}})$${\rm{PO}}_4^{3-}$-P通量$({F_{{\rm{PO}}_4^{3 - } - {\rm{P}}}})$ 4个响应值的二次多项式回归模型分别如式(4)~式(7)所示。

      对4个模型方程分别进行了方差分析和显著性分析,结果如表4所示。4个模型的P值均<0.01,说明显著性高,失拟项也均大于0.5,即所得方程和实际拟合中非正常误差所占比例小,这说明所得模型可信度高,模拟精确。TN、TP、${\rm{NH}}_4^{+}$-N及${\rm{PO}}_4^{3-}$-P模型的变异系数分别为0.077、0.044 4、0.075 4、0.032 7,即数据离散程度较小,表明各方程拟合度好,模型回归性较好。4个模型的校正决定系数$R_{{\rm{Adj}}}^2$分别为0.907 5、0.922 8、0.914 5、0.958 2,即4个模型分别能解释90.75%、92.28%、91.45%、95.82%各自响应值的变化,总变异度分别有9.25%、7.72%、8.55%、4.18%不能用模型解释;且$R_{{\rm{Adj}}}^2$R2值接近。在此基础上,用回归方程进行模拟,图6为氮磷交换通量的实际值和预测值的对比,过原点的斜率为1的直线代表实际值和预测值完全吻合的情况[29]。可以看出,部分点落在直线上,其余点基本均匀地分布在直线两侧且偏离较小,这进一步说明4个模型与实验结果吻合度较高。由此可知,所建立的4个模型稳定性较高,拟合度良好,能够用来分析各环境因子对底泥氮磷交换通量的影响效果。

      表4可知,单因子温度(X1)和溶解氧(X2)对各营养盐通量的影响显著(P<0.01),pH(X3)无显著影响。交互项X1X2${\rm{PO}}_4^{3-}$-P通量影响显著(P<0.01),对TP和${\rm{NH}}_4^{+}$-N通量影响显著(P<0.05),对TN无显著影响;X1X3对TN和${\rm{NH}}_4^{+}$-N通量影响显著(P<0.05),对TP和${\rm{PO}}_4^{3-}$-P通量无显著影响;X2X3对4种营养盐通量均无显著影响。模型${X}_1^{1}$${X}_2^{2}$${X}_3^{3}$对各响应值的均有显著影响(P<0.01),说明环境因子与响应值之间不是简单的线性关系,二次项对响应值有着较大的影响。从整体显著性来看,与TN和TP相比,环境因子对离子态的${\rm{NH}}_4^{+}$-N及${\rm{PO}}_4^{3-}$-P的交换通量影响作用更大。

      3)响应面交互作用分析。为直观反映温度、溶解氧、pH 3个环境因子及其交互作用对各营养盐交换通量的影响,利用回归方程建立对应的等值线和响应面图。一般来说,等高线越接近于椭圆状,响应面坡度越大,说明交互作用越明显[30-31]。根据回归模型,本研究中选取交互作用效果较好的等高线图和曲面图。

      图7反映了温度和pH对TN通量交互作用的影响。可以看出,等高线图呈现明显的椭圆形状,且由图7(b)可以看出,最小响应值在曲面上,与沿着单一因子坐标轴走向的最小值不重合,表明温度和pH对TN通量有交互作用。图8为温度和溶解氧对TP通量交互作用的影响,等高线接近椭圆,响应面存在一定弯曲,表明温度和溶解氧交互作用对TP通量有一定的影响。由等高线图8(a)可知,沿溶解氧轴方向上的等高线比沿温度轴方向上的密集且陡峭,表明相对于温度,溶解氧对TP通量影响更为显著,这与方差分析结果一致(表4)。由图9可知,温度和溶解氧、温度和pH的交互作用均对${\rm{NH}}_4^{+}$-N通量的等高线呈椭圆,且响应面有一定坡度,表明温度和溶解氧交互作用、温度和pH交互作用均对${\rm{NH}}_4^{+}$-N通量影响显著。由图9(a)图9(c)可知,在实验范围内,当温度一定时,${\rm{NH}}_4^{+}$-N通量随pH变化波动很小,而随溶解氧变化很大,表明溶解氧对响应值的贡献更大,即溶解氧对${\rm{NH}}_4^{+}$-N通量的影响作用比pH更显著,这与方差分析结果一致。由图10可知,等高线及响应面形状表明了温度和溶解氧之间存在交互作用。由图10(a)可知,当溶解氧一定时,${\rm{PO}}_4^{3-}$-P通量随温度变化很小;当温度一定时,${\rm{PO}}_4^{3-}$-P通量却随溶解氧的增大而减小,这说明相比于温度,溶解氧对${\rm{PO}}_4^{3-}$-P通量影响更显著。

      4)最不利释放条件的预测与结果验证。在对内源污染进行防控时,可通过对环境因子的控制使得底泥释放量较小,所以对雁鸣湖底泥最不利释放(即释放量最小)条件进行了预测。图7~图10也显示响应面开口向上,即响应值有极小值。以响应面优化得到条件组合,通过响应面分析得出氮磷释放量最小时自变量的数值分别为X1=10 ℃、X2=7.87 mg·L−1X3=7.13。温度为10 ℃、溶解氧为7.87 mg·L−1、pH为7.13时,各响应值取最小值,TN交换通量为3.956 mg·(m2·h)−1,TP交换通量为0.471 mg·(m2·h)−1${\rm{NH}}_4^{+}$-N交换通量为1.469 mg·(m2·h)−1${\rm{PO}}_4^{3-}$-P交换通量为0.146 mg·(m2·h)−1

      为了检验响应面法所得的实验结果,通过实验加以验证,结合实验情况将条件修正为:X1=10 ℃、X2=8.0 mg·L−1X3=7.0,即在温度为10 ℃、溶解氧为8.0 mg·L−1、pH为7.0的条件下,进行了3次平行实验,实验结果如下:TN交换通量为3.779 mg·(m2·h)−1,TP交换通量为0.488 mg·(m2·h)−1${\rm{NH}}_4^{+}$-N交换通量为1.409 mg·(m2·h)−1${\rm{PO}}_4^{3-}$-P交换通量为0.141 mg·(m2·h)−1,与预测值很接近,误差均在5%以内,表明预测结果是可靠的,进一步说明所得模型能较好地预测雁鸣湖界面间营养盐的交换通量。因此,园区应注意对温度、溶解氧及pH进行控制,以防止内源释放对雁鸣湖环境造成更严重的污染。

    • 1)控制单因子变化条件下,雁鸣湖沉积物-上覆水界面间氮磷通量随着温度的升高而增加,30 ℃时达到最大;随溶解氧浓度的增加而减小,在缺氧条件下(DO=3 mg·L−1),氮磷通量取得最大值;pH为中性条件下最小,氮通量在酸性条件下(pH=5)达到最大,磷通量在碱性条件下(pH=9)达到最大。

      2)对环境因子与雁鸣湖沉积物-上覆水界面间营养盐通量的关系进行拟合,结果表明4种营养盐通量与环境因子间均符合二次多项式回归模型。

      3) 方差分析及响应面结果表明,温度和pH交互作用对TN通量影响显著;温度和溶解氧交互作用对TP通量影响显著;温度和溶解氧交互作用、温度和pH交互作用对${\rm{NH}}_4^{+}$-N通量影响显著;温度和溶解氧交互作用对${\rm{PO}}_4^{3-}$-P通量影响极显著。

      4)响应面优化结果表明,温度为10 ℃、溶解氧为7.87 mg·L−1、pH为7.13是雁鸣湖底泥释放的最不利条件,此时TN交换通量为3.956 mg·(m2·h)−1,TP交换通量为0.471 mg·(m2·h)−1${\rm{NH}}_4^{+}$-N交换通量为1.469 mg·(m2·h)−1${\rm{PO}}_4^{3-}$-P交换通量为0.146 mg·(m2·h)−1,各响应值取最小值。

    参考文献 (31)

目录

/

返回文章
返回