王麗艷，楊 樺，唐星林，周 晨，劉光正

(江西省林業(yè)科學院，江西 南昌330013)

贛江源是贛江的最大源頭，是贛鄱大地的生態(tài)屏障，對保障贛江流域、鄱陽湖乃至長江中下游水資源生態(tài)安全具有重要作用。針對源頭區(qū)大面積人工針葉純林及殘次林林相單一、水源涵養(yǎng)能力差等關(guān)鍵生態(tài)環(huán)境問題，開展了人工針葉林結(jié)構(gòu)優(yōu)化、低效殘次林水土流失控制綜合治理關(guān)鍵技術(shù)研究。人工林生態(tài)系統(tǒng)在人為干擾下能否保持穩(wěn)定、變化程度有多大、恢復時間為多長等一系列問題是人們關(guān)注的焦點，也是穩(wěn)定性研究的熱點[1]。這關(guān)系到人們能否合理利用森林資源、保護生態(tài)環(huán)境、促進經(jīng)濟增長和社會穩(wěn)定[2-3]。但大多數(shù)的研究都是針對林地枯落物本身或從不同森林類型枯落物持水特性的比較進行[4-8]，而針對低效林改造后的林分涵養(yǎng)水源性能變化，以及通過空間代時間的方法對低效林改造后期的效果研究較少。

贛江源區(qū)水源涵養(yǎng)功能弱的防護林主要是針葉純林，包括杉木（Cunninghamia lanceolata）林、馬尾松(Pinus massoniana)林、濕地松（P. elliottii）林等，本研究結(jié)合南方丘陵地區(qū)常見闊葉樹，通過林下補植鄉(xiāng)土闊葉樹木荷 (Schima superba)、楓香（Liquidambar formosana）等調(diào)整林分結(jié)構(gòu)。改造后的林分通過群落演替，將向以闊葉樹為優(yōu)勢種群的針闊混交林或闊葉林演替[9-12]。因此，本研究以空間代時間的方法，在贛江源頭區(qū)選擇馬尾松與木荷(Schima superba)的復合林分、杉木與闊葉樹復合林，以及闊葉林為參照，作為改造林分在未來較短時期內(nèi)演替進程中發(fā)育形成的目標林分，分析現(xiàn)改造低效及殘次針葉林和參照林分生態(tài)功能的差異，來判斷針葉林改造模式達到未來10～20 a 后生態(tài)功能狀況。通過對不同優(yōu)化措施林分在各種降雨條件下產(chǎn)流產(chǎn)沙的連續(xù)4 a 動態(tài)監(jiān)測，開展不同人工針葉林結(jié)構(gòu)改造技術(shù)模式的生態(tài)水文效益評價，篩選出適合源頭區(qū)的、低擾動的人工針葉林結(jié)構(gòu)優(yōu)化模式，以期為亞熱帶低效針葉林改造模式的選擇提供可靠方案，進而為長江中下游水源涵養(yǎng)林提質(zhì)增效提供參考。

1 研究地點與方法

1.1 研究地林分概況

研究區(qū)位于江西贛江源國家級自然保護區(qū)（以下簡稱贛江源自然保護區(qū)），該保護區(qū)位于我國華東“屋脊”武夷山脈的中南段，與福建長汀交界，地理坐標為25°56′30″-26°07′42″ N，116°15′01″-116°29′06″ E，是中亞熱帶與南亞熱帶的過渡區(qū)。研究樣地位于石城縣贛江源自然保護區(qū)。低效林和殘次林通過間伐及林下補植鄉(xiāng)土闊葉樹的改造方法，起到林分結(jié)構(gòu)優(yōu)化及防護功能提升的目的，以研究區(qū)馬尾松和木荷的復合林分、杉木闊葉樹復合林，以及闊葉林為參照。不同類型林分如表1。

表1 研究區(qū)不同類型林分概況Tab. 1 General situation of different types of forest stand in study area

1.2 以空間代時間的方法

2012 年對低效和殘次針葉林分進行林下補植闊葉樹、灌木及草本的補植改造。為了進一步分析低效針葉林改造后林分隨時間維度的變化效果，因被改造林分長時間維度監(jiān)測數(shù)據(jù)在短短幾年內(nèi)無法獲取，因此采用空間代替時間法。群落的空間序列研究法，即空間代替時間法，是植被演替研究中常采用的方法，對無長期固定生態(tài)監(jiān)測站點的研究來說是尤為有效的一種研究手段[5]。

本研究以馬尾松和木荷的復合林分、杉木闊葉樹復合林，以及闊葉林為參照，作為低效針葉改造林分在未來較長時間維度演替進程中發(fā)育形成的目標林分（即被改造針葉林分不同演替階段的群落空間），以林分的演替空間代替恢復時間維度的研究方法。

1.3 低效和殘次林的判定方法

根據(jù)國家林業(yè)行業(yè)標準《低效林改造技術(shù)規(guī)程》（LY/T 1690—2007）對低效林分進行判別。

1.4 林分及物種調(diào)查方法

于2014、2016 和2017 年的7-9 月，選擇林分設(shè)置標準樣方。對樣地內(nèi)喬木開展每木檢尺調(diào)查，測定其胸徑、樹高、冠幅等，分別在每個喬木樣方內(nèi)設(shè)置3個5 m×5 m 的灌木樣方，調(diào)查灌木和DBH＜5.0 cm 的喬木樹種的種類、數(shù)量、高度、蓋度等。分別在每個灌木樣方內(nèi)的四角及中心位置設(shè)5 個1 m×1 m 的草本樣方，分別調(diào)查草本層的種類、數(shù)量、高度、蓋度等。同時記錄每塊標準樣地的地理坐標、郁閉度、坡向、坡度、枯枝落葉層厚度及土壤類型等生態(tài)環(huán)境因子。

1.5 徑流小區(qū)布設(shè)和指標觀測

在各類型樣地建有20 m×5 m 標準徑流場，徑流場下方建1 m×1 m×1 m 的集水池，安裝水表和水尺等測流設(shè)施，進行不同改造模式徑流觀測、泥沙觀測和水質(zhì)觀測。在試驗區(qū)內(nèi)建一個小型自動氣象站，自動觀測和記錄各種氣象因子。徑流場觀測內(nèi)容為徑流小區(qū)產(chǎn)生徑流大小和徑流所攜帶的泥沙量，并對水樣進行監(jiān)測。分析確定試驗區(qū)內(nèi)的水土保持、水土涵養(yǎng)和凈化水質(zhì)的作用。

1.6 數(shù)據(jù)計算

降雨量：取自研究區(qū)自動氣象觀測系統(tǒng)和人工觀測的觀測數(shù)據(jù)。徑流系數(shù)：年產(chǎn)地表徑流的量/年降雨量×100%徑流量（t·hm-2）=V×A×ΔT（式中，V為平均水流速；A為過水斷面積；ΔT為一定時段）

使用Excel 2010 和統(tǒng)計軟件DPS7.5 和SPSS 19.0 進行數(shù)據(jù)計算分析。

2 結(jié)果與分析

2.1 不同類型林分生物量和土壤質(zhì)量特征

研究區(qū)改造林分林下植被覆蓋度和生物量大幅提高，土壤理化性質(zhì)得以改善。低效濕地松改造林林地上生物量由57.80 t·hm-2提高到64.17 t·hm-2，林下灌草覆蓋度提高了15.67%，低效馬尾松改造林分林下灌草覆蓋度提高了33.32%、地上部分生物量提高了5.55 t·hm-2。殘次杉木改造林林下灌草覆蓋度提高了26.13%、地上部分生物量提高了8.21 t·hm-2。2017年與2014 年相比，改造低效和殘次林分林地土壤養(yǎng)分含量有所增加。有機質(zhì)含量的增幅為1.04～12.44 g·kg-1，全氮含量增幅為0.02～0.12 g·kg-1，堿解氮含量增幅5.07～7.75 mg·kg-1，速效磷含量增幅為0.04～0.11 mg·kg-1，速效鉀含量增幅為19.77～50.13 mg·kg-1，而土壤中全磷含量增幅僅為0.02～0.04 g·kg-1（表2）。

表2 不同類型林分生物量和土壤性質(zhì)Tab. 2 Stand biomass and soil properties of different types of forest stand

在森林的生長發(fā)育過程中，林分結(jié)構(gòu)是一個重要的驅(qū)動因子。促進森林發(fā)育的干擾，主要表現(xiàn)為森林結(jié)構(gòu)調(diào)整[6]，說明喬、灌、草3 層結(jié)構(gòu)在森林防護效益上有不可忽視的重要性，只有在喬木、灌木、草本和枯落物層具有較大覆蓋度并合理搭配時，林地固持土壤能力的防護功能才是高效的。

2.2 不同類型林分水源涵養(yǎng)能力比較

2.2.1 枯落物持水量

不同改造林分與參照林分枯落物蓄積量如圖1所示。杉木殘次林、濕地松林和馬尾松林改造林分中林下枯落物的蓄積量和厚度都遠小于參照林分。主要原因是喬木層中闊葉樹處于幼林，而杉闊復合林與馬闊復合林中，闊葉樹已成為優(yōu)勢群落，在林分中占有優(yōu)勢地位，隨著林分生長周期林下枯落物不斷增多，其蓄積量和厚度均遠高于人工針葉純林。2016 年枯落物的蓄積量依次為杉闊（8.48 t·hm-2）＞馬闊（8.30 t·hm-2）＞闊葉（6.12 t·hm-2）＞濕地松改造林（4.36 t·hm-2）＞馬尾松改造林（2.75 t·hm-2）＞殘次杉木改造林（2.56 t·hm-2）。枯落物的蓄積量與枯落物的厚度呈正相關(guān)，馬尾松改造林與殘次杉木改造林的林下枯落物厚度分別是1.73 cm 和1.38 cm。

圖1 不同類型林分枯落物蓄積量Fig. 1 Litter accumulation of different types of forest stand

從2016 年枯落物的最大持水量來看（圖2），馬闊復合林最大持水量是23.30 t·hm-2＞闊葉林22.73 t·hm-2＞低效馬尾松改造林（2.53 t·hm-2）。杉闊復合林（26.90 t·hm-2）＞闊葉林（22.73 t·hm-2）＞低效濕地松改造林（5.79 t·hm-2）＞殘次杉木改造林（2.12 t·hm-2）?？萋湮锏挠行r蓄水量依次是杉闊復合林（19.79 t·hm-2）＞馬闊復合林（18.06 t·hm-2）＞闊葉林（17.03 t·hm-2）＞濕地松改造林（4.44 t·hm-2）＞馬尾松改造林（2.02 t·hm-2）＞殘次杉木林（1.61 t·hm-2）。

圖2 不同類型林分枯落物持水量Fig. 2 Interception amount of litter layers of different types of forest stand

2.2.2 土壤容重和土壤持水量

不同類型林分土壤持水特征如表3 所示。2016年測定的闊葉林土壤容重0.98 g·cm-3，杉闊復合林土壤容重1.05 g·cm-3，馬闊復合林土壤容重1.05 g·cm-3，濕地松改造林土壤容重為1.18 g·cm-3，馬尾松改造林土壤容重為1.23 g·cm-3，杉木殘次林的土壤容重是1.27 g·cm-3。土壤容重越小越疏松，土壤蓄水能力越強。在試驗區(qū)的馬尾松改造林分和殘次杉木林地土壤多為砂石土，土壤容重較大，蓄水能力較差。

表3 不同類型林分土壤容重和持水量Tab. 3 Stand soil bulk density and moisture capacity of different types of forest stand

不同林分類型的土壤持水量，與土壤容重的表現(xiàn)狀況相近，土壤儲蓄的水分總量取決于土壤質(zhì)地、土壤層厚度、容重、毛管孔隙度和非毛管孔隙等物理因素。土壤毛管孔隙蓄存的水分提供樹木根系吸收或土壤蒸散，而非毛管孔隙為水分蓄存提供了空間。土壤的最大持水量依次為杉闊復合林(4 721.02 t·hm-2)＞馬闊復合林(4 674.78 t·hm-2)＞闊葉林(4 617.22 t·hm-2)＞濕地松改造林（3 562.27 t·hm-2）＞馬尾松改造林（2 635.67 t·hm-2）＞杉木殘次林改造林（2 482.71 t·hm-2）。不同林分的土壤毛管持水量變化為：杉闊復合林（4 520.44 t·hm-2）＞馬闊復合林（4 430.70 t·hm-2）＞闊葉林（4 349.74 t·hm-2）＞濕地松改造林（3 367.84 t·hm-2）＞馬尾松改造林（2 635.67 t·hm-2）＞杉木殘次林（2 309.50 t·hm-2）。

由此可見，各改造林分的水源涵養(yǎng)能力幾乎都低于參照林分的1/2 或2/3，這表明，改造后林分的水源涵養(yǎng)功能提升空間很大，隨著林分的更新演替，現(xiàn)有林分的水源涵養(yǎng)能力將不斷提高。

2.2.3 地表徑流水質(zhì)

圖3 反映了不同林分類型的地表徑流的pH 值。不同植被狀況對區(qū)域內(nèi)水源的pH 值影響較大。低效濕地松和馬尾松改造林、馬闊復合林的地表徑流水為弱堿性，杉闊復合林、闊葉林、殘次杉木改造林地表徑流水為微酸性。2014—2017 年不同類型改造林分徑流水pH 值差異不顯著（P＜0.05）。

圖3 不同類型林分地表徑流pH 值Fig. 3 pH of surface runoff water of different types ofmodified forest stands

水中溶解氧的含量與空氣中氧的分壓、水的溫度都有密切關(guān)系。水中溶解氧的多少是衡量水體自凈能力的一個指標。水里的溶解氧被消耗，要恢復到初始狀態(tài)所需時間短，說明該水體的自凈能力強，或者水體污染不嚴重。否則說明水體污染嚴重，自凈能力弱，甚至失去自凈能力。因此，溶解氧飽和度越高，說明其自凈能力越強。各林分地表徑流中溶解氧飽和度2017 年變化趨勢為：不同發(fā)育階段林分年際間差異不明顯(圖4)，闊葉林（78.02%）＞杉闊復合林（76.20%）＞馬闊復合林（75.36%）＞濕地松改造林（69.72%）＞馬尾松改造林（64.72%）＞杉木殘次林（63.12%）。根據(jù)國家地表水環(huán)境質(zhì)量標準（GB/T 3838—2002），馬闊復合林地表徑流水質(zhì)達到Ⅱ類水質(zhì)標準，濕地松和馬尾松改造林為Ⅲ級數(shù)值標準。這說明，對照的杉闊和馬闊林地的水自凈能力要高于改造低效和殘次林分，隨著林分的更新演替，各改造林分對水的自凈能力將不斷增強，林分中涵養(yǎng)水的水質(zhì)將不斷改善。

圖4 不同類型林分地表徑流水溶解氧Fig. 4 Dissolved oxygen in surface runoff water of different types of modified forest stands

2.3 不同類型林分水土保持能力比較

不同森林由于結(jié)構(gòu)不同，其林冠層和下木層等地上植被對降水的攔截和阻擋作用差異較大。通過對不同林分類型的地表徑流的長期觀測可以發(fā)現(xiàn)，2014年地表徑流最大的是殘次杉木林，徑流量達到1 490.50 t·hm-2，最小為闊葉林726.50 t·hm-2。通過對不同類型林分不同年份徑流系數(shù)和徑流量方差分析表明，隨著林分生長年限的增長改造林分徑流量減少，并且年際間差異顯著（P＜0.05），2017 年與2014年相比，低效濕地松改造林、馬尾松改造林和杉木殘次改造林徑流量分別減少了6.56%、4.17%、4.76%。

通過觀測和計算得到2014 年闊葉林的土壤侵蝕量為116.90 t·km-2·a-1、杉闊復合林為235.10 t·km-2·a-1、杉木殘次林在585.75～612.15 t·km-2·a-1之間、馬尾松低效改造林為571.00 t·km-2·a-1（表4）。由此可見，隨著林分的發(fā)育更新，改造殘次林和低效林的土壤侵蝕量將有一定幅度的下降，改造林分的土壤侵蝕量減少量在20.95～26.40 t·km-2·a-1之間，減少了3.7%～4.3%，再次證明了對低效針葉林進行闊葉化改造的重要性和必要性，也表明了改造后的低效林和殘次林在未來林內(nèi)水土流失狀況將得到明顯改善。如果采取合理的方案和技術(shù)，可以避免水土流失加劇的情況。

表4 不同類型林分地表徑流與土壤侵蝕量Tab. 4 Stand surface runoff and amount of soil erosion in different types

對不同演替階段的幾種林分類型的地表徑流、土壤侵蝕量等進行了方差分析比較。結(jié)果發(fā)現(xiàn)，馬尾松改造林與殘次杉木改造林和其他更高演替階段的林分在徑流系數(shù)和徑流量上表現(xiàn)出差異顯著性，徑流量更高；更高演替階段的林分與改造的杉木殘次林在土壤侵蝕量與徑流中泥沙含量上表現(xiàn)出差異顯著性，隨著演替階段的上升，林分的土壤侵蝕量與地表徑流中泥沙含量減小。

3 結(jié)論與討論

恢復森林水源涵養(yǎng)功能的核心是調(diào)整林冠組成和結(jié)構(gòu)[8]。植物群落由共存的物種構(gòu)成，組成群落的各物種間的關(guān)系，決定著群落的結(jié)構(gòu)特征和群落的動態(tài)[13-14]。我國亞熱帶廣大地區(qū)的森林次生演替規(guī)律為：馬尾松為優(yōu)勢的針葉林→馬尾松為優(yōu)勢的常綠針闊混交林→闊葉林為優(yōu)勢的常綠針闊混交林→常綠闊葉林[15-16]。王方超[11]等對江西中亞熱帶針闊混交林的物種組成和空間分布特征研究表明，針闊混交林是中亞熱帶先鋒群落馬尾松林相頂級群落常綠闊葉林演替的中間階段。本研究采用了空間代時間的方法，以針闊混交林和闊葉林作為參照林分，研究低效針葉林闊葉化改造林分通過群落演替，將向以闊葉樹為優(yōu)勢種群的針闊混交林或闊葉林演替，從而推斷林分改造后的成效。低效林的改造就是通過調(diào)控林分結(jié)構(gòu)，進而改善土壤生態(tài)功能[12]。補植是低效林改造的有效措施之一，隨著郁閉度增加，植物多樣性增大[16]。許多學者的研究結(jié)果顯示通過補植闊葉樹，低效針葉林林下生物量和物種豐富度增加[17-19]。本研究改造林分林下闊葉樹以鄉(xiāng)土樹種為主，主要是木荷和楓香為主，體現(xiàn)了地帶性針闊混交林的群落類型，以期達到改造目標為喬、灌、草3 層復層結(jié)構(gòu)。改造后的低效和殘次針葉林將加快林上植被更新發(fā)育的進程，未來不斷演替進化為具有喬、灌、草復層結(jié)構(gòu)的林分，其生物多樣性將大大提升。低效和殘次針葉林改造4 a 后，林分地上部分生物量和林下植被的覆蓋度得到了較大幅度的提高，分別提高15.67%和33.32%。研究區(qū)改造林分林下植被覆蓋度和生物量大幅提高，土壤理化性質(zhì)得以改善。

森林的水源涵養(yǎng)功能是森林生態(tài)服務(wù)功能的重要功能之一，不同森林類型由于發(fā)育階段不同、結(jié)構(gòu)存在差異，林分整體水源涵養(yǎng)功能也有一定的差異[20-21]。枯落物的持水能力主要集中在未分解層和半分解層，已分解層補充了土壤有機質(zhì)，增加了土壤通透性和持水能力，從而提高土壤的持水能力，在林分改造和經(jīng)營過程中應(yīng)盡可能地加以保護[22-24]。林分改造由于改變了林地光照、水分等條件，在促進更新樹種及灌草生長的同時，也加劇了地表水土流失和有機物分解，一定程度上導致了土壤抗蝕性降低[25]。林分改造降低了喬木層的林冠截留，有利于林分涵養(yǎng)水源[26]。美國和日本采用林草復合改造增加林地肥力[27]。許多研究發(fā)現(xiàn)改造后的林分凋落物、土壤的持水能力及肥力都有一定程度的提升[28-30]。本研究表明馬尾松改造模式與杉木殘次林改造模式在徑流系數(shù)和徑流量上表現(xiàn)出顯著性差異。改造后的林分可增強水自凈能力，提高了地表徑流水溶解氧飽和度，增幅為4%～12%。本研究再次證明了對低效針葉林進行闊葉化改造的重要性和必要性，也表明了殘次林改造后在未來林內(nèi)水土流失狀況將得到明顯改善。如果采取合理的方案和技術(shù)，可以避免水土流失加劇的情況。林分結(jié)構(gòu)的復雜化與穩(wěn)定需要長時間的發(fā)育，雖然本研究中經(jīng)密度調(diào)控與結(jié)構(gòu)優(yōu)化的林分在物種的豐富性、喬灌草層的生物量、林分立體結(jié)構(gòu)上有明顯的改善，但其水源涵養(yǎng)功能在短短的4 a 項目期中無法得到切實體現(xiàn)。今后，應(yīng)對改造后的林分進行水土流失的長期觀測，進一步深入開展水土流失動態(tài)變化規(guī)律的研究，并且建立針葉林改造不同類型生物量及演替預測模型，為森林生態(tài)系統(tǒng)服務(wù)功能評價提供理論依據(jù)。