-
喀斯特脆弱生态系统的退化是以强烈的人类干扰为驱动力、以植被减少为诱因、以土地生产力退化为本质、以出现类似荒漠化景观为标志的复合过程[1]。我国西南喀斯特在全球三大集中连片分布区域中面积最大、发育类型最全[2],截止2016年底,喀斯特地区石漠化总面积为1007万hm2[3]。该区广泛分布着石灰岩和白云岩,石灰岩区土粒容易聚集在岩体裂隙和地下孔隙中,基岩裸露率高;白云岩中溶蚀残余物质能相对均匀地分布于地表,基岩裸露率低[4]。白云岩的风化过程以物理风化为主,而物理崩解的岩石更有利于化学风化的进行,加上晶间孔隙均匀,推动了溶蚀作用的整体进行,导致白云岩土层厚度往往大于石灰岩区,成土速率较快,土体分布连续[5]。二者在发育形态、岩石裂隙发育程度、持水性和风化产物在地表的堆积和丢失方式都存在差异,这些因素影响着区域水土资源的空间配置,从而影响植被变化,相较于石灰岩地区,白云岩区更有利于草丛和草灌的自然恢复[6]。植被恢复是石漠化地区生态修复的根本手段,必须以蓄水、保土及造林为核心进行综合治理[7];而喀斯特生态系统植被具有喜钙、石生、旱生的特点,一旦破坏很难恢复[8]。土壤水分及其空间变化与植被密切相关,充足的土壤水分是植物正常生长的必要条件,影响着植物的存活和生长特性[9-10]。此外,土壤水分对土壤形成也有重要意义,通过制约溶质的溶解和转移、土壤养分与微生物活动来影响土壤的物理性质,也是评价土壤肥力的重要环境因素[11]。虽然我国西南喀斯特地区降雨充沛,但因该地区特殊的地上地下二元结构,降雨大量入渗,加之土壤总量少、土层浅薄且分布不连续,储水能力低,导致土壤水分成为喀斯特石漠化地区植被恢复的主要限制因子。因此,展开喀斯特地区土壤水分的相关研究十分必要[12]。
地统计学已被证明是分析土壤水分空间特征及其变异规律的有效方法之一,不仅能有效地揭示属性变量在空间上分布和变异的相关特性,还可有效地解释空间格局对于生态过程和功能的影响[13],而传统统计学方法只能描述土壤的全局性变化特征,不能描述土壤局部的变化特征,同时也忽略了每一个观测值的空间位置,不能准确的描述空间变异特性[14]。大量学者应用地统计学方法对土壤水分动态变化进行了研究,发现土壤水分均存在不同程度的空间自相关,如李春茂等[15]对峰丛洼地不同土地利用方式坡面土壤水分的研究发现,雨季土壤水分空间相关性较旱季更好,坡上未扰动的自然植被区明显高于坡下人为改造区;徐慧芳等[16]研究发现,喀斯特峰丛洼地土壤含水量均较高,呈现中等或强烈的空间自相关性。Castrignanò等[17]结合地统计学方法提出了一种评估意大利南部森林地区表土水分耗竭水平的概率方法,允许在土壤评估中剔除明显的无关因素,可有效指导林区土壤水分敏感区的划分;Landrum等[18]采用多因子协同克里金法(Kriging)确定了土壤性质和地形属性影响土壤水分格局的空间尺度,利用遥感与多元地统计学结合可以更好地制定土壤水分管理策略,但研究对象较为单一,多为喀斯特石灰岩质地区,且缺少乔灌草系统性研究。因此,本研究以喀斯特白云岩区草地、灌丛、人工林和天然次生林4种典型植被类型为研究对象,采用标准网格法,利用经典统计学和地统计学相结合的方法研究不同植被类型坡地土壤水分的空间变化及分布规律,结果有助于深入了解喀斯特白云岩石漠化地区坡地土壤水分变化过程,为退化喀斯特生态系统植被恢复和生态重建提供科学依据。
HTML
-
研究区位于贵州省安顺市普定县,其地理坐标为105º7′49″~105º58′51″ E,26º9′36″~26º31′42″ N,东西宽51.4 km,南北长40.0 km,总面积1091.6 km2。全县地势南北高、中部低,海拔高度为1044.0~1846.0 m,属北亚热带季风湿润气候区,气候温和,年平均气温15.1℃,年平均降水量1396.9 mm,且降雨多集中在5—9月,占全年降水量的75%以上,年平均日照时数1202 h,无霜期289 d,相对湿度79%,土壤主要以石灰土为主。
-
在普定县喀斯特白云岩区选择4种典型植被类型——草地(Ⅰ)、灌丛(Ⅱ)、人工林(Ⅲ)和天然次生林(IV)作为研究对象(表1),选取坡位和坡向相近的区域设立标准样地,草地和灌丛样地面积为200 m2,人工林和天然次生林样地面积为1200 m2。草地植物主要以白茅(Imperata cylindrical (L.) Beauv.)、荩草(Arthraxon hispidus (Trin.) Makino)、委陵菜(Potentilla chinensis Ser.)为主,样地中部有基岩出露(基岩裸露率约为12%),偶有火棘(Pyracantha fortuneana (Maxim.) Li)分布;灌丛样地内主要植被有山鸡椒(Pyracantha fortuneana (Lour.) Pers.)、滇鼠刺(Itea yunnanensis Franch)、算盘子(Glochidion puberum (L.)Hutch.)、红泡刺藤(Rubus niveus Thunb.)等,沿坡面中部向下有基岩出露(裸露率约为30%),间生草本和藤蔓植物;人工林样地植被主要为30 a林龄滇柏(Cuprrssus duclouxiana Hichel),其平均胸径为14.47 cm,平均树高为15.42 m;天然次生林样地中乔木主要有光叶海桐(Pittosporum glabratum Lindl.)、润楠(Machilus nanmu (Oliver) Hemsley)、青冈栎(Cyclobalanopsis glauca (Thunberg) Oersted)和盐肤木(Rhus chinensis Mill.)等,草本植物稀疏。
植被类型
Vegetation types经度(E)
Longitude纬度(N)
Latitude坡位
Slope position坡度
Slope/(°)海拔
Altitude/m土壤类型
Soil types干扰程度
DisturbⅠ 草地 105º45′38″ 26º16′44″ 中 14 1356 石灰土 中度 Ⅱ 灌丛 105º50′2.91″ 26º23′3.3″ 中 19 1340 石灰土 轻度 Ⅲ 人工林 105º48′22.77″ 26º16′40.43″ 中 21 1417 石灰土 轻度 Ⅳ 天然次生林 105º45′41.70″ 26º16′40.05″ 中 25 1393 石灰土 轻度 Table 1. Basic information of different plots
于2019年11月利用网格法将土壤分为0~10 cm(浅层)和10~20 cm(深层)进行取样。草地和灌丛采用2 m × 2 m网格,人工林和天然次生林采用5 m × 5 m网格,在每个网格中心点附近用梅花型布点法布置3个点进行混合采样,每个样点采集土样约500 g,用四分法取舍后分为两部分装袋带回实验室,一部分用来测定土壤水分含量,另一部分在阴凉处风干用来测定土壤养分。同时利用便携式GPS记录取样点中心坐标,由于取样点附近地表有岩石出露而无法取样,共计采样188个点(草地、灌丛各47个,人工林48个,天然次生林46个)。为保持土壤水分处于相对稳定水平,保证每次取样前3天均无降雨,每个样地内同样随机选择四点利用环刀取土来测定土壤密度及孔隙度。
-
用烘干法测定土壤质量含水量(%),环刀法测定土壤密度和孔隙度,设环刀体积为V,在取样前先记录每个环刀的质量M0,取样后土壤和环刀的质量为M1,然后将其去掉顶盖、保留底网放在平底容器内,向容器内加水至环刀上沿(水不漫过),待环刀内土壤充分吸水至饱和后从容器中拿出,迅速擦干环刀外面的水分后加顶盖称质量为M2,称量完毕后将环刀放置在支架上静置12 h后称质量为M3,最后打开顶盖,将环刀放入105℃烘箱中烘干12 h至恒质量,干燥后待温度降到室温加顶盖后称质量为M4,分别按照下列公式计算土壤密度、土壤总孔隙度、土壤毛管孔隙度和非毛管孔隙度。
土壤有机质用重铬酸钾容量法测定,土壤pH值用pH计(土水比为1:2.5)测定,全氮(TN)采用浓硫酸消煮-凯氏定氮法测定,全磷(TP)采用HClO4-H2SO4氧化钼锑抗比色法测定,全钾(TK)采用NaOH熔融-火焰光度法测定,水溶性有机碳用TOC分析仪测定(表2)。
土层
Soil layer/cm植被类型
Vegetation types有机质
SOM /(g·kg−1)水溶性有机碳
DOC/(mg·kg−1)pH 全氮
TN/(g·kg−1)全磷
TP/(g·kg−1)全钾
TK/(g·kg−1)土壤质地
Soil texture0~10 Ⅰ 112.00 ± 7.07 323.25 ± 11.73 7.75 ± 0.06 5.20 ± 0.77 1.12 ± 0.08 3.67 ± 1.00 轻壤土 Ⅱ 67.15 ± 6.76 201.50 ± 51.03 7.73 ± 0.42 2.91 ± 0.31 0.32 ± 0.02 8.09 ± 2.46 轻壤土 Ⅲ 68.10 ± 6.35 217.25 ± 27.00 6.63 ± 0.15 2.63 ± 0.76 0.32 ± 0.04 10.62 ± 0.93 重壤土 Ⅳ 141.00 ± 18.46 229.75 ± 51.19 6.99 ± 0.36 5.78 ± 0.87 0.43 ± 0.06 19.08 ± 1.35 中壤土 10~20 Ⅰ 96.90 ± 11.85 320.50 ± 68.56 7.78 ± 0.06 4.92 ± 0.53 1.07 ± 0.08 4.06 ± 0.54 轻壤土 Ⅱ 49.35 ± 2.78 180.50 ± 24.24 7.64 ± 0.09 2.24 ± 0.20 0.27 ± 0.03 8.98 ± 2.30 轻壤土 Ⅲ 55.25 ± 8.98 212.75 ± 60.00 6.63 ± 0.21 2.55 ± 0.33 0.30 ± 0.05 10.49 ± 1.32 重壤土 Ⅳ 113.43 ± 17.59 232.50 ± 56.45 7.04 ± 0.51 4.65 ± 0.77 0.40 ± 0.06 19.15 ± 0.99 中壤土 Table 2. The physical and chemical soil properties of different vegetation types
-
实验数据均采用域法[5]进行处理,即利用样本均值加减3倍标准差来识别特异值,在此区间外的数据均看成是特异值,分别用正常的最大和最小值代替,后续的计算均采用处理后的数据。利用Excel 2016进行数据的录入和初步整理,利用DPS 18.0进行单因素方差分析、相关性分析和半方差函数最优模型的拟合[19],利用Sigma Plot 14.0绘制直方图,运用ArcGIS 10.6的地统计模块进行克里金插值的绘制。
1.1. 研究区概况
1.2. 野外试验设计
1.3. 测定方法
1.4. 数据处理
-
不同植被类型土壤水分含量均值分布在18.59%~40.22%之间(表3),灌丛浅层土壤水分含量最低,天然次生林深层土壤水分含量最高,是灌丛浅层土壤水分的2.16倍;变异系数范围为5.34%~29.42%,其中,人工林和天然次生林呈现弱变异特征(变异系数 < 10.00%),草地和灌丛呈现中等程度变异(10.00% < 变异系数 < 100.00%)[20],表现出较强的异质性。整体来看,同种植被类型深层土壤水分含量高于浅层,但仅有天然次生林呈极显著差异(p < 0.01);在0~10 cm和10~20 cm土层,除草地和人工林外,不同植被类型下土壤水分含量呈极显著差异(p < 0.01),其含量由高到低依次为:天然次生林 > 草地 > 人工林 > 灌丛。通过对偏度和峰度的观察,结合Kolmogorov-Smirnov法进行正态分布检验,结果表明:除灌丛浅层土壤水分不服从正态分布以外,其余均服从正态分布,将灌丛0~10 cm和10~20 cm土层土壤水分进行对数转化后服从正态分布,可进行地统计分析。
土层
Soil layer/cm植被类型
Vegetation types最大值
Max/%最小值
Min/%平均值
Mean/%标准差
SD变异系数
CV/%偏度
Skew峰度
Kurtosis分布类型
Distribution types0~10 Ⅰ 29.79 13.62 23.08B 3.37 14.59 −0.751 1.368 N Ⅱ 28.18 12.25 18.59C 4.26 22.93 0.406 −1.134 N* Ⅲ 24.75 16.53 21.01B 2.08 9.94 −0.368 −0.507 N Ⅳ 43.92 28.20 35.22A 3.37 9.58 0.344 0.503 N 10~20 Ⅰ 31.23 17.08 23.55B 2.99 12.69 0.164 0.032 N Ⅱ 37.81 11.79 19.11C 5.62 29.42 0.887 0.788 N Ⅲ 25.53 19.99 22.37B 1.18 5.34 0.328 0.052 N Ⅳ 47.60 31.79 40.22A 3.36 8.35 −0.121 0.326 N 注:不同大写字母表示同一土层不同植被类型土壤水分含量差异极显著(p < 0.01),N表示服从正态分布,N*表示经过对数转化后服从正态分布。
Notes: Different capital letters indicate significant difference (p < 0.01), and N means normal distribution, N* means data is normally distributed after the logarithmic transformation.Table 3. Descriptive statistics and normal test of soil moisture in different vegetation types
-
采用普通克里金法对各土壤水分进行半变异函数拟合,按照一定标准[13]分别对指数、高斯和球面模型进行验证选出最优半变异函数模型,得到各水分空间变异特征(表4)。灌丛和草地浅层土壤水分为高斯模型,人工林和天然次生林土壤水分为指数模型,且各模型决定系数均较高(0.559~0.977),这说明模型有较好的拟合精度,可以很好地反映土壤水分的空间结构特征。草地到天然次生林浅层土壤块金值依次为0.84、0.29、0.82和1.27,深层土壤分别为0.10、0.35、0.28和0.98,4种植被类型的块金值均较小,这表明随机变异在系统中所占的比例小。灌丛土壤浅层基台值最大(14.42),人工林土壤深层基台值最小(2.05),这表明土壤水分在研究范围内灌丛变异大而人工林变异小。块金系数(块基比)是块金值和基台值的比值,能够确定空间自相关性,表示随机部分引起的空间变异占系统总变异的比例,一般用>25%、25%~75%、<75%这3个区间来表示系统具有强烈、中等和较弱的空间相关性[13]。4种植被类型块金系数变化范围为1.14%~16.73%,最大值出现在人工林土壤浅层,是最小值草地浅层土壤的14.68倍,表明样地内土壤水分呈现强烈的空间自相关关系,随机变异占总变异较小。变程表示具有相似性质斑块的空间连续性的范围,在变程以内的取样点具有空间自相关性,其余则不存在相关性[15]。从变程来看,变程随着植被的正向演替呈现先增加后减小的趋势,草地最小(4.24 m),人工林最大(16.48 m),人工林土壤水分变程明显大于其它3种类型,表明人工林内部空间自相关尺度较大,土壤水分空间连续性较好。
土层
Soil layer/cm植被类型
Vegetation types模型
Model块金值
Nuggets基台值
Sill value块金系数
Nugget sillRatio/%变程
Range/m决定系数
R20~10 Ⅰ 高斯 0.84 6.27 13.40 4.24 0.842 Ⅱ 高斯 0.29 14.42 2.01 6.37 0.977 Ⅲ 指数 0.82 4.90 16.73 16.48 0.781 Ⅳ 指数 1.27 10.81 11.75 8.82 0.906 10~20 Ⅰ 指数 0.10 8.77 1.14 5.66 0.890 Ⅱ 指数 0.35 7.83 4.47 5.79 0.801 Ⅲ 指数 0.28 2.05 13.66 15.06 0.559 Ⅳ 指数 0.98 10.64 9.21 8.18 0.930 Table 4. Related parameters of semi-variogram theoretical model of soil moisture
-
为更直观、全面地了解不同植被类型坡面土壤水分含量的分布特征,利用普通克里金法插值得出的预测图可以直接反映土壤含水量的空间分布格局,由图1可知:不同植被类型土壤水分空间分布不同,但均呈不规则的斑块或条带状分布,同一植被类型浅层、深层土壤水分具有较一致的空间分布格局。样地内浅层土壤水分空间结构复杂,土壤水分大都呈斑块分布,且草地和灌丛样地破碎化程度较高,微地貌结构复杂,导致空间连续性较差,也形成了较小的变程;样地内深层土壤水分呈条带状分布,表明深层土壤水分结构变异小,空间连续性较好。此外,草地土壤水分含量自左向右呈现先减少后增加的趋势,可能是因为样地中部基岩出露造成,同时零星分布的灌木也对土壤水分有一定的影响;灌丛土壤水分由上向下则呈递减趋势,样地上半部分土壤水分含量明显高于下半部分,这是因为样地内下半部相对于上半部有较多的裸岩分布,土壤水分蒸发较快所致;人工林土壤水分含量自上向下呈现先减少后增加的趋势,天然次生林坡面土壤水分含量呈现两边低中间高的特征,二者均随着样地内地势高低和植被分布情况发生波动。综上分析,草地、灌丛和天然次生林土壤水分含量大体沿着坡面自上向下呈递减趋势,人工林土壤水分含量自上而下沿坡面呈先减小后上升趋势,4种植被类型土壤水分最大值均出现在样地上半部,草地、灌丛和天然次生林土壤水分最小值出现在样地下半部,而人工林则出现在样地中部,这可能与样地坡度、植被和裸岩分布相关。
-
土壤密度和孔隙度反映了土壤密实程度,密度低、孔隙度越大表明土壤越疏松,土壤持水能力越强[21]。不同植被类型土壤物理性质差异明显(图2),土壤密度变化范围为0.78~1.30 g·cm−3,由大到小依次为:灌丛 > 人工林 > 草地 > 天然次生林,且灌丛明显高于其它类型土壤,与草地和天然次生林土壤密度差异显著(p < 0.05),同时浅层土壤密度小于深层,但差异不显著。土壤总孔隙度与土壤密度呈相反的规律,浅层大于深层,整体来看,土壤非毛管孔隙度随着植被正向演替呈现下降趋势,草地深层土壤最高(2.41%),灌丛深层土壤最小(0.99%),且灌丛土壤深层、浅层非毛管孔隙度差异显著(p < 0.05),深层土壤非毛管孔隙度仅为浅层土壤的59.2%,其余植被类型土壤浅层均小于深层。土壤总孔隙度和毛管孔隙度呈现相同的变化趋势,均随着植被正向演替方向而增加,即草地 < 灌丛 < 人工林 < 天然次生林,人工林和天然次生林显著高于另外2种植被类型(p < 0.05),同时其土壤总孔隙度和毛管孔隙度浅层均大于深层,且差异显著(p < 0.05)。由土壤水分与土壤孔隙度的相关性分析(图3)可知:土壤水分含量和土壤总孔隙度总体呈正比例趋势变化,土壤总孔隙度越大,土壤水分含量越高,其中,土壤浅层相关系数为0.731,深层相关系数为0.737,均表现出极显著的正相关关系(p < 0.01)。土壤水分含量和土壤密度的相关性分析(图4)表明:二者呈负相关关系,浅层土壤水分与密度的相关系数为 ‒0.820,深层土壤水分与密度的相关系数为 ‒0.725,且均为极显著负相关(p < 0.01)。