劉 崴， 魏天興,， 朱清科

（1.北京林業(yè)大學(xué) 水土保持學(xué)院，北京100083；2.北京林業(yè)大學(xué) 水土保持國家林業(yè)局重點(diǎn)實(shí)驗(yàn)室，北京100083；3.北京林業(yè)大學(xué) 林業(yè)生態(tài)工程教育部工程研究中心，北京 100083）

土壤水分是黃土丘陵溝壑區(qū)植被恢復(fù)的主要限制因子，水分虧缺會(huì)導(dǎo)致植物生長衰退甚至死亡等變化［1－5］。盡管通常認(rèn)為該區(qū)植物是耐旱的，但是黃土丘陵溝壑區(qū)自然條件嚴(yán)酷，生態(tài)系統(tǒng)脆弱，加上極端干旱事件和全區(qū)域干旱天氣的增加，必然會(huì)引起一系列嚴(yán)重的生態(tài)環(huán)境問題［5－6］。因此，基于當(dāng)前水資源利用的矛盾［7］，需要全面了解主要造林樹種的耗水特性。河北楊Populus hopeiensis是陜北黃土丘陵區(qū)的鄉(xiāng)土樹種，是楊屬Populus中相對(duì)耐旱、耐瘠薄，能 “上山”的樹種之一，在當(dāng)?shù)鼐G化造林及維持生態(tài)系統(tǒng)的穩(wěn)定性中發(fā)揮著重要的作用［8－9］。油松Pinus tabulaeformis是黃土丘陵區(qū)的常見造林樹種，被廣泛應(yīng)用于 “三北”（東北、華北、西北）防護(hù)林、退耕還林工程等造林過程中。但陜北地區(qū)水資源供需不平衡，存在大量生長緩慢的油松［7,10］。所以，有必要對(duì)河北楊和油松在黃土丘陵區(qū)的環(huán)境適應(yīng)性及蒸騰耗水特性展開研究。國內(nèi)外相關(guān)研究表明，通過測(cè)定樹干液流可以較準(zhǔn)確地估算植物蒸騰耗水量，而熱技術(shù)法是測(cè)量喬木個(gè)體蒸騰耗水規(guī)律的主要方法，其中Granier熱擴(kuò)散探針法以其簡(jiǎn)單高效的特點(diǎn)得到廣泛的應(yīng)用［3,7,9,11］。周海光等［9］利用熱擴(kuò)散法研究了黃土丘陵溝壑區(qū)河北楊，小葉楊Populus simonii和山杏Armeniaca sibirica的樹干液流變化規(guī)律，萬艷芳等［12］利用熱擴(kuò)散法研究青海云杉Picea crassifolia樹干液流密度變化特征及其影響因素。上述研究發(fā)現(xiàn)：植物為滿足自身水分平衡，因根壓作用會(huì)在夜間產(chǎn)生液流活動(dòng)，張涵丹等［7］和溫杰等［13］的研究也表明：油松普遍存在夜間液流活動(dòng)。夜間液流廣泛存在植物生命活動(dòng)中，具有促進(jìn)早晨碳固定、營養(yǎng)供給、補(bǔ)充水分虧缺等重要的生理及生態(tài)意義，可能是揭示樹木抗旱機(jī)制的一種重要現(xiàn)象［14］。相關(guān)研究表明：樹干液流除了能直觀監(jiān)測(cè)植物的蒸騰耗水過程，還能揭示樹木生長變化［1,3］。觀測(cè)發(fā)現(xiàn)：樹干液流活動(dòng)和生長變化的耦合關(guān)系是植物在水分虧缺條件下的綜合反應(yīng)［3,15］。 如對(duì) 挪 威 云 杉Picea abies［16－17］，胡 楊Populus euphratica［18］， 華 北 落 葉 松Larix principis-rupprechtii［19］的研究表明：樹木蒸騰作用強(qiáng)度的日變化引起樹干組織水分的減少和增加，會(huì)導(dǎo)致樹干收縮和膨脹的變化。鑒于此，本研究選取陜北半干旱區(qū)典型樹種河北楊和油松為研究對(duì)象，采用熱擴(kuò)散式探針法，研究河北楊和油松樹干液流動(dòng)態(tài)變化特征及其對(duì)環(huán)境因子的響應(yīng)關(guān)系，并分析樹干液流密度與樹干莖周長的關(guān)系，旨在揭示液流活動(dòng)對(duì)外界環(huán)境因子的響應(yīng)規(guī)律，增加對(duì)黃土丘陵區(qū)人工林水分利用策略的理解，為該地區(qū)水資源承載力研究和樹種選擇提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

本研究地位于陜西省延安市吳起縣退耕還林森林公園大吉溝北京林業(yè)大學(xué)試驗(yàn)基地（36°54′20″N，108°10′28″E）。吳起縣屬于黃土高原丘陵溝壑區(qū)，海拔為1 233～1 809 m，氣候?yàn)榕瘻貛Т箨懶愿珊导撅L(fēng)氣候，日照充足，晝夜溫差大。年平均降水量和氣溫（1957－2013年）分別為464.6 mm和7.8℃，降水主要集中在7－9月，無霜期為120～155 d，年均蒸發(fā)量為891.23 mm。土壤類型主要為黃綿土。現(xiàn)有林分以人工灌喬木林為主，主要喬木樹種有河北楊，油松，刺槐Robinia pseudoacacia和小葉楊等，灌木主要有沙棘Hippophae rhamnoides和檸條Caragana korshinskii等。

試驗(yàn)地為2008年?duì)I造的人工油松、河北楊混交林，初植密度為1 400株·hm－2。河北楊是當(dāng)?shù)剜l(xiāng)土樹種，油松是主要造林樹種，林下植被主要有胡枝子Lespedeza davurica，冰草Agropyron cristatum，草木犀Melilotus officinalis等。試驗(yàn)在林內(nèi)1塊大小為20 m×20 m，坡向西偏南10°，坡度約為10°，樹種立地條件較一致的樣地內(nèi)進(jìn)行。

1.2 樹干液流及莖周長測(cè)定

利用熱擴(kuò)散探針（TDP-10，TDP-30，Dynamax，美國）于2016年6－10月對(duì)河北楊和油松人工林進(jìn)行樹干液流的觀測(cè)。該探針傳感器是基于Granier熱擴(kuò)散理論［20－21］。選擇生長良好的河北楊2株、油松3株作為觀測(cè)樣木（表1），將5個(gè)探針分別安裝在樣樹離地面1.3 m處。使用DT80數(shù)據(jù)采集器（Datataker，澳大利亞）設(shè)置采樣程序記錄并存儲(chǔ)數(shù)據(jù),數(shù)據(jù)采樣間隔為30 min。使用Granier公式計(jì)算樹干液流密度。Fd=0.000 119×［（ΔTmax–ΔT）/ΔT］1.231×3 600。 其中：Fd是樹干液流密度（m3·m－2·h－1）， ΔTmax為 24 h 內(nèi)無液流時(shí)上下探針之間的最大溫差值，ΔT為瞬時(shí)溫差值。依據(jù)樹干液流密度推算單株單位時(shí)間液流量Q（kg·h－1）的計(jì)算公式為：Q=Fd×As×103。 其中：As為邊材面積（m2）， 根據(jù)時(shí)間尺度可以換算為日總液流量（kg·d－1）。將1 d中太陽輻射為0時(shí)段內(nèi)的樹干液流定義為夜間液流。

在各樣樹上緊鄰熱擴(kuò)散探針處安裝D6樹木生長儀（UMS，德國，精度5 μm，工作溫度－30～50℃，單通道4線橋式應(yīng)變傳感器）。該儀器安裝簡(jiǎn)單，無需破壞樹皮或影響樹木的生長，可連續(xù)、高精度測(cè)量樹木的莖周長變化。D6樹木生長儀由記錄器和樹木生長應(yīng)變傳感器組成，樹木莖周長的變化能直接傳遞到傳感器，瞬時(shí)記錄樹木生長變化對(duì)環(huán)境影響的反應(yīng)等。傳感器受溫度影響小，其與樹皮之間有一個(gè)特富龍層，可減少傳感器和樹皮的摩擦，同時(shí)也能不受結(jié)冰、樹脂或結(jié)疤的影響。數(shù)據(jù)記錄在另一臺(tái)DT80數(shù)據(jù)采集器中。

表1 樣樹基本特征Table 1 Characteristics of the monitored trees

1.3 邊材面積計(jì)算及環(huán)境因子

利用生長錐取樣計(jì)算邊材面積。為避免生長錐取樣對(duì)樹體造成破壞，影響液流觀測(cè)結(jié)果，本研究在樣樹附近選擇其他植株進(jìn)行取樣。河北楊由于數(shù)量較少，故選擇與樣樹徑級(jí)較一致的各3株河北楊鉆取木芯，取平均值來代替樣樹的邊材厚度。對(duì)油松取樣發(fā)現(xiàn)，其邊材厚度與胸徑（DBH）之間存在極顯著的二次函數(shù)關(guān)系：y=0.062x2－0.707x+4.482（R2=0.880 3，n=10），其中y為邊材厚度（cm）；x為油松胸徑（cm）。在此基礎(chǔ)上計(jì)算油松樣樹的邊材面積。距離樣地附近的HOBO自動(dòng)氣象觀測(cè)站記錄太陽輻射（W·m－2）， 氣溫（℃）， 降水量（mm）， 風(fēng)速（m·s－1）和空氣相對(duì)濕度（%）等。 在各樣樹下使用 HOBO 土壤水分傳感器（S-SMC-M005，美國）長期測(cè)定土壤體積含水量，觀測(cè)深度為30 cm，數(shù)據(jù)采樣間隔為30 min。

1.4 數(shù)據(jù)分析

使用Excel軟件對(duì)樹干液流等數(shù)據(jù)進(jìn)行整理和計(jì)算，所有統(tǒng)計(jì)分析均采用SPSS 18.0軟件完成，作圖采用Origin 9.0軟件。

2 結(jié)果與分析

2.1 樹干液流密度動(dòng)態(tài)變化

目標(biāo)樹種典型晴天的樹干液流密度變化（圖1）顯示：河北楊液流密度日變化特征以單峰型為主，在黎明前保持極低水平，6:00－8:30樹干液流啟動(dòng)，在11:00－15:00達(dá)到峰值。7月觀測(cè)數(shù)據(jù)顯示，液流密度在午間出現(xiàn)短暫下降，導(dǎo)致呈現(xiàn)雙峰現(xiàn)象。之后樹干液流密度迅速減小，在17:00－21:30達(dá)夜間低值水平。各月夜間均存在微弱的液流。圖2顯示：生長季中期的6，7和8月液流密度明顯高于9和10月，7月白天液流密度明顯高于其他月份，且曲線幅度比其余月份寬。日均液流量大?。?.69 kg·d－1（7月）＞2.90 kg·d－1（6 月）＞2.63 kg·d－1（8 月）＞1.47 kg·d－1（9 月）＞1.16 kg·d－1（10 月）。 油松的液流密度在8:30迅速上升，在11:00－15:00達(dá)到峰值，結(jié)束時(shí)間在18:00－22:00，比河北楊晚1 h左右。7月也出現(xiàn)雙峰現(xiàn)象，各月存在微弱的夜間液流活動(dòng)。日均液流量大小分別為5.74 kg·d－1（7 月）＞3.30 kg·d－1（6 月）＞3.29 kg·d－1（8 月）＞2.06 kg·d－1（10 月）＞1.76 kg·d－1（9 月）。

在晴天（7月30日）和雨天（8月1日）的天氣下，2種樹種的液流變化曲線相似（圖3）。河北楊在晴天和雨天條件下液流密度峰值分別為（0.126 0±0.046 0）和（0.112 8±0.052 0） m3·m－2·h－1； 油松液流密度峰值分別為（0.120 3±0.044 0）和（0.099 8±0.016 0） m3·m－2·h－1。 晴天白天液流密度明顯高于清晨和夜晚，變化幅度大；雨天變化曲線平緩，雨天液流密度峰值出現(xiàn)明顯下降。河北楊雨天比晴天下降41.67%，油松下降10.48%。但雨天夜間液流密度較高，可能與白天太陽輻射減少且降雨補(bǔ)充了林下土壤水分有關(guān)。

圖1 河北楊與油松典型晴天的液流密度變化Figure 1 Daily course of sap flow density of each species

圖2 河北楊與油松樹干液流量變化（平均值±標(biāo)準(zhǔn)誤）Figure 2 Seasonal means of daily sap flow

2.2 夜間液流

觀測(cè)期間發(fā)現(xiàn)2種樹種的夜間液流明顯，只有在極少數(shù)的情況下液流密度為0。6－8月河北楊液流密度低且平穩(wěn)，凌晨1:00開始逐漸上升，且夜間液流主要發(fā)生在后半夜（0:00－6:00），9－10月河北楊夜間液流活動(dòng)較強(qiáng)。河北楊6－10月夜間液流量占全天總液流量的比例分別為：9.7%，4.8%，14.2%，8.0%和23.8%。觀測(cè)日平均比例為12.1%。在季節(jié)之間沒有顯著差異。6－8月油松夜間都保持較低液流密度，9和10月波動(dòng)較大，夜間液流活動(dòng)明顯增強(qiáng)，可能與土壤水分增加有關(guān)。6－10月夜間液流量占全天總液流量的比例分別為：9.5%，7.2%，11.4%，5.5%和16.6%。觀測(cè)日平均比例為10.0%（圖4）。

圖3 河北楊與油松晴天、雨天樹干液流密度變化Figure 3 Daily course of sap flow density in a sunny day （30 July） and a rainy day（1 August） of monitored trees

2.3 液流密度與環(huán)境因子的關(guān)系

環(huán)境因子與河北楊和油松樹干液流密度之間的相關(guān)分析（表2）表明：河北楊液流密度與太陽輻射和氣溫極顯著正相關(guān)（P＜0.01），與風(fēng)速顯著正相關(guān)（P＜0.05），而與相對(duì)濕度極顯著負(fù)相關(guān)（P＜0.01），與其他環(huán)境因子的相關(guān)關(guān)系不顯著（P＞0.05）。對(duì)于油松，液流密度與太陽輻射、風(fēng)度、氣溫和土壤含水量顯著正相關(guān)（P＜0.01），而與相對(duì)濕度和氣壓顯著負(fù)相關(guān)（P＜0.01）。圖5比較了8月典型晴天的液流密度均值在氣象因子較一致的情況下對(duì)土壤含水量下的響應(yīng)，表現(xiàn)為土壤含水量從10.7%下降到4.9%導(dǎo)致河北楊液流密度均值減少27.0%，油松減少33.0%。

圖4 河北楊與油松夜間液流量占日總液流量的比例Figure 4 Ratio of nocturnal sap flow to total daily sap flow in different month

圖5 河北楊與油松在晴天不同土壤含水量的樹干液流密度Figure 5 Sap flow density of monitored trees during the different soil moisture

表2 油松和河北楊液流密度與環(huán)境因子的相關(guān)性Table 2 Pearson correlation coefficient of the monitored trees sap flow density to environmental factors

由于樹體自身存在著巨大的儲(chǔ)水性，樹體內(nèi)部的水分消耗能通過樹干莖周長變化反映出來。為了分析2種樹種莖周長日變化與液流活動(dòng)的關(guān)系，繪制液流密度與莖周長的日變化過程圖（圖6）。結(jié)果表明：日出前，2種樹種存在較弱的夜間液流活動(dòng)，莖周長表現(xiàn)出微弱的下降。日出之后即太陽輻射開始增強(qiáng)，植物耗水量加大，2種樹種液流密度迅速升高，而莖周長卻迅速減小。觀測(cè)發(fā)現(xiàn)莖周長變化總的趨勢(shì)表現(xiàn)為先降低再升高后降低，與液流密度變化規(guī)律相反。其次，從表3發(fā)現(xiàn)：2種樹種的樹干收縮期正好處在液流活動(dòng)期間內(nèi)，在整個(gè)研究期間均存在這一現(xiàn)象。相關(guān)分析表明：2種樹種的液流密度與莖周長日變化呈顯著負(fù)相關(guān)關(guān)系（河北楊：R2=-0.50，P＜0.01；油松：R2=-0.36，P＜0.01），說明樹木莖周長的日變化與樹干液流活動(dòng)緊密相關(guān)。

圖6 河北楊與油松液流密度與莖周長日變化及太陽輻射的日變化Figure 6 Diurnal variation of daily sap flow density and daily perimeter variation,solar radiation for Populus hopeiensis and Pinus tabulaeformis

表3 河北楊與油松液流密度與莖周長日變化特征Table 3 Diurnal variation of daily sap flow density and daily perimeter variation for Populus hopeiensi and Pinus tabulaeformis

3 討論

本研究中河北楊和油松的樹干液流密度晝夜變化相似，主要呈單峰型曲線，與其他樹種一致［12-15,22-29］，差異表現(xiàn)在啟動(dòng)時(shí)間、峰值和下降時(shí)間的不同。晝夜變化規(guī)律是因?yàn)闃涓梢毫鞯漠a(chǎn)生主要是由于蒸騰拉力的作用，影響蒸騰作用強(qiáng)度的因子就會(huì)影響樹干液流密度，在土壤水分不受限的條件下，液流密度就隨著蒸騰強(qiáng)度的增加而升高［22-26］。本研究結(jié)果也證實(shí)了這一點(diǎn)。6-8月液流密度均值與峰值均較高，與太陽輻射和氣溫具有較強(qiáng)正相關(guān)性。其次，王華田等［24］認(rèn)為：樹木日液流密度最大值可能僅與土壤水分狀況及樹木自身蒸騰耗水的生物學(xué)、生理學(xué)特性有關(guān)。而KANALAS等［27］的研究發(fā)現(xiàn)：土壤水分的減少及干旱脅迫的加劇會(huì)令無梗花櫟Quercus petraea產(chǎn)生獨(dú)特的氣孔調(diào)節(jié)，使最大液流密度從7-8月下降81%。在本研究中，土壤水分的減少，河北楊液流密度均值減少27%，油松減少33%，觀測(cè)結(jié)果基本符合上述規(guī)律。另外，本研究發(fā)現(xiàn)10月油松最大液流密度不是出現(xiàn)在胸徑最大的樣樹上，而是出現(xiàn)在胸徑較小的樣樹上。趙平等［29］認(rèn)為對(duì)于相同種類的樹木，由于生長過程中所處的環(huán)境不同，如：不同的方位、林內(nèi)風(fēng)向與風(fēng)速、接受太陽輻射的強(qiáng)弱等，均會(huì)引起液流密度的差異，也可能與樹種自身水分利用效率有關(guān)［29-30］。

液流密度的日變化主要還是受環(huán)境因子的影響，例如土壤水分、太陽輻射、氣溫、風(fēng)速和相對(duì)濕度等因素［3,7,9,12-13,25-31］。徐世琴等［32］認(rèn)為：液流與風(fēng)速的相關(guān)性較大主要是由于風(fēng)速增大能夠增加植物冠層邊界層導(dǎo)度，從而加速葉片蒸騰。本研究也表明，河北楊液流密度與太陽輻射和氣溫極顯著正相關(guān)，與風(fēng)速顯著正相關(guān)，而與相對(duì)濕度顯著負(fù)相關(guān)。油松液流密度與太陽輻射、風(fēng)速、氣溫和土壤含水量顯著正相關(guān)，而與相對(duì)濕度和氣壓顯著負(fù)相關(guān)。相關(guān)程度絕對(duì)值前3位的依次為太陽輻射＞氣溫＞相對(duì)濕度，說明太陽輻射對(duì)液流活動(dòng)的驅(qū)動(dòng)作用最強(qiáng)。這與溫杰等［13］的研究較一致。有研究表明：夜間液流在一定程度上反映了植物體內(nèi)水分對(duì)由蒸騰作用的響應(yīng)［33-34］。本研究發(fā)現(xiàn)：降雨過后，河北楊和油松的夜間液流密度迅速增加。NAKAI等［35］認(rèn)為：這是干旱區(qū)植物適應(yīng)干旱的方式，植物增加夜間液流活動(dòng)補(bǔ)充水分來應(yīng)對(duì)干旱環(huán)境，且能反映出植物木質(zhì)部的導(dǎo)水效率。其次，本研究發(fā)現(xiàn)2種樹種夜間液流量對(duì)日總液流量的比例為4.8%～23.8%，王艷兵等［34］認(rèn)為，夜間液流量所占日總液流量的比例是由土壤水分狀況和樹木生長狀況共同決定的，最大值出現(xiàn)在10月。本研究觀測(cè)最大值也出現(xiàn)在10月，造成這種現(xiàn)象的原因可能與7-9月雨季過后土壤水分得到較大補(bǔ)充有關(guān)。

本研究在對(duì)液流密度與莖周長變化的相關(guān)分析中，發(fā)現(xiàn)2種樹種莖周長變化與液流密度變化呈顯著負(fù)相關(guān)關(guān)系，與前人研究結(jié)果一致［16-19］。ROMà等［1］發(fā)現(xiàn)樹干收縮變化與樹干液流有關(guān)，認(rèn)為蒸騰作用消耗樹干儲(chǔ)水會(huì)導(dǎo)致樹干徑向尺寸減小，日落后蒸騰作用減少，樹干開始儲(chǔ)水會(huì)使樹干膨脹，這種雙向變化反映植物對(duì)環(huán)境條件變化的響應(yīng)。而造成這種現(xiàn)象的原因可能由植物體內(nèi)水分利用策略決定的，當(dāng)相關(guān)性大時(shí)，樹干水分調(diào)節(jié)能力強(qiáng)，而當(dāng)相關(guān)性小時(shí)，水分調(diào)節(jié)能力弱［14-15,19］。FORNER等［3］發(fā)現(xiàn)植物在缺水的情況下，植物對(duì)蒸騰耗水過程和生長發(fā)育過程的響應(yīng)是高度結(jié)合的，關(guān)于兩者的關(guān)系需要在后續(xù)的試驗(yàn)中進(jìn)一步研究。

4 結(jié)論

對(duì)半干旱黃土丘陵區(qū)河北楊和油松生長季樹干液流的研究發(fā)現(xiàn)，2種樹種液流密度的日變化過程主要呈 “晝高夜低”單峰或雙峰型曲線，樹干液流密度主要受太陽輻射、氣溫和相對(duì)濕度影響。河北楊日間液流密度峰值較高，夜間液流活動(dòng)明顯，意味著河北楊具備高效的導(dǎo)水結(jié)構(gòu)。油松在觀測(cè)期間能保持較高的液流量，與氣象因子的相關(guān)性較高。2種樹種的液流活動(dòng)能夠指示樹木莖周長白天收縮和夜間膨脹的日變化。