1. 引言
湖泊现代过程研究有助于了解湖泊在自然和人类活动影响下生态环境变化历程 [1] 。高原湖泊水体理化参数的季节性垂直分布特征研究对于深入理解不同指标的环境指示意义、探讨高原湖泊生态系统的变化十分重要，因此有必要开展长期连续的高原湖泊现代过程监测与研究。
湖水的热量分配会影响湖水的分层与混合，进而影响到湖水中的其他溶解性物质及其理化参数变化 [1] [2] 。通过对水温的测定表明，湖泊温度具有很强的季节性变化，分层也呈现出季节性格局，因此根据水温就能对湖泊进行识别和分类。湖泊水温不仅决定了湖泊的分类和理化状况，还决定了生物群落的结构和水生生态系统的生产力 [2] 。在深水型湖库中，存在季节性分层现象，并且水温分层温差很大，如天山天池夏季出现明显水温分层，并且影响水体溶解氧浓度、电导率、pH值、叶绿素-a浓度和蓝绿藻细胞密度在水体中的分布 [3] 。这说明季节性温度成层现象会影响紊流的上下混合运动，从而影响到湖泊的生态模型变动。与深水湖库相比，高原浅水湖泊一般并不形成水体分层，由于监测数据缺乏也不系统，研究程度有待提高。
近年来，高原湖泊现代基本理化特征研究与综合考察对湖沼学的发展起到巨大推动作用 [4] [5] 。就青藏高原地区而言，王君波 [5] 等对西藏纳木错水深分布及现代湖沼学特征进行分析，研究表明：纳木错湖水的湖沼学参数垂向变化趋势分明，表层变温层、中部温跃层与底部均温层的理化性质差异明显；黄磊 [6] 等针对纳木错水温变化及热力学分层进行深入研究，发现纳木错为典型双季对流、完全混合型湖泊，不同湖盆形状及水深分布是造成其中部与东部两个湖盆热力学特征差异的主要原因。然而在云贵高原地区对于长期监测湖泊基本理化特征的研究并不多见也不系统，赵雪枫 [7] 、文新宇等 [8] 对泸沽湖温度分层季节变化及其环境效应的研究，其结果揭示了泸沽湖季节变化的过程，及其对溶解氧、pH、叶绿素-a和总磷的时空分布影响。
以往对于洱海的研究多集中在营养化的评估、水质及沉积物微体分析等领域 [9] [10] [11] [12] [13] ，本文通过对洱海水体2015年1月~2016年4月的监测数据的分析，探讨近年洱海水体理化参数的垂直分布特征与空间分布特征和季节性变化过程，为洱海富营养化和污染治理提供科学依据。
2. 材料与方法
2.1. 研究区概况
洱海位于云南大理白族自治州中部(100˚05'E~100˚17'E, 25˚36'N~25˚58'N)，为一构造断陷形成的高原浅水湖(图1)。洱海属澜沧江流域，发源于洱源县的芘碧湖，源头出自罢谷山，北有弥直河注人，东南收波罗江水，西纳苍山18溪水，西出洱河，湖水依赖地表径流和湖面降水补给，年入湖水量13.78 × 108 m3，湖面面积256.5 km2，蓄水量30亿立方米，平均水深约10.5 m，最深可达20.9 m，湖水停留时间约2.75年。洱海属亚热带高原季风气候，四季温和，平均温度较低，日较差较大，光照充足，干湿季节分明，据1951~1980年气象资料，年平均年降水量约为1078.9 mm；5~10月雨季降水量占全年降水量的85%~96% [14] (图2)。
2.2. 数据采集与研究方法
为较全面、系统的了解洱海水质现状与季节性变化，在洱海从北到南共设置3个点进行定点监测，分别位于洱海的北部湖区、中部湖区与南部湖区。在采样现场以1 m为间隔，运用YSI(预先校准，测量结果可靠性及对比分析另有专文讨论)测定水样的叶绿素a(Chla)、温度(Tem)、溶解氧(DO)、pH、浊度(Turbid+)，并于2016年5月19日采集柱状水样，低温封装带回，进行总磷(TP)、总氮(TN)实验测试。水质评价主要采用相崎守弘的湖泊营养程度评分标准对各采样点进行评分(保留两位小数)，然后用卡森指数法(见表1)对各监测点进行分级 [15] 。
3. 结果分析与讨论
3.1. 湖水基本理化特征
3.1.1. 水体中叶绿素垂直分布特征
叶绿素a是表征浮游植物现存量的主要参数，也是湖泊水质的重要指标之一 [16] ，叶绿素a浓度的高低变化可以反映水体的营养状况。以前的研究结果认为，洱海的叶绿素a呈现明显的季节变化，即4月整体含量开始逐渐上升，7月大幅增大，并在10月达到峰值，在次年1月骤降 [16] 。但我们的监测数据显示出以下几个特点：1) 2015年1月的叶绿素在中部和南部水域表层呈现明显的高值，但在水深0~7 m之间迅速递减，其叶绿素a浓度从36.3 μg/L递减到3.7 μg/L，随着水深增加稳定在3 μg/L左右。2) 随
Figure 1. Location and selected monitoring stations in Erhai Lake
图1. 洱海地理位置与选取的监测点
Figure 2. Average monthly temperature and precipitation in 1951-1980 in Dali
图2. 大理1951~1980年月均气温与降水量
Table 1. Grading and classification standard
表1. 评分与分级标准
着季节变化，叶绿素a的浓度在4月稍有增长但整个湖泊从南到北十分均一，为3.5~5.6 μg/L，10月增长为8.8~17.7 μg/L，在2016年1月又下降为0.6~5.5 μg/L，在2016年4月又增长为2.1~6.6 μg/L (图3)。
从空间分布来看，2015年1月洱海中部和南部湖区的叶绿素浓度在表层均呈现由表层高值大幅递减的分布趋势，而北部湖区则呈现稳定的低值，这与中部和南部人口密度大、湖泊较浅，特别是南部湖泊中间水位反而较浅的特征相关(可能与湖流和湖底地形有关)。在2015年的10月份明显表现出较浅的南部湖区叶绿素a浓度整体偏高，北部湖区次之，中部湖区的叶绿素a浓度相对较低的特点。其他季节三个测区除表层有所差异，湖水较深处叶绿素a浓度变化不明显。
3.1.2. 水温分布
就我们目前监测的水温季节性变化而言，洱海1月平均水温最低(图4)，稳定在11℃左右，但在2015年1月表层水温达12.5℃。4月水温上升至16℃~17℃之间，7月洱海平均水温最高，且表面水温相对较高，随湖泊水体的深度增加迅速降温稳定在23℃~24℃，时至10月，水体温度下降至20℃左右。整体来讲，洱海水温在垂直断面上十分一致，不同部位之间差异很小，这可能与水体温度传输均匀且稳定、并主要受环境大气温度影响有关。另外，从2015年7月不同湖区水温差异可以看出，南部湖区从水温从表层到下部维持在23℃左右，而北部与中部湖区尽管下层水温维持在23℃左右，但表层水温则达26℃~27℃，可能与南部湖区风力扰动作用所致 [17] 。
图5展示了洱海在2015年1月到2016年4月的周降水量与周平均气温，结合全湖各监测点的表层水平均温度，可以发现表层水温、周平均气温与周降水的季节变化显著，水温与周气温变化存在很好的同步协调变化特征，并且表层水温总是略高于周平均气温。
3.1.3. 溶解氧垂直分布特征
溶解氧是指溶解在水中的分解态氧，是维持水体生态环境动态平衡的重要因子和维持水生生物生存的必要条件，同时溶解氧还参与部分物质转化 [18] 。有关研究表明，较清洁的湖泊溶解氧一般在7.5 mg/L以上，当湖水溶解氧含量低于5 mg/L以下时，浮游生物则无法生存 [19] 。实测数据显示，洱海的溶解氧季节变化分明，1月整体浓度最高稳定在8~9 mg/L，7月整体最低但变化幅度较大，从4~8.6 mg/L，4月
Figure 3. From January 2015-April 2016 water temperature changes in Erhai Lake
图3. 洱海2015年1月~2016年4月水温变化
Figure 4. From January 2015-April 2016 vertical distribution of chlorophyll-a in Erhai Lake
图4. 洱海2015年1月~2016年4月叶绿素垂直分布
Figure 5. Erhai lake surface water temperature and weeks precipitation, weeks average temperature
图5. 洱海表层水温与周降水、周平均气温图
和10月则处于转换过渡时期(图6)。
湖区不同部位水体DO在1月和4月差异不大，但在7月和10月出现明显分异。7月的溶解氧含量北部湖区的波动 > 中部湖区 > 南部湖区，而10月则表现为南部湖区溶解氧含量 > 北部湖区 > 中部湖区。
3.1.4. pH垂直分布特征
洱海湖水的pH值表现为以弱碱性为特征(图7)，整体来看4月和10月各湖区的pH值高于7月和1月，其中4月pH值在8.7左右变化，且深层湖水的pH相比于表层小且稳定。7月和10月不同湖区的pH差异比较明显，且以7月北部湖区表层pH值最高，随后随水深增加迅速降低；南部湖区次之，底层略高于表层；其中中部湖区的pH值最低，并且中部湖区在水深2~9 m 处出现pH值突然降低然后又升高的过程；10月3个湖区的水体剖面pH值虽然垂向变化不大，却呈现明显空间异质性，出现南部湖区pH
Figure 6. From January 2015-April 2016 dissolved oxygen (DO) vertical distribution in Erhai Lake
图6. 洱海2015年1月~2016年4月溶解氧(DO)垂直分布
Figure 7. From January 2015-April 2016 pH vertical distribution in Erhai Lake
图7. 洱海2015年1月~2016年4月pH垂直分布
值 > 北部湖区 > 中部湖区，相对而言，最高值的南部湖区与最低值的中部湖区pH值相差0.4左右。
3.1.5. 浊度垂直分布特征
洱海水体浊度的季节差异较大(图8)，7月和10月的浊度及其波动均大于1月。其中1月份的浊度在1 NTU左右波动，7月份的浊度在0.3~3.6 NTU之间波动，10月在1~3 NTU之间波动。南部湖区浊度在4月和7月高于北部湖区和中部湖区，10月北部湖区的浊度高于南部湖区和中部湖区。
3.1.6. 总磷(TP)总氮(TN)含量与水质评价
评价水体富营养化的方法主要有：营养状态指数法、营养度指数法和评分法。营养状态指数法中根
Figure 8. From January 2015-April 2016 turbidity vertical distribution in Erhai Lake
图8. 洱海2015年1月~2016年4月浊度垂直分布
据水体透明度制定的卡森指数法是最常用的评价水体富营养化的方法之一。后来，日本的相崎守弘等人提出了修正的营养状态指数(TSIM)，即以叶绿素a浓度为基准的营养状态指数 [20] 。其中总氮(TN)通过碱性过硫酸钾消解–紫外分光光度法，总磷(TP)通过钼酸铵分光光度法测得，样品的取样、保存及分析监测方法均依照《水和废水监测分析方法》进行 [21] ，并结合洱海水体中不同深度的叶绿素结果，经过计算分析得到表2。
从表2中可以看出洱海3个湖区的营养级均为中营养，整体而言，较深的中部湖区和北部湖区相比于南部湖区评分较低。总氮相对总磷而言对营养级评价的贡献率更高，但相差不大，且都在各湖区底部出现高值。
3.2. 讨论
王立前 [22] 等人在1999~2003年对云南省主要湖库水的叶绿素a (μg/L)浓度进行了调查，统计结果显示洱海的叶绿素a浓度范围为0.5~6.0 μg/L，显然如今洱海的叶绿素a浓度相比之前的统计结果有很大的增加。从2015年4月~2016年4月监测数据分析可以看出，可能由于洱海夏季大量藻类的爆发，而秋季大量浮游植物死亡，10月出现了较高的叶绿素a浓度，这与洱海浮游植物数量的季节变化 [23] 是一致的。但在2015年1月，中部湖区表层水中叶绿素a浓度异常偏高，但随深度增加迅速下降，在水深7米处形成稳定低值，在南部湖区也呈现相似的特点，只有北部湖区没有体现。吴艳龙等 [24] 曾对滇池水华束丝藻–微囊藻演替过程及其机制进行研究，发现水温是驱动两者复苏、生长和演替的重要环境因子。水华束丝藻休眠孢子在3月份(水温14℃左右)开始复苏，在无竞争压力下持续生长。从洱海2015年1月份的水温数据可以推断，很可能由于较暖的气温以及人类活动在南部和中部湖区注入大量的营养物质，短期内催生出大量浮游植物有关，但受冬季太阳辐射强度影响，水体热量停留在表层，随深度增加迅速降低，这在2015年的水温中可以看出水深0~2 m比下层水温高出1℃~2℃，这些浮游植物也受热量影响停留在表层，但在较远的北部湖区尚未受到影响。这种情况给我们一个警示，虽然洱海目前保存为中营养水平，但即使再气温最低的1月份，突然的气温升高同样会导致藻类的暂时爆发，形成生态危机，2017年初的情况就是这种情况的演示。另一方面，从全年来看，南部湖区叶绿素a浓度普遍高于北部与中部湖区，
Table 2. Water quality evaluation in Erhai Lake
表2. 洱海水质评价表
注：叶绿素a评价见评分1，TN评价见评分2，TP评价见评分3，EH-1表示洱海北部采样点，EH-2代表洱海中部采样点，EH-3代表洱海南部采样点。
这可能是因为湖水较浅，混合程度相对较好，光照充足，且南部湖区风速较大 [16] ，风力强烈扰动的水团带动着湖水充分搅动，使水温、藻类与营养物质垂向混合均匀，整个湖区垂向分布的浮游植物在相似的条件下繁殖更快有关。
水温是叶绿素a浓度变化的重要影响因素，也是浮游植物生长的关键因子 [25] [26] [27] 。与湖面附近气温相似，湖泊的水温也存在明显的季节性变化，相应的也出现不同的分层格局。气温可以直接影响到湖水的表层水温，但洱海作为大型浅水湖泊，往往比热容较大，因而出现了表层水温略高于周平均气温的现象。从不同的时间序列来看整个洱海水温垂向剖面，从表层到底层都受到气温变化的强烈影响，这与湖下层水温常年稳定的深水湖库 [7] 截然不同。就空间分布而言，洱海的主风向为西南风，杨威 [16] 等人通过监测发现，洱海的风速在空间上存在明显差异，南部湖区 > 中部湖区 > 北部湖区，况且南部湖区最浅，风力扰动强度大，可以从图2中展示的水温变化明显看出，南部湖区从表层到底层的水温完全受到季节性气温变化的影响，水体垂直剖面温度保持一致，这表明南部湖区在风力影响下的湖水混合彻底，不存在水温垂向变化。而北部湖区与中部湖区由于风速较低且水深大于南部湖区，因此在7月出现了暂时性的分层现象，但由于洱海风力很强，风生流使得冷却的上部水体密度差异变小，抗拒风力的能力下降，温跃层的厚度就会发生变化；研究表明 [28] ，水色重且较为浑浊的湖泊表层水获取的太阳能更多，表层水和深层水的温差增加，因此往往温跃层较浅；此时表层冷却的温度与温跃层温度相同，在7月温跃层很容易与湖上层融合，而形成洱海北部湖区与中部湖区在7月的分层模式。
洱海的溶解氧显示了鲜明的季节变化，形成整个1月值最高、7月最低的势态，这可能是因为1月洱海受到冬季低温、大风和较弱太阳辐射，使水体从表层到底层垂向充分接触大气，实现氧气交换，复氧过程频繁。而7月温度高、降水量较大，藻类死亡分解作用消耗大量氧气，加上混合层以下水中微生物分解水中有机物作用的加强，三个湖区的平均溶解氧浓度都有显著降低。10月南部湖区的溶解氧高于北部和中部除了受到湖区深度和风力作用的影响，与叶绿素a也有关系(10月的叶绿素a同样有此趋势)，通过SPSS相关性分析，洱海10月的三个湖区溶解氧浓度与叶绿素浓度显著相关(r = 0.6, p < 0.01)。当湖水中的叶绿素a的浓度偏高时，藻类的生物量也相对较高，产生的光合作用大量消耗CO2，并释放O2，使得溶解氧含量升高，呈现出南部 > 北部 > 中部的势态。由此推测，在洱海的7月可能由于大量浮游植物生存于水深2m左右，因此7月的北部湖区和中部湖区都出现在水深3 m以上溶解氧含量偏高，而在水深3 m后由于缺少光照，浮游植物难以产生光合作用，加上微生物分解耗痒及水生生物呼气耗痒出现了大幅度下降。而南部湖区则由于风速较大藻类不易汇聚成层，并且风力带动了湖水的垂向充分混合，溶解氧很低且垂向分布较为均匀。
以湖泊整体来看，pH的季节性变化主要是由于水生植物光合作用消耗大量CO2，从而使pH值在1月很低，从4月开始增大，7月垂向波动较大且空间分异开始出现，10月由于气温开始逐渐降低，因缺乏生物作用消耗而使pH值有所下降，并且有明显空间分异。7月份在水深2~9 m处达到pH最低值，之后又缓慢增加。这与西藏纳木错 [5] 以及三峡水库大宁河库湾 [29] 的研究结果中显示的变化趋势相似。出现这样的现象可能与7月份在水深1~3 m层位大量出现的浮游植物有关，因夏季湖泊表层藻类的大量爆发给湖中的大量生物提供了丰富的食物来源，但是正因为表层藻类的覆盖3 m以下的浮游动植物、微生物受到光照的影响较小，反而这些生物的呼吸作用释放的大量CO2使得pH有所降低。因位于湖泊中部，水深约20 m，湖泊底层的生物量远低于表层与中层，因此湖泊下层pH未出现显著变化。10月份3个湖区pH值空间分异明显，结合叶绿素a浓度数据与溶解氧浓度数据，可以明显看出藻类活动对于湖水pH值变化影响很大，进而可能影响到湖泊中其他生物生存。
浊度是湖泊水环境特征之一，其季节性变化可能与水中的生物量有关。经过相关性分析发现，湖水浊度只在中部湖区与叶绿素a浓度显著相关(p < 0.01)，说明湖泊水深是影响浊度受到湖泊内部生物量驱动的关键性因子。深水环境下，湖水的流动性较差，水流不易混合，较多的生物量(如藻类)在水体内垂向分布不均匀，会显著体现在水体浊度数据中。浮游植物春季开始逐渐增加，夏季大量产生，秋季消亡，使得湖水中的浊度也随之发生变化。
目前洱海虽然仍处于中营养化阶段，但是通过评分可以看出已经存在向中富营养化发展的趋势。从总磷、总氮的空间分布上分析，洱海底部氮磷的增加，可能是由于底泥和沉积物中氮和磷的释放。并且受人类活动强烈影响下，氮和磷的双重驱动很可能加速湖泊的富营养化进程。
4. 结论
湖泊调查是湖泊学研究的基本工作，基于对洱海跨时一年多的监测数据，对其现代湖沼学参数进行了初步的分析，并得到了对洱海的一些基本认识，获得以下结论：
洱海作为高原浅水湖泊，不存在典型的湖泊温度分层，仅在7月份在中部湖区与北部湖区出现温跃层和湖上层的混合现象。洱海较浅的南部湖区混合度 > 北部湖区 > 中部湖区。水温的季节变化会影响叶绿素a的变化，但人类活动排放的营养物质同样会在叶绿素a的增长中起重要作用；较浅的南部湖区更易受到人类活动与自然因素(如风力、吹程、湖流)等多重影响。
洱海10月的数据在叶绿素a浓度、溶解氧浓度和pH值之间均表现出显著的相关性。洱海水体的溶解氧含量存在季节变化，在10月的空间分布显著；溶解氧不仅受到湖泊中藻类的影响也受到入湖河流的影响。洱海的pH值季节性差异明显，10月同样出现明显空间差异；pH值易受水生动植物的光合作用和呼吸作用的影响，7月在湖泊水体中波动较大。洱海水体中的浊度在季节性变化中深水湖区受浮游动植物量的影响很大。
致谢
赵帅营老师、张云鹰、文新宇、廖君、王教元、卢同平同学在野外监测和室内数据分析处理过程中给予了帮助，在此一并表示感谢。
基金项目
云南省领军人才项目“云贵高原湖泊演化与水安全”(2015HA024) 和高端人才引进项目“云南(云贵高原)湖泊记录与生态环境及可持续发展研究”(2010CI111)资助。