我国耕地水土流失面积达4.54×10⁷ hm²，占耕地总面积的35%，每年因土壤侵蚀而造成的水土流失中磷素损失总量达7.65×10⁶ t^[1-2]。磷素随径流进入水体，使得磷素养分流失成为最大的农业面源污染之一^[3]，同时也是造成水体富营养化的主要限制因子^[4-5]。作物秸秆占其生物总量的50%左右，是一种极为丰富且可直接利用的可再生有机资源^[6-7]。秸秆还田可减少因侵蚀、径流或淋溶而造成的磷素损失^[8-9]，改善土壤养分状况，改良土壤理化性质，对作物的产量及品质有较好的促进作用^[10]；同时可避免秸秆焚烧造成的资源浪费和环境污染问题^[11]。因此，选择适宜的秸秆堆沤还田方式及对施肥量进行控制，可在较大程度上改善土壤养分状况，有效控制农业面源污染。

目前，有关秸秆还田对磷素的影响已有较多研究，主要集中在秸秆还田、秸秆还田与不同覆盖方式、减量化施肥与秸秆还田对径流泥沙和磷素流失的影响。王静等^[12]研究表明，与传统耕作方式相比，稻季秸秆还田使磷素流失分别减少了5.13%和18.98%；侯红波等^[13]研究表明，洞庭湖区红壤坡耕地常规施肥+秸秆覆盖(3 000 kg·hm⁻²)处理后，磷素总流失量较CK降低了15.4%；王志荣等^[14]研究表明，减量化施肥与秸秆还田施肥可降低地表径流磷素流失。但针对堆沤方式与秸秆还田密度对坡耕地产流产沙和磷素流失的影响研究较少。因此，本试验在自然降水条件下，设置0.75、1.5 kg·m⁻²秸秆密度，1、5 cm秸秆粒度，水、水与尿素堆沤方式，研究其对滇中二龙潭小流域烟草坡耕地径流泥沙及磷素流失的影响，拟分析不同秸秆堆沤方式与秸秆还田密度结合下径流泥沙及磷素流失特征，以期得到该流域烟草坡耕地秸秆资源最佳的利用方式以及为控制农业面源污染提供科学的理论依据。

1.   材料与方法

1.1   研究区概况

研究区位于云南省玉溪市红塔区二龙潭小流域，地理坐标：E102°34′12.30″、N24°17′32.33″，地处玉溪市中心城区东南方向，流域控制面积为53 km²。雨季为5—10月，雨季降水量占全年降水量的85%，其中降水在8月达到峰值，年平均暴雨天数为15 d。年平均降水量为909.1 mm，最大年降水量为1 413.7 mm，年均蒸发量为1 624.9 mm。该流域地貌属波状起伏线切割中山区，土壤属山地红壤，肥力较低，介于强酸至微酸性之间。土层浅薄，岩石裸露，属水土流失高发区域，流域内主要树种为云南松，烟草和玉米为当地主要农作物。

1.2   试验设计

本试验地为15°烟草坡耕地，2019年5月7日进行烟草移植，烟草种植密度为每公顷16 500株。氮肥施用量为105~135 kg·hm⁻²，N、P₂O₅、K₂O质量比为1∶0.5∶2.5~3；肥料在烟草移栽后25 d内全部施完，共施肥3次：底肥、提苗肥、追肥。肥料施用量如表1所示。

表  1  肥料施用情况

Table  1.  Fertilizer details

种类 Species	N∶P∶K1)	施用量/(kg·hm−2) Fertilization dosage
		总量 Total	N	P
底肥 Base fertilizer	12∶6∶24	247.5	29.7	14.85
提苗肥 Seedling fertilizer	28∶0∶5	45.0	12.6	0
追肥 Additional fertilizer	12∶6∶24	772.5	92.7	46.35
1)复合肥中的氮磷钾质量比 1)The mass ratio of nitrogen, phosphorus and potassium in compound fertilizer

下载: 导出CSV 

| 显示表格

小区规格为1 m×1 m，坡度、坡向基本保持一致，呈2行对称排列，试验田外围设有试验保护区，小区田埂筑高20 cm，小区周围用塑料薄膜包被，以减少侧渗。每个小区外设置1个集流装置，用10 L的大塑料桶埋入土中收集，塑料桶上盖有盖子，以防止雨水进入塑料桶，每次取样后，测量集流装置中的雨水量并将水清除干净，方便后续观察。用自动雨量计观测试验期间的降水量与降水强度。

本试验将当地种植的玉米秸秆堆沤处理后施用于烟草坡耕地。试验设置9个处理，空白对照即不添加秸秆(T₁)、0.75 kg·m⁻²秸秆密度+ 5 cm秸秆粒度+水堆沤(T₂)、0.75 kg·m⁻²秸秆密度+ 5 cm秸秆粒度+水与尿素堆沤(T₃)、0.75 kg·m⁻²秸秆密度+ 1 cm秸秆粒度+水堆沤(T₄)、0.75 kg·m⁻²秸秆密度+ 1 cm秸秆粒度+水与尿素堆沤(T₅)、1.5 kg·m⁻²秸秆密度+ 5 cm秸秆粒度+水堆沤(T₆)、1.5 kg·m⁻²秸秆密度+ 5 cm秸秆粒度+水与尿素堆沤(T₇)、1.5 kg·m⁻²秸秆密度+ 1 cm秸秆粒度+水堆沤(T₈)、1.5 kg·m⁻²秸秆密度+ 1 cm秸秆粒度+水与尿素堆沤(T₉)。秸秆还田后及时翻耕，把秸秆翻入0~20 cm土层，使秸秆分散均匀且与土壤充分混合。不同处理见表2。

表  2  各样方处理方式及施肥量

Table  2.  Treatment and fertilization dosage of different plots

处理 Treatment	施肥量/(g·m−2) Fertilization dosage					秸秆还田量/(g·m−2) Straw-returned amount		秸秆粉碎 程度/cm Degree of straw crushing	秸秆堆沤添加 尿素量/(g·m−2) Straw compost urea addition amount
	基肥 Base fertilizer	提苗肥 Seedling fertilizer	追肥 Additional fertilizer	折合磷 Pure phosphorus		还田量 Straw-returned amount	折合磷 Pure phosphorus		施用量 Fertilization dosage	折合磷 Pure phosphorus
T1	70	13	214	17.04		0	0	0	0	0
T2	70	13	214	17.04		750	2.025	5	0	0
T3	70	13	214	17.04		750	2.025	5	3.75	0.01
T4	70	13	214	17.04		750	2.025	1	0	0
T5	70	13	214	17.04		750	2.025	1	3.75	0.01
T6	70	13	214	17.04		1500	4.050	5	0	0
T7	70	13	214	17.04		1500	4.050	5	3.75	0.01
T8	70	13	214	17.04		1500	4.050	1	0	0
T9	70	13	214	17.04		1500	4.050	1	3.75	0.01

下载: 导出CSV 

| 显示表格

1.3   样品采集与测定

用对角线法采集0~20 cm耕层土壤，取上、中、下坡0~5 cm、5~10 cm、10~20 cm 3个土层的土样，调查土壤养分背景值。根据每场降水历时，在每次降水结束后，分别采集1次收集装置中的水样，并根据降水次数增加适当增加采集频率，每个重复取200 mL，再将这些样品混合后取250 mL作为本次降水采样点测定水样，4 ℃条件下保存，在24 h内过滤，测定水样中的总磷、磷酸(PO₄⁻-P)质量浓度。取完径流水样后，将其静置、晾干，测定泥沙中总磷含量。体积法测定径流量，711型便携式悬浮物/界面分析仪测定泥沙含量。水样总磷用过硫酸钾氧化–钼锑抗比色法(GB 11893—1989)^[15]测定；PO₄⁻-P用钼锑抗比色法(NY/T 2421—2013)^[16]测定；泥沙和土样总磷用自动定氮仪(NY/T 1121.24—2012)^[17]测定。各个指标均做3组平行试验，取平均值。

用自动雨量计观测降水量与降水强度，本试验降水集中在6月至9月，在记录的降水中，选取了4场典型的降水事件，即7月28日、8月7日、8月15日和9月8日，降水量分别为6.2、22.0、12.4、17.2 mm，降水强度分别为2.6、23.6、3.1、21.0 mm·h⁻¹。其中最小降水量与最小降水强度发生在7月28日(6.2 mm、2.6 mm·h⁻¹)，8月7日的降水量与降水强度达到峰值(22.0 mm、23.6 mm·h⁻¹)。

1.4   数据处理与分析

在本试验中，通过烘干法确定各小区产沙量，各个小区的平均径流流量由径流收集器中水样的质量转换得出。径流和泥沙中磷素流失量由以下公式计算得到：

数据整理和图表制作采用Microsoft Excel 2010进行，相关性分析采用SPSS 23.0软件进行，方差分析及显著性检验采用LSD法。

2.   结果与分析

2.1   不同堆沤方式与秸秆还田密度对坡耕地产流产沙的影响

从产流产沙量(图1A、1B)来看，7月28日的径流量与泥沙量均为各次降水的最小值，8月7日的径流量与泥沙量多达到各次降雨的最大值，8月7日各样方平均土壤侵蚀量为5.51 g·m⁻²，7月28日则只有2.75 g·m⁻²，说明产流量与产沙量随着降水强度的增大而增加，且降水量与径流量、降水强度与产沙量相关性显著(P<0.05)。

图  1  4场降水中不同处理产流、产沙特征

各小图中柱子上不同的小写字母表示同一降水不同处理之间差异显著，不同的大写字母表示同一处理不同降水之间差异显著(P<0.05，LSD法)

Figure  1.  Runoff and sediment characteristics under four rainfalls in different treatments

Different lowercase letters on the columns in each figure indicate significant differences among different treatments under the same rainfall, different capital letters indicate significant differences among different rainfalls under the same treatment (P<0.05, LSD method)

下载: 全尺寸图片 幻灯片

从图1A可以看出，在４场降水中，T₂~T₉样方产流量较T₁平均减少了15.15%~40.58%；不同堆沤方式下各样方产流量不同：施用1.5 kg/m²秸秆密度(T₆、T₇、T₈、T₉)，5 cm粗秸秆粒度(T₂、T₃、T₆、T₇)，加水堆沤(T₂、T₄、T₆、T₈)分别较施用0.75 kg/m²秸秆密度(T₂、T₃、T₄、T₅)，1 cm细秸秆粒度(T₄、T₅、T₈、T₉)，加水与尿素堆沤(T₃、T₅、T₇、T₉)产流量减少了2.01%~30.38%，5.04%~43.26%，0.67%~37.13%，说明施用1.5 kg·m⁻²秸秆还田密度、5 cm秸秆粒度、加水堆沤均可较好地减少坡面产流量。从产沙量(图1B)来看，T₂~T₉样方产沙量较T₁平均减少了0.97%~19.68%；施用0.75 kg/m²秸秆密度(T₂、T₃、T₄、T₅)，5 cm粗秸秆粒度(T₂、T₃、T₆、T₇)，加水堆沤样方(T₂、T₄、T₆、T₈)分别较1.5 kg·m⁻²秸秆密度(T₆、T₇、T₈、T₉)，1 cm细秸秆粒度(T₄、T₅、T₈、T₉)，加水与尿素堆沤(T₃、T₅、T₇、T₉)产沙量减少了0.22%~20.82%，0.58%~24.64%，0.63%~36.10%。说明施用0.75 kg·m⁻²秸秆密度、5 cm秸秆粒度、加水堆沤可有效地减少坡面产沙量。

2.2   不同堆沤方式与秸秆还田密度对坡耕地地表径流磷素质量浓度的影响

不同处理磷素流失量见图2，4场降水中，各取样时期的磷素浓度变化差异较显著，4次降水径流中总磷浓度呈现出8月7日>8月15日>9月8日>7月28日的规律，8月7日径流总磷浓度达到峰值，且与其他降水差异显著(P<0.05)。相同降水条件下，处理T₂~T₉总磷流失浓度较T₁高17.13%~71.15%；不同堆沤方式径流总磷流失浓度呈现出施用0.75 kg·m⁻²秸秆密度(T₂、T₃、T_4、T₅)，5 cm秸秆粒度(T₂、T₃、T₆、T₇)，加水堆沤(T₂、T₄、T₆、T₈)分别较1.5 kg·m⁻²秸秆密度(T₆、T₇、T₈、T₉)，1 cm秸秆粒度(T₄、T₅、T₈、T₉)，加水与尿素堆沤(T₃、T₅、T₇、T₉)降低了7.79%~66.67%，3.33%~66.51%，2.82%~61.92%(图2A)。处理T₂~T₉ PO₄⁻-P流失浓度较T₁高16.29%~71.14%；不同堆沤方式径流PO₄⁻-P浓度变化呈现出施用0.75 kg·m⁻²秸秆密度(T₂、T₃、T_4、T₅)，5 cm秸秆粒度(T₂、T₃、T₆、T₇)，加水堆沤(T₂、T₄、T₆、T₈)分别较1.5 kg·m⁻²秸秆密度(T₆、T₇、T₈、T₉)，1 cm秸秆粒度(T₄、T₅、T₈、T₉)，加水与尿素堆沤(T₃、T₅、T₇、T₉)降低了1.97%~64.47%，3.23%~64.36%，3.39%~64.73%(图2B)。PO₄⁻-P/总磷变化幅度为35.30%~77.59%，9月8日占7月28日的91%，在施肥后期溶解态磷素流失的风险增加，但总体上呈现随时间的推移而逐渐减小的趋势(图2C)。

图  2  4场降水中不同处理对径流总磷、PO₄⁻-P流失量及其比值与泥沙总磷流失量的影响

各小图中柱子上不同的小写字母表示同一降水不同处理之间差异显著，不同的大写字母表示同一处理不同降水之间差异显著(P<0.05，LSD法)

Figure  2.  Effect of different treatments on the loss of total phosphorus、PO₄⁻-P and loss ratio of PO₄⁻-P to total phosphorus in runoff and total phosphorus loss in sediment under four rainfours

Different lowercase letters on the columns in each figure indicate significant differences among different treatments under the same rainfall, different capital letters indicate significant differences among different rainfalls under the same treatment (P<0.05, LSD method)

下载: 全尺寸图片 幻灯片

2.3   不同堆沤方式与秸秆还田密度对坡耕地泥沙中磷素含量的影响

由图2D得出，在4场典型降水中，泥沙总磷浓度呈现出8月7日>8月15日>9月8日>7月28日的规律，且差异显著(P<0.05)。处理T₂、T₃、T₄、T₆泥沙总磷浓度较T₁平均减少38.61%、12.11%、24.68%、20.99%，处理T₅、T₇、T₈、T₉泥沙总磷浓度较T₁平均增加1.13%、3.02%、19.02%、26.92%；不同堆沤方式泥沙总磷流失浓度表现为施用0.75 kg·m⁻²秸秆密度(T₂、T₃、T_4、T₅)，5 cm秸秆粒度(T₂、T₃、T₆、T₇)，加水堆沤(T₂、T₄、T₆、T₈)分别较1.5 kg·m⁻²秸秆密度(T₆、T₇、T₈、T₉)，1 cm秸秆粒度(T₄、T₅、T₈、T₉)，加水与尿素堆沤(T₃、T₅、T₇、T₉)降低了8.38%~57.53%，0.54%~55.84%，0.38%~55.55%。与8月7日相比，8月15日、9月8日泥沙总磷流失浓度分别降低了11.27%~42.97%和9.33%~39.26%，说明随着烟草生长时间的推移，各处理泥沙总磷流失浓度呈现逐渐降低的趋势。

2.4   不同堆沤方式与秸秆还田密度对坡耕地磷素流失量的特征

4次降水中，不同堆沤方式和秸秆还田密度下径流总磷流失量分别为1.97~5.35、3.51~11.24、2.25~5.74和2.54~8.52 mg·m⁻²，泥沙总磷流失量分别为0.61~4.45、1.47~4.16、0.91~2.16和1.05~2.70 mg·m⁻²。不同处理径流和泥沙总磷流失量表现为8月7日>8月15日>9月8日>7月28日的规律。降水量增加，总磷流失量也增加，变化显著(P<0.05)。与8月7日相比，8月15日和9月8日总磷流失量分别下降了43.81%和22.10%，且随着时间的推移逐渐下降。不同堆沤方式、秸秆还田密度下降低了磷素流失量。除T₂外，其余处理径流总磷流失量较T₁均增加16.17%~62.17%，除T₈、T₉外，其余处理泥沙总磷流失量较T₁均减少8.99%~38.44%；同样的降水条件下，各处理径流和泥沙总磷流失量表现出施用0.75 kg·m⁻²秸秆密度(T₂、T₃、T_4、T₅)，5 cm秸秆粒度(T₂、T₃、T₆、T₇)，加水堆沤(T₂、T₄、T₆、T₈)分别较1.5 kg·m^-2秸秆密度(T₆、T₇、T₈、T₉)，1 cm秸秆粒度(T₄、T₅、T₈、T₉)，加水与尿素堆沤(T₃、T₅、T₇、T₉)最高降低了63.64%和64.74%，63.89%和57.87%，63.89%和64.74%。说明随秸秆密度降低、粉碎粒度变粗以及仅加水堆沤径流和泥沙总磷流失量呈下降趋势。在4次典型降水中，径流输出是坡耕地总磷流失的主要方式。占总磷总流失量的68.78%~74.44%。不同处理径流和泥沙总磷流失量见图3。

图  3  4场降水中不同处理对坡耕地径流和泥沙总磷流失量的影响

Figure  3.  Effect of different treatments on total phosphorus loss in runoff and sediment of sloping farmland under four rainfours

下载: 全尺寸图片 幻灯片

3.   讨论与结论

3.1   耕地径流泥沙及磷素流失形态的变化特征

降雨量、降雨强度、施肥方式、施肥量、土壤性质等是磷素在坡耕地上迁移的影响因素^[18-20]。本研究表明，与T₁相比，各处理可显著减少坡耕地的产流产沙量0.97%~40.58%。这是由于秸秆均匀翻入土表，可避免雨滴直接击打地表，增加地表的粗糙程度，减缓径流流速，从而减轻坡耕地地表溅蚀，减少坡面径流和泥沙的产生；同时，秸秆本身结构比较疏松，具有较强的持水能力，使土壤空隙得到改善，大量降水储存在土壤中，从而降低了坡耕地的产流产沙量^[21]。其中，1.5 kg·m⁻²秸秆密度、5 cm秸秆粒度、加水堆沤的处理能较好降低坡耕地产流产沙量。这可能是由于较高秸秆密度、粗颗粒秸秆及加水堆沤的秸秆翻入土表，相比之下地表更为粗糙，截留的降水较多，易于减缓雨水对坡面的击溅，从而较大程度上减少径流泥沙的产生。

本试验中，与T₁相比，各处理径流总磷浓度增加，这可能是由于秸秆本身含有丰富的氮、磷、钾、等营养元素^[22]，随秸秆腐解，堆肥成肥增加了氮磷钾含量，致使径流泥沙磷素浓度上升。在施肥前期，由于产流量较大，径流中携带着大量磷素，导致磷素极易随径流流失。到施肥后期，泥沙总磷流失浓度逐渐降低9.33%~42.97%，其原因可能是，秸秆矿化分解为土壤提供丰富的有机碳及氮、磷、钾等矿质营养元素来改善土壤理化性状和生物学性状，同时秸秆腐解时可富集大量微生物，吸收水土环境中的磷素^[23]，降低磷素流失风险。

径流中磷素损失包括土壤全磷和土壤中可溶性磷的损失^[24]。本研究表明，径流PO₄⁻-P浓度变化规律同径流总磷浓度变化规律一致。在4场降水中，PO₄⁻-P/总磷变化幅度为35.30%~77.59%，且9月8日占7月28日的91%。这可能是因为溶解态磷素为各处理径流总磷流失的主要方式。在本试验中，各处理PO₄⁻-P流失浓度较T₁高16.29%~71.14%，这可能是由于秸秆在腐解时产生有机酸，有机酸与磷酸根之间竞争吸附，降低土壤矿物仅对磷酸根的吸附，导致磷素在土壤中迁移变得相对容易^[25]，徐泰平等^[26]研究表明，秸秆还田较化肥氮磷钾对照处理渗漏径流增加了30%~52%，可溶性磷素释放与土体内水流大小有很大关系，秸秆堆沤还田可改善土壤结构，增加土壤渗漏径流，故增加了土壤磷素淋失风险，因此在从事农事活动时要严格控制秸秆还田量及密度并且减少尿素施用量。

3.2   耕地磷素流失量的变化特征

本试验中，径流输出是二龙潭小流域坡耕地磷素流失的主要途径，径流中磷素流失量为1.97~11.24 mg·m⁻²，占磷素流失总量的68.78%~74.44%。由于二龙潭小流域坡耕地为红壤且呈酸性，土壤养分较少，加之雨季多发暴雨，土壤淋溶作用较强，土壤孔隙相对较大，坡面易于产生径流^[27]。Li等^[28]研究表明，地表径流是土壤养分流失的关键因素。降水后在坡耕地上从事农业生产活动会引发土壤磷素流失，其流失量与坡耕地利用方式、土壤磷含量、施肥量密切相关，故在降水频发季节进行耕作、施肥、种植等农事活动是导致土壤磷流失的重要影响因素^[29]。在本研究的4场典型降水中，施用0.75 kg·m⁻²秸秆密度、5 cm秸秆粒度、加水堆沤均能有效减少径流和泥沙中总磷流失量。各处理中，T₂较其他处理可以显著降低径流泥沙总磷流失量，这是由于粗颗粒秸秆在还田后降解速度相对较缓慢，可防止更多的土壤颗粒被雨水冲刷，降雨能被更好地截留，从而减少径流的侧向流动^[30-31]，最终降低坡面的产流产沙量，达到降低磷素流失的目的；也可能是由于在秸秆还田过程中，秸秆能吸附一定量的磷，改变耕地土壤的生态环境，微生物能吸收环境中的磷^[32]，从而阻止磷和一些养分的流失和溶淋，降低坡耕地径流和泥沙中总磷含量。同时，在4场降雨中，可以看出当降水量达到最高时，产流产沙量也达到最大且携带的总磷含量也达到最高，这是由于在较大的降水强度下坡耕地上水土流失严重，增加了磷素流失风险。

综上所述，二龙潭小流域烟草坡耕地农田生态系统可采用0.75 kg·m⁻²秸秆密度、5 cm秸秆粒度以及加水堆沤处理降低该流域坡耕地径流泥沙及磷素的流失。为了控制该流域磷素流失的发生，应尽量减少或避免在雨季进行农事活动，减少对土壤的扰动，以利于降低磷素流失含量，并提高基肥、追肥和秸秆的利用效率；在进行农事活动时考虑秸秆还田与减量施肥相配合，在减小坡耕地产流产沙量的同时提高化肥利用率，降低坡耕地面源污染，有效避免磷素淋溶流失对该流域坡耕地带来的不利影响。

3.3   结论

1)施用0.75 kg·m⁻²秸秆密度、5 cm秸秆粒度、加水堆沤，径流总磷与PO₄⁻-P浓度分别降低了2.82%~66.67%和1.97%~64.73%，溶解态磷素占总磷浓度的35.30%~77.59%，为径流磷素流失的主要方式；

2)施用0.75 kg·m⁻²秸秆密度、5 cm秸秆粒度以及加水堆沤泥沙总磷流失浓度呈规律性下降，幅度为0.38%~57.53%；随着烟草作物生长时间的推移，各处理泥沙总磷流失浓度呈现逐渐降低的趋势；

3)施用0.75 kg·m⁻²秸秆密度、5 cm秸秆粒度及加水堆沤可降低二龙潭小流域坡耕地磷素的流失风险，径流和泥沙总磷流失量分别最高降低63.64%和64.74%、63.89%和57.87%、63.89%和64.74%；径流中磷素流失量为1.97~11.24 mg·m⁻²，占磷素总流失量的68.78%~74.44%。说明造成该流域农业面源污染的主要影响因素是降雨产生的径流所携带的磷素流失。