研究区位于四川盐亭农田生态系统国家野外科学观测研究站(105°27′E，31°16′N)所在的小流域—四川盆地中北部的盐亭县林山乡截流村. 该区地形为中丘，海拔400~600 m，由水平砂泥岩互层形成多级梯地，山顶为园丘、长岗状. 沟谷切割较深，冲沟发育，相对高差100~200 m，谷底宽50~150 m，比降1/100~1/350，两侧山坡较陡，平均坡度为8.7%~36.4%^[10]. 该区属中亚热带季风气候，年均气温17.3 ℃，多年平均降水量为826 mm，春季降水量占全年降水量的5.9%，夏季占65.5%，秋季占19.7%，冬季占8.9%. 土壤类型以水稻土和石灰性紫色土为主，主要农作物有水稻、玉米、小麦、油菜等，林地植被以桤木和柏木为主^[11].

小区(4 m×2 m)试验地从2005年开始自然撂荒，至2018年，已有13年时间. 撂荒前小区内主要种植作物为玉米，目前小区内主要植被物种有：白茅、马桑、小蓬草(表 1)，其中白茅分布最多，覆盖度最大. 撂荒开始时(2005年)采集了5个试验小区0~20 cm土层的土壤混合样，测试了各理化指标. 2018年8月2日从试验地中5个小区分别采集0~20 cm土层的土壤混合样，然后，用撂荒开始时同样的方法测试以下理化指标：土壤颗粒组成、土壤有机碳、全氮、全磷、全钾、碱解氮、有效磷、速效钾. 由于耕作深度和撂荒后植物根系作用深度基本都在0~20 cm，因此采集0~20 cm土壤样品，能较好地反映撂荒前后土壤理化性质的变化.

测试所得数据，用Excel 2010做基本处理，用IBM SPSS Statistics 22进行撂荒前后土壤各养分的配对样本t检验，用OriginLab origin 2019作图.

设参与评价土壤肥力的因子有n个，用U={U₁，U₂，…，U_n}表示，集合中U_i是n个土壤养分因子的检测数值集合；设有m种不同的评价等级，它们组成评价集V={V₁，V₂，…，V_m}，集合中V_j是与U_i中指标相对应的评价标准集合. 用r_ij表示第i种因子的数值可以被评价为第j级的可能性(即i对j的隶属度，它们的关系即为隶属函数)，构成U，V之间模糊关系矩阵R：

r_ij通过以下隶属函数来确定：

第1级，即j=1，其隶属函数为：

第2~(m-1)级，即j=2，3，…，m-1时，其隶属度函数为：

第m级，即j=m时，其隶属度函数为：

式中：X_i为第i种因子实测值；S_ij为第i种因子的第j级评价标准值.

因子集U中各参评指标的权重系数w_i构成模糊子集W=(w₁，w₂，…，w_n)，其中：$\mathop \sum \limits_{i = 1}^n {w_i} = 1$. 用OriginLab origin 2019对各指标检测数值进行主成分分析，根据分析结果中主成分贡献率和因子载荷确定各参评指标的权重.

已知权重矩阵W及模糊关系矩阵R，则n个评价指标的多指标模糊综合评价矩阵可以表示为：

模糊综合指数可以表示为：

式中：G为土壤养分等级标准向量，G^T=(1，2，…，m)；FCI为土壤养分的综合等级.

用MATLAB R2013b进行矩阵运算.

据图 1所示，撂荒前土壤砂粒(0.05~2 mm)比例为15.8%，撂荒后为28.4%，撂荒后增加12.6%；撂荒前土壤粉粒(0.002~0.05 mm)比例为65.4%，撂荒后为51.9%，撂荒后减少13.5%；撂荒前后土壤黏粒(＜0.002 mm)比例分别为18.8%和19.7%，变化不大. 根据美国土壤质地分类标准，撂荒前后的土壤质地均为粉壤土，不过由于砂粒的增加和粉粒的减少，撂荒后的土壤质地有由粉壤土变为壤土的趋势，即撂荒后土壤有粗骨化的趋势.

撂荒前，由于耕作熟化，土壤的颗粒组成偏细，砂粒比例低. 在耕时及撂荒初期，地表植被盖度低，土地裸露面积较大，降雨时雨滴、径流动能和地表径流量大，无论是雨滴击打对土壤的分散作用，还是径流对泥沙的携带能力都很强，从而土壤侵蚀强烈，带走大量泥沙，且泥沙颗粒分选性较差^[12-13]，土壤无论是粗颗粒还是细颗粒均容易流失. 相比在耕地，撂荒后人为干扰减少，随着撂荒年限的增加，植被盖度增加，植被对降雨的拦截率增加，坡面粗糙度增大，径流速率和地表径流量减小，植被根系发达，能疏松土壤，促进团聚体形成，从而入渗及壤中流增加^[14-16]，并且白茅低矮且盖度高，可以有效地削弱雨滴对土壤的击溅动能，减少土粒的分散^[16]，因此土壤侵蚀减少，且侵蚀带走的泥沙颗粒具有较高的分选性：较粗的砂粒因质量高不容易被较小的径流搬运而得以保留，黏粒由于其比表面积及电荷密度大的特性，容易胶结土粒形成团聚体也不容易流失，相比之下粉粒流失的比例就较多. 另一方面，由于紫色土母岩极易风化且物理风化强烈，成土快速^[17]，植被恢复过程中植物根系的增加会加强母岩母质的破碎风化. 采样时可见小区内土壤中含有大量的植物根系，由于小区内土层较浅，植物根系很容易穿透土壤到达母质母岩层，促进母岩裂隙的发育，加快母岩破碎风化，从而增大了土壤中粗颗粒的比例. 并且，和在耕时不同，撂荒后不存在植物收割和土地翻耕，植物根系对母岩的破碎风化作用是持续的.

由图 2可知，撂荒前后土壤有机碳、全钾、碱解氮质量分数差异均无统计学意义(p＞0.05)，但是撂荒后这几种养分的均值都大于撂荒前；撂荒后土壤全氮质量分数显著增加(p＜0.05)，且增幅可达61.5%；撂荒后土壤全磷、有效磷、速效钾质量分数显著减少(p＜0.05)，有效磷减少幅度最大，可达61.4%，全磷、速效钾分别减少14.0%，22.6%.

土壤有机质是土壤非常重要的组成部分，在土壤的形成和肥力发展过程中起着极其重要的作用，自然植被条件下植物残体和根系分泌物是土壤有机质的主要来源. 耕地撂荒后，无人为收割，植物残体归还到土壤中，可以增加土壤中有机质的比例，但是，紫色土区气候温暖湿润，微生物活动强烈，不利于有机质的积累，因此撂荒后有机碳比例增加并不显著. 土壤中的氮素主要来源于生物固氮和氮沉降、施肥等，耕地撂荒后不再向土壤中投入肥料，然而撂荒13年后小区土壤中的全氮比例增加了61.5%，一个很重要的原因可能是氮沉降的增加. 况福虹^[18]的研究显示：氮沉降量在四川盆地的增长非常迅速，2009-2014年该区域湿沉降平均增至22.5 kg N/(hm²·a)，比起2002-2004年，几乎翻了一倍. 磷在土壤中移动性小，土壤中的磷素比例与母质有很大关系^[14]. 撂荒后土壤中的全磷显著降低，可能是由于相比在耕地，撂荒后少了磷肥的投入，一方面磷被植物吸收却不能及时得到补充，另一方面磷随土壤流失流出土体. 土壤中的钾主要来自成土母质及含钾矿物岩石的风化和分解^[19]，紫色母岩母质含钾丰富，紫色土的全钾比例也比较高^[20]，虽然撂荒后不再施肥，减少了钾肥的投入，且湿热条件下钾易流失，但由于紫色母岩快速风化成土的特性，紫色土中的钾素很容易得到补充，因此撂荒后土壤全钾比例变化不大.

土壤中的碱解氮、有效磷、速效钾是植物可以直接吸收利用或者容易被分解转化以供植物利用的氮磷钾素，能很好地表征土壤的供肥能力. 撂荒后土壤全氮比例增加61.5%，然而碱解氮增加并不显著，可能是因为采样时正值植被生长旺盛的季节，可利用的氮素被植物吸收较多，从而在一定程度上减少了土壤中的碱解氮比例. 紫色土中的有效磷比例很低，常常成为提高作物产量的限制性因子^[20]；虽然紫色母岩成土快，钾素比例多，但含钾矿物通过风化直接转换成速效钾的比例却是极低的^[20]. 因此，施肥对于提高土壤中有效磷和速效钾至关重要，撂荒后有效磷和速效钾显著降低，和施肥停止有很大关系. 另外，采样时茂盛生长的植被对养分的强烈吸收也是土壤中有效磷和速效钾降低的原因之一.

图 3显示，撂荒后，土壤的碳氮比显著降低(p＜0.05)，由撂荒前的13.88变为撂荒后的8.84，降低36.3%；碳磷比极显著升高(p＜0.01)，撂荒前后分别为12.02和14.79，升高23.0%；氮磷比极显著升高(p＜0.01)，升高89.8%，并且比例由小于1变化至大于1，撂荒前后分别为0.88和1.67.

土壤C∶N∶P可以反映土壤内部碳氮磷循环，可以指示碳氮磷的矿化或者固持作用，可作为养分限制和碳氮磷饱和诊断的指标^[21]. 一般来说，土壤碳氮比越低，分解速度越快，元素周转速率越快；磷的有效性由土壤有机质的分解速率确定，较低的碳磷比是磷有效性高的一个指标；氮磷比可以指示土壤养分的限制元素^[21]. 撂荒后，土壤碳氮比降低，说明土壤中微生物活性增强，有机质分解加快，这可能与撂荒后人类干扰少，小区内逐渐恢复为良好的自然生态环境有关. 碳磷比和氮磷比升高，指示着撂荒过程中氮元素的积累及磷元素的缺乏，随着撂荒的进行，植被生长会受到磷元素的限制.

1) 确定评价因子及评价标准

选取有机质、全氮、全磷、全钾、碱解氮、有效磷、速效钾这7种养分作为评价因子，根据全国第二次土壤普查养分分级标准(表 2)确定评价标准，则：撂荒前土壤各养分检测值组成集合U₁={12.24，0.52，0.59，19.45，41.30，5.28，109.46}；撂荒后土壤各养分检测值组成集合：U₂={12.91，0.85，0.51，20.20，47.13，2.04，84.75}. 集合U₁和U₂第1至7个元素依次为有机质、全氮、全磷、全钾、碱解氮、有效磷、速效钾的检测平均值.

评价标准集构成矩阵S：

2) 建立模糊关系矩阵

根据公式(1)(2)(3)，可以得出U，V之间的模糊关系矩阵R，即每一种评价因子对于每一等级的隶属度矩阵.

撂荒前：

撂荒后：

3) 确定各项评价因子的权重

利用主成分分析，得到主成分的贡献率及每个评价因子对每个主成分的因子载荷，以此计算出每个评价因子的权重，如表 3所示.

4) 根据公式(4)，计算多指标综合评价矩阵.

撂荒前土壤养分的多指标综合评价矩阵为：

撂荒后土壤养分的多指标综合评价矩阵为：

5) 根据公式(5)，计算模糊综合指数FCI，量化等级

根据表 2，G^T=(1，2，3，4，5，6)，因此：

撂荒前土壤养分综合指数为：

撂荒后土壤养分综合指数为：

FCI代表各评价因子根据表 2的评价标准所算出来的综合等级，与表 2的等级一样，FCI越小，综合来看土壤的养分越好. FCI₁＜FCI₂，说明撂荒前土壤养分的综合等级高于撂荒后，即撂荒后整体来说土壤的肥力下降了，但是下降得并不多，还是处于同一等级区间，都位于三、四级之间.

撂荒13年后，紫色土坡耕地：

砂粒比例增加，粉粒比例减少，黏粒比例变化不大，依据美国土壤质地分类，撂荒前后的土壤质地均为粉壤土，但撂荒后的土壤质地接近于壤土，土壤有粗骨化的趋势.

土壤全氮比例显著增加，有机碳、全钾、碱解氮比例略有增加，但是变化并无统计学意义，全磷、有效磷、速效钾比例显著降低. 用模糊数学法对土壤肥力进行综合评价，结果显示撂荒后土壤肥力略有下降，但变化不大.

土壤碳氮比显著升高，或可指示有机质分解速率的加快，碳磷比、氮磷比极显著降低，均表明撂荒过程中氮元素的积累和磷元素的缺乏，随着撂荒的进行，磷会成为植被生长的限制元素.

关于耕地撂荒的环境效应，到底是更积极的还是更负面的，一直有争论^[3]. 本研究显示，紫色土坡耕地撂荒13年后，演替为荒草地，同时生长有马桑等小型灌木. 撂荒后砂粒比例增加，粉粒减少，黏粒变化不大，土壤有粗骨化的趋势，这与廖超林等^[22]的研究结果类似，但是撂荒后土壤质地并没有发生改变. 撂荒后除土壤全氮比例以外，其他养分比例并没有显著增加，甚至部分养分特别是有效态养分比例显著降低，这与梁楚涛^[23]的研究结果相似，因为都以在耕地做对比；与乔文静等^[24]的研究结果不同，这是因为他们的研究是选择已经撂荒的不同时间序列的撂荒地做对比，这两种情况下肥料投入的影响是不同的，而且紫色土区由于其温暖湿润的气候，土壤中的养分本就容易流失. 撂荒后土壤碳氮比显著降低，碳磷比、氮磷比极显著升高，但是土壤全氮比例的变化太大，从一定程度上削弱了碳氮比指示有机质分解速率的可靠性，因为氮沉降对于土壤有机碳比例变化的影响是复杂的，受生态系统类型、土壤有机碳组成、施氮类型和剂量等多种因素的影响^[25]. 撂荒13年以来，一直没有肥料投入，采样时间又正是植被生长茂盛、对土壤养分吸收较多的季节，但是综合来看土壤肥力变化却并不大，因此我们认为撂荒是有利于土壤养分的固持与积累的.

关于撂荒对土壤的影响研究，多数都是以空间代替时间的方式来对比^[22-24]，这种方法，即使选择的地块立地条件相似，也很难保证土壤理化背景一样，可能会对分析结果产生影响，而本研究是对比同一地块撂荒13年前后土壤理化性质的变化，相比之下分析结果会更准确. 不过，撂荒前测定的指标也限制了能对比分析的指标种类，另外，撂荒13年过程中间没有再采样，对于撂荒过程中土壤理化性质的动态变化难以掌握.

因为缺乏撂荒早些年的土壤养分数据，我们无法判断撂荒早期土壤养分是如何变化的. 不过梁楚涛^[23]对黄土丘陵区不同年限撂荒地的土壤养分研究显示，随着撂荒的进行，土壤碳氮磷均表现出先降低后升高的趋势，并均在撂荒17年左右达到最低值，土壤有机质直到撂荒演替20年才恢复升高. 参考此研究结果，对于该撂荒地块后续仍有必要进行观测、采样与测定，也可以再加测一些指标，有利于进一步判断紫色土坡耕地撂荒的环境效应随时间的变化情况.