王夢(mèng)蕾,袁紅娟,白 宇,孫照斌*
(1.河北農(nóng)業(yè)大學(xué) 林學(xué)院,河北 保定 071000;2.河北省文安縣自然資源和規(guī)劃局,河北 文安 065800)
樹(shù)干液流是植物進(jìn)行蒸騰作用時(shí),因蒸騰拉力所產(chǎn)生的樹(shù)木水分流動(dòng)過(guò)程[1-3]。樹(shù)干液流的空間差異主要表現(xiàn)在徑向深度、方位和垂直高度3個(gè)方面,目前研究主要集中在不同徑向深度和不同方位,并發(fā)現(xiàn)大部分樹(shù)種邊材徑向不同深度處和不同方位處的樹(shù)干液流存在變化[4-5]。
由于林木間水分傳輸結(jié)構(gòu)差異,距形成層不同深度的樹(shù)干液流速率存在較大差異,如二白楊[6](Populusgansuensis)邊材液流速率最大值在距形成層30 mm深處;尾巨桉(Eucalyptusurophylla×E.grandis)[7]、青海云杉(Piceacrassifolia)[8]樹(shù)干液流速率在距形成層0~20 mm深處液流最大;四川山礬[9](Symplocossumuntia)在形成層10~30 mm處樹(shù)干液流速率最大,杉木(Cunninghamialanceolata)[9]在距形成層30~50 mm處樹(shù)干液流速率最大;黃土區(qū)刺槐(Robiniapseudoacacia)和遼東櫟[10](Quercuswutaishanica)樹(shù)干液流速率隨徑向深度增加而遞減。
研究發(fā)現(xiàn),不同方位間的樹(shù)干液流速率也存在很大差異,如西伯利亞落葉松(Larixgmelinii)[11]和樟子松(Pinussylvestrisvar.mongolica)[12]液流速率在南側(cè)較高,但棗樹(shù)(Ziziphusjujuba)液流速率在西側(cè)較高[13],胡楊(P.euphratica)的液流速率在東側(cè)最高[14]。旱柳(Salixmatsudana)與小葉楊(P.simonii)的方位間差異隨季節(jié)而變化等[15]。黨宏忠等[16]研究發(fā)現(xiàn),蘋(píng)果(Maluspumila)不同方位液流速率大小排序?yàn)楸眰?cè)>南側(cè)>東側(cè)>西側(cè);韓輝等[17]研究發(fā)現(xiàn),沙地赤松(P.densiflora)的不同方位液流速率大小一般為南側(cè)>西側(cè)>北側(cè)>東側(cè)??梢?jiàn),在樹(shù)干液流速率測(cè)定過(guò)程中必須考慮到徑向深度和方位差異,以免造成較大測(cè)算誤差。
目前部分研究認(rèn)為,樹(shù)干液流的空間差異可能與樹(shù)木邊材水力結(jié)構(gòu)的非均質(zhì)性、樹(shù)冠幾何形狀或根系分布的空間不對(duì)稱(chēng)性等有關(guān)[18]。對(duì)于影響樹(shù)干液流的研究主要集中在外界因素(環(huán)境因子),其中溫度、風(fēng)速、濕度、光合有效輻射、水汽壓虧缺等因素均會(huì)造成樹(shù)干液流的變化[19-21],而樹(shù)木本身的一些生物特征(木材解剖)對(duì)樹(shù)干液流的影響研究并不多見(jiàn)[22]。
油松(P.tabuliformis)是中國(guó)北方地區(qū)的典型樹(shù)種之一,廣泛分布于東北、華北等地區(qū),在中國(guó)北方森林群落中具有一定代表性。在外界環(huán)境大致相同的條件下,本研究選取河北省冀北山區(qū)大力推廣的油松,通過(guò)測(cè)定不同徑向深度和不同方位樹(shù)干液流與木材解剖特征,分析樹(shù)木構(gòu)造對(duì)水分輸送功能的影響,以期為不同樹(shù)種樹(shù)干液流的比較、耗水評(píng)價(jià)、林木水分傳輸機(jī)理提供理論依據(jù)。
試驗(yàn)地位于河北省承德市隆化縣徐八屋林場(chǎng)(117°62′56″E,41°27′70″N),土壤類(lèi)型為棕壤土,坡向?yàn)闁|坡,南坡度約13°。地勢(shì)北高南低,海拔650~1 000 m。屬暖溫帶半濕潤(rùn)大陸性季風(fēng)型山地氣候。全年日照時(shí)數(shù)為2 777 h,年平均溫度5.1 ℃。日照時(shí)間隨季節(jié)的變化與太陽(yáng)輻射的變化規(guī)律一致,日照分布的季節(jié)以春季為最高,日照時(shí)間在1月和12月最少。全年降水量為480 mm左右,其中3-5月占全年降水量的14%,6-8月占70%左右,9-12月占12.5%,冬季占3.5%。年均有風(fēng)日數(shù)110 d,最多130 d,風(fēng)沙日數(shù)62 d,大風(fēng)日數(shù)38~40 d,風(fēng)向以西北風(fēng)為主,春、秋兩季多風(fēng),空氣干燥,蒸發(fā)量比較大,多年相對(duì)濕度年平均為55%,年蒸發(fā)量1 400 mm左右。
該林場(chǎng)主要喬木有油松、落葉松、樺樹(shù)(Betula)、山楊(P.davidiana)、椴樹(shù)(Tiliatuan)、蒙古櫟(Q.mongolica)、五角楓(Acermono)等。按照樹(shù)冠圓滿(mǎn)、樹(shù)干通直、長(zhǎng)勢(shì)良好的原則,選擇3株近熟油松個(gè)體作為研究對(duì)象,樣樹(shù)基本情況見(jiàn)表1。
表1 樣樹(shù)基本情況Table 1 Basic information of sample trees
測(cè)定方法按照馬長(zhǎng)明等[23]所采用的熱擴(kuò)散樹(shù)干液流測(cè)定試驗(yàn)方法進(jìn)行。2021年8-9月使用Granier熱擴(kuò)散式莖流計(jì)連續(xù)監(jiān)測(cè)。樹(shù)干液流速率按不同徑向深度與不同方位分別測(cè)定。
測(cè)定不同徑向深度液流速率時(shí),選取探針直徑2 mm,長(zhǎng)度規(guī)格為10、20、30、50 mm,在樣樹(shù)胸高位置的北側(cè)1/3圓周內(nèi),平行安裝4組探針,探針長(zhǎng)度從左至右分別為10、20、30 mm和50 mm,4組探針之間水平相隔約5 cm,保證探針之間互不干擾。
測(cè)定不同方位液流速度時(shí),使用手持羅盤(pán)儀,確定東西南北4個(gè)方位,在4個(gè)方位上均插上直徑2 mm、長(zhǎng)度為30 mm的探針,保證探針之間互不干擾。
待樹(shù)干液流測(cè)定完成后,將3株油松樣樹(shù)伐倒,在樹(shù)干液流測(cè)定的1.3 m左右,截取3個(gè)厚度為50 mm圓盤(pán)(其中2個(gè)分別用于不同徑向深度和不同方位解剖特性測(cè)定,另1個(gè)備用)。
測(cè)定不同徑向深度解剖特性時(shí),將1個(gè)圓盤(pán)通過(guò)髓心沿南北方向鋸成與圓盤(pán)厚度相同的2 cm寬長(zhǎng)條,在該長(zhǎng)條上從形成層到髓心每隔10 mm取樣,即將長(zhǎng)條分為0~10、10~20、20~30、30~40、40~50 mm 5個(gè)部分,在各部分各取火柴桿粗細(xì)試樣4個(gè)用于測(cè)定,代表各個(gè)徑向深度的解剖特性。
測(cè)定不同方位解剖特性時(shí),另取1個(gè)圓盤(pán),在東、南、西、北4個(gè)方位,距形成層0~10、10~20、20~30 mm處各取火柴桿粗細(xì)試樣2個(gè),共6個(gè)試樣,混合后用于測(cè)定,代表各個(gè)方位的解剖特性。
將上述各木材試樣采用富蘭克林離析法解離、制片[24],用3%番紅水溶液染色1 min,將制好的片放置在數(shù)碼顯微圖像電腦分析系統(tǒng)上對(duì)木材管胞形態(tài)測(cè)定。在10×物鏡下觀(guān)察木質(zhì)部解剖特征,選取30個(gè)視野拍照。使用Image-J軟件測(cè)量管胞寬度、管胞長(zhǎng)度和管胞長(zhǎng)寬比等數(shù)據(jù)后取平均值。
運(yùn)用Execl2013進(jìn)行數(shù)據(jù)整理、計(jì)算,運(yùn)用SPSS20.0進(jìn)行曲線(xiàn)擬合及相關(guān)性分析,運(yùn)用Origin2021做圖。
由圖1可知,隨著木質(zhì)部徑向深度的增加,油松的液流速率徑向變化規(guī)律呈先降低后升高的趨勢(shì),其中距形成層30~50 mm處液流速率最大,液流速率峰值為0.009 7 g·cm-2·s-1;其次為在20~30 mm和0~10 mm處,液流峰值分別為0.008 4、0.005 0 g·cm-2·s-1;在距形成層10~20 mm處最小,峰值僅為0.002 5 g·cm-2·s-1。其中30~50 mm處液流速率分別為20~30、10~20、0~10 mm處液流速率的1.15、3.88、1.94倍。油松樹(shù)干液流徑向分布格局屬于不規(guī)則型[25],液流密度由外向內(nèi)的分布存在2個(gè)或者多個(gè)峰值,變化規(guī)律不明顯。
圖1 不同徑向深度樹(shù)干液流速率Fig.1 Sap flow rate of trunkat different radial depths
根據(jù)油松不同徑向深度液流變化規(guī)律,建立不同深度處液流速率的數(shù)量關(guān)系,將單一深度的樹(shù)干液流速率轉(zhuǎn)化為多深度處樹(shù)干液流速率。由表2可知,油松V0~10 mm與V10~20 mm液流速率趨勢(shì)擬合程度較高,線(xiàn)性擬合的決定系數(shù)(R2)為0.945;油松V0~10 mm、V10~20 mm與V30~50 mm的線(xiàn)性擬合直線(xiàn)斜率較大,分別為4.999和4.265,說(shuō)明在0~10、10~20 mm與30~50 mm處液流速率差異較大,V0~10 mm與V10~20 mm的線(xiàn)性擬合直線(xiàn)斜率最小,V0~10 mm是V10~20 mm液流速率的1.153倍,說(shuō)明在0~10 mm與10~20 mm處液流速率差異較小。
表2 不同徑向深度樹(shù)干液流速率間線(xiàn)性回歸方程Table 2 Linear regression equation of sap flowat different radial depths
由圖2可知,樹(shù)干液流速率值最高的來(lái)自于西側(cè)(0.011 1 g·cm-2·s-1),其次為東側(cè)(0.008 6 g·cm-2·s-1)和北側(cè)(0.004 6 g·cm-2·s-1),南側(cè)最低(0.003 0 g·cm-2·s-1)。東、南、西、北各個(gè)方位液流速率的峰值時(shí)間出現(xiàn)在12:00左右,12:00-16:00開(kāi)始樹(shù)干液流速率逐漸降低,16:00后樹(shù)干液流速率急速下降。
圖2 不同方位樹(shù)干液流速率Fig.2 Sap flow rates in different directions
為了明確不同方位間液流速率間的定量關(guān)系,以北側(cè)(通常選擇在樹(shù)干北側(cè)測(cè)量,減少陽(yáng)光直射帶來(lái)的測(cè)量誤差)為因變量,建立其他方位與北側(cè)液流速率間的回歸關(guān)系(圖3)。由圖3可知,各方位液流速率與北側(cè)液流速率均存在極顯著的線(xiàn)性相關(guān)(P<0.01)。西側(cè)與北側(cè)間的線(xiàn)性擬合直線(xiàn)斜率為1.50,線(xiàn)性擬合R2為0.983;南側(cè)與北側(cè)間的線(xiàn)性擬合直線(xiàn)斜率為1.25,線(xiàn)性擬合R2為0.956;東側(cè)與北側(cè)間的線(xiàn)性擬合直線(xiàn)斜率為0.81、線(xiàn)性擬合R2為0.973。
圖3 北側(cè)液流速率與其他3個(gè)方位液流速率間的相關(guān)性Fig.3 Correlationof the flow rate among the north side and the other three sides
樹(shù)干液流特征反映了木質(zhì)部不同深度和不同方位導(dǎo)水性的差異,而組成木質(zhì)部管胞的結(jié)構(gòu)特征(如管胞直徑、長(zhǎng)度和長(zhǎng)寬比)是導(dǎo)致其水分運(yùn)輸能力差異的直接原因,進(jìn)而反映了不同空間樹(shù)干液流對(duì)木材構(gòu)造的響應(yīng)。
由圖4可知,管胞長(zhǎng)度和長(zhǎng)寬比在距形成層40 mm處最大(3 196.256 μm),在20 mm處最小(2 644.93 μm)。管胞寬度在距形成層40 mm處最小(35.56 μm),在30 mm處最大(42.36 μm)。不同徑向深度樹(shù)干液流速率與管胞長(zhǎng)度、管胞長(zhǎng)寬比變化趨勢(shì)一致,液流速率與管胞長(zhǎng)度徑向深度在10~20 mm處先降低,在20~50 mm處逐漸升高;管胞寬度與不同徑向深度樹(shù)干液流速率變化規(guī)律相反,在徑向深度50 mm處樹(shù)干液流速率最高,但管胞寬度最小。
a、b、ab、bc表示不同方位管胞指標(biāo)之間差異性。下同。圖4 油松不同徑向深度管胞解剖特征Fig.4 Tracheid length and tracheid width in different radial depths
由圖5可知,管胞長(zhǎng)度數(shù)值大小順序?yàn)?西側(cè)>南側(cè)>北側(cè)>東側(cè),西側(cè)平均管胞長(zhǎng)度(3 190.53 μm)分別為南側(cè)(2 958.27 μm)、北側(cè)(2 913.10 μm)、東側(cè)(2 875.79 μm)的1.09、1.10、1.11倍。管胞寬度數(shù)值大小順序?yàn)?東側(cè)>北側(cè)>南側(cè)>西側(cè),東側(cè)管胞寬度(42.00 μm)分別為北側(cè)(40.56 μm)、南側(cè)(37.09 μm)、西側(cè)(32.09 μm)的1.04、1.13、1.31倍。管胞長(zhǎng)寬比大小順序?yàn)?西側(cè)>南側(cè)>北側(cè)>東側(cè),不同方位管胞長(zhǎng)寬比與管胞長(zhǎng)度規(guī)律一致,西側(cè)管胞長(zhǎng)寬比(99.44)分別為南側(cè)(79.78)、北側(cè)(71.82)、東側(cè)(68.47)的1.24、1.39、1.45倍。結(jié)合圖1和圖2的樹(shù)干液流速率可知,樹(shù)干液流變化特征與解剖構(gòu)造之間響應(yīng)程度較高。不同方位樹(shù)干液流速率與管胞解剖構(gòu)造變化趨勢(shì)一致,西側(cè)液流速率與管胞長(zhǎng)度、管胞長(zhǎng)寬比均高于其他3個(gè)方位,東側(cè)液流速率與管胞長(zhǎng)度、管胞長(zhǎng)寬比均為最低值;管胞寬度與不同方位深度樹(shù)干液流速率變化規(guī)律相反,管胞寬度越小,液流速率越大。
圖5 油松不同方位管胞解剖特征Fig.5 Tracheid length and tracheid width in different orientations
隨著木質(zhì)部徑向深度增加,徑向液流速率變化規(guī)律為V30~50 mm>V20~30 mm>V0~10 mm>V10~20 mm,最大液流速率為0.009 7 g·cm-2·s-1,各深度液流速率之間呈顯著線(xiàn)性相關(guān)關(guān)系;不同方位液流變化規(guī)律為西側(cè)>東側(cè)>北側(cè)>南側(cè)。樹(shù)干北側(cè)方位的液流速率與其他方位間均具有緊密的線(xiàn)性關(guān)系。
油松的管胞長(zhǎng)度、寬度和管胞長(zhǎng)寬比與樹(shù)干液流的徑向和方位變化規(guī)律趨勢(shì)大致相同,即管胞長(zhǎng)度越大,管胞寬度越小;管胞長(zhǎng)寬比越大,樹(shù)干液流速率越高;樹(shù)干液流變化特征與解剖構(gòu)造之間響應(yīng)程度較高。
本研究中油松樹(shù)干液流徑向變化規(guī)律與張璇等[8]研究杉木在形成層30~50 mm處樹(shù)干液流最大的結(jié)果一致。此外,孫守家等[26]、劉超等[5]、孫慧珍等[27]、王華等[28]、張小由等[29]大量學(xué)者都探索了樹(shù)干液流密度沿徑向分布的不均勻性,由于不同學(xué)者所研究的樹(shù)木大小和樹(shù)木種類(lèi)存在很大差異,其心、邊材的劃分也不同,進(jìn)而導(dǎo)致不同樹(shù)種木質(zhì)部徑向分布規(guī)律不一致,所以液流速率最大值出現(xiàn)的位置也存在很大差異。由于本試驗(yàn)條件限制,最大徑向深度僅做到50 mm處,并未進(jìn)行更深層次樹(shù)干液流的測(cè)定,有可能在距離形成層50 mm至心材部分液流速率還有減小,因此并未得出油松樹(shù)干液流徑向變化的完整變化規(guī)律,還有待于深入研究。
本研究中樹(shù)干北側(cè)方位的液流速率與其他方位間均具有緊密的線(xiàn)性關(guān)系,油松西側(cè)液流量最大,與黃土區(qū)棗樹(shù)[13]研究結(jié)果類(lèi)似,但與其他學(xué)者對(duì)于西伯利亞松[11]、樟子松[12]等樹(shù)種研究結(jié)果不一致。造成不同方位樹(shù)干液流差異的原因可能與不同方位邊材寬度、相應(yīng)方位的冠幅大小、所處地理位置和太陽(yáng)輻射等因素有關(guān),本研究所測(cè)量樹(shù)干液流的油松西側(cè)冠幅較大,也可能是造成西側(cè)液流量最大的原因之一。
木材中的水分輸送通道是影響樹(shù)干液流的自身因素,油松的主要水分運(yùn)輸通道是管胞,管胞的數(shù)量和直徑都會(huì)對(duì)木質(zhì)部水分運(yùn)輸產(chǎn)生一定的影響,隨著管胞直徑的增大,最大導(dǎo)水率和運(yùn)輸水分的效率也隨之提高[29]。這是由于不同的木質(zhì)部直徑和長(zhǎng)度影響木質(zhì)部運(yùn)輸水分的阻力,木質(zhì)部管胞的直徑較小,長(zhǎng)度較短,則會(huì)增加水分輸送的阻力,使樹(shù)木水分運(yùn)輸效率降低。這與本研究中樹(shù)干液流速率與管胞寬度呈正相關(guān)關(guān)系結(jié)果一致。木質(zhì)部解剖結(jié)構(gòu)是栓塞抗性和水力導(dǎo)度的基礎(chǔ),與葉琳峰[30]研究的管胞與水力直徑對(duì)樹(shù)木輸水效率的相關(guān)性較高,且管胞直徑越大管胞密度越小的結(jié)果一致。
由于油松主要由管胞組成,結(jié)構(gòu)較為簡(jiǎn)單,管胞的主要功能為運(yùn)輸水分和支持作用[32-33]。有學(xué)者認(rèn)為對(duì)針葉樹(shù)而言,具有輸導(dǎo)作用的主要是軸向管胞,管胞直徑對(duì)導(dǎo)水性的影響小于管胞長(zhǎng)度對(duì)導(dǎo)水性的影響,因?yàn)楣馨麅?nèi)水分運(yùn)輸?shù)闹饕枇?lái)自于紋孔面上的紋孔,短的管胞比長(zhǎng)的管胞阻力大[34],這與本文中管胞長(zhǎng)度較管胞寬度與樹(shù)干液流響應(yīng)程度更高的結(jié)果一致。管胞長(zhǎng)度和寬度主要是受內(nèi)部因素(如植物激素等)和外部因素(如溫度和降雨量)共同作用的結(jié)果[35],單位面積內(nèi)管胞數(shù)量的多少也是影響樹(shù)木水分運(yùn)輸?shù)闹匾蛩?。油松為早晚材區(qū)別明顯的樹(shù)種,與早材相比,晚材管胞直徑小,細(xì)胞壁較厚,管胞腔較窄,晚材水分流速較早材大。
樹(shù)干液流的差異還可能受解剖構(gòu)造的其他方面影響(如紋孔)。水分在木質(zhì)部管孔中的傳輸阻力主要從管胞中產(chǎn)生的,其中紋孔口面積和紋孔膜面積越小,樹(shù)干液流速率越大,樹(shù)干液流率與紋孔結(jié)構(gòu)之間存在著比較復(fù)雜的聯(lián)系,紋孔面積假說(shuō)中[36],當(dāng)其上紋孔膜的總面積和數(shù)目增大時(shí),在管孔內(nèi)發(fā)生栓塞的可能性增大,本研究只對(duì)解剖構(gòu)造進(jìn)行了初步探索,還未涉及紋孔膜總面積、管胞上的紋孔密度、紋孔膜厚度等其余指標(biāo),尚不足以解釋紋孔結(jié)構(gòu)與樹(shù)干液流速率的關(guān)系,尚需深入探究。
樹(shù)種間樹(shù)干液流徑向分布格局的差異,不僅取決于木質(zhì)部結(jié)構(gòu)特征,還受氣象因素、立地條件等外界的影響,未來(lái)對(duì)樹(shù)干液流規(guī)律的研究需將內(nèi)部與外部影響因素相結(jié)合,多角度、全面、綜合地開(kāi)展。