土壤侵蚀是近年来人类社会面临的最严峻的全球性生态环境问题，已严重威胁到人类生存和土地可持续发展^[1-2]，体现为土地资源破坏、江河湖库淤积、洪涝干旱加剧等，对农业生产和环境都产生巨大的威胁^[3-4]。坡耕地作为主要侵蚀地类^[5]，其水土流失严重、土层变薄、土壤肥力下降^[6-7]，同时由于其侵蚀产沙量相对较高，已经成为河流泥沙的主要策源地^[8]。贵州省作为西南山地丘陵区最具代表性的省份之一，山地和丘陵面积高达92.5%^[9]，而不合理的开垦导致大量坡耕地出现^[10]，据资料显示，2014年贵州省坡耕地面积占全省耕地面积的82.9%^[11]，该区域坡耕地土壤侵蚀情况需密切关注；因此明晰坡耕地土壤侵蚀特征，对于区域水土流失治理具有重要意义^[12]。

土壤可蚀性是评价土壤质量的指标之一，也是衡量土壤对侵蚀的敏感程度及土壤被外营力作用搬运和分散的难易程度的参数^[13]，国际上通常使用土壤可蚀性因子K值来衡量这一指标^[14]。近年来，众多学者都对不同地区的土壤可蚀性进行不同层面的研究。研究结果表明植被类型^[15]、气候^[16]、种植模式^[17]和冻融作用^[18]等都对土壤可蚀性存在一定的影响。现有研究大多数从单一影响因素出发，缺乏人为因素和自然因素相结合的研究。而岩性作为影响土壤发育的重要自然因素，制约着土壤的理化性质^[19]，不同岩性发育形成的土壤在土层厚度^[20]和土壤矿物学特征^[21]上均有所差异，因而母岩的性质控制着土壤侵蚀的空间布局和发展程度^[22]；耕犁扰动则是影响土壤性质的主要人为因素，长期耕犁会改变土壤的理化性质，破坏土壤结构，造成土壤养分流失等^[23-24]。

目前关于贵州省水土流失研究多集中于喀斯特地区，由碳酸盐岩发育形成的土壤上，但贵州地质情况复杂，还分布着很多非碳酸盐岩形成的土壤。同时关于岩性和耕犁对可蚀性影响的研究较少，特别是不同岩性区土壤由于耕犁扰动导致可蚀性差异变化尚未见报道；因此，本研究以贵州省贵阳市花溪区不同岩性发育形成的坡耕地为研究对象，通过研究不同岩性岩石发育形成的土壤在耕犁扰动下的土壤理化性质变化，进而分析岩性和耕犁扰动对坡耕地土壤可蚀性的影响，旨在为区域坡耕地水土流失治理、提高土壤抗侵蚀能力提供参考。

研究区为贵州省贵阳市花溪区(E 106°27′~106°52′，N 26°11′~26°34′)，平均海拔1 226 m，年均降雨量1 178.3 mm。全区地貌以山地和丘陵为主，以广泛出露的二叠系、三叠系石灰岩、白云岩、砂页岩及侏罗系紫色砂页岩为其主要特征，全区地带性土壤为黄壤，土壤类型主要为酸性黄壤土和石灰土，小部分紫色土和水稻土等^[25]。

根据研究区地层分布特征，选择由石灰岩、紫色砂岩和砂页岩发育形成的坡耕地为取样区域，其中石灰岩坡耕地位于花溪区小围寨，砂页岩坡耕地位于花溪区孟关乡上板村，紫色砂岩坡耕地位于花溪区孟关乡翁岩村。研究区坡耕地因常年耕作，土壤受到耕犁这一物理过程扰动影响，导致耕犁层和犁底层表现出土壤性质、颜色和土层硬度的差异；其中耕犁层表现出受耕犁扰动后的土壤性质，厚度通常介于20~40 cm之间，犁底层因未受耕犁扰动，表现出土壤原有性质，厚度不定；因此分别采集耕犁层和犁底层的土壤样品进行测定。

从2018年10月—2019年7月，共调查土壤剖面60个(石灰岩区、砂页岩区和紫色砂岩区坡耕地各20个)，采集土壤样品120个(石灰岩区、砂页岩区和紫色砂岩区坡耕地的耕犁层和犁底层各采集20个)。样地海拔为1 047~1 198 m，样点均位于坡中下部位，样点坡度为8°~25°，坡向多为南坡、东坡和西南坡，所种植的作物均为当地的典型农作物，种植方式多为连作和轮作。

土壤颗粒组成采用吸管法测定，采用美国制粒径分级标准；土壤有机质采用重铬酸钾外加热法测定，每个土壤样品的各项指标重复测定3次，最终取平均值记录结果。

土壤可蚀性因子K值计算采用Williams等在1990年建立的EPIC模型，该模型在土壤侵蚀预测模块中采用土壤有机碳和颗粒组成数据来估算土壤可蚀性因子K值^[26]。K值单位为美国制单位(short tont·acre·h/(100 ft·acre·short ton·in))，为便于与同类研究对比，将其乘以0.131 7后转为国际制单位(t·hm²·h/(hm²·MJ·mm))。

采用Excel和SPSS 19.0对试验数据进行统计分析、显著性分析和相关性分析，用Sigmaplot软件绘制图形。

由3种岩石发育形成的土壤颗粒组成情况如图 1所示。可见，3个岩性区土壤砂粒比例依次表现为：紫色砂岩>砂页岩>石灰岩；粉粒比例依次表现为：紫色砂岩>砂页岩>石灰岩；黏粒比例依次表现为：石灰岩>砂页岩>紫色砂岩。紫色砂岩区土壤和砂页岩区土壤在3种粒径比例上均无明显差异，而石灰岩区土壤粉粒和黏粒比例显著高于紫色砂岩区和砂页岩区。

图 1(Fig. 1)

不同小写字母表示相同指标不同岩性之间差异显著(P＜0.05)。下同。 Different lowercase letters indicate the same index difference significant among different lithology(P < 0.05). The same below. 图 1 土壤颗粒组成分布特征 Fig. 1 Distribution characteristics of soil particle composition

由表 1可知，研究区土壤有机质质量分数介于3.81~68.16 g/kg之间，其中紫色砂岩区、砂页岩区和石灰岩区土壤有机质质量分数均值分别为21.23、32.89和37.72 g/kg；3个岩性区土壤有机质质量分数均属于中等变异。经数据分析后得出，紫色砂岩区土壤有机质质量分数显著低于砂页岩区和石灰岩区，砂页岩区和石灰岩区无显著差异。

表 1(Tab. 1)

表 1 不同岩性区有机质质量分数的描述性统计 Tab. 1 Descriptive statistics of organic matter contents in different lithologic regions 岩性类型 Lithology category 样本数 Number of samples 平均值±标准差 Average±Standard deviation/(g·kg-1) 最大值 Maximum/ (g·kg-1) 最小值 Minimum/ (g·kg-1) 变异系数 Coefficient of variation/% 紫色砂岩Purple sandstone 40 21.23±7.01b 44.64 3.81 33.00 砂页岩Sand shale 40 32.89±14.79a 68.16 12.19 44.98 石灰岩Limestone 40 37.72±8.97a 59.25 18.89 23.79

表 1 不同岩性区有机质质量分数的描述性统计 Tab. 1 Descriptive statistics of organic matter contents in different lithologic regions

表 2描述了不同岩性发育形成的土壤，耕犁层和犁底层的土壤颗粒组成特征，可知，紫色砂岩区和石灰岩区的耕犁层和犁底层土壤各粒级比例均无显著差异，但砂页岩区耕犁层土壤黏粒比例显著小于犁底层。石灰岩区耕犁层黏粒比例显著高于砂页岩和紫色砂岩，紫色砂岩和砂页岩耕犁层黏粒比例无显著差异，但3个岩性区犁底层土壤黏粒比例两两差异显著；3个岩性区耕犁层粉粒比例无显著差异，但砂页岩区犁底层粉粒比例显著低于石灰岩区。

表 2(Tab. 2)

表 2 不同岩性区不同耕层的土壤颗粒组成特征 Tab. 2 Characteristics of soil particle composition in different plough layers in different lithologic regions 岩性类型Lithology category 耕层Soil layer 砂粒Sand/% 粉粒Silt/% 黏粒Clay/% 紫色砂岩Purple sandstone 耕犁层Plough layer 34.92±12.72Aa 37.60±10.96Aa 27.48±5.04Ba 犁底层Plowpan layer 35.36±14.03Aa 38.84±9.90ABa 25.81±5.56Ca 砂页岩Sand shale 耕犁层Plough layer 33.08±7.31Aa 39.89±5.85Aa 27.04±4.82Bb 犁底层Plowpan layer 32.51±8.62Aa 36.14±6.88Ba 31.35±10.02Ba 石灰岩Limestone 耕犁层Plough layer 13.39±6.64Ba 44.70±4.89Aa 41.91±6.01Aa 犁底层Plowpan layer 10.95±6.49Ba 44.32±6.30Aa 44.73±6.99Aa 注：不同小写字母表示相同岩性不同耕层之间差异显著，不同大写字母表示不同岩性相同耕层之间差异显著(P < 0.05)。下同。Notes: Different lowercase letters indicate the same lithology, and the difference between different soil layers is significant (P < 0.05). The same below.

表 2 不同岩性区不同耕层的土壤颗粒组成特征 Tab. 2 Characteristics of soil particle composition in different plough layers in different lithologic regions

由图 2可知，3个岩性区土壤有机质质量分数均表现出耕犁层显著大于犁底层。其中耕犁层表现出：砂页岩>石灰岩>紫色砂岩，而犁底层表现出：石灰岩>砂页岩>紫色砂岩；同时石灰岩区犁底层有机质质量分数显著高于其他2个岩性区，但耕犁层却表现出石灰岩和砂页岩区无显著差异，且显著高于紫色砂岩区。

图 2(Fig. 2)

图 2 不同岩性区不同耕层的土壤有机质质量分数特征 Fig. 2 Characteristics of organic matter contents in different soil layers in different lithologic regions

图 3是根据研究区土壤颗粒组成绘制的土壤质地三角图，可知，研究区土壤质地主要为壤土、黏壤、粉黏壤和粉黏土，其中紫色砂岩和砂页岩发育形成的土壤以砂粒和粉粒为主，因此多为黏壤和壤土，石灰岩则以粉粒和黏粒为主，所以多为粉黏土、黏土和粉黏壤。

图 3(Fig. 3)

图 3 不同岩性区土壤质地三角图 Fig. 3 Triangle of soil texture in different lithologic regions

研究区各土壤质地K值特征见表 3，可知研究区土壤一共分为8个类别，其K值大小依次表现为：粉壤>壤土>砂壤>粉黏壤>黏壤>粉黏土>砂黏壤>黏土，其中黏土K值显著低于其他土壤质地类型，粉壤K值显著高于其他土壤质地类型。

表 3(Tab. 3)

表 3 不同土壤质地可蚀性因子特征 Tab. 3 Erodibility characteristics of different soil textures 土壤质地Soil texture 样本数 Number of samples 土壤可蚀性因子K值Soil erodibility K/(t·hm2·h·hm-2·MJ-1·mm-1) 均值Average 差异性Difference 标准差Standard deviation 粉壤Silt loam 3 0.058 6 a 0.000 64 壤土Loam 24 0.056 8 b 0.001 62 砂壤Sandy loam 4 0.056 1 cb 0.000 79 粉黏壤Silty clay loam 12 0.056 0 cb 0.000 65 黏壤Clay loam 42 0.054 8 cd 0.001 50 粉黏土Silty clay 20 0.053 3 ed 0.001 27 砂黏壤Sandy clay loam 7 0.052 8 e 0.001 82 黏土Clay 8 0.050 9 f 0.003 25 总计Total 120 注: 不同小写字母表示不同质地土壤可蚀性因子K值差异显著(P < 0.05)。Notes: Different lowercase letters indicated that the K value of soil erodibility of different texture had significant difference (P < 0.05).

表 3 不同土壤质地可蚀性因子特征 Tab. 3 Erodibility characteristics of different soil textures

研究区K值描述性统计分析见表 4。由表可知，研究区K值变化范围介于0.041 9~0.059 3之间，3个岩性区K值大小依次表现为：紫色砂岩>砂页岩>石灰岩，其中最大值和最小值均出现在砂页岩区，同时砂页岩区的变异系数也是3个区中最大的，但3个岩性区土壤可蚀性因子K值变异系数均在10%以下，说明均只存在低等程度的变异性。

表 4(Tab. 4)

表 4 不同岩性区的土壤可蚀性因子K值描述性统计分析 Tab. 4 Descriptive statistical analysis of soil erodibility K in different lithologic regions 岩性类型 Lithology category 样本数 Number of samples K/(t·hm2·h·hm-2·MJ-1·mm-1) 变异系数 Coefficient of variation/% 平均值±标准差 Average±Standard deviation 最大值 Maximum 最小值 Minimum 紫色砂岩Purple sandstone 40 0.055 6±0.002 1 0.059 1 0.050 5 3.78 砂页岩Sand shale 40 0.055 1±0.003 0 0.059 3 0.041 9 5.44 石灰岩Limestone 40 0.054 8±0.001 8 0.056 6 0.049 6 3.28

表 4 不同岩性区的土壤可蚀性因子K值描述性统计分析 Tab. 4 Descriptive statistical analysis of soil erodibility K in different lithologic regions

图 4可知，3个岩性区中，耕犁层土壤K值按大小变化依次表现为：砂页岩>紫色砂岩>石灰岩，犁底层依次表现为：紫色砂岩>砂页岩>石灰岩；其中紫色砂岩区和石灰岩区的耕犁层和犁底层土壤可蚀性因子K值未表现出明显差异，但砂页岩区的耕犁层K值显著高于犁底层。从犁底层来看，紫色砂岩区土壤可蚀性因子K值显著高于砂页岩区和石灰岩区，而砂页岩和石灰岩区则无显著差异，但从耕犁层来看，紫色砂岩区土壤可蚀性因子K值与砂页岩区无显著性差异，且和石灰岩区也无显著性差异，但砂页岩区K值显著高于石灰岩区。

图 4(Fig. 4)

图 4 不同岩性不同耕层的土壤可蚀性因子K值分布特征 Fig. 4 Distribution characteristics of soil erodibility K value of different lithology and different soil layer

本研究区中，土壤共分为8种质地类型，而土壤质地是土壤的一种十分稳定的自然属性，可以反映母质来源及成土过程某些特征，因此3种岩性区的土壤质地类型存在差异。研究结果表明，基于EPIC模型下，研究区8种土壤质地中，黏土可蚀性最小，粉壤可蚀性最大，且黏粒比例越高的土壤可蚀性越小，因为黏粒具有强吸附力，土壤胶结作用较强，因此抗侵蚀能力较强。

研究区中由紫色砂岩发育形成的坡耕地土壤K值最大，最易被侵蚀，砂页岩发育形成的坡耕地土壤次之，石灰岩发育形成的坡耕地土壤可蚀性最弱。是由于紫色砂岩是岩石经风化、剥蚀、搬运在盆地中堆积形成的陆相碎屑沉积岩，其中含大量难以风化的石英颗粒，发育形成的土壤砂粒比例高，黏粒比例较少，因此多为壤土和黏壤，所以抗蚀能力弱；石灰岩主要是在浅海的环境条件下，由化学沉积作用和生物化学沉积形成的，以方解石为主的碳酸盐类岩石，石灰岩发育的土壤黏粒比例高，土壤黏重，多为粉黏土、粉黏壤和黏土，因此抗蚀性相对较强 ；砂页岩是在静水的环境中，泥沙经过长时间沉积的黏土岩类砂质页岩，发育形成的土壤各粒级颗粒分布情况与紫色砂岩相似，多为黏壤和壤土^[27-28]。因此不同岩性发育形成的坡耕地土壤由于成土过程差异可能导致土壤理化性质差异，土壤可蚀性作为土壤固有属性之一也受到岩性影响。

此外，研究结果表明，砂页岩发育形成的坡耕地土壤耕犁层和犁底层K值存在显著差异；紫色砂岩区犁底层和砂页岩区及石灰岩区犁底层土壤K值均存在显著差异，但耕犁层却表现出无显著差异；砂页岩区和石灰岩区犁底层土壤K值未表现出显著差异，但耕犁层却出现显著差异。这些结果均可以说明耕犁扰动对于土壤的影响是明显存在的，一方面耕犁层长期受到耕犁扰动，导致同一岩性发育形成的耕犁层和犁底层土壤可蚀性因子K值出现显著差异，另一方面由于耕犁扰动，耕犁层土壤养分、细粒流失，导致土壤性质发生改变，使得由不同岩性发育形成的土壤原本无显著差异，经过耕犁扰动以后出现显著差异，或者原本存在显著差异的土壤经过耕犁扰动后反而无显著差异2种情况出现。已有研究表明耕地的土壤侵蚀是所有地类中最剧烈的^[29]，所以可以认为耕犁扰动对坡耕地土壤侵蚀存在影响；已有研究表明耕犁扰动会通过影响土壤含水量^[30]、土壤内聚力^[31]、土壤紧实度、土壤密度及抗剪强度^[32]等来影响土壤侵蚀的强度和程度。总的来说，耕犁扰动对土壤性质的改变是确定的，但耕犁扰动和土壤侵蚀之间怎样相互作用还有待进一步研究。

研究区土壤可蚀性因子介于0.041 9~0.059 3之间。其中粉壤最易被侵蚀，黏土抗蚀性最强；而3种岩性区中，石灰岩发育形成的坡耕地土壤抗蚀能力较强；紫色砂岩发育形成的坡耕地土壤抗蚀能力较弱；从犁底层来看，可蚀性因子K值大小依次表现为：紫色砂岩>砂页岩>石灰岩，但从耕犁层来看，依次表现为：砂页岩>紫色砂岩>石灰岩，说明耕犁扰动对坡耕地土壤可蚀性因子存在一定的影响，而砂页岩区表现尤为明显。