试验小区设在陕西地建土地工程技术研究院富平中试基地(108°57'—109°26'E,34°42'—35°6'N),海拔375.8～1 420.7 m,半干旱大陆性气候,年平均降水量527.2 mm(1960—1995年),年内降雨分布极不均匀,多集中在7—9月,占全年降雨量的49%,其他季节较干旱。年总辐射量5 187.4 MJ /m²,年平均日照时数约2 389.6 h,年均气温为13.1℃,降水年际变率大,年降水变异系数(CV)达21.2%。该区自然条件基本上满足小麦、玉米、棉花等作物的生长需求^[12]。

试验处理设计为砒砂岩与沙的混合体积比例分别为1:1,1:2,1:5,每个处理3次重复,共计9个试验小区,每个小区面积为2 m *2 m。小区排列顺序采用随机排列方法,均在表层30 cm(耕作层)覆盖不同比例的复配土,30 cm以下为当地土壤。添加的砒砂岩处理过筛小于4 cm,确保砒砂岩与沙按比例均匀混合,砒砂岩与沙均拉运于毛乌素沙地榆林市榆阳区大纪汗村。小区采用小麦—玉米轮作模式,其中小麦品种为小堰22,玉米品种为伟科702。2010年6月—2018年10月,各处理均采用常规耕作和当地传统水肥管理措施,施肥为有机肥。不同复配比例的复配土及原始风沙土的主要理化性质见表1。

表1 土壤主要理化性质

每年于作物种植前以及收获后采集土壤样品。每个试验小区按照“S”形5点取样法采集0—30 cm耕层土样。除去粗根及小石块,将土样自然条件下风干。采用四分法(对角线法)分土,然后所有风干土壤过2.0 mm筛子。用马尔文激光粒度分析仪(Mastersizer 2000)分别测定其颗粒组成,每种土样重复3次。依据美国农部制土壤质地三角图确定质地。

土壤颗粒分形维数通常可分为质量分形维数和体积分形维数。土壤颗粒质量分形维数在计算时需要假设不同粒级的土壤具有相同的密度,存在着一定的不合理性。随着激光衍射技术的深入研究和发展应用,可准确获得土壤颗粒体积分布的大小和数量,因此本文采用Tyler等^[13]提出的体积分形模型计算土壤颗粒体积分形维数来表示土壤颗粒分形维数,公式为:

(V(di))/(V_t)=((d_i)/(d_max))^3-D_v

式中:d为粒径; d_i为粒径分区中的第i级粒级,其值一般为该粒级区间上下限的算术平均值; d_max为所有土壤颗粒中的最大粒径,本研究中d_max=2000; V为粒径小于d_i的土壤颗粒累积体积; V_t为全部土壤颗粒总体积。

计算步骤为:在各分级区间,计算d_i的算术平均值及累积体积,画出对数图,拟合的直线斜率为3-D_v,根据斜率从而求出土壤颗粒体积分形维数D_v值。

依据中国制分级标准(1987年)对土壤粒径进行分级:细黏粒(<0.001 mm),粗黏粒(0.001～0.002 mm),细粉粒(0.002～0.005 mm),中粉粒(0.01～0.005 mm),粗粉粒(0.01～0.05 mm),细砂粒(0.05～0.25 mm),粗砂粒(0.25～1 mm)。砒砂岩与沙按不同比例复配后土壤颗粒粒径分布见图1。1:1复配土的中粉粒(0.01～0.005 mm)含量最高,体积分数为26.82%,次之为粗砂粒(0.01～0.05 mm)25.55%,最低为细黏粒(<0.001 mm)5.13%; 1:2复配土中细砂粒(0.05～0.25 mm)最高为38.65%,次之为粗粉粒(0.01～0.05 mm)21.30%,最低为细黏粒(<0.001 mm)3.94%; 1:5复配土的中粗砂粒(0.25～1 mm)体积高达55.48%,并从粗砂粒到细黏粒依次递减,细黏粒(<0.001 mm)含量仅为2.40%。从复配土黏粒(<0.002 mm)、粉粒(0.05～0.002 m)、砂粒(1～0.05 mm)的累积粒径来看,黏粒(<0.002 mm)和粉粒(0.05～0.002 m)的累积粒径含量均表现为1:1>1:2>1:5,砂粒(1～0.05 mm)的累积粒径含量表现为1:5>1:2>1:1,表明复配土中砒砂岩含量越高,颗粒粒径越细。

图1 复配土颗粒粒径分布

用最小二乘法对复配土粒径分布的对数图中回归曲线进行拟合(表2),由表2可知,不同复配比例的拟合方程决定系数(R²)均高于0.8,最大值达到1,平均值均在0.9以上,说明拟合方程效果极好,达到了显著和极显著水平。从同一复配比土壤颗粒体积分形维数来看,土壤颗粒分形维数变化差异不大,变异系数均低于0.5,说明样品的颗粒粒径组成较为稳定。从不同复配比的土壤颗粒体积分形维数来看,土壤颗粒分形维数表现为1:1>1:2>1:5>0:1。

表2 不同配比下土壤颗粒分形维数及拟合方程决定系数

由土壤颗粒分形维数计算公式可知,土壤颗粒分形维数与土壤各粒级粒径的累积量密切相关,按美国制土壤粒径分级标准,将土壤粒级分为黏粒(<0.002 mm)、粉粒(0.05～0.002 m)、砂粒(1～0.05 mm),对各粒级体积百分比与土壤颗粒分形维数进行回归分析(图2—4)。

土壤颗粒分形与黏粒体积百分含量直接呈显著的对数关系(图2),p=0.001^**,决定系数(R²=0.9982)较高。表现为随着黏粒含量的增加,分形维数显著增加。复配土粉粒相互之间呈显著相关,p=0.001^**,但拟合方程的R²(0.089 8)较小,说明土壤颗粒分形与粉粒质量百分含量直接的相关性很小(图3)。由土壤颗粒分形与砂粒体积百分含量的拟合与粉粒体百分含量拟合类似,p=0.002^**,说明复配土砂粒体积百分含量显著相关,但砂粒的拟合方程R²(0.030 2)更小,说明土壤颗粒分形维数与砂粒体积百分含量相关性更小(图4)。

质地是土壤中大小各异的颗粒组合状况,是土壤结构体形成的基础,与土壤通气、保水、保肥性及耕作难易程度有着密切的关系。不同配比下砒砂岩改良风沙土质地变化见图5。随着种植年限的增加,风沙土结构体中赋存的关键粒级含量(黏粒和粉粒)逐渐提高,在种植6 a内,1:1,1:2,1:5黏粒平均含量分别为(4.60±2.15)%,(3.34±1.62)%,(3.08±1.40)%,较沙地而言,1:1黏粒含量提高了62.13%,1:2,1:5也提高了40%以上; 1:1,1:2,1:5粉粒平均含量分别为(31.14±13.72)%,(21.72±7.38)%,(19.80±6.20)%,较沙地而言,不同配比下的粉粒含量均提高了70%以上,其中1:1粉粒含量提高了81.93%。

图2 黏粒含量对土壤颗粒体积分形维数的影响

图3 粉粒含量对土壤颗粒体积分形维数的影响

图4 砂粒含量对土壤颗粒体积分形维数的影响

图5 不同配比下砒砂岩改良风沙土质地变化

从颗粒组成来看,砒砂岩混合比例越高,复配土中的黏粉粒含量越高; 此外,1:1复配土土壤颗粒主要集中在粉粒段(0.002～0.05 mm),粉粒含量占49.37%,而1:2,1:5复配土土壤颗粒主要集中在砂粒段(0.05～2 mm),含量分别为52.74%,75.55%。这可能是由于砒砂岩黏粒含量高达10.3%～30.3%^[14],而沙土90%为砂粒,且颗粒孔隙大,当砒砂岩与风沙土混合后,砒砂岩的黏粒填充到沙土孔隙中,因此砒砂岩混合越多,复配土中的黏粒、粉粒含量越高。

从土壤颗粒体积分形维数来看,分形维数表现为1:1>1:2>1:5,这是由于砒砂岩在风沙土中的含量越高,黏粉粒含量就越高,土壤颗粒就越细粒化,分形维数越高,表征土壤结构逐渐变好的变化趋势。陈小红等^[15]在研究沙漠化的逆转过程时发现,随着沙漠化程度的减小,土壤颗粒分形维数逐渐变大,说明砂粒含量越小,分形维数越大,与本研究结果一致。罗雅曦等在研究腾格里沙漠草方格固沙林地土壤颗粒组成及分形维数的关系时也得到了一致的研究结果。

从土壤粒级与颗粒体积分形维数拟合来看,土壤体积分形维数和土壤黏粒体积含量呈显著的正相关关系,而与粉粒和砂粒的体积含量相关性较小(图4),这与王国梁等^[17]研究土壤颗粒的体积分形维数及其在土地利用中的应用结果一致。杨金玲等^[11]在对比土壤颗粒粒径分布重量分形维数和体积分形维数也发现,体积质量分形维数与黏粒关系更加密切。

复配土中粉粒含量基本随着种植年限的增加而改变,尤其是1:1复配土中粉粒含量明显增加,质地也随之改变,总的来说质地是一个变化较为缓慢的过程。这可能是由于随着种植年限的增加,有机肥不断施入,作物根系分泌物增加,土体中微生物种类变多,代谢产物也随之增加^[18-19],促进了复配土颗粒的优化,但粒径的改变需要达到一定数量后才能引起质地的改变。