崔 健, 党晓宏, 汪 季, 张 超, 李婉娇, 靳灵娜

(内蒙古农业大学 沙漠治理学院, 呼和浩特 010018)

PLA沙障是一种新型的可降解沙障。对于它的研究目前主要涉及到防风固沙效益、降解性能、铺设方式等方面。党晓宏等[1]通过对比PLA沙障和麦草沙障后发现，PLA沙障防风效能显著大于麦草沙障。在PLA沙障降解性能方面，风沙吹蚀活动是影响PLA沙障材料降解的主要因素[2-4]。PLA沙障影响近地表风速的同时，进而影响沙丘表层土壤理化性质，目前针对PLA沙障内土壤粒度特征方面的研究，党晓宏等[5]对PLA沙障障内0—11 cm的土壤采样并试验，结果发现铺设PLA沙障后表层土壤粒径增大，坡顶处土壤粒径变化最大，表层土壤粒径变化最大；丁延龙等[6]通过研究发现铺设PLA沙障后表层土壤出现粗粒化现象。PLA沙障使用寿命一般在8～10 a，前人只是对刚铺设后或铺设1～2 a后障格内的土壤粒度进行了研究，而铺设5 a后的土壤粒度呈现一种什么样的特征，目前还尚未明确。同时在铺设5 a以后，随着风沙活动沙物质沉积、搬运、再沉积、再搬运这一循环的过程中，对于有机质是否发生积累或者流失还没有明确研究，与土壤粒度存在一种怎样的关系到目前为止也尚未清晰。

基于此，本文以库布齐沙漠东南缘铺设5 a后的PLA沙障为研究对象，以未铺设PLA沙障的裸沙丘为对照，通过分析不同规格PLA沙障障格内表层沉积物粒度及有机质含量，明晰铺设5 a后不同规格PLA沙障内土壤理化性质如何变化，对土壤粒度和有机质特征进行深入分析的同时对二者相关性进行分析，探究二者存在一种怎样的联系，根据最后得出的结论为区域土壤质量评价与管理以及风蚀防治提供基础数据和科学依据。

1 研究区概况

研究区位于库布齐沙漠东南缘，地理坐标106°55′16″—109°16′08″E，39°22′22″—40°52′47″N，属于温带大陆性半干旱气候，多年平均降水311.8 mm，全年盛行西风和西北风，年平均风速2.8 m/s，土壤类型主要有栗钙土、棕钙土和灰漠土；天然植被以沙生植物为主，主要有沙拐枣(Calligonumarborescens)、黑沙蒿(Artemisiaordosica)、羊柴(Hedysarummongolicum)、沙蓬(Agriophyllumsquarrosum)。

2 材料与方法

2.1 样地选择与样品采集

选取库布齐沙漠东南缘铺设5 a的PLA沙障作为研究对象，以未铺设PLA沙障的流动沙丘作为对照，于2020年7月采用5点取样法，分别采集3种规格PLA沙障(0.5 m×0.5 m，1 m×1 m，2 m×2 m)、裸沙丘4种样地、3个坡位(迎风坡底、坡中、坡顶)0—3 cm近地表土样，每种规格内各个坡位设置3个障格为1组重复,重复取样3次，将所取土样装入塑封袋带回实验室进行后处理。

将所有土样去除根系和枯落物等杂质后过2 mm筛，加体积分数为30%的H2O2去除土壤中的有机质；加体积分数为10%的HCL溶液煮沸去除土壤中的碳酸盐；加去离子水稀释并静置12 h，去除上清液，并反复进行稀释，直至pH值为6.5～7.0；加六偏磷酸钠，在超声波环境下进行30 s处理，用激光粒度仪测量土壤粒径的体积分数此试验土壤粒径分级标准采用美国制分级标准[7]。有机质含量将待测土壤过0.15 mm筛后采用常规的重铬酸钾容量法-稀释热法测定[8]。

2.2 粒度参数

采用伍登-温德华粒级标准，根据克鲁宾对数转化法，分别将先前输出的各土壤颗粒累积体积分数对应的颗粒直径进行转换，有利于计算Φ值[9]，见公式(1)。

Φ=-L=-log2D

(1)

通过克伦拜因和福克计算累计频率曲线图上累计百分含量5%，10%，16%，25%，50%，75%，84%，95%所对应粒度对数值Φ5，Φ16，Φ25，Φ50，Φ75，Φ84，Φ95，据此计算粒度特征参数[9]，见公式(2)—(6)。

平均粒径(d0)反映沙物质粒度平均状况的参数，d0越大，表示细粒物质越多。计算公式为：

(2)

标准偏差(σ0)反映沙物质粒径分布的分散程度，σ0越小，颗粒越均匀，粒径分配的分散程度就越小。根据取值范围划分7个分选级别(从集中至分散)：分选极好(σ0≤0.35)、分选好(0.35<σ0≤0.50)、分选较好(0.50<σ0≤0.71)、分选中等(0.71<σ0≤1.00)、分选较差(1.00<σ0≤2.00)、分选差(2.00<σ0≤4.00)、分选极差(σ0>4.00)。计算公式为：

(3)

偏度(S0)反映了沙物质粒度粗细分配的对称性，根据取值范围分5个偏度等级，从细粒物质占比大过渡到粗粒物质占比较大：极负偏度(-1.0≤S0≤-0.3)、负偏度(-0.3

(4)

峰态值(K0)反映颗粒粒度的集中程度，根据取值范围分6个峰度等级，从比较分散过渡到比较集中：很宽平(K0≤0.67)、宽平(0.673.00)

(5)

2.3 数据处理与分析

采用Excel 2010进行前期的数据统计和处理工作，运用Origin 2019b和IBM SPSS Statistics 25软件对不同规格的PLA沙障障格内表层沉积物粒度组成及有机质含量进行单因素方差分析(One-way ANOVO)，采用最小显著差异法(LSD)进行差异显著性检验，置信区间为95%。

3 结果与分析

3.1 不同规格PLA沙障表层沉积物极配特征

由表1可以看出，土壤颗粒组成主要以细砂和粗砂为主。首先，裸沙丘土壤颗粒以细砂和中砂为主，黏粒和粉砂含量很少接近0，细砂含量为33.95%～41.85%，中砂含量为34.41%～39.83%。3种规格沙障障格内土壤颗粒变化情况相比于裸沙丘均表现为细砂、中砂含量降低，中砂含量依次为0.02%～0.56%，0.1%～1.23%，0.37%～3.78%，细砂含量为10.50%～23.39%，18.83%～31.73%，13.28%～41.17%，细砂含量随着铺设规格的增加而增加；粗砂和极粗砂含量随着铺设规格的增加而减少，粗砂含量分别为33.71%～44.91%，25.50%～36.42%，17.53%～42.88%；极粗砂含量依次为35.72%～38.93%，30.21%～35.61%，28.31%～38.24%。其次，每种障格内不同坡位土壤颗粒含量也存在一定差异，3种规格沙障内坡顶处细砂含量均为最高依次为23.39%，31.73%，41.17%，且在1 m×1 m和2 m×2 m规格障格内细砂含量均呈现出坡顶>坡底>坡中的规律，坡中粗砂和极粗砂含量均表现为最高，粗砂含量依次为44.91%，36.42%和42.88%，极粗砂含量依次为38.93%，35.61%，38.24%。

表1 不同规格PLA沙障表层沉积机械组成 %

3.2 不同规格PLA沙障沉积物频率分布曲线

从图1可以看出，4种样地3个坡位的土壤颗粒频率分布曲线总体走向一致，但个别存在一定的差异。3种规格PLA沙障内和裸沙丘沉积物颗粒频率分布曲线均为多峰型曲线，3种规格分布曲线在<100 μm，100～1 000 μm和>1 000 μm这3种范围内均出现波峰。3种坡位的颗粒频率分布曲线图中在10 μm附近均有明显的凸起，且土壤体积百分含量均高于裸沙丘，由表1可知，在此范围内是由于黏粒和粉砂含量增加；在100～1 000 μm范围内，3种坡位均呈现出障格内土壤颗粒粒径相比裸沙丘土壤颗粒粒径都具有明显减小的趋势，但其含量却存在一定的差异，图1A和图1B中3种障格内土壤体积百分含量显著低于裸沙丘土壤体积百分含量，在坡底土壤体积百分含量由高到低依次为2 m×2 m>1 m×1 m>0.5 m×0.5 m，在坡中3种障格内表现为1 m×1 m土壤体积百分含量最高，其余2种规格土壤体积百分含量基本一致无明显差异，但在图1C中3种障格内土壤体积百分含量与裸沙丘土壤体积百分含量并无差异；在>1 000 μm范围内土壤体积百分含量变化明显，处于3种坡位内不同规格障格中土壤体积百分含量均显著高于裸沙丘土壤体积百分含量，且在此范围内0.5 m×0.5 m规格的土壤体积百分含量均为最多且高于其余2种规格。

图1 PLA沙障内表层土壤颗粒频率分布曲线

累积频率分布曲线能反映土壤颗粒的分布情况，通常曲线越陡，颗粒分布越均匀。图2粒径<100 μm时，沙障规格越小，障格内土壤颗粒分布越均匀，由图2C可以看出，位于坡顶时土壤颗粒累积情况发生明显的变化，随着沙障规格的增加障格内土壤颗粒分布越均匀，土壤分布均匀程度由高到低依次为2 m×2 m，1 m×1 m，0.5 m×0.5 m，CK；当粒径>100 μm时，3种规格内土壤颗粒累积分布均低于裸沙丘，相反，在坡底和坡中是1 m×1 m规格内土壤颗粒分布最均匀，在坡顶2 m×2 m规格内土壤颗粒分布最均匀，在3种坡位中0.5 m×0.5 m规格内土壤颗粒分布情况是最差的。

图2 PLA沙障内表层土壤颗粒累积频率曲线

3.3 不同规格PLA沙障沉积物粒度参数

按Folk-Ward图解法的划分标准，CK表层沉积物平均粒径为1.293～1.458 Φ，均值为1.402 Φ，0.5 m×0.5 m，1 m×1 m，2 m×2 m平均粒径为0.474～0.881 Φ，0.745～2.028 Φ，0.613～1.451 Φ，均值分别为0.707 Φ，1.197 Φ，0.971 Φ；裸沙丘沉积物分选系数为0.96～1.047，均值为0.990，0.5 m×0.5 m，1 m×1 m，2 m×2 m规格沙障内沉积物分选系数分别为1.096～2.194，1.868～3.131，1.227～2.121，均值分别为1.697，2.310，1.835；裸沙丘沉积物偏度-0.404～-0.293，均值为-0.363，为负偏度。0.5 m×0.5 m，1 m×1 m，2 m×2 m规格障内沉积物偏度分别为0.769～0.871，0.819～0.884，0.016～0.841，均值分别为0.825，0.849，0.549，裸沙丘表层沉积物峰态值为0.810～1.273，均值为1.072，颗粒集中程度为中等。0.5 m×0.5 m，1 m×1 m，2 m×2 m规格障格内沉积物峰态分别为1.025～1.424，1.195～1.358，1.344～1.384；均值分别为1.205，1.288，1.364。

由图3可知，3种规格PLA沙障内表层土壤平均粒径均显著低于裸沙丘土壤平均粒径，仅坡顶2 m×2 m和坡底1m×1m规格内土壤平均粒径高于裸沙丘；在坡底和坡中位置3种规格分选性都高于裸沙丘，且都呈现出1 m×1 m规格内分选性最高；偏度则表现出均高于裸沙丘，且在坡底和坡中处3种规格偏度变化情况基本保持一致，无明显的起伏变化，坡顶出现较大的差异；峰态在坡中呈现剧烈波动的现象，在坡底处0.5 m×0.5 m规格内峰态显著高于其余2种规格，且1 m×1 m， 2 m×2 m规格与裸沙丘相比没有明显变化，在坡中处0.5 m×0.5 m规格均显著低于其余2种规格及裸沙丘，在坡顶处3种规格均显著高于裸沙丘，2 m×2 m规格峰态最高。

注：MZ表示平均粒径；σ表示标准偏差；SK表示偏度值；Kg表示峰态值。

3.4 土壤有机质特征

研究区土壤有机质的测试分析结果(表2)显示，3种规格PLA沙障内表层土壤有机质含量为0.7～1.14 g/kg。0.5 m×0.5 m规格内呈现出坡底(1.05 g/kg)>坡中(0.78 g/kg)>坡顶(0.75 g/kg)，其中，坡底和坡顶有机质含量均有显著差异；1 m×1 m规格内呈现出坡底(1.14 g/kg)>坡中(0.77 g/kg)>坡顶(0.7 g/kg)且坡顶与坡底、坡中与坡底有机质含量均具有显著差异；2 m×2 m规格内呈现出坡底(1.04 g/kg)>坡中(0.75 g/kg)>坡顶(0.73 g/kg)，其中，坡顶与坡底、坡中与坡底均具有显著差异；CK各坡位呈现出坡顶(0.92 g/kg)>坡中(0.85 g/kg)>坡底(0.8 g/kg)且有机质含量均无显著差异。由此可见，整体上4种不同类型PLA沙障的表层土壤有机质含量表现出1 m×1 m>0.5 m×0.5 m>CK>2 m×2 m的规律。

3.5 土壤粒度组分与有机质的相关性

从表3可以看出，不同规格PLA沙障表层土壤有机质与土壤粒级含量的相关性有着显著差异。其中，0.5 m×0.5 m规格PLA沙障内有机质与极细砂相关系数为0.966，呈现显著正相关关系，与中砂相关系数为-0.877，呈现显著负相关；1 m×1 m规格PLA沙障内有机质与极细砂相关系数为0.999，呈现显著正相关，与中砂相关系数为-0.98，呈现显著负相关；2 m×2 m规格PLA沙障内有机质与极细砂相关系数为0.975，呈现显著正相关，与中砂相关系数为-0.707，呈现显著负相关。由此可知，研究区中极细砂和中砂是3种规格PLA沙障内表层土壤有机质积累与否的关键粒级。

表2 不同规格PLA沙障表层土壤有机质含量

表3 不同规格PLA沙障表层土壤有机质与粒度组分的相关性

4 讨 论

4.1 不同规格PLA沙障防沙效果

分析PLA沙障铺设后表层沉积物易蚀颗粒的变化可有助于细化PLA沙障阻滞风沙流效果，从图4可以看出，3种规格PLA沙障在不同坡位下均表现为细砂和中砂含量相较裸沙丘减少，粗砂含量较裸沙丘升高，粉砂和黏粒含量较裸沙丘升高。铺设5 a后的PLA沙障土壤粒径均有增大趋势，0.5 m×0.5 m规格沙障由于在铺设时设计铺设的规格较小，在铺设5 a后易蚀颗粒几乎被完全风蚀，细砂和中砂含量远低于裸沙，粗砂和极粗砂含量远高于裸沙，大部分的障体已经完全被风沙所掩埋，导致0.5 m×0.5 m规格的沙障防护效益大大降低，这与屈建军等[10]的研究结果发现0.5 m×0.5 m规格沙障障体很快被沙埋后失去防护作用结论相一致。沙丘各部位的1 m×1 m规格PLA沙障内极细砂、细砂含量均降低，而粗砂和极粗砂颗粒相较于裸沙丘含量均升高；通过分析可知，铺设沙障后，沙障能改变地表粗糙度，使近地表风速降低，当饱和风沙流在遇到沙障后，由于沙障的阻挡导致其风速降低，同时在障格内形成小范围紊流，从而造成风沙流携砂能力降低，导致质量较大的砂砾在沿着主风向继续向上运动过程中动能减小从而掉落于障格内，由于风沙流要继续向上运动，而此时的风沙流结构变为不饱和风沙流，沿着主风向的过程中携带走障格内细小砂砾，此时的沙障内部以风蚀为主。2 m×2 m规格沙障通过实地观测发现由于障格面积过大并且铺设时间过长，部分障体受到风力作用障体发生扭曲，防护功能大大降低。风沙流由坡底向上运动到坡中的过程中，部分障体受到破坏阻挡作用降低，携带走大量的细小物质，所以使其障格内细砂含量降低，粗砂和极粗砂含量相较于裸沙丘含量升高，而此时风沙流要继续向上运动，粗砂继续掉落，风沙流趋于不饱和，当风沙流到达坡顶时发生蚀积转换，中砂物质被风蚀，虽然土壤粒径整体变粗，但细砂含量却有所增加。综上所述，0.5 m×0.5 m规格防护效果最差，1 m×1 m规格防护效果次之，2 m×2 m规格虽然较1 m×1 m规格能够阻拦更多细粒物质，但由于部分障体受损，防护功能有所下降。因此，对于本研究区域1 m×1 m规格防护效益最佳，这与丁延龙等[6]研究发现2 m×2 m规格防护效果最佳结果略有不同。

4.2 不同规格PLA沙障对土壤粒度特征的影响机制

在PLA沙障障格内，细砂和中砂属于障格内的易蚀颗粒，受风沙活动影响程度较大。铺设5 a后0.5 m×0.5 m障格内细砂含量为45.15%，易蚀颗粒相较裸沙含量低，其原因是单一障格障体所能承受的积沙量较少，因此极易容易被沙埋，而在实际布设并观察后发现阻沙作用可能会有所降低，有些沙障甚至会失去阻沙作用[11]。1 m×1 m和2 m×2 m规格沙障内细砂含量分别为75.8%，80.66%，这两种规格障格内都有不同面积的侵蚀，但两者易蚀颗粒含量大致相同，由此可以看出，2种规格沙障对风沙流结构中沙粒拦截效果基本相同，且能控制过境饱和风沙流中大部分易蚀颗粒不被带走；从其分选系数上看，在铺设PLA沙障5 a后，3种规格障格内土壤颗粒分选性均变差，PLA沙障部分障体防护功能可能受到削弱，障格内颗粒分布不均匀，大量易蚀颗粒被风吹蚀，而使得大量的粗砂粒留在障格内；在风沙流中，0～10 cm高度内的沙砾约占整体的75%[12]，在这一高度下，粒径在2.32～2.74 Φ里的沙砾极易被风吹蚀[13]，但本研究样地里裸沙丘平均粒径为1.402 Φ，远大于0～10 cm高度内沙砾粒径。从土壤频率分布曲线图来看，3种规格PLA沙障每个坡位内土壤颗粒粒径均比裸沙地的颗粒粒径要细但是含量却远低于裸沙地，从分选系数上看，在铺设5 a后每种规格PLA沙障障格内分选性都变差，偏度变大。

图4 PLA沙障表层沉积物颗粒与裸沙丘差值曲线

4.3 不同规格PLA沙障粒级组分与有机质分析

土壤粒度和有机质含量是构成土壤综合体的重要组成部分，二者含量的多寡直接影响着土壤理化性质和结构的稳定性；同时，粒度组分与有机质存在着相互关联性，在我国西北沙质荒漠化地区，通常黏粒含量较多的土壤，其有机质含量也相对较高；而有机质含量较高的土壤，其物理性质和结构也相对稳定[14]；土壤有机质易与细颗粒结合形成有机-无机复合体，而且一旦与有机质结合便难以被微生物分解而且极易发生积累[15]。因此，不同规格PLA沙障铺设5 a后其土壤组分也存在差异，而有机质随着土壤组成的变化随之变化；研究区不同规格PLA沙障表层土壤粒度组成中粒级含量的变化情况大致变现为：铺设沙障前土壤粒级主要以细砂和中砂为主，随着铺设PLA沙障5 a后土壤粒级主要表现为细砂和粗砂为主。土壤有机质能改善土壤的物理和化学性质，有机质含量的增加，能提高土粒间微结构的胶结力和土壤结构体间的抵抗离散力，从而增加土壤的抗风蚀能力[16]。一般而言，土壤细颗粒对有机质的吸附和有机质的积累起着关键作用，二者存在着关联性，但不同地类的土壤粒度组成与有机质的相关性也存在差异。首先，研究区3种规格PLA沙障内表层土壤有机质表现为1 m×1 m>0.5 m×0.5 m>CK>2 m×2 m这与土壤沉积物颗粒分布曲线变化情况大致相同，虽然细砂含量相较于CK出现降低现象，但在PLA沙障内土壤粒级表现出明显变小的状态，这与有机质含量变化情况基本保持一致。从而不同规格PLA沙障内有机质含量与极细砂呈现显著正相关关系，与中砂呈现显著负相关关系。这与杨梅焕等[17]对毛乌素沙地东南缘榆阳区、靖边、横山等地土壤表层(0—5 cm)粒度与有机质含量关系分析得出的黏粒含量与有机质存在显著相关性结果有所不同。

5 结 论

(1) 研究区不同规格PLA沙障内表层土壤粒度组分及参数特征存在显著差异。PLA沙障内主要以细砂和粗砂为主，细砂含量分别为2 m×2 m(80.66%)>1 m×1 m(75.8%)>0.5 m×0.5 m(45.15%)，每种规格坡顶处细砂含量最高分别为23.39%，31.73%，41.17%，坡中粗砂含量最高分别为44.91%，36.42%，42.88%。

(2) 研究区障格内土壤分选性均变差，峰态偏离正态分布，沉积物颗粒频率分布曲线部分波段变窄，表层颗粒出现细粒化。

(3) 研究区不同规格PLA沙障内表层土壤有机质含量存在显著差异。整体上表现出1 m×1 m>0.5 m×0.5 m>CK>2 m×2 m的规律，3种规格PLA沙障内有机质含量均与极细砂含量呈现显著正相关，相关系数依次为0.996，0.999，0.975，与中砂呈现显著负相关，相关系数依次为-0.877，-0.980，-0.707；且极细砂和中砂是有机质积累的关键粒级。