西苕溪流域主要经济林土壤磷素流失风险研究*

王肖君，王季丰，侯 琼，冷明珠，倪吾钟†

（1. 浙江大学环境与资源学院，浙江省农业资源与环境重点实验室，杭州310058；2. 浙江省安吉县农业农村局，浙江安吉313300）

水体富营养化是当今世界面临的最主要的水污染问题，严重影响了人们的正常生活，制约了区域社会经济的可持续发展。西苕溪作为太湖上游的重要支流，入湖水量大，其携入的氮磷等污染物均会对太湖水环境质量产生最直接的影响[1]。据统计，西苕溪流域林地面积占了流域总面积的65%以上[2]，每年可产生超过650 t 的总氮流失量及270 t 的总磷流失量。其中，经济林面积占了林地面积的 56%，土壤侵蚀及其产生的养分流失量占比77%[3]。

大量研究表明，由于化肥及有机肥的不合理施用，土壤磷素积累现象明显。以钱塘江附近农田为例，超过土壤磷素流失阈值60 mg·kg–1的土样占比可达11.60%[4]。当土壤中磷素积累到一定水平时，其释放潜力会急剧增加，土壤磷释放能力突变点所对应的土壤有效磷含量即土壤磷素流失阈值[5]。综合Hesketh 和Brookes[6]、Zhao[7]等的研究结果可以发现，不同类型的土壤磷素流失阈值差异很大，最大值甚至可达最小值的16 倍。自从提出土壤磷素流失阈值概念后，土壤磷素流失阈值的研究引起了国内外学者的广泛关注。但是大多数研究集中在农田土壤[8-11]上，对经济林土壤磷素流失阈值的研究相对较少。

相关数据显示，浙江省酸雨覆盖面积可达90%以上，其中重酸雨区约占20%，酸雨污染严重[12]。西苕溪流域位于浙江省西北部的湖州市境内，2012—2017 年年均降水pH 在4.60～4.98 之间，6年pH 中位数为4.75，酸雨率超过了90%。刘旭阳等[13]发现，酸雨显著提高了土壤磷的含量。大量研究表明，长期的酸雨作用会导致土壤磷素流失，且土壤磷流失总量与酸雨pH 存在显著的非线性关系。随着pH 的降低，土壤磷流失总量先增加后降低[14]。Liang 等[15]研究发现，pH1.6～6.0 范围内，pH 越低，胶体态磷流失风险就越高。

综上，土壤磷素流失阈值研究以及酸雨对土壤磷素流失的影响一直是学者们关注的重点。虽然也有少量土壤磷素流失阈值与土壤 pH 关系的研究[7]，但是关于酸雨对土壤磷素流失阈值的影响仍未见报道。因此，本文拟通过模拟酸雨浸提实验，根据土壤水溶性磷、模拟酸雨浸提磷与有效磷之间的相关关系，对西苕溪流域主要经济林土壤磷素流失风险进行研究，旨在为酸雨频发区土壤磷素流失阈值研究提供新思路，对流域内磷素的管理、农业面源污染及水体富营养化治理也具有十分重要的意义。

1 材料与方法

1.1 土样的采集与预处理

浙江省安吉县境内西苕溪流域面积为1 806 km2，约占西苕溪流域总面积的79.42%。根据安吉县主要经济林的分布情况，采集毛竹林、白茶园、板栗林土壤样品（0～20 cm）105 个，采样点分布如图1所示。采集的土样在室温下自然风干后，研磨过筛（20 目、100 目），保存备用。

1.2 测定项目与方法

水溶性磷（WEP）采用蒸馏水浸提，水土比为10∶1，振荡时间30 min，0.45 μm 滤膜过滤后用磷钼蓝比色法测定[16]；模拟酸雨浸提磷（SARP）测定方法参照水溶性磷，其中浸提液换成pH 为4.75 的模拟酸雨；土壤有效磷含量采用Bray 1 法，盐酸-氟化铵溶液提取，钼蓝比色法测定[17]，其他土壤理化性质如全氮（TN）、全磷（TP）、全钾（TK）、碱解氮（AN）、速效钾（AK）、pH、有机碳（SOC）等测定方法参照《土壤农业化学分析方法》[17]。

1.3 数据处理与分析

经济林土壤采样点的分布采用ArcMap 10.2 作图。试验数据用Excel 2010 软件整理、作图，用SPSS 20.0 做相应指标间的相关性分析。土壤活性磷与Bray 1-P 之间的定量关系利用分段线性回归模型拟合，拟合方程如下：

式中，y 为土壤活性磷，mg·kg–1；x 为Bray 1-P，mg·kg–1；a 为纵坐标上的截距，b、c 为拟合回归直线方程的斜率；α 为分段线性回归方程转折点所对应的Bray 1-P 含量，即阈值，mg·kg–1。

2 结 果

2.1 西苕溪流域经济林土壤基本理化性质与活性磷含量

调查区域安吉县105 个经济林土壤的基本理化性质汇总情况如表1 所示，土壤SOC 含量9.00～26.35 g·kg–1，平均为16.59±3.65 g·kg–1；土壤pH3.55～5.92，平均为4.69±0.45；TN、TP、TK 含量分别在0.72～2.07、0.22～0.73、3.57～22.80 g·kg–1之间，平均值分别为1.38±0.33、0.42±0.11、9.81± 3.56 g·kg–1；土壤AN、AK 含量分别在74.10～265.7、48.00～379.5 mg·kg–1之间，平均值分别为152.2±39.90、138.5±69.02 mg·kg–1。与其他指标相比，Bray 1-P 含量变幅较大，范围在0.93～313.2 mg·kg–1之间，平均为30.87±63.18 mg·kg–1，变异系数可达204.7%。统计数据表明，约17.14%的土壤Bray 1-P 含量超过40 mg·kg–1，但也有38.01%的土壤Bray 1-P 含量低于5 mg·kg–1。

土壤活性磷的提取分别采用蒸馏水和模拟酸雨浸提，结果如图2 所示：土壤水溶性磷（WEP）含量0.03～38.15 mg·kg–1，平均为2.64±7.06 mg·kg–1，变异系数为267.5%；模拟酸雨浸提磷（SARP）含量0.03～42.91 mg·kg–1，平均为2.86±7.69 mg·kg–1，变异系数为268.6%。SARP 含量总体高于WEP，但未达到显著性水平（P＞0.05）。

2.2 土壤磷素流失风险与土壤有效磷的关系

分段线性回归分析得出，土壤水溶性磷含量与有效磷含量的回归方程为：

统计检验结果表明，回归方程的显著性达到极显著水平（P＜0.01）。从图3 可以看出，拟合曲线上有一明显的转折点，相对应的Bray 1-P 含量为93.63 mg·kg–1。当Bray 1-P低于93.63 mg·kg–1时，WEP含量随着Bray 1-P 变化很小，拟合线段斜率为0.056 9；当Bray 1-P超过93.63 mg·kg–1，WEP 含量迅速增加，具体表现为：Bray 1-P 每增加10 mg·kg–1，WEP 增加1.48 mg·kg–1，拟合线段斜率增加了1.61 倍。因此，Bray 1-P 含量93.63 mg·kg–1可作为研究区域主要经济林土壤磷素流失阈值。统计数据表明，西苕溪流域主要经济林土壤Bray 1-P 超过阈值的样品占比12.38%。

表1 调查区域105 个经济林土壤基本理化性质Table1 Basic soil physico-chemical properties of the 105 sampling sites in the main non-wood forests under survey

注：SOC、TN、TP、TK、Bray 1-P、AN、AK 分别表示土壤有机碳、全氮、全磷、全钾、有效磷、碱解氮和速效钾。下同。Note：SOC，TN，TP，TK，Bray 1-P，AN and AK stands for soil organic carbon，total nitrogen，total phosphorus，total potassium，available phosphorus，alkalyzable nitrogen，readily available potassium，respectively. The same below.

图2 土壤水溶性磷和模拟酸雨浸提磷含量Fig. 2 Water-extractable phosphorus（WEP）and simulated-acidrain-extractable phosphorus（SARP）content in the soils

图3 土壤水溶性磷与有效磷的关系Fig. 3 Relationship between WEP and Bray 1-P

酸雨会影响土壤磷的溶出，土壤模拟酸雨浸提磷与有效磷含量的关系见图4，分段线性回归方程为：

显著性检验结果表明，回归方程达到极显著水平（P＜0.01）。SARP 含量普遍高于WEP，且其上限较WEP 上限高4.76 mg·kg–1。在模拟酸雨的作用下，拟合方程的两个斜率均有所增加。此外，SARP 所对应的土壤磷素流失阈值为87.68 mg·kg–1，较WEP所对应的阈值降低了5.95 mg·kg–1。

图4 土壤模拟酸雨浸提磷与有效磷的关系Fig. 4 Relationship between SARP and Bray 1-P

2.3 土壤测试磷与土壤基本理化性质的关系

本文对土壤全磷、有效磷、活性磷与土壤基本理化性质做了相关性分析，结果如表2 所示：土壤SOC、TN、AN、AK 与所有土壤测试磷指标均呈极显著正相关（P＜0.01）；土壤pH 与所有土壤测试磷指标呈极显著负相关；土壤TK 与Bray 1-P、WEP、SARP 呈极显著负相关，但决定系数R2较小。

表2 土壤全磷、有效磷、活性磷与土壤基本理化性质的关系Table2 Relationships of total phosphorus，available phosphorus and active phosphorus with basic physico-chemical properties of the soils

3 讨 论

3.1 西苕溪流域主要经济林土壤磷素积累状况

土壤全磷和有效磷作为评价土壤肥力、磷素水平的指标，常被用来反映土壤磷素水平及磷素积累情况。西苕溪流域经济林土壤TP 含量0.22～0.73 g·kg–1，平均为0.42 g·kg–1。Bray 1-P 含量变幅极大，范围在0.93～313.2 mg·kg–1之间，平均为30.87 mg·kg–1，变异系数可达204.7%（表1）。其中，Bray 1-P 含量超过40 mg·kg–1、低于5 mg·kg–1的土壤样品分别达到了17.14%、38.01%，且二者占比超过了55%。可见，本研究地经济林磷缺乏与磷过度积累土壤并存。这可能与山区林农的粗放式管理有关。通常认为，Bray 1-P 含量在20～40 mg·kg–1之间即可基本满足作物的需求，过多土壤磷的积累将导致土壤中磷素吸附位点减少，无法继续容纳外源磷素，大大增加土壤磷素的流失风险[19]，而Bray 1-P低于5 mg·kg–1又将影响作物的正常生长[17]。因此，有必要对该地经济林土壤磷素水平及其分布进行调查，根据相应的土壤有效磷含量确定适宜的施肥量，同时加强对土壤磷关键源区的识别[20]，严格控制磷肥施用。

3.2 土壤基本理化性质与土壤磷素流失风险的关系

水溶性磷（包括文中的模拟酸雨浸提磷）与径流水中磷的浓度密切相关，因此可用来分析土壤磷素随径流及酸雨径流流失的风险。大量研究表明，土壤磷素的流失受土壤有效磷、全磷、pH、有机质等诸多因素的影响。其中，土壤有效磷对土壤磷素流失影响最大，可作为土壤磷的环境指标[21]。这也解释了本文中Bray 1-P 的决定系数R2远高于TP、pH、SOC 等因素的结果（表2）。本研究还发现，随着土壤pH 下降，土壤Bray 1-P、WEP、SARP 在统计学水平上有显著的增加（表2），与黄敏等[22]的研究结果一致。土壤TN、AN、AK 与WEP、SARP均呈极显著正相关（表2），这可能与农户普遍施用NPK 复合肥有关。土壤SOC 与WEP、SARP 以及Bray 1-P 的相关关系均可达到极显著水平（表2），与 Huang 等[23]的研究结果相一致。Heredia 和Cirelli[24]及Meng 等[25]认为，土壤有机质易与磷酸根离子竞争土壤固相表面的专性吸附点位，从而减少了土壤对磷的吸附，且土壤有机质的存在有助于土壤磷向土壤溶液中释放。夏文建等[26]和Yan 等[27]的研究发现，土壤有机质含量高的土壤磷吸持饱和度（DPS）较大，土壤吸附磷的能力较弱。因此，土壤有机质含量的增加将导致土壤吸附磷的能力下降，土壤磷更容易向液相中转移，从而使土壤磷素流失阈值降低。Fortune 等[28]研究发现，土壤磷素流失阈值与土壤有机碳呈显著负相关，与本文推断结果一致。但是也有学者认为，有机质的存在可以增加土壤团聚体的稳定性[29]。团聚体作为土壤磷素释放与流失的重要载体，对土壤磷素流失具有重大影响，且土壤磷素流失阈值随着土壤有机碳的增加而增大[7]。关于有机质与土壤磷素流失风险的关系有待于进一步的研究。

3.3 模拟酸雨对经济林土壤磷素流失风险的影响

西苕溪流域年均降水pH 在4.60～4.98 之间，酸雨率超过 90%，处于较重酸雨区与重酸雨区之间[30]，酸雨污染严重。袁宇志等[31]研究表明，相对于其他类型耕地土壤而言，红壤旱地更易受到酸雨的影响，导致土壤酸化。同时磷素作为水体富营养化的关键限制性因子，酸雨对土壤磷素流失的影响不容忽视。土壤磷素流失阈值的研究一直是国内外学者研究的热点。阈值所对应的有效磷Olsen-P、Mehlich-3 P 含量得到了广泛研究，相比之下，Bray-P在土壤磷素流失风险研究中的应用较少。仅有零星报道显示，美国俄亥俄州农业土壤磷素淋失阈值为Bray 1-P 122 mg·kg–1[32]。本研究中，WEP、SARP与Bray 1-P 的相关分析结果显示，西苕溪流域经济林土壤磷素流失阈值分别为93.63、87.68 mg·kg–1（图3，图4）。当土壤Bray 1-P 含量超过土壤磷素流失阈值时，WEP 及SARP 急剧增加，土壤磷素流失风险加大。关于酸雨对土壤磷素流失的影响主要体现在以下两个方面。一方面，一定强度的酸雨有助于土壤无机磷组分的溶解，可以促进土壤微生物生物量磷和酸性磷酸酶活性提高，从而提高土壤磷的有效性[33-34]。本研究也发现，土壤SARP 含量普遍高于WEP，且分段线性拟合曲线的两段斜率均有了一定程度的提高（图3，图4）。可见pH4.75 的模拟酸雨有助于土壤磷向土壤溶液中释放。另一方面，酸雨的长期作用将导致土壤磷素流失阈值降低。土壤磷素流失阈值越低，土壤磷素流失风险就越大。Zhao 等[7]的研究结果显示，当土壤pH＜6.0 时，土壤磷素流失阈值随着 pH 的升高而升高；而当土壤pH＞6.0 时，随着pH 的升高土壤磷素流失阈值将降低。本研究发现，模拟酸雨作用下土壤磷素流失阈值下降了5.95 mg·kg–1。这表明，相同背景下将会有更多经济林土壤处于土壤磷高流失风险范畴，土壤磷素的潜在流失风险将进一步加大。而常规的根据水溶性磷与土壤有效磷的相关关系得出的磷素流失阈值可能会低估酸雨频发区土壤磷素流失的风险。

4 结 论

西苕溪流域安吉段经济林土壤磷素状况变异较大，Bray 1-P 含量超过40 mg·kg–1、低于5 mg·kg–1的土壤样品占比分别可达17.14%、38.01%，缺磷和磷过度积累现象并存。土壤磷素积累导致水溶性磷、模拟酸雨浸提磷大幅度升高，存在较大的磷素流失潜在风险。影响土壤磷素流失的因子主要有土壤有效磷、全磷、pH、有机质等，其中有效磷是最重要的影响因子。水溶性磷、模拟酸雨浸提磷与土壤Bray 1-P 的分段线性回归分析得出，西苕溪流域经济林土壤磷素流失阈值分别为93.63 mg·kg–1（水溶性磷）、87.68 mg·kg–1（模拟酸雨浸提磷），酸雨会进一步增大磷素的流失风险。鉴于西苕溪流域为酸雨频发区，在土壤有效磷水平的实际限控中应充分考虑酸雨增大磷素流失风险的作用。