湖泊富营养化分析.docx

湖泊富营养化分析湖泊富营养化导致的藻类暴发一直是我国最为突出的水环境问题之一藻类过度生长是 多种因素共同作用的结果，既包括水温、光照、风速等自然因素，也包括氮(N)、磷(P)、 铁(Fe)等营养物过量排放的人类活动因素.在诸多因素中，全球气候变暖背景下的水温变化 与高强度人类活动所引起的N、P排放增加被认为是导致湖泊富营养化最关键因素，因此， 同时考虑水温、N、P因子变化的湖泊富营养化相关研究在逐渐增多，但温度与营养物对 湖泊藻类生物量的交互作用等还需要深入研究［7］，比如水温、N、P促进藻类生物量增长 的相对重要度的长期变化规律和季节性特征.富营养化湖泊的藻类生长是自然界中一个非期望或非平均的现象,藻类生物量数据异质 性很强，水华期间的藻类数据会呈“高峰厚尾”的分布，或存在显著的异方差等情况.近年 来在环境科学和生态学领域受到重视的分位数回归(quan tile regression)方法特别适合处 理这种波动性大、异质性很强的环境数据.该方法可针对回归变量任何一个分位点进行回 归分析，且在存在极端值或重尾情况时仍能保持较好的稳健性，适宜处理应变量对自变量的 极端响应，而不只是平均水平的响应，因此能更加全面地反映藻类生物量对水温、N、P 等环境指标的响应特征.本研究基于云南洱海长时间尺度(1990-2013年)的水质观测数据， 运用分位数回归方法，按不同年份区间和不同季节分别分析洱海藻类生物量［以叶绿素 a(Chl-a)表征］对N、P、水温的定量响应关系，探讨营养物因子与水温因子相对重要性的 长时间尺度演变规律和季节性变化规律，对制定洱海富营养化控制策略提供科学依据1材料与方法1.1研究区域洱海是云南省第二大高原淡水湖泊，为滇西最大的断陷湖，跨洱源、大理两县市，处 于东经 100° 06’ -100° 17',北纬 25° 36’ -25° 55’ 之间.水面面积 249.80 km2，汇水 面积2 565.0 km2，最大水深21.0 m,平均水深10.5 m,库容28.8亿m3(图1).洱海是沿 湖人民生活、灌溉、工业用水的主要水源地，是整个流域社会经济可持续发展的基础［14］. 洱海湖面多年平均海拔1 965.8 m光照充足，辐射强，气温温和，为浮游藻类的大量繁殖 提供了有利条件.区年均气温15°C左右，年均降雨量1 055 mm,年均蒸发量1 970 mm.流 域水系发达，入湖河流大小共117条.图1洱海流域及水质监测点分布示意1.2数据来源从云南省环境监测中心站、中国大理洱海湖泊研究中心等环境监测和研究部门收集了 洱海1990-2013年的水环境常规监测数据，全湖水质监测点为13个(图1),经纬度依次为: 25° 51 36〃 N 100° 10’ 12〃 E,25° 51’ 18〃 N 100° 11’ 24〃 E,25° 51 00〃 N 100° 12’ 36〃 E,25° 41’ 42〃 N 100° 12’ 54〃 E,25° 42’ 00〃 N 100° 13’ 48〃 E,25° 42’ 18〃 N 100° 15’ 00〃 E,25° 37’ 48〃 N 100° 13’ 48〃 E,25° 36’ 0〃 N 100° 15’ 0〃 E,25° 37’ 60〃 N 100° 15’ 50〃 E,25° 55’ 30〃 N 100° 6’ 54〃 E,25° 54’ 54〃 N 100° 8’ 42〃E,25° 54’ 36〃 N 100° 10’ 48〃 E,25° 47’ 29〃 N 100° 11’ 43〃 E.水质指标主要包括总 氮(TN)、总磷(TP)、水温、藻类叶绿素a(Chl-a),采样频次为一年6次或12次.各水 质指标的主要测试方法如下：TN为过硫酸钾氧化紫外分光光度法，TP为钼锑抗分光光度法, 水温为温度计现场实时测定，Chl-a为四波段分光光度法[15].数据主要统计特征见表1, TP的变异系数最大(239.7%)，其次是Chl-a和TN,这可能与各指标的年份间差异和季节性 差异有关.各指标数据的Kolmogorov-Smirnov正态分布检验，除了水温数据呈正态分布以 外，所有指标的数据都呈非正态分布.水质指标单位 平均値最小值最大值中位数标准差变异浆数/品———概计P水温r175D.827174.525.89.0950.94TNmg-L'1D.41D.012.3D.40D.2253.29.7370TPmg-L'1D.D23D.0012.000D.020D.056239.7211770Chl-amgm-3S.20.1S5 35.SS.6103.96.0660TN/7P无里纲235SO211564.34.5340表1洱海水温、总氮、总磷及藻类叶绿素a的统计值(n=1 419)1.3分位数回归方法分位数回归模型(quan tile regression)是依据因变量的条件分位数对自变量进行回 归，得到所有分位水平的回归模型.假设随机变量的分布函数为F(y)=Prob(YWy),Y的T分位数为满足F(y)三T的最小 y值：Q(F) = P 書匚丁亡(° J)F(y)的T分位点Q(t )由最小化关于§的目标函数得到:Q(r) = arg min J | pT ( y其中P T(u)称为检验函数:T — f(uTU假设因变量Y和自变量X在T分位的线性函数关系为Y=XZ B +£ 条件分布函数为FY (y | x),贝叮分位数为：，给定X=x时，Y的Q(t\ x)=int y: A\(y I x) 3：丁｝ , t e(0,1)线性条件分位数通常表示为:I X) = Xfp(T)分位数回归能在不同的分位数T得到不同的分位函数，残差计算方法不同于最小二乘 法，具体如下：Residual = min \ I 兀 一 y\ I T1.4数据处理与分析考虑长时间尺度变化，将1990-2013年的监测数据按照相同年限(6 a)分成4组： 1990-1995年、1996-2001年、2002-2007年、2008-2013年，分别对这4组进行三元分 位数回归分析.为分析季节变化，将所有监测数据按春、夏、秋、冬四季分成4组，按不同季节分别进行三元分位数回归分析，比较各因子的斜率值变化，并计算95%置信区间 (CI), CI覆盖0(包含0),表示相关性“不显著”，而CI未覆盖0(不包含0),说明相关性 “显著”.为便于比较3个因子的斜率值变化，在进行多元分位数回归之前，所有原始数据 取对数(lg)处理.使用统计软件STATA V.12.0进行分位数回归分析和参数估计.2结果与讨论2.1水质指标的年际变化与季节变化2.1.1年际变化从TN变化趋势图看，2001年之前，TN浓度呈缓慢上升趋势，但总体处于较好的II类水 质.在2001-2003年期间有一次跳跃式增长，从II类区跨入III类区，近10年(2003-2013 年)TN浓度维持在III类水平，达到2001年之前浓度水平的近2倍[图2(a)]. TP浓度年均值 的变化趋势表现为有升有降，基本在II类水平线上下浮动，浓度高峰出现在2002-2004年区 间和2013年.近10余年(2002-2013年)的TP浓度水平比前10年(1990-2001年)的浓度水 平总体高出15%左右[图2(b)].洱海的年平均水温在16-19°C之间周期性波动，未见明显的 上升或下降趋势[图2(c)].iiiiuii "iitniu图2洱海TN、TP、水温和藻类Chl-a的逐年变化(1990-2013年)从藻类Chl-a变化趋势看，2000年之前洱海全湖Chl-a维持在低浓度水平，变化范围 为1.0-3.5 mg・m-3. 2000-2002的3年期间，Chl-a浓度呈直线上升，藻类生物量从低水 平跨越至较高水平.2003年之后的Chl-a平均浓度达到2000年之前的10余倍[图2(d)]. Pearson相关性分析显示藻类Chl-a与TN以及水温均呈显著正相关性(P<0.01)，而与TP的 相关性不显著(P〉0.05).2.1.2季节变化4项水质指标具有明显的季节性特征，均表现为夏、秋两季较高，且两季数据接近， 而春、冬季数值较低.8-10月TN浓度最高，5月是一年中TN上升的起点，11月是TN开 始下降的拐点[图3(a)]. 一年中TP浓度在4-7月是洱海上升速率最快的时段[图3(b)]. 水温在夏季和秋季最高，分别为23.3C和20.2C，很适宜藻类生长[图3(c)].藻类Chl-a 浓度在6-7月出现跳跃式上升，平均值由10 mg • m-3跃升至20 mg • m-3，在7-11月维持 在最高水平[图3(d)].Utt图3洱海TN、TP、水温和藻类Chl-a的季节变化2.2藻类对氮、磷及水温的长时间尺度响应特征通过三元分位数回归方法分析洱海N、P和水温对藻类Chl-a的耦合效应随时间变化特 征4个年份区间的lg(TN)、lg(TP)和lg(水温)的斜率值随分位点的变化谱图见图4,斜 率值及其95%置信区间的统计特征见表2.年饴区间斜率蜃数统计重 1990-1995 1996-2001 2002-2007 2003-2013平均值0.060.200.160.03变化范圉-0.09-0.16-0.12-0.39-0.00-0.35-0.52-0.4 -95%^信区间陛頗盖口0.32

平均值2 701.Q41.010.44切(水温)变化范国2.05-4.780.53-1.240.43-1.300.15-0.75眄滋置信区间陛床覆盖口陛;未覆盖0〔全;未覆盖口0.45