消落带是河流、湖泊、水库特有的一种现象，主要由于湖库河流水位的季节性消涨和周期性蓄水，形成的分布在最高水位线与最低水位线之间的区域^[1].消落带是衔接水-陆生态系统的生态交错带，同时也是水陆生态系统间物质能量转换最活跃、最重要的区域^[2-3].消落带由于阶段性的淹没-落干过程，使其在营养盐的分布和贮存方面与长期处于淹水状态的沉积物存在明显差异，消落带的淹水-落干的频率和时间对消落带营养盐的形态转化与水界面的交换过程有重要影响.周期性的淹水对消落带的影响主要表现，导致土壤中的氮磷营养盐向水中转移，造成水体的富营养化，影响母岩养分的释放^[4-5].有研究表明，消落带是库区流域富营养化的主要来源^[6].因此研究消落带营养盐的储量和分布情况，对进一步管理和防护水体富营养化有显著作用和意义^[7].

目前，国内外针对消落带的研究侧重各有不同，国外的研究对象大多针对河流^[8]、湖泊^[9-11]以及湿地^[12-13]自然水位的涨落或者强降雨事件下消落带与水界面进行的物质交换及营养盐的形态转化，而国内对消落带的研究则集中在三峡地区^[14-17]，由人为对水资源的调控造成的水位落差形成的消落带.

然而，密云水库作为北京市重要的饮用水源地，同样值得关注.密云水库建于1960年，水库库容40×10⁸ m³，多年蓄水量为6.5×10⁸~12×10⁸ m³，历年水位在131~137 m之间.由于受到降雨和人为调控水量的影响，其水位呈现出季节性的波动^[18]，每年都会形成1~4 m左右的消落带.除此之外，考虑到“南水北调”工程作为一项浩大的国家性工程，为流经城市带来珍贵的水资源之外，可能也会对受水区的生态环境带来影响.由于密云水库是“南水”调入北京的储存库，一方面水库的水源得到补给，另一方面由于储量的增加，势必在原有基础上水位有所抬升，最初的消落带可能会成为永久的淹没带，库岸带部分也会有新的消落带形成.这对密云水库的水质安全可能造成一定的风险.因此，本研究以密云水库内湖为研究对象，考虑到“南水北调”来水后，不同时间段对水量有不同调度要求，蓄水量可能存在阶段性变化，将蓄水量分12×10⁸ m³(对应水位137 m)、17×10⁸ m³(对应水位142 m)、22.3×10⁸ m³(对应水位146 m)和30×10⁸ m³(对应水位152 m)4个阶段，研究分析4个不同高程下(133~137、137~142、142~146、146~152 m)消落带碳、氮、磷的分布特征，为评估消落带氮、磷入库风险负荷量，维护密云水库的水质安全提供依据.

1 材料与方法 1.1 采样点设置与样品采集

消落带样品于2012年3月进行采集，采样点主要分布在内湖区域.采集不同高程梯度(133~137、137~142、142~146、146~152 m)、土层深度(0~20 cm，现场5 cm分层)和土地利用类型(农田、草地、林地)，共54个采样点(图 2)，112个消落带样品.采集的样品立即放入封口袋，4℃保存带回实验室，用于各指标分析.

1.2 样品分析方法

消落带样品总磷和有机磷分析方法采用欧洲标准测试委员会框架下发展的SMT(the Standards Measurements and Testing Program of the European Commission)法，即称取0.2 g样品，加20 ml 1 mol/L HCl，振荡16 h后提取无机磷；称取0.2 g样品，于450℃马弗炉中灰化3 h，加20 ml 3.5 mol/L HCl，振荡16 h后提取总磷，总磷和无机磷提取液采用钼锑抗分光法进行测定^[19-20].全氮含量采用凯氏定氮法测定，土壤有机质含量采用重铬酸钾氧化-外加热法测定^[21]，氨氮采用KCl浸提-蒸馏法测定，硝态氮采用酚二磺酸比色法测定^[22]，土壤含水率采用铝盒烘干差减法测定，容重采用环刀法进行^[23]测定.

1.3 数据处理和统计方法

研究数据中Duncan多重比较以及相关性分析均用SPSS 16.0软件进行统计分析，图表由Excel和Origin软件进行绘制.消落带有机质、总氮和总磷含量计算公式为：

$ F = \sum\limits_{i = 1}^n {\frac{{Co{n_{ik}} \cdot {\rho _{ik}}}}{{100}} \cdot Are{a_{ik}}\left( {1 - {\varphi _{ik}}/100} \right) \cdot 10 \cdot L} $

(1)

式中，F为研究区域内消落带有机质、营养盐(碳、氮、磷)储量(t)，n为土层划分层数，L为常数，5 cm，代表土层厚度(cm)，Con_ik、ρ_ik、Area_ik和φ_ik分别代表第i层第k种土地利用类型的土壤全氮、全磷和有机质含量(g/kg)、土壤容重(g/cm³)、面积(m²)和含水率(%)，10为单位转换系数.

2 结果与讨论 2.1 133~146 m高程范围内消落带有机质和营养盐含量

密云水库内湖133~146 m高程范围内农田和林地消落带(0~20 cm均值)有机质、总磷、总氮、氨氮、硝态氮、亚硝态氮、无机磷含量分别为23.15±13.65 mg/g、0.17±0.09 mg/g、1.44±0.81 mg/g、10.86±3.54 μg/g、8.07±2.73 μg/g、0.41±0.71 μg/g、9.09±4.18 μg/g.土壤湿容重和含水率分别为0.94±0.08 g/cm³和12.54 % ±3.49 % (表 1).

2.2 消落带有机质和营养盐分布特征影响因素分析 2.2.1 高程对消落带(0~20 cm)有机质和营养盐含量的影响

同一高程，不同土地利用类型各营养盐含量有明显变化(表 1).林地的有机质、总氮和硝态氮含量较高，分别是均值的1.66、1.54、1.22倍.而农田的营养盐含量皆处于较低的水平，或与均值持平(无机磷、氨氮、硝态氮)或低于均值水平(总磷、有机质、总氮为均值的0.7倍).因此分析高程对营养盐分布的影响，分别就同一土地利用类型进行讨论.对林地而言，有机质和硝态氮受高程影响，主要表现在133~137 m和142~146 m有机质与硝态氮含量差异显著(P < 0.05)，而137~142 m与133~137 m和142~146 m有机质和硝态氮含量差异不显著(P>0.05).对农田而言，133~137与137~142 m高程范围为总氮和氨氮含量无显著差异(P>0.05)，而与142~146 m高程范围内总氮和氨氮含量具有显著差异(P < 0.05).分析主要原因可能是目前密云水库的水位在132~138 m范围之间变动，对此范围的消落带影响更加频繁，因此133~137和137~138 m范围内消落带的总氮、氨氮无显著差异.研究表明，一定的淹没不利于氮磷元素的累积，同时也会加速氮的矿化^[24-25]，因此在高程梯度上，总氮、氨氮、硝态氮含量总体表现是随着高程的增加而增加，这与水位的消长、氮元素的累积和矿化产生的影响不无关系.而高程对总磷和磷酸盐含量无显著影响(P>0.05).王业春等^[26]对三峡消落带的研究同样得到类似结果，这可能是因为不同高程的土壤母质具有相似性.

2.2.2 土地利用类型对消落带(0~20 cm)有机质和营养盐含量的影响

土地利用类型对各营养盐含量的影响在不同高程情况下有所不同(表 2).在133~137、137~142 m高程范围内，土地利用类型对有机质、总氮和氨氮含量有显著影响(P < 0.05)，在142~146 m高程范围内，土地利用类型对有机质、总氮和硝态氮含量有显著影响(P < 0.05).在各高程范围内土地利用类型对总磷、磷酸盐含量影响不显著(P>0.05).有相关研究表明，不同土地植被覆盖下，土壤碳氮含量存在很大差异^[27-29]，这可能是因为林地的枯落物丰富，且根系较多，交错复杂，其土壤具有良好透性，能够促进土壤微生物活动，对土壤氮素循环和土壤矿化过程有积极的促进作用^[30]，所以林地总氮和有机质含量比农田高(表 1).而农田和林地的总磷和磷酸盐含量差别较小，表明土地利用类型不是影响密云水库消落带总磷和磷酸盐含量的主要因素.

2.2.3 土壤深度对消落带总氮、总磷、有机质的影响

有机质、总磷和总氮含量随着深度的增加而递减(图 4).在133~137、137~142、142~146 m高程范围内，有机质随土层深度增加而降低，总磷、总氮含量变化趋势相似，先增高后降低，在5~10 cm土壤深度达最大值.说明有机质、总磷、总氮含量分布有较强的表聚性，有随土层增加而降低的趋势^[31-32].此种现象可能与地表植被根系吸收土壤营养物质进行生长，再以枯落物的形式归还土壤，营养元素的循环过程有关^[33].有机质、总磷和总氮含量在各土壤深度的变化范围也有明显趋势，其中0~5和5~10 cm的变异系数较大，在10~15、15~20 cm土壤深度变异系数较小，表明在表层0~10 cm土壤各营养物交换较为频繁，可能受地表覆盖物、表层水分及光照的影响较大.

2.3 消落带土壤基本理化特征及其与各营养盐间的相互关系

有机质与总氮、氨氮含量呈极显著正相关(P < 0.01)(表 3)，表明有机质与总氮、氨氮含量存在相似的消长关系，这也与之前的研究结论一致^[34-35].总磷与磷酸盐、总磷与亚硝态氮含量呈显著正相关(P < 0.05)，这与孙斌^[36]得到的结果一致.硝态氮与含水率、磷酸盐与湿容重呈显著正相关(P < 0.05)，由于土壤含水率与土壤容重的比值是常数，表明土壤水分可能是影响土壤硝态氮和磷酸盐的因素.因土壤水分影响微生物的活性，从而影响氮磷的形态转化.

2.4 消落带各营养盐储量估算 2.4.1 不同土壤深度各营养盐储量分布特征

内湖消落带0~5、5~10、10~15、15~20 cm土壤深度有机质储量分别为2630.27、1610.58、1268.07、974.57 t, 储量由表层到深层有逐渐降低的趋势；总磷储量分别为18.75、18.40、11.92、10.47 t，表层10 cm磷储量差别不大，随后有缓慢降低的趋势；总氮储量分别为142.20、117.38、89.73、64.70 t，储量也随着土层深度增加而降低.

2.4.2 不同高程和土地利用类型各营养盐储量分布特征

在133~137和137~142 m高程范围内，林地和农田中有机质、总磷和总氮储量有显著变化.其中，在133~137 m高程内，农田的营养盐储量显著高于林地，有机质、总磷和总氮储量分别为林地储量的4.8、10.38和9.38倍.而在137~142 m高程范围内，林地有机质、总磷和总氮储量分别为农田的1.7、2.8和4.1倍. 142~146 m高程范围内，林地与农田各营养盐含量差别不大(表 5).

3 结论

1) 密云水库内湖消落带有机质、总磷、总氮、氨氮、硝态氮、亚硝态氮和无机磷含量分别为23.15±13.65 mg/g、0.17±0.09 mg/g、1.44±0.81 mg/g、10.86±3.54 μg/g、8.07±2.73 μg/g、0.41±0.71 μg/g和9.09±4.18 μg/g.

2) 土地利用类型和高程是影响消落带有机质和氮含量的重要因素，而对于磷的作用不显著.

3) 内湖消落带在133~137 m高程范围内有机质、总磷和总氮储量分别为566.32、20.93和100.97 t；137~142 m高程范围内有机质、总磷和总氮储量分别为2407.62、14.00和116.03 t，有机质储量与133~137 m高程范围内差别较大，总磷和总氮储量相近.在133~146 m高程范围内有机质、总磷和总氮的储量分别为5324.07、59.56和414.02 t.

4) 因此，需要加强对133~142和142~146 m高程范围的消落带进行防护，以防止水位的波动使大量有机质和总氮进入水体.

5) 密云水库内湖消落带是巨大的有机质、氮磷等元素的储藏库，阶段性的淹水-落干过程对水库水质有重要的影响，需要进一步研究.