溶解氧对湖库热分层和富营养化的响应%26mdash%3B%26mdash%3B以枣庄周村水库为例.pdf

中国环境科学 2016,36(5)：15471553 china environmental science 溶解氧对湖库热分层和富营养化的响应 以枣庄周村水库为例 邱晓鹏,黄廷林*,曾明正 (西安建筑科技大学环境与市政工程学院,陕西 西安 710055) 摘要：为研究热分层和富营养化对湖库溶解氧变化特征的影响,于 2014 年 1 月12 月对周村水库水温、溶解氧、叶绿素以及初级生产力 的季节变化及垂向分布进行了监测.结果表明:水温和溶解氧的分层期均为 411 月份;分层期叶绿素在 2050g/l 之间,初级生产力为 2.162.23go2/(m3d),光补偿点在 13m 之间;恒温层在 5 月中旬进入厌氧状态;由于光补偿点位置较高,58 月份氧跃层位置为 16m,高于温 跃层上界面;而氧跃层位置偏高造成溶解氧在垂向上的极值一般在表层,且变温层溶解氧浓度梯度较大;911月份温跃层的下移使得氧跃层 和厌氧区界面同时下移,厌氧区界面与温跃层上界面的位置变化始终同步,而氧跃层受水体耗氧作用的影响在热分层结构相对稳定时会再 次上移.热分层和水体富营养化均对溶解氧的浓度和分布有重要的影响. 关键词：温跃层；光补偿点；初级生产力；氧跃层 中图分类号：x524,tv697.2 文献标识码：a 文章编号：1000-6923(2016)05-1547-07 responses of dissolved oxygen on thermal stratification and eutrophication in lakes and reservoirsan example in zhoucun reservoir in zaozhuang city. qiu xiao-peng, huang ting-lin*, zeng ming-zheng (school of environmental and municipal engineering, xian university of architecture and technology, xian 710055, china). china environmental science, 2016,36(5)：15471553 abstract：to investigate the effects of thermal stratification and eutrophication on dissolved oxygen, seasonal variations and vertical distribution of water temperature, dissolved oxygen, chlorophyll a and primary production were monitored from january 2014 to december 2014 in the zhoucun reservoir. the results showed that water temperature and dissolved oxygen both experienced stratification from april to november. the chlorophyll a was 2050g/l, the primary productivity was 2.162.23g o2/(md) and the position of light compensation point was between 1m and 3m during stratification. the hypolimnion has been in anaerobic condition since mid-may. because of the high position of light compensation point, the oxyclines position was between 1m and 6m during stratification, higher than thermoclines. the maximum value of dissolved oxygen was usually at the surface since the position of oxycline was high and the dissolved oxygen gradient was large in the empilimnion. the vertical migration of thermocline made oxycline and anaerobic zone interface both moved down between september and november. the anaerobic zone interface position was always the same as thermoclines, while the oxycline would move up again when the thermal structure was relatively stable. the thermal stratification and eutrophication have significant effects on dissolved oxygen concentration and distribution. key words：thermocline；light compensation point；primary production；oxycline 溶解氧在地表水体的生物化学循环中起着 关键作用,是湖库水生态系统生化条件改变的敏 感指标1 -2.当水体溶解氧浓度低于 6.5mg/l 时, 鱼类的生理学过程会受到影响3,水体进入厌氧 状态时,水生生物死亡4.沉积物会释放还原性污 染物如氨氮、正磷酸盐、铁、锰、硫化物等5 -6, 使水质恶化.研究湖库溶解氧的季节变化和垂向 分布,对于维护良好的湖库生态环境、保证城市 供水安全有重要的指导意义. 热分层是深水湖库所共有的湖沼学特性,分 层期间湖库溶解氧的垂向分布也有明显的分层 现象,热分层是影响溶解氧垂向分布的关键因 素 7 -8.此外,我国湖库正面临越来越严重的富营 收稿日期：2015-10-16 基金项目：国家科技支撑计划(2012bac04b02) * 责任作者, 教授, 1548 中 国 环 境 科 学 36 卷 养化问题9,富营养水体中藻类对溶解氧的浓度 和分布也有显著的影响10. 目前,湖库溶解氧的时空变化已受到越来越 多的关注.zhang等7,11探讨了千岛湖水库热力结 构对长期气候变化的响应以及溶解氧分层对热 分层的响应,张军毅等12研究了富营养水体梅梁 湖水域溶解氧的时空分布及与 ph 值的关系,曾 春芬等13分析了天目湖溶解氧的水平和垂向分 布特征,成晓奕等14采用三维水动力-水质模型 模拟了天目湖水温和溶解氧的分层过程,段燕 等 15探讨了千岛湖溶解氧的动态变化特征及其 影响因子.然而,以往的研究结合热分层结构变化 和藻类代谢活动分析溶解氧动态和垂向分布特 征的研究较少.此外,对湖库溶解氧的原位监测频 率较低(多为1个月1次)12 -13,不利于揭示不同季 节溶解氧的变化特征.本研究以我国枣庄周村水 库为研究对象,通过高频率的水温、溶解氧的原 位监测,结合水体叶绿素和初级生产力季节变化 和垂向分布情况,探讨了水库热分层和富营养化 对溶解氧分布的耦合作用,以期为富营养型深水 湖库的管理提供理论指导. 1 材料与方法 1.1 周村水库概况 1173930e 1174030e1174130e 345630n 34570n 345730n n 西伽河 徐洼河 下十河 监测点 出水口大坝 图 1 周村水库及主库区监测点位置 fig.1 location of the zhoucun reservoir and sampling site 周村水库始建于1960年,位于山东省枣庄市 东北角,是枣庄市重要的水源地之一.周村水库总 库容 8404 万 m3,防洪库容 3304 万 m3,设计控制 流域面积为 121km2,水库面积 6.5km2,平均水深 13m,最大水深 18m.周村水库属于温带大陆性季 风气候,年平均气温 13.9,夏季炎热多雨,冬季 寒冷干燥.20 世纪 90 年代,周村水库实行大规模 的网箱养鱼,虽然在 2008年 7月已完成了库区内 的网箱清理工作,但是周村水库沉积物污染仍然 严重16.20142015 年枣庄市供水总公司对周村 水库表层水体的监测结果表明,水库总氮的年平 均浓度为 1.63mg/l,总磷为 0.056mg/l,水体处于 富营养化状态. 1.2 采样与监测 本实验于2014年112月对周村水库主库区 的水温、溶解氧、叶绿素 a 和初级生产力进行了 监测.根据周村水库的地形特征,将监测点设置在 坝前深水区(图1),可以反映周村水库主库区水质 的变化特征.对不同水层的水温和溶解氧采用哈 希 hq30d 单路输入多参数数字化分析仪进行原 位监测,垂向间隔为 1m,监测频率每周 12 次.叶 绿素 a 的监测水深为 0.5,2.5,5,7.5,10,13,16m,每 个水深取水样 500ml 保存于聚乙烯瓶中,水样运 回实验室后立即用 0.45m 的醋酸纤维膜进行过 滤,并采用丙酮萃取分光光度法进行检测17,叶 绿素 a 的监测频率为每周 1 次.初级生产力的测 定采用黑白瓶测氧法17,每月中旬选取晴天进行, 培养时间为 24h,监测水深为 1,2,3,4,5,7.5,10m (当某水层的溶解氧小于 2mg/l 时,不监测),每个 水深设置 2 个对照. 1.3 数据处理与分析 本研究所指各水层初级生产力为培养 24h 后总初级生产力,即 24h 后的白瓶溶解氧浓度与 初始溶解氧浓度之差.目前,国内外对温跃层和 氧跃层仍没有统一的定义.因温度梯度大于 1 /m 的水层对湖库中溶解性物质以及浮游生物 的垂向分布影响显著18 -19,本研究定义温跃层 为温度梯度大于 1/m 的水层.采用 2 个溶解氧 垂向分布特征参数讨论水体溶解氧的垂向变化: 定义氧跃层为溶解氧浓度梯度最大的水层;厌 氧区为溶解氧浓度小于 0.2mg/l 的水层.温跃层 深度和氧跃层深度分别为温跃层和氧跃层上界 面处的水深. 5 期 邱晓鹏等：溶解氧对湖库热分层和富营养化的响应以枣庄周村水库为例 1549 数据统计和分析采用 spss19.0.由于氧跃层 和厌氧区深度在周年内不存在正态分布,其与环 境因子的相关性分析采用spearman秩相关系数. 2 结果与讨论 2.1 水温、溶解氧和叶绿素的季节变化 由图 2a 可知,周村水库存在明显的温度分 层.2014年表层(0.5m水深处)水温和底层(16m水 深处)水温的变化范围分别为 3.830.1和 3.8 13.4,表层水温受气温影响较大,呈现夏季高冬 季低的变化特征,底层水温变化缓慢,表底层温差 在夏季最大,冬季最小.周村水库热分层期为 4 11 月份,混合期为 1、2、3 和 12 月份. 水深(m) 0 2 4 6 8 10 12 14 16 0 4 8 12 16 20 24 28 32 0 2 4 6 8 10 12 14 16 0 2 4 6 8 10 12 14 16 水温() 溶解氧(mg/l) 月份 2014-01 2014-03 2014-05 2014-07 2014-09 2014-11 2015-01 0 2 4 6 8 10 12 14 16 0 5 10 15 20 25 30 35 40 45 50 叶绿素(g/l) a b c 水深(m) 水深(m) 图 2 2014 年周村水库水温、溶解氧及叶绿素 a 的季节 变化 fig.2 seasonal variations of water temperature, dissolved oxygen and chlorophyll a in zhoucun reservoir in 2014 表层和底层溶解氧的变化范围分别为 5.02 16.29mg/l和0.0012.70mg/l(图2b).冬季水库完 全混合,大气复氧充分,造成冬季水库溶解氧水平 整体较高(1113mg/l).3 月份,水体仍处于混合状 态,但是随着水温的升高,水中溶解氧的饱和度随 之降低,所以水库整体溶解氧浓度降低.进入分层 期,水库表层水温较高,表层水体受光合作用的影 响,在晴天无大风时常达到超饱和状态.周村水库 表层水体溶解氧最高为 16.29mg/l,大于千岛湖 的 14.91mg/l20和泽雅水库的 12mg/l21,这表明 周村水库表层水体藻类光合作用较强,水体富营 养化相对严重.由于温跃层的存在阻碍了溶解氧 的垂向传递,加之受水体耗氧和沉积物耗氧作用 影响,造成恒温层在分层期一直处于缺氧状态或 厌氧状态.分层期结束后,随着水库水体的完全混 合,整个水体在 11 月中旬再次进入好氧状态,并 随着气温的降低,溶解氧整体呈升高趋势. 表层叶绿素浓度为 3.254.4g/l,底层叶绿 浓度为 1.619.1g/l (图 2c).叶绿素呈现夏季高 冬季低的季节变化规律,最大值和最小值分别出 现在 5 月 28 日和 3 月 23 日.夏秋季叶绿素在垂 向上也存在分层现象,叶绿素的极值一般出现在 表层,这主要是由于周村水库光补偿点位置较高, 光补偿点以下的藻类因呼吸作用强于光合作用, 藻类生物量减少.分层期间表层水体叶绿素在 2050g/l, 高 于 湖 库 富 营 养 化 的 衡 量 标 准 10g/l22,周村水库处于富营养状态.在混合期, 即 13 月份和 12 月份,周村水库水温、 溶解氧和 叶绿素在垂向上无差异. 2.2 分层期水温、溶解氧和初级生产力的垂向 分布 如图 3 所示,水库表层水体(1m 处)初级生产 力在夏季(68 月份)较高,高达 2.162.23go2/ (m3d),春秋季较低,为 0.451.78go2/(m3d).周村 水库初级生产力水平与富营养型水体三峡水库 香溪河(水深 1m 处初级生产力为(2.191.88)go2/ (m3d)相近23. 周村水库4月份光补偿点位置最深,为56m, 初级生产力沿水深逐渐降低.周村水库属温带湖 库,藻类季节演替模式符合典型的 peg 模型24, 在春季,由于浮游动物的牧食作用而进入“清水 期”,此时水体的叶绿素处于全年最低水平(图 2c), 所以水体光衰减系数较低,光补偿点位置也较 低.510 月份,周村水库光补偿点在 13m 之间, 1550 中 国 环 境 科 学 36 卷 这主要是由于该时期藻类大量繁殖,水体光衰减系数较高,藻类光合作用随水深的增加迅速降低. 18 16 14 12 10 8 6 4 2 0 0 5 10 15 2025 水深(m) 水温() 0 4 8 1216 溶解氧(mg/l) -2 -1 0 1 23 初级生产力go2/(m3d) a. 4 月 17 日 18 16 14 12 10 8 6 4 2 0 05 10152025 水深(m) 水温() 03 691215 溶解氧(mg/l) -2-1 0123 b. 5 月 15 日 初级生产力go2/(m3d) 18 16 14 12 10 8 6 4 2 0 51015 20 25 30 水深(m) 水温() 036 9 12 15 溶解氧(mg/l) -2-10 1 2 3 c. 6 月 10 日 初级生产力go2/(m3d) 18 16 14 12 10 8 6 4 2 0 5 10 15 20 2530 水深(m) 水温() 0 2 4 6 810 溶解氧(mg/l) -2 -1 0 1 23 d. 7 月 17 日 初级生产力go2/(m3d) 18 16 14 12 10 8 6 4 2 0 510 15202530 水深(m) 水温() 02 4 681012 溶解氧(mg/l) -2-1 0123 e. 8 月 17 日 初级生产力go2/(m3d) 18 16 14 12 10 8 6 4 2 0 0510 15 20 25 水深(m) 水温() 024 6 8 10 溶解氧(mg/l) -2-10 1 2 3 f. 9 月 20 日 初级生产力go2/(m3d) 18 16 14 12 10 8 6 4 2 0 0 51015 2025 水深(m) 水温() 0 246 810 溶解氧(mg/l) -2 -101 23 g. 10 月 13 日 初级生产力go2/(m3d) 18 16 14 12 10 8 6 4 2 0 05101520 水深(m) 水温() 02468 溶解氧(mg/l) -2-10123 h. 11 月 8 日 初级生产力go2/(m3d) 图 3 分层期水温、溶解氧和初级生产力的垂向分布 fig.3 vertical distribution of water temperature, dissolved oxygen and primary production during stratification 水温溶解氧初级生产力 在热分层期,溶解氧浓度的极值均出现在表 层(水深 0.5m 处),高于其他深水湖库溶解氧极值 出现的位置20,25,这主要是水库光补偿点位置较 高所致.4、10 和 11 月份水库溶解氧和水温的垂 向分布特征相似,而 58 月份溶解氧在变温层的 变化幅度明显高于水温.这主要是由于 58 月份 周村水库光补偿点在 13m 之间,2m 以下水体一 般处于耗氧状态,如 6 月 10 日,2m 水深处初级生 5 期 邱晓鹏等：溶解氧对湖库热分层和富营养化的响应以枣庄周村水库为例 1551 产力为-1.71go2/(m3d).同时,由于夏季水库受风 力影响较小,大气复氧作用不充分,2m 以下水体 溶解氧迅速降低.在 4、 10 和 11 月份,水温和溶解 氧在变温层的垂向差异性较小,这表明春秋季水 库变温层的混合强度大于夏季. 2.3 溶解氧分层的季节变化及影响因素 0 4 8 12 16 03-204-2106-1007-3009-1811-712-27 日期 水深(m) 温跃层(上) 温跃层(下) 氧跃层 a 0 4 8 12 16 03-204-2106-1007-3009-1811-712-27 日期 水深(m) 温跃层(上) 温跃层(下) 厌氧区界面 b 图 4 分层期温跃层、氧跃层和厌氧区界面季节变化 fig.4 seasonal variations of thermocline, oxycline and the anaerobic zone interface during stratification 周村水库主库区温跃层深度和氧跃层深度的 季节变化明显,而温跃层下界面的位置则相对稳定 (图 4a).4 月份周村水库温跃层形成,氧跃层的深度 为 47m,低于温跃层深度,这主要是由于 4 月份周 村水库光补偿点在 56m,5m 以上溶解氧变化较 小.58 月份,温跃层深度会随着气温的波动上下浮 动,而氧跃层受藻类光合作用的影响其深度的变化 范围大于温跃层.特别是在风力较小,大气复氧作 用较弱时,周村水库光补偿点较低,氧跃层深度只 有 1m(如 6 月 25 日,8 月 13 日).911 月份,气温降 低,温跃层上界面逐渐下移,氧跃层深度也随之增 加.值得注意的是,当温跃层上界面发生下移时,氧 跃层深度随之增加;而当温跃层界面保持相对稳定 时,氧跃层深度又减小.这主要是由于当温跃层上 界面下移时,变温层水体混合强度较大,变温层的 耗氧作用会被大气复氧抵消;而当温跃层界面相对 稳定时,变温层水体也相对稳定,光补偿点以下的 水体耗氧不能得到及时的补充,造成氧跃层上移. 如图 4b所示,周村水库在 5月中旬开始形成 厌氧区,且厌氧区的起始深度与温跃层下界面相 近,这表明周村水库分层形成 1 个半月后恒温层 进入厌氧状态.厌氧区形成后,周村水库温跃层内 的水体耗氧作用使得厌氧区界面逐渐上移.至 6 月份,水体厌氧区界面与温跃层上界面重合,且此 后二者的变化基本同步,这表明水体热分层结构 决定了周村水库厌氧区的范围. 表1 周村水库溶解氧分层与环境因子的相关性分析 table 1 correlation analysis of oxygen stratification and environmental factors in zhoucun reservoir 参数 氧跃层深度(m) 厌氧区深度(m) 日平均气温() -0.734* -0.869* 温跃层深度(m) 0.686* 0.842* 叶绿素(g/l) -0.459* -0.511* 注:*在置信度(双侧)为0.05的条件下显著相关;*在置信度(双侧) 为0.01的条件下显著相关. 对溶解氧垂向分布特征参数与环境因子做 spearman 相关性分析可知,氧跃层深度和缺氧区 深度均与日平均气温和温跃层深度在 p0.01 的 条件下显著相关,与叶绿素在p0.05的条件下显 著相关.相关性分析表明,周村水库热分层和浮游 植物均对水库溶解氧的分布影响显著. 2.4 讨论 热分层是影响 do 垂向分布的关键因素8. 热分层期间,水体温跃层会阻碍变温层和恒温层 的能量交换和物质迁移,温跃层的垂向迁移,会改 变水库溶解氧的垂向分布形式.春季温跃层的形 成促进了水体溶解氧分层;秋季周村水库温跃层 的下移使得厌氧区范围随之减小.已有研究表明, 秋季 “翻库” 现象还会将恒温层中的沉积物和还 原性污染物带入到表层水体,从而增加水体浊度 降低表层水的溶解氧浓度26 -27.但由于周村水库 秋季降温相对缓慢,温跃层的并没有在短时间内 1552 中 国 环 境 科 学 36 卷 快速下移,这使得翻库期间单位时间内被带入表 层水体的污染物量较少,大气复氧可以及时补充 被污染物消耗的氧,所以周村水库变温层水体并 未因翻库出现缺氧现象.热分层期间恒温层中水 体和沉积物消耗的氧不能得到补充,造成恒温层 溶解氧持续降低,最后进入厌氧状态28.在厌氧 状态下,湖库沉积物又会向水体中释放大量还原 性污染物,造成恒温层水体和温跃层水体中还原 性污染物超标,化学需氧量增大5,这进一步降低 了温跃层溶解氧的浓度.而变温层水体在热分层 期间始终可以获得充足的大气复氧和光合作用 复氧,从而维持了较高的溶解氧浓度. 此外,热分 层还会促进水体富营养化.热分层可使藻类获得 更高的水温、 稳定的水环境、 充足的太阳辐射以 及更多的营养盐(如氨氮、 正磷酸盐)29,促进表层 水体藻类的增长,从而加剧富营养化对溶解氧的 影响. 水体富营养化对深水湖库溶解氧浓度变化 和垂向分布均有较大的影响.在浓度变化方面,富 营养水体表层有较强的光合作用,使得表层水体 溶解氧常常达到超饱和状态;此外,富营养水体会 给湖库沉积物和底层水体带来更多的有机物从 而导致需氧量的上升,恒温层水体迅速进入缺氧 或厌氧状态4,7.王晓东等21研究表明,与周村水 库水深条件相近、 处于轻度富营养化的泽雅水库, 其恒温层在分层形成后 5 个月才进入厌氧状态. 而周村水库恒温层形成厌氧状态仅用了 1 个半 月,这表明水库富营养化促进了恒温层水体的氧 的消耗.在垂向分布方面,由于叶绿素浓度与水体 透明度呈显著负相关30,富营养水体中光衰减系 数的较大.这造成周村水库光补偿点位置较低,使 得变温层水体也存在较强的水体耗氧作用.当大 气复氧不充分时,变温层水体溶解氧浓度梯度迅 速增加,氧跃层发生上移.因此,周村水库垂向溶 解氧的极值一般在表层且变温层在个别时期也 出现了厌氧区(如 6 月 25 日,9 月 28 日). 3 结论 3.1 周村水库溶解氧分层与热分层均为 4 月份 出现,11 月份消失,溶解氧与水温在 4、10、11 月 份的垂向分布形式相似,58 月份氧跃层的位置 普遍高于温跃层. 3.2 热分层结构对水库溶解氧的垂向分布起决 定性作用,热分层和溶解氧分层生消时间一致,厌 氧区界面与温跃层上界面的位置变化同步,温跃 层的垂向迁移对氧跃层的位置也有很大影响. 3.3 周村水库处于富营养状态,藻类的过量繁 殖造成热分层期表层水体溶解氧常常过饱和,恒 温层溶解氧迅速减少,光补偿点在 13m 之间,水 库变温层溶解氧浓度梯度较大,溶解氧在垂向上 的极值出现在表层水体. 参考文献： 1 jones i d, winfield i j, carse f. assessment of long-term changes in habitat availability for arctic charr (salvelinus alpinus) in a temperate lake using oxygen profiles and hydroacoustic surveys j. freshwater biology, 2008,53(2):393-402. 2 孙寓姣,陈 程,丁爱中,等.官厅水库水质特征及水体微生物多 样性的响应 j. 中国环境科学, 2015,35(5):1547-1553. 3 pollock m s, clarke l m j, dube m g. the effects of hypoxia on fishes: from ecological relevance to physiological effects j. environmental reviews, 2007,15:1-14. 4 welch e b, cooke g d, jones j r, et al. do-temperature habitat loss due to eutrophication in tenkiller reservoir, oklahoma, usa j. lake and reservoir management, 2011,27(3):271-285. 5 mueller b, bryant l d, matzinger a, et al. hypolimnetic oxygen depletion in eutrophic lakes j. environmental science & technology, 2012,46(18):9964-9971. 6 曾 康,黄廷林,马卫星,等.金盆水库汛期高浊水径流的潜入及 热分层水体水质响应 j. 中国环境科学, 2015,35(9):2778- 2786. 7 zhang y, wu z, liu m, et al. dissolved oxygen stratification and response to thermal structure and long-term climate change in a large and deep subtropical reservoir (lake qiandaohu, china) j. water research, 2015,75(2):49-58. 8 jones j r, knowlton m f, obrecht d v, et al. temperature and oxygen in missouri reservoirs j. lake and reservoir management, 2011,27(2):173-182. 9 陈小锋,揣小明,杨柳燕.中国典型湖区湖泊富营养化现状、 历史 演变趋势及成因分析 j. 生态与农村环境学报, 2014,30(4): 438-443. 10 foley b, jones i d, maberly s c, et al. long-term changes in oxygen depletion in a small temperate lake: effects of climate change and eutrophication j. freshwater biology, 2012,57(2): 278-289. 5 期 邱晓鹏等：溶解氧对湖库热分层和富营养化的响应以枣庄周村水库为例 1553 11 zhang y, wu z, liu m, et al. thermal structure and response to long-term climatic changes in lake qiandaohu, a deep subtropical reservoir in china j. limnology and oceanography, 2014,59(4):1193-1202. 12 张军毅,黄 君,严 飞,等.梅梁湖水体溶解氧特征及其与ph的 关系分析 j. 复旦学报(自然科学版), 2009,48(5):623-627. 13 曾春芬,黄文钰,王伟霞,等.天目湖溶解氧分布特征及环境影响 因子 j. 长江流域资源与环境, 2010,19(4):445-451. 14 成晓奕,李慧赟,戴淑君.天目湖沙河水库溶解氧分层的季节变 化及其对水环境影响的模拟 j. 湖泊科学, 2013,25(6):818- 826. 15 殷 燕,吴志旭,刘明亮,等.千岛湖溶解氧的动态分布特征及其 影响因素分析 j. 环境科学, 2014,35(7):2539-2546. 16 曹占辉,黄廷林,邱晓鹏,等.周村水库沉积物污染物释放潜力模 拟 j. 环境科学与技术, 2013,36(9):41-44. 17 国家环境保护总局水和废水监测分析方法编委会.水和废 水监测分析方法 m. 北京:中国环境科学出版社, 2002. 18 fee e j, hecky r e, kasian s e m, et al. effects of lake size, water clarity, and climatic variability on mixing depths in canadian shield lakes j. limnology and oceanography, 1996, 41(5):912-920. 19 gauthier j,