劉澤勇，馬長明，劉春鵬，李向軍，徐振華，楊 飛

（1.河北省林業(yè)科學(xué)研究院 河北省林木良種工程技術(shù)研究中心，河北 石家莊 050061；2.河北豐寧沙地生態(tài)系統(tǒng)國家定位觀測研究站，河北 承德 068357；3.河北農(nóng)業(yè)大學(xué) 林學(xué)院，河北 保定 071000；4.柏鄉(xiāng)縣林業(yè)局，河北 邢臺 055450）

干旱缺水一直是我國北方干旱半干旱地區(qū)植被恢復(fù)的主要限制因子[1-2]，尤其是石質(zhì)山區(qū)，由于土層淺薄、氣候干旱、保水能力低下，一直是北方植被恢復(fù)的難點，經(jīng)常出現(xiàn)“年年造林不成林”的現(xiàn)象[3]。對于干旱區(qū)的植被恢復(fù)，節(jié)水型造林樹種的選擇以及造林密度的配置尤為重要，而這兩方面均與植物自身的耗水特性有直接關(guān)系，因此研究植物的耗水特性對干旱半干旱區(qū)植被的快速恢復(fù)具有極其重要的意義。

目前，有關(guān)植物蒸騰耗水的研究很多，研究方法各有不同，其中熱擴散探針法（TDP）由于具有操作簡便、精確度高、能夠連續(xù)監(jiān)測且對植物傷害較小等眾多優(yōu)點，成為國內(nèi)外學(xué)者廣泛采用的方法之一[4-8]。縱觀植物蒸騰耗水的相關(guān)研究，研究方向大多集中于樹干液流的變化規(guī)律、樹干液流與環(huán)境因子的關(guān)系等方面[9-11]。而對于構(gòu)建耗水模型的研究較少，有些研究雖然構(gòu)建了耗水模型，但對耗水模型的可靠性缺乏驗證[12-13]，因此植物耗水模型的構(gòu)建與驗證是其研究發(fā)展的一個新方向。

山杏Prunus sibirica作為我國北方的鄉(xiāng)土樹種，具有抗旱、耐瘠薄、少有病蟲害等優(yōu)點，并且兼有生態(tài)和經(jīng)濟兩大效益，一直是我國北方干旱地區(qū)造林的首選樹種之一[14]。研究其耗水特性，構(gòu)建其耗水預(yù)測模型，對華北石質(zhì)山區(qū)造林具有現(xiàn)實指導(dǎo)意義。因此，本研究以山杏為研究對象，利用莖流儀和小型自動氣象站，對其樹干液流和環(huán)境因子進行同步監(jiān)測，采用多元逐步回歸的方法，構(gòu)建其耗水預(yù)測模型，并對其可靠性進行檢驗，以期為石質(zhì)山區(qū)山杏的造林提供理論依據(jù)和技術(shù)支持。

1 試驗區(qū)概況

本試驗于河北平山縣（113°31′～114°15′E，38°9′～38°45′N）寺家溝村進行。該地區(qū)位于河北省西部、太行山東麓，屬于低山丘陵區(qū)，海拔一般在500 m以下。氣候?qū)儆谂瘻貛О霛駶櫦撅L(fēng)型大陸性氣候，立體氣候明顯，季節(jié)分明，光照充足。年均氣溫12.7 ℃，歷史最高溫度達(dá)40.8℃，歷史最低溫度達(dá)-16.4 ℃；全年平均太陽輻射量為131～136千卡/cm2，日照時數(shù)為2 600～2 750 h；年均降水量約為500 mm，蒸發(fā)量為1 815.4 mm，干燥度為1.38；年均風(fēng)速為2 m/s，無霜期為185 d；土壤類型為石灰性褐土。

試驗地設(shè)置于陽坡坡中位置，周圍主要植被有喬木樹種山杏Prunus sibirica、黃連木Pistacia chinensis、刺槐Robinia pseudoacacia、君遷子Diospyros lotus等；灌木主要樹種有荊條Vitex negundovar.heterophylla、酸 棗Zizyphus jujuba等；草本主要有黃背草Themeda triandra、白羊草Bothriochloa ischaemum、蒺藜Tribulus terrestris等。

2 材料與方法

2.1 試驗材料的選擇

由于受莖流儀數(shù)據(jù)采集線的限制，要求在試驗點方圓40 m的范圍內(nèi)，選擇胸徑大小相近、無病蟲害、生長良好的山杏3株作為試驗材料，經(jīng)篩選試驗材料的胸徑大小分別為14.1、12.6、11.1 cm，冠幅分別為4.9、4.7、4.0 m，枝下高分別為1.2、1.3、1.6 m。

2.2 樹干液流的測定

采用熱擴散式液流探針（TDP）法，利用RR-8210莖流儀（國內(nèi)北京雨根科技有限公司設(shè)計生產(chǎn)）對選定的山杏進行樹干液流的測定，監(jiān)測時間為2009年4月21日—11月8日。安裝前，于山杏樹干正南方向胸徑（距地面1.3 m）處，上下垂直，削去長度為10～15 cm的一小片樹皮，削割深度剛好削去韌皮部而不傷木質(zhì)部為宜；然后，用直徑為2.5 mm的鉆頭依靠打孔板打2個孔，打孔深度略大于10 mm，將10 mm的探針插入孔中，用泡沫盒子罩上，外部用反光紙包上，進而防止自然環(huán)境對測量數(shù)據(jù)的影響；將探針通過數(shù)據(jù)采集線與數(shù)據(jù)采集器相連，連接電源，設(shè)置每隔2 s測讀1次，每隔10 min自動平均并儲存數(shù)據(jù)。

2.3 環(huán)境因子的測定

在試驗地范圍內(nèi)，采用美國SPECTRUM設(shè)計生產(chǎn)的watchdog900ET小型全自動氣象站，對林內(nèi)環(huán)境因子進行同步監(jiān)測，即同樣設(shè)置每2 s測讀一次，每10 min進行平均并儲存數(shù)據(jù)。監(jiān)測的指標(biāo)有太陽輻射（wat/m2）、大氣溫度（℃）、空氣濕度（%）、20 cm土壤溫度（℃）、20 cm土壤濕度（%）、風(fēng)速（km/h），其中太陽輻射、大氣溫度、空氣濕度以及風(fēng)速監(jiān)測裝置距地面高度為1.5 m。

2.4 邊材面積的測定

試驗結(jié)束后，于監(jiān)測樹木探針安裝處，利用生長錐鉆取木芯，觀察量取去皮直徑和邊材厚度，并計算邊材面積。經(jīng)計算，監(jiān)測樹木的邊材面積分別為97.2、86.9、76.6 cm2。

2.5 數(shù)據(jù)處理

2.5.1 不同月份耗水特征分析

由于樹干液流受天氣變化的影響很大，陰雨天氣往往會造成樹干液流的不規(guī)律變化，因此，本研究將各月所有晴天的數(shù)據(jù)進行平均后，分析討論山杏不同月份的耗水特征。選擇的晴天日期為：4月，25—27日、29—30日；5月，1—8日、12—13日、17—20日、22—26日、29日、31日；6月，1—2日、4—5日、9—12日、20—25日、28—30日；7月，1—4日、7日、9—10日、12日、14日、18日、20日、22日、25—26日、28日；8月，2日、10—15日、23日、26—28日、30—31日；9月，9—13日、15—17日；10月，1—5日、12—14日、28日；11月，1—6日。

2.5.2 不同月份耗水與環(huán)境因子關(guān)系及耗水模型的構(gòu)建

于各月選擇典型連續(xù)的三個晴天（天空晴朗無云、無大風(fēng)），將樹干液流與各環(huán)境因子進行相關(guān)性分析，并以樹干液流速度為因變量，以各環(huán)境因子為變量，進行逐步回歸分析，以5 %和10 %作為變量入選和剔除臨界值，進而得到山杏不同月份的耗水預(yù)測模型。選擇的日期分別為：4月，25—27日；5月，5—7日；6月，9—11日；7月，2—4日；8月，11—13日；9月，11—13日；10月，3—5日；11月，4—6日。

2.5.3 不同月份耗水模型的驗證

利用未參與建模的3個典型晴天（盡量連續(xù)）的環(huán)境因子數(shù)據(jù)，利用預(yù)測模型計算樹干液流速度以及耗水量，并于實測值進行比較，進而評價預(yù)測模型的可靠性。選擇的用于驗證耗水模型的晴天為：4月，22、29—30日；5月，8、12—13日；6月，12、20—21日；7月，7—10日；8月，10、14—15日；9月，15—17日；10月，12—14日；11月，1—3日。

將獲取的數(shù)據(jù)利用EXCEL 2003對觀測數(shù)據(jù)進行整理、計算和作圖。利用SPSS 13.0中的相關(guān)分析處理程序，進行樹干液流與各環(huán)境因子間的簡單相關(guān)分析；利用多元回歸分析處理程序，進行以樹干液流速度為因變量，以各環(huán)境因子為變量的多元逐步回歸分析。

3 結(jié)果與分析

3.1 不同月份耗水特征分析

如圖1和表1所示，各月份樹干液流均呈現(xiàn)“晝高夜低”的變化規(guī)律，白天樹干液流較大，夜間樹干液流較小，但依然存在。除6月樹干液流呈現(xiàn)“雙峰”外，其他各月份樹干液流均呈現(xiàn)“單峰”。液流啟動時間呈現(xiàn)“晚-早-晚”的趨勢，4—6月逐漸變早，7—11月逐漸變晚；液流達(dá)到峰值時間4月份較早，5—8月較晚，9—11月更晚；而低谷出現(xiàn)時間并未表現(xiàn)出明顯的規(guī)律。各月蒸騰耗水呈現(xiàn)“低-高-低”的趨勢，其中以5月最大，之后逐月降低，至11月份，其樹干液流僅為5月的10%。在整個生長季中，山杏耗水集中于4—7月，其耗水量占整個生長季的78.2%。

圖1 不同月份的樹干液流變化Fig.1 The fluctuation of sap flow velocity in different months

表1 不同月份的蒸騰耗水進程?Table 1 The process of transpiration water consumption in different months

3.2 不同月份耗水與環(huán)境因子的關(guān)系

表2匯總了山杏各月份樹干液流與各環(huán)境因子的相關(guān)性。山杏各月份樹干液流均與太陽輻射、大氣溫度、空氣濕度三者呈現(xiàn)極顯著相關(guān)性，表現(xiàn)出與太陽輻射和大氣溫度呈正相關(guān)、與空氣濕度呈負(fù)相關(guān)。就相關(guān)性大小而言，除5月外，其他月份的樹干液流與太陽輻射、大氣溫度和空氣濕度的相關(guān)性最大，而5月樹干液流與大氣溫度、空氣濕度和風(fēng)速的相關(guān)性最大。樹干液流與土壤溫度、土壤濕度和風(fēng)速并沒有表現(xiàn)出穩(wěn)定的相關(guān)性，但就其與風(fēng)速的相關(guān)性而言，除了8月份以外，其他月份樹干液流均與風(fēng)速表現(xiàn)出或大或小的極顯著正相關(guān)性。

3.3 不同月份耗水模型的構(gòu)建

表3列出了山杏各月樹干液流與環(huán)境因子的逐步回歸方程。如表3所示，各月份山杏耗水預(yù)測模型的調(diào)整R2有所差異，4—7月調(diào)整R2較高，均達(dá)到了0.85以上，說明4—7月山杏耗水預(yù)測模型擬合度較高，其中以5月擬合度最高，調(diào)整R2達(dá)到0.935。而8—11月耗水預(yù)測模型的調(diào)整R2一般在0.63～0.78之間，說明4—7月山杏耗水預(yù)測模型擬合度較低，其中以10月份最低，其調(diào)整R2僅為0.635。另外，山杏各月份耗水預(yù)測模型經(jīng)F檢驗均達(dá)到了顯著或極顯著水平。

表2 不同月份樹干液流與環(huán)境因子的相關(guān)系數(shù)?Table 2 Correlation coefficients between sap flow velocity and environmental factors in different months

表3 不同月份耗水預(yù)測模型?Table 3 The water-consumption prediction models in different months

3.4 不同月份耗水模型的驗證

利用山杏各月份耗水預(yù)測模型，選擇未參與建模的3個典型晴天，通過監(jiān)測的各環(huán)境因子指標(biāo)，計算山杏各月份各時間點的樹干液流速度，并與實測值進行比較。圖2所示，各月份樹干液流速度的實測值與預(yù)測值雖然存在一些偏差，但變化趨勢基本重合。通過預(yù)測模型計算山杏各月份典型晴天的日均耗水量與實測值進行t檢驗，發(fā)現(xiàn)二者不存在顯著性差異（P＞0.05）。統(tǒng)計各月日均耗水量預(yù)測值與實測值，計算相對誤差率，山杏各月份典型晴天日均耗水量預(yù)測值與實測值的相對誤差率均在10%以下（表4），按生物統(tǒng)計的要求，10%以內(nèi)的誤差都是可接受的[15]，因此，說明本研究構(gòu)建的山杏各月份耗水預(yù)測模型具有一定的可靠性。

4 討 論

整個生長季節(jié)，山杏樹干液流均呈現(xiàn)“晝高夜低”的規(guī)律，夜間液流較小，但仍然存在。這與眾多[11,16-17]研究結(jié)果一致。眾多研究一致認(rèn)為夜間液流存在的原因在于“根壓”。除6月份山杏樹干液流呈現(xiàn)“雙峰”外，其他各月份樹干液流均呈現(xiàn)“單峰”。筆者認(rèn)為，6月份山杏處于果實成熟期，水分需求較大，而此時期降水較少，土壤水分很難滿足其需求，因此植物在蒸騰耗水最大的中午，采用閉合部分氣孔的措施，進行暫時休眠，減少自身的蒸騰，進而達(dá)到節(jié)水的目的。這與劉龍等[18]以及劉瀟瀟等[19]的觀點相近。山杏樹干液流啟動4月份較晚，5—6月份逐漸變早，7—11月份逐漸變晚；樹干液流的啟動與氣象因子的變化有關(guān)，太陽輻射、大氣溫度與空氣濕度均伴隨日出而發(fā)生變化，日出早其各氣象因子開始變化的時間就早，樹干液流啟動就早。4—6月日出時間逐漸變早，7月份（夏至之后）日出時間逐漸變晚，進而導(dǎo)致樹干液流啟動出現(xiàn)“晚-早-晚”的變化趨勢。這與任啟文等[20]的觀點一致。山杏各月蒸騰耗水主要集中在4—7月。這主要是因為4—7月是山杏開花、快速生長、結(jié)果的時期，生長旺盛，水分需求大，耗水量隨之增大；之后山杏生長逐漸進入衰退期，水分需求變小，其耗水量隨之減少，這與高峻等[21]的觀點一致。

圖2 樹干液流預(yù)測值和實測值比較Fig.2 Comparison between predicted and measured value of sap flow velocity

表4 耗水預(yù)測模型誤差評價Table 4 Error evaluation of water-consumption prediction models

山杏各月份樹干液流主要與太陽輻射、大氣溫度、空氣濕度三者極顯著相關(guān)，表現(xiàn)出與太陽輻射和大氣溫度呈正相關(guān)、與空氣濕度呈負(fù)相關(guān)。這與陳婧等[13]的觀點一致。另外，除了8月份以外，山杏樹干液流均與風(fēng)速也表現(xiàn)出或大或小的極顯著正相關(guān)性。有關(guān)樹干液流與風(fēng)速間相關(guān)性的研究存在兩種不同的觀點，許文豪等[10]與Kitaya Y等[22]認(rèn)為二者存在顯著相關(guān)性，而Fisher J B等[23]和李廣德等[24]的觀點恰恰相反。本研究結(jié)論恰處于兩種觀點之間，筆者認(rèn)為與太陽輻射、大氣溫度和空氣濕度等氣象因素相比，風(fēng)速變化規(guī)律性較弱，且經(jīng)常沒有連續(xù)性，因此，其與樹干液流的相關(guān)性不穩(wěn)定。

本研究對各月山杏的樹干液流和各環(huán)境因子進行了逐步回歸分析，得到了山杏各月耗水預(yù)測模型，其中4—7月耗水模型調(diào)整R2達(dá)到了0.85以上，說明擬合度較好。而8—11月耗水預(yù)測模型調(diào)整R2較低，擬合度稍弱。筆者認(rèn)為，在生長旺盛的時期（4—7月），山杏的蒸騰耗水對環(huán)境因子的反應(yīng)更靈敏，因而模型擬合度較高；而生長衰退時期，其蒸騰耗水對環(huán)境因子的響應(yīng)遲鈍，因此擬合度較低。

5 結(jié) 論

整個生長季節(jié)，山杏樹干液流表現(xiàn)出“晝高夜低”的變化規(guī)律，除6月呈現(xiàn)“雙峰”外，其他各月份均呈現(xiàn)“單峰”。液流啟動時間呈現(xiàn)“晚-早-晚”的趨勢，達(dá)到峰值時間4月份較早，5—8月較晚，9—11月更晚。各月蒸騰耗水呈現(xiàn)“低-高-低”的趨勢，且多集中于4—7月。

山杏各月份樹干液流與太陽輻射、大氣溫度、空氣濕度均呈現(xiàn)極顯著相關(guān)性，且相關(guān)性最大，表現(xiàn)出與太陽輻射和大氣溫度呈正相關(guān)、與空氣濕度呈負(fù)相關(guān)。另外，除8月份外，山杏樹干液流與風(fēng)速也表現(xiàn)出或大或小的極顯著正相關(guān)性。

通過多元逐步回歸分析，本研究得到了山杏各月的耗水預(yù)測模型，經(jīng)F檢驗均達(dá)到極顯著水平，調(diào)整R2達(dá)到了0.85以上，擬合良好。

通過檢驗，山杏各月份日均耗水量預(yù)測值與實測值的誤差率均在10%以下，說明山杏各月耗水模型均具有一定的可靠性。

綜上所述，本研究所構(gòu)建的山杏耗水模型具有一定的實際應(yīng)用性。在整個生長季節(jié)，可以根據(jù)不同時間段內(nèi)的各環(huán)境因子，利用耗水預(yù)測模型，對山杏耗水情況進行預(yù)測，并結(jié)合降雨情況，以水循環(huán)和平衡原理為基礎(chǔ)，計算出該地區(qū)山杏最佳的造林密度，進而為山杏的造林實踐提供理論依據(jù)和技術(shù)支持。