-
喀斯特地区长时间强烈的岩溶作用造成了水土资源不协调的地表、地下双层空间结构,导致地表水易流失、地下水深埋,土层浅薄、土壤持水性能差,致使该区域生态系统十分脆弱,石漠化问题突出,植被恢复困难[1-2]。断陷盆地作为喀斯特石漠化综合治理的重要类型区,主要分布在滇东-攀西,受西南印度洋季风影响,该区季节性干旱严重,有明显的雨季和旱季,即使在降水较多的雨季,由于土层浅薄、土壤持水性能差等原因,也只能在每次降雨后维持短时间的水分充足状态,水资源矛盾极为突出[3]。目前,关于喀斯特区土壤水分动态变化的研究主要集中在岩溶高原、峰丛洼地等区域。如李菲[4]研究了贵州普定生态示范区不同植被类型土壤水分的变化特征,发现原生林与灌丛地的土壤含水量高于次生林地与荒草地,而土壤密度变化与之相反。张川等[5]通过比较桂西北喀斯特地区典型灌丛与草灌坡地旱季表层土壤水分含量空间变异特征,发现灌丛坡地的土壤水分含量与变异系数明显高于草灌坡地,两种坡地土壤水分与降雨量波动变化趋势相同。颜蒙蒙等[6]分析了贵阳市花溪区2013年和2015年土壤水分的动态变化规律,认为土壤水分呈季节性动态变化,土壤剖面水分的垂直变化具有明显层次性。
特殊的地质结构和水文地质条件,导致喀斯特断陷盆地降雨、蒸发、水分分配等均与岩溶高原、峰丛洼地等其它喀斯特类型区存在显著差异[7-8]。目前,对该区不同植被恢复模式土壤水分时空变化规律的研究鲜有报道。本研究通过对喀斯特断陷盆地6种植被恢复模式不同土层深度土壤体积含水量进行全年连续监测,研究林地土壤水分时空变异规律,分析降雨前后各林地土壤水分变化及持续干旱条件下土壤水分衰减特征,以期为该地区植被恢复与石漠化治理工作提供参考。
HTML
-
研究区位于西南喀斯特断陷盆地区云南省建水县九标小流域范围内(地理坐标102°54′00″~102°54′55″ E,23°36′50″~23°37′30″ N),该流域为东西走向,主沟道长约1.32 km,流域面积约为1.12 km2(图 1),小流域内地貌为典型的岩溶断陷盆地地貌。海拔高度1 350~1 700 m,处于低纬度地区,属南亚热带季风气候。受西南印度洋季风影响,形成明显干湿季节,季节性干旱严重,旱季(11月—次年4月)空气干燥,降雨稀少。雨季(5—10月)湿度大,降雨较多。该地区年均气温19.8℃,年均地温20.8℃,年均相对湿度72%,年均日照时数2 322 h,年均降水量805 mm,全年无霜期307 d。
Figure 1. Spatial distribution map of six vegetation restoration patterns in Jiubiao small watershed
选取小流域中6种植被恢复模式作为研究对象,除天然次生灌丛为1970年开始自然恢复植被外,其余5种皆为1985年人工造林恢复植被,恢复时间均已30年以上,土壤结构稳定且受干扰较少。其中3种为乔木林,即桉树林(直杆蓝桉(Eucalyptus maideni F. V. Muell.))、冲天柏(Cupressus duclouxiana Hickel)林和马尾松(Pinus massoniana Lamb.)林,林地土层深厚,分别设置40 m×40 m的标准样地。桉树林具有乔-灌-草层次结构,乔木均高14.20 m,灌木均高1.22 m,草本均高0.30 m,主要植被有直杆蓝桉、车桑子(Dodonaea viscosa L.)和紫茎泽兰(Eupatorium adenophora Spreng.)等。冲天柏林地岩石裸露率25%,具有乔-灌-草层次结构,乔木均高11.78 m,灌木均高0.84 m,草本均高0.21 m,主要植被有冲天柏、小石积(Osteomeles anthyllidifolia Lindl.)和铁仔(Myrsine africana Linn.)等。马尾松林具有乔-灌-草层次结构,乔木均高13.42 m,灌木均高0.45 m,草本均高0.18 m,主要植被有马尾松和刺芒野古草(Arundinella setosa Trin.)等。另外3种为灌丛(天然次生灌丛、低盖度车桑子灌丛和高盖度车桑子灌丛),土层较薄,仅30 cm左右,分别设置20 m×20 m的标准样地。天然次生灌丛具有灌-草层次结构,灌木均高2.71 m,草本均高0.22 m,主要植物有假虎刺(Carissa spinarum L.)、薄叶鼠李(Rhamnus leptophylla Schneid.)、铁仔(Myrsine africana Linn.)和石山羊蹄甲(Bauhinia comosa Craib)等。低盖度车桑子灌丛岩石裸露率70%,具有灌-草层次结构,灌木均高2.20 m,草本均高0.20 m,主要植被有车桑子(Dodonaea viscosa L.)和硬秆子草(Sclerochloa dura Steud.)等。高盖度车桑子灌丛岩石裸露率60%,具有灌-草层次结构,灌木均高2.34 m,草本均高0.27 m,主要植被有车桑子(Dodonaea viscosa L.)和荩草(Arthraxon hispidus Thunb.)等。样地其他信息如表 1所示。
样地Plot 地理坐标
Geographical coordinates坡度
Slope/(°)坡向
Aspect海拔
Elevation/m造林时间
Afforestation time树高
Height/m植被盖度
Vegetation cover/%植株密度
Plant density/(tree·hm-2 )冠幅
Crown diameter/(m×m)桉树林
Eucalyptus maideni102°54′47″ E
23°37′19″ N2 北偏东North east 1 424 1985 14.20 50 1 388 2.06×1.96 冲天柏林
Cupressus duclouxiana102°54′34″ E
23°37′23″ N5 北North 1 466 1985 11.78 35 2 731 1.60×1.54 马尾松林
Pinus massoniana102°54′23″ E
23°37′11″ N5 西West 1 398 1985 13.42 65 2 128 2.63×2.66 天然次生灌丛
Secondary shrubs102°54′38″ E
23°37′12″ N22 北偏东North east 1 475 1970 2.71 90 4 550 1.86×1.81 低盖度车桑子灌丛
low-coverage Dodonaea viscosa102°54′46″ E
23°37′08″ N15 西West 1 509 1985 2.20 45 1 575 1.48×1.44 高盖度车桑子灌丛
high-coverage Dodonaea viscosa102°54′47″ E
23°37′03″ N12 西偏南South west 1 573 1985 2.34 60 1 850 1.53×1.55 Table 1. Basic profiles of six vegetation restoration patterns
-
采用土壤水分传感器CS616(Campbell Scientific Inc., USA)记录土壤水分数据。3种乔木林地分4层剖面(0~10, 10~30, 30~50, 50~70 cm)对土壤体积含水量进行实时监测。3种灌丛林地分3层剖面(0~10, 10~20, 20~30 cm)对土壤体积含水量进行实时监测。测定前对CS616测定的数据进行了相应校准,仪器测定结果可代表研究区土壤水分变化情况。采用CR800(Campbell Scientific Inc., USA)进行数据采集,采集频率设置为10 min,日记录数据的平均值作为本研究的分析数据,数据记录为期1年(2016年5月—2017年4月)。
-
采用自动气象观测系统(WatchDog 2900ET, Spectrum Technologies Inc., USA)同步测定降雨量与蒸散发(ET0),2016年5月—2017年4月期间降雨总量为876.3 mm,蒸散发全年为1 397.88 mm,降雨量及蒸散发分布情况见图 2。
-
变异系数cv表征随机变量的离散程度,其定义为:
式中:s为标准差,x为土层样本均值。cv≤10%属于弱变异,10%≤cv≤100%属于中等变异,cv≥100%属于强变异[9]。本文中的变异系数为土壤水分在统计周期内的波动情况。
所有数据统计在Excel 2016中进行,采用普通克里金插值方法生成土壤水分随时间消退过程的空间分布图,数据分析与作图在R 3.3.3 for Windows和Golden Software Surfer 8.0软件中进行。
1.1. 研究区概况
1.2. 研究方法
1.2.1. 土壤含水量的测定
1.2.2. 降雨量及蒸散发测定
1.3. 数据处理
-
根据研究区的气候特点,分两个时间段来研究不同植被恢复模式土壤水分垂直变化特征:雨季(2016年5月—2016年10月),该时段降雨量为729.7 mm,占全年降雨量的83.27%,旱季(2016年11月—2017年4月),该时段降雨量为146.6 mm,占全年降雨量的16.73%。综合分析不同植被恢复模式在雨季和旱季各土层间土壤水分含量(表 2和表 3),结果表明:同一植被类型不同土层深度雨季和旱季土壤体积含水量垂直变化趋势相同,且雨季土壤水分含量显著高于旱季(P < 0.05)。随土层深度的增加,雨季和旱季土壤体积含水量的差值有减小的趋势,以桉树林和高盖度车桑子灌丛为例,桉树林0~10、10~30、30~50和50~70 cm土层深度对应的差值分别为:7.09%、6.90%、3.39%和1.31%。高盖度车桑子灌丛0~10、10~20和20~30 cm土层深度对应的差值分别为:9.09%、4.64%和2.15%。不同植被恢复模式土壤体积含水量变化趋势各异,其中3种乔木林地雨季和旱季在0~70 cm剖面范围内的垂直分布规律为:桉树林土壤体积含水量随土层深度的增加而减小,冲天柏林和马尾松林则表现出波动变化,冲天柏林土壤水分先增大后减小,马尾松林则先减小后增大再减小。3种灌丛雨季和旱季在0~30 cm剖面范围内的垂直分布规律为:天然次生灌丛土壤体积含水量先增大后减小,低盖度和高盖度车桑子灌丛土壤体积含水量随土层深度的增加而减小。3种乔木林(0~70 cm土层)和3种灌丛(0~30 cm土层)在雨季和旱季的平均土壤体积含水量均表现为桉树林>天然次生灌丛>冲天柏林>马尾松林>高盖度车桑子灌丛>低盖度车桑子灌丛。另外,从表 2和表 3得出,无论雨季还是旱季,6种植被恢复模式各土层土壤水分均为中等变异(Cv:12.88%~28.66%)。
土层深度
Soil depth/cm0~10 10~30 30~50 50~70 土壤水分
Soil moisture/%变异系数
Variation coefficient/%土壤水分
Soil moisture/%变异系数
Variation coefficient/%土壤水分
Soil moisture/%变异系数
Variation coefficient/%土壤水分
Soil moisture/%变异系数
Variation coefficient/%桉树林
Eucalyptus maideni雨季
Rainy season37.26±7.91 21.23 34.71±7.92 22.82 26.30±4.65 17.68 20.20±3.70 18.32 旱季
Dry season30.17±6.54 21.67 27.81±5.81 20.89 22.91±5.19 22.65 18.89±3.90 20.65 冲天柏林
Cupressus duclouxiana雨季
Rainy season28.99±6.27 21.63 32.55±6.30 19.35 24.98±4.08 16.33 21.16±4.25 20.09 旱季
Dry season21.50±4.83 22.47 24.28±5.72 23.56 20.40±5.37 26.32 19.75±5.28 26.73 马尾松林
Pinus massoniana雨季
Rainy season25.86±4.55 17.59 21.75±3.26 14.99 22.98±3.91 17.01 19.50±4.87 24.97 旱季
Dry season20.02±4.28 21.38 18.22±2.44 13.39 18.93±3.33 17.59 16.33±2.57 15.74 注:表中数据形式为均值±标准差,下同。
Note: The data form in the table is the mean ± standard deviation. The same below.Table 2. Soil moisture content and variation coefficient of three forest land
土层深度
Soil depth/cm0~10 10~20 20~30 土壤水分
Soil moisture/%变异系数
Variation coefficient/%土壤水分
Soil moisture/%变异系数
Variation coefficient/%土壤水分
Soil moisture/%变异系数
Variation coefficient/%天然次生灌丛
Secondary shrubs雨季Rainy season 25.05±5.22 20.84 30.70±4.63 15.08 25.87±4.39 16.97 旱季Dry season 22.63±5.67 25.06 27.64±4.39 15.88 23.91±5.25 19.96 低盖度车桑子
Low coverage Dodonaea viscosa雨季Rainy season 21.68±5.09 23.48 21.21±4.26 20.08 17.51±3.63 20.73 旱季Dry season 18.77±5.38 28.66 18.10±5.32 27.85 14.58±3.74 25.65 高盖度车桑子
High coverage Dodonaea viscosa雨季Rainy season 29.72±7.02 23.62 23.57±3.60 15.27 13.12±1.69 12.88 旱季Dry season 20.63±5.72 27.73 18.93±3.23 17.06 10.97±1.82 16.59 Table 3. Soil moisture content and variation coefficient of three shrubs
一般采用变异系数(Cv)和标准差(S)将剖面土壤水分垂直变化划分为4层:速变层(Cv≥30%且S≥4%)、活跃层(Cv:20%~30%且S:3%~4%)、次活跃层(Cv:10%~20%且S:2%~3%)和相对稳定层(Cv≤10%且S≤2%)[10]。但在实际应用过程中,标准差和变异系数很难同时满足以上要求,此时可依据变异系数来划分土壤水分的垂直变化层次[11]。综合分析表 2、表 3得出:雨季桉树林0~30 cm土层为活跃层,30~70 cm土层为次活跃层,冲天柏林0~10 cm、50~70 cm土层为活跃层, 10~50 cm土层为次活跃层,马尾松林0~50 cm土层为次活跃层,50~70 cm土层为活跃层。旱季表现为:桉树林和冲天柏林0~70 cm土层均为活跃层,马尾松林0~10 cm土层为活跃层, 10~70 cm土层为次活跃层。在雨季和旱季时段内,天然次生灌丛和高盖度车桑子灌丛0~10 cm土层均为活跃层,10~30 cm土层均为次活跃层。低盖度车桑子灌丛0~30 cm土层均为活跃层。
-
不同植被恢复模式各土层月均土壤体积含水量如图 3所示。结果表明:由于降雨的季节性分布,土壤水分波动范围较大,但这种波动随土层深度的增加而减小,高盖度车桑子灌丛20~30 cm土层砾石含量较多,故其全年土壤体积含水量波动范围较小。林地各土层间土壤体积含水量差异显著(P < 0.05)。以2016年10月为分界线,6种植被类型土壤水分含量在5—10月呈上升趋势,10月至次年4月呈下降趋势,各土层土壤体积含水量均在8—9月达到最大值,其变化范围为26.25%~44.08%,3—4月达到最小值,土壤体积含水量变化范围在9.48%~17.47%之间。3—4月各植被类型表层土壤体积含水量较高,但3种乔木林地30~70 cm土层和3种灌丛20~30 cm土层深度土壤水分含量接近于植物萎蔫系数(13%~16.6%)[3, 12],当雨季刚来临时,少量的降雨主要补充表层土壤,这对于浅根性草本植物种子萌发尤为重要。另外,不同植被类型间月均土壤体积含水量差异显著(P < 0.05),表现为桉树林(27.30%)>天然次生灌丛(26%)>冲天柏林(24.23%)>马尾松林(20.47%)>高盖度车桑子灌丛(19.52%)>低盖度车桑子灌丛(18.84%)。
-
喀斯特断陷盆地具有土层浅薄、石灰岩的高渗透性以及较高的蒸散速率等特点,研究充足降雨后持续干旱条件下土壤水分的衰减规律,对于该地区植被恢复工作意义重大。此次降雨时间为2016年11月8日—9日,总降雨量83.2 mm,此次降雨前9天(10月30日—11月7日)和降雨后19天(11月10日—11月28日)均无有效降雨。
分析图 4可知,3种乔木林0~30 cm土层、3种灌丛0~10 cm土层降雨后约10天土壤体积含水量均恢复为降雨前1天水平,主要原因可能是旱季来临,部分草本植物叶片萎蔫,地表植被覆盖度有所降低,土壤表层水分蒸发强烈。3种乔木林30~70 cm土层土壤水分波动较小,土壤水分降雨后19天均高于植物萎蔫系数,主要原因可能是此次降雨充足,另外植物在旱季耗水量较低,且土层较深,受蒸发影响较小。天然次生灌丛10~30 cm土层土壤水分表现为降雨后约15天土壤体积含水量恢复为降雨前1天水平,降雨后19天土壤体积含水量仍高于植物萎蔫系数,主要原因可能是该林地高达90%的植被覆盖度,生物多样性丰富,样地植株冠幅较大(表 1),起到遮阴、保水的作用,从而使10~30 cm土层的内部环境较稳定,土壤保水性较强。低盖度和高盖度车桑子灌丛10~30 cm土层土壤体积含水量均表现为降雨后约15天恢复为降雨前1天水平,并且接近植物萎蔫系数,主要原因可能是样地植被少,覆盖度低,土层薄,受植被盖度与土壤蒸发等因素共同影响。