滇中亚高山林地土壤团聚体结构及其对优先路径的影响

卢华兴，赵洋毅,2＊，段 旭,2，朱梦雪，涂晓云，杨文艳

(1.西南林业大学生态与环境学院，云南 昆明 650224；2.国家林业和草原局云南玉溪森林生态系统国家定位观测研究站，云南 昆明 650224)

土壤水分和溶质运移是土壤学研究的难点和热点[1]，优先流是一种常见的非均匀性土壤水分运动，是土壤水和溶质绕过土壤基质优先通过导水能力强的优先路径，快速向土壤深层和地下水运移的现象[2]。优先流在地表水、地下水和土壤水的运移和转化过程中起着关键作用，其发生可以减少地表径流和增加入渗[3]，提高土壤含水率、碳储量和养分含量，促进林木生长[4]，同时提高水资源利用效率与森林涵养水源能力[5]，对地下水补给和森林水文有显著影响[6]。土壤团聚体结构作为评价土壤特征和土壤质量的重要指标[7]，在维持土壤孔隙和通气性、提高水分渗透性和持水能力[8]、协调养分的消耗和积累、改善土壤肥力[9]、促进植物根系伸展[10]等方面起着重要作用。开展土壤优先路径和团聚体结构研究，在生态文明建设新形势下对提高水源涵养林质量、林区生态功能重建和生态环境保护具有深远意义。

土壤结构是土壤水分运动和优先流及优先路径形成的影响因素之一，牛健植等[11]根据“空间异质性=空间变异=土壤结构”的观点，认为土壤优先流的异质性主要来源于土壤结构。优先流可以发生在任何非结构性土壤层，而结构性土壤首先出现在表层[12]；吴庆华[13]证实，结构性好的土壤，优先流现象和染色深度更明显。程金花等[14]发现，土壤结构差异导致土壤含水量不同，对优先流和优先路径的形态特征和表现形式起重要作用。通常情况下，透水性不好的土层容易发生侧向流，结构性好的粉砂和黏质土壤层，优先路径以大孔隙流为主[15]。团聚体结构是最常见的土壤结构，然而，国内外学者对于土壤团聚体结构与优先流及优先路径的关系持有不同的观点，Jarvis等[16]研究发现，土壤团聚体结构在水平方向上促进土壤水分、溶质与土壤基质接触，延长部分土壤水分快速向下运移的时间，从而在垂直方向上抑制优先流及优先路径的发育；Bouma[17]提出，大团聚体结构具有毛细屏蔽作用，使土壤水分和溶质仅能沿着大孔隙通道运动而不能进入土壤基质，对优先流和基质流起分离作用；而Wuest[18]发现，优先流区域的土壤团聚体含量显著高于基质流区域。由此可见，土壤结构（尤其是团聚体结构）是土壤水分运动不可或缺的因素，对优先流和优先路径的形成、发生有直接作用，但土壤团聚体结构对优先流及优先路径如何响应仍然具有不确定性。

滇中亚高山磨盘山地处我国云贵高原、横断山地和青藏高原南缘的地理结合部，区位特征明显，是云南亚热带北部气候与南部气候的复合型生态过渡带，兼具低纬气候、季风气候和山地气候的特点[19]。植被类型多以中山半湿性常绿阔叶林为主[20]，由于长年受强风、低温和土壤脊薄影响，海拔2900 m左右地带形成罕见的高山矮曲林（高山栎林），保留着原始的森林群落面貌。此外，磨盘山作为云南新平县重要的水源涵养林区，其森林生态系统发挥着涵养水源、保持水土和调节气候等重要的生态功能[20]。目前，针对云南低纬度高海拔森林区域土壤优先路径和团聚体结构的研究鲜有报道。因此，本研究以该地区云南松(Pinus yunnanensisFranch.)林、华山松(Pinus armandiiFranch.)林、高山栎(Quercus semecarpifoliaSmith.)林和滇油杉(Keteleeria evelynianaMast.)林4种林地土壤为研究对象，通过野外染色示踪、室内图像处理和土壤团聚体测定，分析4种林地土壤团聚体结构及优先路径特征，探讨以下问题：（1）研究区4种林地类型对土壤团聚体结构、土壤优先路径是否有影响？（2）土壤团聚体结构与优先路径存在何种关系？通过以上研究，揭示滇中林地土壤优先路径形成机理，深化对磨盘山森林生态系统土壤水文过程的认识，以期为滇中亚高山森林管理、林区水源涵养林建设、天然林植被恢复与重建提供科学依据。

1 研究区概况

研究区位于云南省玉溪市新平县磨盘山国家林业和草原局森林生态系统定位观测研究站和长期科研基地内（23°46′18″～23°54′34″ N，101°12'06″～101°12'12″ E），气候类型为中亚热带高原性气候，海拔1 260.0～2 914.4 m，相对高差大；最高气温33 ℃，最低气温-2.2 ℃，年平均气温15 ℃，年日照时数为2 380 h；降雨主要集中在6—8月，年均降水量1 050 mm，干湿分明、雨热同季。土壤基岩以砂岩类岩石为主，土壤类型主要是第三纪古红土发育的山地红壤、玄武岩红壤和棕黄壤[19]，局部地区有黄壤分布，土壤整体偏酸性，pH值3.74～4.96，土层厚度1 m左右，以中厚土壤层为主。滇中亚高山磨盘山是以云南特有中山半湿性常绿阔叶林为主的重要林区，森林资源丰富，植被类型多样，且随海拔上升呈垂直分布特征，有亚热带常绿阔叶林、亚热带中山针阔混交林、针叶林（以华山松、滇油杉、云南松纯林为主）和高山栎林等森林类型[21]，是云南森林群落的重要组成部分。

2 研究方法

2.1 样地布设

在研究区4种林地内设置面积为100 m×100 m的固定样地各1块，分别在各个固定样地内选择植被均匀、林分结构一致、地形相似、土壤形成母质一致的地段布设4块20 m×20 m的样方，在每个样方内选择3个试验点开展染色示踪试验，同时，在各个样方内利用5点取样法分层采集原状混合土样，采样深度为0～50 cm，以土层深度10 cm为间隔，混合土样带回实验室风干后，用于测定土壤团聚体结构，样地基本状况见表1。

表1 研究区4种林地基本状况Table 1 Basic conditions of four forest lands in the study area

2.2 染色示踪及优先路径提取

于2020年7月1日至9月10日开展染色示踪试验，以浓度4 g·L-1的亮蓝溶液20 L（模拟当地24 h累计降雨量40 mm）作为染色示踪剂，试验过程参考卢华兴等[22]前期的研究报道。示踪试验所获取的染色图像通过Photoshop CS5、Image Pro Plus 6.0和Surfer 15等图像处理软件进行解译，以获取不同影响半径范围下优先路径的位置分布特征和数量特征，优先路径的提取过程参考卢华兴等[23]在滇南地区的研究方法。

2.3 土壤团聚体结构测定与参数计算

原状混合土样自然风干后除去粗根及石砾等杂物，将大土块沿着自然断裂面剥离为1 cm大小的土块，参照林业行业标准LY/T 1 227—1999，用机械筛分法测定＞5.0、2.0～5.0、1.0～2.0、0.5～1、0.25～0.5、＜0.25 mm 6个粒级水稳定团聚体含量，然后计算粒级＞0.25 mm水稳定性团聚体含量(WSA)、平均质量直径(MWD)、几何平均直径(GMD)、平均质量比表面积(MWSSA)、结构体破坏率(PAD)、分形维数(D)等团聚体结构特征参数[24-27]。

2.4 数据分析

所有数据用Excel 2016软件进行初步整理和分析；用SPSS 20.0统计软件进行Duncan差异性分析、双因素方差分析；用Origin 2021软件进行主成分分析和Pearson相关分析；用DPS软件进行灰色关联分析，然后计算灰色关联度[28]和耦合度[29]。

3 结果与分析

3.1 土壤团聚体结构特征

3.1.1 不同林地及粒级土壤团聚体特征变化 由图1可知：6个粒级水稳定性团聚体中，高山栎林地＞5.0 mm水稳定性团聚体含量最大，含量为60.65%；华山松林地1.0～2.0、0.5～1.0 mm水稳定性团聚体含量最大，而其它2.0～5.0、0.25～0.5、＜0.25 mm 3个粒级水稳定性团聚体含量最大的是云南松林地。比较同一林地不同粒级水稳定性团聚体含量，研究区4种林地土壤水稳定性团聚体多分布于＞5.0、2.0～5.0 mm粒级中，大粒级的水稳定性团聚体占主导地位。总体看，研究区4种林地土壤水稳定性团聚体含量变化趋势表现为随粒级的减小基本呈减小的趋势，且4种林地不同粒级土壤团聚体含量差异显著 (P＜0.05)，在同一粒级范围内，林地类型对团聚体含量的影响差异显著(P＜0.05)。

图1 研究区4种林地土壤水稳定性团聚体含量Fig.1 Content of soil water stable aggregates in four forest lands in the study area

由表2可知：不同林地类型对土壤团聚体结构影响显著(P＜0.05)，且不同林地类型对各个土壤团聚体结构特征指标的影响作用不同。研究区4种林地土壤团聚体结构差异大，较其他3种林地，高山栎林地MWD、GMD、D值最大，而WSA、MWSSA、PAD等特征指标占有优势的林地依次分别是滇油杉林、云南松林、华山松林，其中，4种林地WSA平均值均大于90%。WSA值越高，土壤结构越稳定，土壤抗蚀性越好[24]，表明研究区4种林地土壤具有稳定的土壤结构和良好的抗蚀性。MWD、GMD可反映不同粒级团聚体的分布状况，MWD和GMD值越大，土壤平均粒径团聚度和稳定性越好[25]，由此说明，高山栎林地有助于增加土壤团聚度，提高土壤稳定性。分形维数是土壤结构几何形状的参数，可以反映水稳定性团聚体含量对土壤结构与稳定性的影响趋势，即土壤分形维数越小，土壤团聚体结构的稳定性越好[27]。根据分形维数大小可认为，4种林地土壤团聚体结构的稳定性依次为：华山松林＞滇油杉林＞云南松林＞高山栎林。综上所述，林地类型对土壤团聚体结构稳定性的影响程度较大。

表2 研究区4种林地土壤团聚体结构特征Table 2 Structural characteristics of soil aggregates in four forest lands in the study area

3.1.2 林地类型与土层深度对团聚体结构的影响

由表3可知：林地类型和土层深度对土壤团聚体结构产生的影响作用不同，林地类型和土层深度对WSA、MWD、GMD、MWSSA、PAD和D均影响极显著(P＜0.01)；林地类型与土层深度的交互效应对WSA、MWD、GMD、MWSSA和PAD的影响极显著(P＜0.01)，而对D的影响显著(P＜0.05)。

表3 林地类型、土层深度与土壤团聚体结构的双因素方差分析Table 3 Two factor ANOVA of forest land type,soil depth and soil aggregate structure

3.1.3 土壤团聚体结构与各粒级的相关分析 由表4可知：各个粒级水稳定性团聚体含量与团聚体结构特征指标存在相关关系。粒级＞5.0 mm水稳定性团聚体含量与MWD、GMD呈极显著正相关(P＜0.01)，1.0～2.0 mm水稳定性团聚体含量与MWD呈极显著负相关(P＜0.01)，与GMD、D呈显著负相关(P＜0.05)，但与PAD呈显著正相关(P＜0.05)。0.5～1.0 mm水稳定性团聚体含量与WSA、PAD呈极显著正相关(P＜0.01)，与MWD、D呈极显著或显著负相关。0.25～0.5 mm水稳定性团聚体含量与WSA、MWD、GMD呈极显著负相关(P＜0.01)，但与MWSSA、D呈极显著或显著正相关。粒级＜0.25 mm水稳定性团聚体含量与WSA、PAD呈极显著负相关(P＜0.01)，与MWSSA、D呈极显著正相关(P＜0.01)。用相关系数来评价不同粒级水稳定性团聚体含量对团聚体结构特征指标的贡献，粒级2.0～5.0 mm水稳定性团聚体含量对所有土壤团聚体结构影响不大；粒级＞5.0、1.0～2.0、0.25～0.5 mm水稳定性团聚体含量对MWD、GMD的贡献较大，相关系数的绝对值大于0.56；粒级0.5～1.0、＜0.25 mm水稳定性团聚体含量对PAD、D的贡献较大，相关系数的绝对值大于0.64。由此说明，不同粒级水稳定性团聚体含量是团聚体结构稳定性的重要影响因素。此外，各个团聚体结构特征指标之间也存在相关性，WSA与MWSSA、D、PAD均呈极显著相关(P＜0.01)，MWD与GMD、MWSSA，GMD与MWSSA以及PAD与D均呈极显著相关(P＜0.01)，说明上述各个指标适用于评价团聚体结构特征。

表4 不同粒级水稳定性团聚体含量与土壤团聚体结构的相关分析Table 4 Correlation analysis between water stable aggregate content of different particle sizes and soil aggregate structure

3.2 土壤优先路径特征

3.2.1 不同林地土壤优先路径的位置分布 由图2可看出：云南松林地染色半径＞10 mm的优先路径在0～50 cm土层中的连通性最好，染色半径≤1 mm的优先路径连通性最差；华山松林地染色半径1～2.5、2.5～5、5～10和＞10 mm的优先路径连通性相似，相比以上4个染色半径，染色半径≤1 mm优先路径的连通性最差。滇油杉林地和高山栎林地在0～50 cm土壤优先路径的连通程度具有一致性，即染色半径5～10、＞10 mm的优先路径连通性较好，染色半径≤1、1～2.5 mm的优先路径连通性相对较差。研究区内4种林地土壤优先路径的连通性存在一定差异性，主要原因是林木根系[23]和土壤结构[16-18]显著影响其对应染色半径范围内优先路径的连通性。

图2 研究区4种林地土壤优先路径的位置特征Fig.2 Location characteristics of soil priority paths of four forest lands in the study area

比较4种林地同一染色半径下不同土层深度优先路径的连通性发现：华山松林地染色半径≤1、1～2.5 mm的优先路径在20～30 cm的土层中连通性较好，而染色半径5～10、＞10 mm的优先路径在30～40、40～50 cm的土层中连通性最好。滇油杉林地染色半径≤1、1～2.5、2.5～5、5～10 mm的优先路径在0～20 cm的表层土壤中连通性较好，染色半径＞10 mm的优先路径在30～40 cm土层连通性最好。高山栎林地土壤优先路径的连通性在30～40、40～50 cm土层中最佳，与云南松林地基本一致。整体而言，以华山松林地5个染色半径的优先路径连通性最好，其次是云南松林地和高山栎林地，滇油杉林地的优先路径连通性最差。

3.2.2 不同染色半径土壤优先路径的数量特征 由表5可知：研究区4种林地中，高山栎林地优先路径数量最多（3974），其次是云南松林地和滇油杉林地，华山松林地的土壤优先路径数量最少（2770）。水平方向上，4种林地同一土层不同染色半径的优先路径数量由大到小均表现为：≤1、1～2.5、2.5～5、5～10、＞10 mm，即随着染色半径的增大，同一土层优先路径的数量逐渐减小，说明林地土壤中小径级的土壤孔隙在数量上占有一定的优势，大径级土壤孔隙可能是小径级土壤孔隙演化形成的。竖直方向上，4种林地5个染色半径的优先路径数量在土层深度方向的变化存在差异，4种林地土壤优先路径总数以及滇油杉林地、高山栎林地不同染色半径的优先路径数量均随土壤深度的增加而降低，表层土壤优先路径数量显著大于深层，说明土壤优先路径数量随土壤深度增加而降低[23]。由于林地土壤空间异质性，云南松林地和华山松林地不同染色半径的优先路径数量均随土壤深度的变化无明显的规律性。

3.3 土壤团聚体结构与优先路径的关系

3.3.1 土壤优先路径数量与团聚体结构的相关分析 由表6可知：各个染色半径范围内（≤1、1～2.5、2.5～5、5～10和＞10 mm）优先路径数量之间及与优先路径总数均呈极显著正相关(P＜0.01)。土壤优先路径总数与MWD、GMD、PAD均呈显著正相关(P＜0.05)，与MWSSA呈显著负相关(P＜0.05)，而与D呈极显著正相关(P＜0.01)。土壤优先路径数量与MWD、GMD、MWSSA、PAD、D等多个土壤团聚体结构特征指标具有相关性，说明土壤优先路径的形成在一定程度上受土壤团聚体结构的影响。

表6 不同影响半径优先路径数量与土壤团聚体结构的相关分析Table 6 Correlation analysis between the number of priority paths with different influence radii and soil aggregate structure

3.3.2 优先路径影响因素的主成分分析 通过主成分分析对6个土壤团聚体结构特征指标进行主成分提取，由表7可知：主成分1、2、3的初始特征值均＞1，累积方差贡献率为94.643%，说明这3个主成分基本能反映不同林地土壤团聚体结构特征的有效信息占94.643%，仅有5.357%的信息损失。在主成分1中，MWD、GMD、MWSSA载荷值的绝对值较大，均＞0.87，其中，MWSSA的载荷值为负值，其它指标的载荷值均为正值，说明MWSSA越小，土壤团聚体结构越好。在主成分2中，WSA、PAD、D载荷值的绝对值＞0.79，其中，WSA、PAD的载荷值为负值。主成分3中，PAD、D载荷值的绝对值较大。比较3个主成分所有指标载荷值的绝对值，＞0.8的指标依次为GMD(0.955) ＞MWSSA (-0.930) ＞MWD (0.875) ＞PAD (-0.871)，说明这6个土壤团聚体结构特征指标中，GMD、MWSSA、MWD、PAD能较好地评价土壤团聚体结构特征，对优先路径形成的贡献度较大。

表7 土壤团聚体结构特征的主成分分析Table 7 Principal component analysis of structural characteristics of soil aggregates

3.3.3 土壤团聚体结构和优先路径的灰色关联分析 将土壤优先路径数量作为母序列，土壤团聚体结构特征作为子序列，即子序列为：土壤团聚体结构特征 ={WAS，MWD，GMD，MWSSA，PAD，D}，进行灰色关联分析，由图3可知：各个土壤团聚体结构特征指标与优先路径数量的关联系数介于0.601～0.782，均值为0.672，属于中等或较强关联，表明不同林地土壤优先路径数量与土壤团聚体结构特征指标间有一定的关联作用，且联系和反馈较为紧密。云南松、华山松、滇油杉林地PAD的关联系数最大，分别为0.782、0.701、0.742，属于较强关联，说明PAD对土壤优先路径的影响和反馈作用较明显；高山栎林地MWSSA的关联系数达0.771，说明高山栎林地MWSSA是影响土壤优先路径形成的关键土壤团聚体结构特征指标。将4种林地土壤优先路径数量与各个团聚体结构特征指标的平均关联度按大小排序依次为0.742(PAD) ＞0.683 (GMD) ＞0.675 (MWD) ＞0.666 (D) ＞0.660 (WSA) ＞0.630 (MWSSA)，均＞0.62，关联度为中度，说明这类指标是影响土壤优先路径的重要土壤团聚体结构指标，与上文主成分分析的结论基本一致。

图3 研究区4种林地土壤团聚体结构的灰色关联度Fig.3 Grey relational degree of soil aggregate structure of four forest lands in the study area

3.3.4 土壤团聚体结构和优先路径的耦合度分析

考虑土壤团聚体结构与优先路径的相互作用和耦合程度，利用耦合度模型计算不同林地土壤团聚体结构特征指标和优先路径的耦合协调程度。4种林地土壤团聚体结构和优先路径的耦合度介于0.662～0.728，平均耦合度为0.683，系统整体水平属于弱协调，说明土壤团聚体结构与优先路径的形成并未达到最佳的耦合协调状态，可能是4种林地的植被类型存在差异，植被在生长过程中根系向四周延伸，破坏土壤结构的稳定性，不利于形成良好的土壤团聚体结构。4种林地类型的系统耦合协调程度按耦合度C值排序为：0.728（云南松林）＞0.681（华山松林）＞0.663（高山栎林）＞0.662（滇油杉林），其中，云南松林地耦合度最大，属于中度协调，其粒级2.0～5.0 mm水稳定性团聚体和MWSSA显著大于其他3种林地，且云南松林地优先路径总数量显著大于华山松、滇油杉林地，说明大粒级的团聚体提高土壤结构的复杂和变异程度，增加土壤孔隙数量，有利于优先流现象的发生和优先路径数量的增多，从而使得优先路径和土壤团聚体结构性能稳步改善，协调发展。滇油杉林地土壤层中存在不同粒级的石砾，破坏土壤团聚体结构，阻碍优先路径的形成，导致系统耦合协调程度下降。

4 讨论

4.1 林地类型对土壤团聚体结构的影响

团聚体结构是土壤结构的基本反映，团聚体的粒径分布和稳定性对土壤肥力状况、结构稳定性及抗蚀性有重要影响[8]。研究表明，研究区4种林地土壤水稳定性团聚体差异明显，大粒级团聚体含量占一定优势，且随粒级减小团聚体含量呈减小的趋势，这可能与研究区林地土壤缺少人为活动有关。王轶浩等[30]认为，土壤团聚体结构在林木生长和植被恢复过程中依靠自然演变，小粒级团聚体逐渐向大粒级团聚体转变；吴秀坤等[31]发现，林地表层土壤地被植物和枯落物的覆盖，增加表层土壤入渗率，有效防止雨滴对土壤大粒级团聚体的冲击，有利于大团聚体的形成和保持。4种林地粒径＞0.25 mm水稳定性团聚体含量均大于90%，说明研究区林地具有良好的土壤结构和抗蚀性[24]。林地土壤表层枯落物被分解后，归还土壤转化为有机质，有效促进土壤颗粒之间的胶结作用[25]，增加土壤团聚体的数量和稳定性[32]。聂浩亮等[33]发现，林分类型对土壤有机质含量影响显著，林分类型影响凋落物产量和分解速率，进而造成土壤有机质含量存在差异[34]。研究区4种林分中，华山松和滇油杉林地土壤有机质含量较高，且这2种林分土壤有机质含量差异规律与团聚体结构稳定性一致。孙金兵等[35]证实，土壤有机质含量与水稳定性团聚体含量显著相关，土壤颗粒与有机质的胶结作用可以改善土壤结构，有利于土壤团聚体的形成。土壤团聚体结构稳定性受土层深度和植被类型的影响，魏霞等[26]研究发现，土层深度变化对土壤团聚体结构特征存在显著影响，随土层深度的增加，土壤结构稳定性逐渐减弱，且植被类型对土壤结构的影响比土层深度更大，本研究结果与之基本一致。

4.2 林地类型对土壤优先路径的影响

研究区4种林地土壤优先流现象明显，但由于受土壤空间异质性[12]和优先流自身特征[11]（环绕性、快速穿透性和侧向入渗等）影响，优先路径的发育程度存在差异。从优先路径的位置分布看，华山松林地优先路径的连通性最好，其次是云南松林地和高山栎林地，滇油杉林地连通性最差；从优先路径数量看，数量大小依次为：高山栎林（3 974）＞云南松林（3513）＞滇油杉林（3171）＞华山松林（2770），优先路径连通性较好的林地土壤在数量上不一定占有优势，相反，优先路径数量最大的林地土壤连通性并不是最佳，这可能是染色图像处理过程中造成的误差，导致同一林地土壤优先路径的连通性和数量不一致。研究区4种林地优先路径发育均较好，程竞萱等[36]研究发现，在重庆四面山土壤优先流发育程度为天然针阔混交林地＞天然阔叶林地＞农地；Zhang等[37]提出，林地土壤比果园土壤有更多的优先流动路径，均说明林地土壤优先流具有优势性。研究区4种林地表层土壤优先路径数量较多，显著多于深层土壤，其原因主要还是表层土壤团聚体结构稳定性好，枯落物的分解增加表层土壤有机质含量和团聚体数量[38]。随着土层深度增加，土壤团聚体结构稳定性和水分下渗能力减弱，优先路径数量逐渐减少，这与卢华兴等[23]在滇南地区的研究结果一致。华山松林地土质疏松，大孔隙数量多，有利于水分在土壤层中垂直运移，且其土壤表层枯落物和腐殖质较厚，增加土壤有机质含量，形成更多的团粒结构，提高土壤孔隙的通透性和优先路径的连通性[39]。挖掘染色剖面时发现，滇油杉林地土壤中含有不同粒级的石砾，这可能是导致滇油杉林地土壤优先路径连通性较差的原因。

4.3 土壤团聚体结构对优先路径的影响

通过主成分分析、相关分析和灰色关联分析，表明土壤团聚体结构影响土壤优先路径的形成。华山松林地表层土壤土质疏松，云南松林地深层土壤黏性大、透水性差，土壤结构可能是导致云南松和华山松林地优先路径数量随土壤深度分布不规律的主要原因。程金花等[14]认为，土壤结构是形成优先路径的基础，土壤结构通过改变土壤含水率间接影响优先路径的形成；Wuest[18]研究发现，基质流区域的土壤团聚体结构含量明显比优先路径附近区域低。土壤优先路径数量与各个土壤团聚体结构特征指标的关联水平属于中等或较强关联，其中，PAD的关联度平均值达0.742，MWSSA的关联度最小，说明PAD最能反映土壤团聚体结构对土壤优先路径的影响程度。魏霞等[26]认为，PAD是反映土壤团聚体受到侵蚀后崩解、破坏的程度，而MWSSA反映土壤团聚体质量的准确性不高，不利于土壤团聚体结构特征的表达。4种林地土壤团聚体结构和优先路径的平均耦合度为0.683，系统整体水平属于弱协调，说明土壤团聚体结构并不是优先路径形成的限制因素，优先路径的形成可能受其他因素（如林木根系）的影响。程金花等[40]认为，林地土壤中的林木根系能够产生孔隙，促进优先路径的形成，同时也能使土壤团聚体破碎，堵塞土壤孔隙而不利于优先路径的形成。

5 结论

(1) 4种林地土壤具有良好的团聚体结构，其稳定性依次为：华山松林＞滇油杉林＞云南松林＞高山栎林，土壤团聚体结构受林地类型、土层深度的影响显著(P＜0.05)，且与不同粒级水稳定性团聚体有相关性。

(2) 林地土壤优先流具有优势性，华山松林优先路径的连通性优于云南松林、高山栎林和滇油杉林，0～10 cm表层土壤土质疏松，优先路径数量较多，随着土层深度增加，优先路径数量逐渐减少。

(3) WSA、MWD、GMD、MWSSA、PAD、D能较好地评价土壤团聚体结构特征，与优先路径数量的关联度大于0.62，属于中度关联，土壤团聚体结构是影响优先路径形成的重要因素。

(4) 不同林地类型土壤优先路径与土壤团聚体结构的耦合度为0.683，系统整体水平属于弱协调，土壤团聚体结构不是优先路径形成的限制因素，4种林地耦合协调程度依次为：0.728（云南松林）＞0.681（华山松林）＞0.663（高山栎林）＞0.662（滇油杉林）。