李 雙，杜建括

(陜西理工大學(xué)歷史文化與旅游學(xué)院，陜西 漢中 723001)

我國荒漠區(qū)植被以耐旱草本和灌木為主。了解干旱荒漠區(qū)灌木個(gè)體和群叢對有限水資源的消耗量，對于研究荒漠植被的生長狀況、個(gè)體和群落尺度的水量平衡至關(guān)重要。目前，基于植物個(gè)體和林分水平研究植物耗水過程的方法很多。測定林分水平耗水的方法大致可以分為兩類，一類是直接測定的方法，主要包括空氣動(dòng)力學(xué)法、渦動(dòng)相關(guān)法、Penman-Monteith方程法、水分平衡法、能量平衡-波文比法和紅外遙感法[1]；另一類方法是借助一定的尺度轉(zhuǎn)換參數(shù)，由單株/木耗水推求林分水平的耗水信息[2]?；诔叨绒D(zhuǎn)換的林分蒸騰耗水量通常可以依據(jù)液流與胸徑[3,4]、莖干橫截面積[5,6]、邊材面積[7]、葉面積[1,2]的關(guān)系以及植被種植密度[8]和單株木占地面積[9]等進(jìn)行尺度擴(kuò)展后求得。選擇何種尺度擴(kuò)展方法推求林分蒸騰主要取決于植被群落特征[10]。在個(gè)體尺度水平上，測定植物耗水的方法主要有整樹容器法、風(fēng)調(diào)室法、盆栽法、蒸滲儀法、液流測定法(同位素示蹤法和熱力學(xué)方法)。熱力學(xué)技術(shù)實(shí)現(xiàn)了對植物莖干液流實(shí)時(shí)、準(zhǔn)確、高時(shí)間分辨率的自動(dòng)監(jiān)測[11,12]，已在草本、農(nóng)作物、喬木和灌木耗水研究方面得到廣泛應(yīng)用[13]。

灌木沒有明顯主干或主干很短，熱平衡液流儀通常監(jiān)測的是灌木單枝耗水量，需要借助合適的尺度轉(zhuǎn)換參數(shù)來推導(dǎo)和獲取單株、灌叢群落的耗水量。已有學(xué)者通過莖流與莖干橫截面積[6]、葉面積[1,2,14]的函數(shù)關(guān)系推算了灌木單株和群落的蒸騰耗水量，取得了較好的實(shí)驗(yàn)結(jié)果。在干旱荒漠區(qū)，灌木等植被的分布相對稀疏，陽光照射時(shí)冠層葉片之間相互遮擋的程度不高，灌叢群落耗水量與葉片總面積之間存在著極強(qiáng)的相關(guān)性[2]。但與灌叢葉片總面積相比，葉片總干重更容易直接獲取，若荒漠植被蒸騰耗水與葉片干重顯著相關(guān)，則葉干重可作為擴(kuò)展參數(shù)推求單株、灌叢尺度的耗水量。

近些年來，在我國西北荒漠植物耗水與生態(tài)需水量計(jì)算方面已取得一些成果[15-19]，但基于莖干液流量與植物形態(tài)參數(shù)估算小尺度范圍檉柳耗水的研究很少。本文選擇甘肅省景泰縣寺灘鄉(xiāng)典型荒漠灌木多枝檉柳作為研究對象，在單枝莖流測定基礎(chǔ)上，以葉干重作為擴(kuò)展變量來推算單株蒸騰耗水量，并對樣地檉柳總耗水量作估算，旨在尋求一種適合小尺度范圍或種植布局和長勢較統(tǒng)一的灌叢單株/群落蒸騰耗水量的測量方法和理論，為科學(xué)評價(jià)荒漠區(qū)水分平衡及人工植被管理提供一定依據(jù)。

1 研究區(qū)概況

甘肅省景泰縣寺灘鄉(xiāng)位于黃土高原與騰格里沙漠過渡地帶，屬于溫帶大陸性干旱氣候，海拔1 828 m(37°14′N,103°48′E)，年均氣溫8.6 ℃，極端最高氣溫為38.6 ℃，極端最低氣溫-27.3 ℃，≥0 ℃的年活動(dòng)積溫3 614.8 ℃，≥10 ℃的年有效積溫3 038 ℃，無霜期120 d左右，年日照時(shí)數(shù)2 723.9 h。年平均降雨量180 mm，降水分配不均，夏秋多雨，冬春干旱，約60%的降水集中在7-9月，年潛在蒸發(fā)量3 038 mm。風(fēng)沙日數(shù)較多，平均風(fēng)速2.0～3.1 m/s，瞬時(shí)最大風(fēng)速可達(dá)21.7 m/s。土壤類型以灰鈣土和灰棕荒漠土為主。

為了保護(hù)脆弱的生態(tài)環(huán)境，研究區(qū)為2003年退耕后種植的多枝檉柳(以下簡稱檉柳)生態(tài)建設(shè)林用地。檉柳行株距為4×2 m。林地內(nèi)草本植物主要有堿蓬(Suaedaglauca)和灰藜(Chenopodium album Linn.)等。檉柳的平均株高170 cm，莖干基徑平均為3.2 cm。本區(qū)地下水埋藏較深，無人工灌溉檉柳生長的水源主要來自于大氣降水和土壤水。

2 研究方法

2.1 莖干液流測定

在研究區(qū)選擇一塊50 m×50 m的樣地進(jìn)行植被調(diào)查和莖流監(jiān)測。選擇長勢良好、健康的檉柳作為樣本，使用熱平衡莖流儀(Flow 32, Dynamax Inc., Houston, TX, USA)監(jiān)測它們的液流速率變化。使用游標(biāo)卡尺測量莖干直徑范圍以選取合適的莖流儀探頭型號。2013年6-9月、2014年6月和8月均在研究區(qū)樣地開展了檉柳莖干液流監(jiān)測。數(shù)據(jù)采集器(CR1000, Campbell Scientific, Logan, UT, USA)每10 s掃描一次，每1 min平均一次，每6 min記錄一次數(shù)據(jù)。實(shí)驗(yàn)期間用到的探頭型號有SGA10、SGA13、SGB16、SGB19和SGB25，每個(gè)探頭有2～3個(gè)重復(fù)。

2.2 地上生物參數(shù)的測定

2014年6月和8月分別選取檉柳個(gè)體各6株，共計(jì)12株作為標(biāo)準(zhǔn)木，在其上安裝熱平衡液流儀傳感器。每一塊樣地內(nèi)地形平坦，同一立地條件下絕大多數(shù)檉柳個(gè)體差異不顯著。測量所選樣株的冠幅大小(東西、南北長)、株高、所有枝條的基徑和長度。每株檉柳樣本中選取3個(gè)枝條，其中至少有1個(gè)來自包裹著傳感器的枝條，共計(jì)36枝，采集其上的所有葉片，每枝的葉片單獨(dú)裝入樣品袋，帶回實(shí)驗(yàn)室，于80 ℃鼓風(fēng)恒溫烘箱中烘至恒重，獲得葉片干物質(zhì)量。其中24枝樣品的葉干重及枝條基徑、長度用于建立單枝葉干重與枝條基徑、長度的關(guān)系模型，其余12枝用于模型檢驗(yàn)。

(1)檉柳單枝葉干重估算模型。參考董道瑞等[20]關(guān)于塔里木河下游檉柳地上生物量估算方法，根據(jù)單枝檉柳的基徑和長度，用線性、冪、對數(shù)和指數(shù)4種函數(shù)形式來建立檉柳單枝葉干物質(zhì)量的估算模型，并用決定系數(shù)R2和F檢驗(yàn)對擬合的回歸方程進(jìn)行評判，從中選出最優(yōu)估算模型。

w=a(d2h)+b

(1)

w=a(d2h)b

(2)

w=aln (d2h)+b

(3)

w=aebd2h

(4)

式中：w為單枝葉干重，g；d為單枝檉柳的基徑，cm；h為單枝檉柳的長度，cm；a、b均為系數(shù)。

(2)檉柳單株葉干重估算模型。根據(jù)12株樣本所有枝條的基徑和長度數(shù)據(jù)，利用從公式(1)～(4)中選出的最優(yōu)單枝葉干重估算模型獲得每株檉柳葉干總干重，其表達(dá)式為：

Wi=niwi(i=1,2,…,12)

(5)

利用公式(5)計(jì)算的單株葉片總干重，結(jié)合每株樣本的冠幅及株高，建立單株葉片總干重與冠幅周長和株高之間的函數(shù)關(guān)系式，其表達(dá)式為：

W=a+b(P2H)c

(6)

P=2 πs+4 (l-s)

(7)

式中：W為單株檉柳葉片總干重，g；P為冠幅周長，m；H為株高，m；s為冠幅的短半徑，m；l為冠幅的長半徑，m；a、b和c均為系數(shù)。

2.3 快速稱重法測蒸騰

試驗(yàn)地檉柳快速離體稱重法參考楊劼等[21]離體稱重測定人工灌木蒸騰強(qiáng)度的方法，在包裹液流儀探頭的檉柳植物附近選擇一株長勢和大小較一致的檉柳，在植株冠層中部四周剪取長5～10 cm的枝條6～8個(gè)，剪下后立即用精度為0.001 g的電子天平稱重，記為W1，并記錄稱重時(shí)間；稱重后將枝條迅速放回原處，使離體枝條處在原來的環(huán)境下蒸騰，3 min后取下迅速重新稱重W2，兩次重量之差即為蒸騰耗水量。本實(shí)驗(yàn)設(shè)計(jì)6次重復(fù)以減少蒸騰失水誤差，每間隔1 h做一組實(shí)驗(yàn)。將第二次稱重后的枝條摘下葉片(包括嫩梢)，并裝入已知重量的信封，帶回實(shí)驗(yàn)室烘干樣品，稱重，凈葉片干重用Wd表示。利用以下公式即可求出檉柳蒸騰速率 。

Tr=60 (W1-W2)/Wd/t

(8)

式中：Tr為蒸騰速率，g/(g·h)；W1為離體枝葉的初始重量，g；W2為離體枝葉的t分鐘后的重量，g；t為兩次稱重的蒸騰時(shí)間間隔，min；Wd為葉片凈干重,g。

3 結(jié)果與分析

3.1 檉柳莖干液流與形態(tài)參數(shù)的關(guān)系

為了同一尺度水平上比較不同單位表示的液流速率變化的差異性，單位橫截面積和單位葉干重表示的液流速率分別乘以系數(shù)2和100，這樣3種單位表達(dá)的液流速率數(shù)值均介于0～300之間(圖1)。無論用何種單位表示液流速率，3個(gè)莖干的液流變化趨勢相同，但3種單位表示的不同徑級枝條液流速率間的差異明顯不同，例如，用單位時(shí)間通過橫截面的液流量表示的檉柳不同徑級間的液流速率相差顯著[圖1(a)]；用單位時(shí)間通過2個(gè)單位橫截面積的液流量表示的各莖干間液流值差異略微縮小，尤其是枝條S1與S2的液流速率相近[圖1(b)]；用單位時(shí)間每100 g葉干重的液流量表示的不同莖干液流變化的曲線彼此間幾乎相互重疊[圖1(c)]，說明檉柳莖干液流與枝條葉干重密切相關(guān)，液流量主要是由枝條葉片量決定的，經(jīng)葉干重標(biāo)準(zhǔn)化處理后可以很好地消除不同枝莖干液流的差異。

圖1 不同單位表示的檉柳不同枝條液流速率變化的日過程

為進(jìn)一步明確檉柳液流量與生態(tài)學(xué)參數(shù)的相關(guān)性，對其進(jìn)行了線性回歸分析。檉柳7∶00至19∶00累計(jì)液流量與莖干邊材面積、橫截面周長和面積的相關(guān)度R2分別為0.739、0.669和0.679(P<0.001)[圖2(a)～圖2(c)]，而與葉干重存在極其顯著的正相關(guān)(R2=0.965，P<0.001)[圖2(d)]。雖然檉柳液流量與枝條斷面參數(shù)整體上呈正相關(guān)，但也可能出現(xiàn)例外，譬如，直徑12.09 mm莖干的液流速率始終大于直徑12.16 mm莖干的液流速率，表明液流與莖干直徑?jīng)]有很好的一致性，類似結(jié)果也存在于液流與莖干橫截面積之間。圖1和圖2揭示了與莖干直徑和橫截面積相比，葉干重作為一個(gè)基礎(chǔ)參數(shù)用于植物耗水研究將更加準(zhǔn)確和有效。因此，本文將葉干重作為擴(kuò)展參數(shù)來獲取整株或灌叢尺度的耗水量是可行的。

3.2 檉柳單枝/株葉片總干重的估算

3.2.1 檉柳單枝葉片總干重回歸模型

選取12株檉柳中的24個(gè)單枝樣本，在SPSS17.0統(tǒng)計(jì)分析軟件中利用公式(1)～(4)構(gòu)建4種函數(shù)形式的生物量估測回歸模型(見表1)。經(jīng)查F值表，F(xiàn)0.01(1，22)=7.94，結(jié)合表1可知，4種擬合函數(shù)均通過顯著性檢驗(yàn)，且線性函數(shù)的R2值最高，故認(rèn)為估算檉柳單枝葉片總干重的最優(yōu)模型為：

w=0.074 (d2h)+14.538 (9)

表1 4種函數(shù)擬合的檉柳單枝葉片總干重回歸方程統(tǒng)計(jì)分析表

注：**表示顯著性水平P<0.01。

用其余的12個(gè)單枝樣本對線性模型進(jìn)行檢驗(yàn)，均方根誤差RMSE為2.014 g，實(shí)測值與模擬值的斜率為0.937，接近1(見圖3)。所以，利用基徑平方與長度乘積建立的枝條葉片總干重線性模型是可信的。

圖2 檉柳日液流量與莖干周長、橫截面積、邊材面積及葉干重的關(guān)系

圖3 檉柳單枝葉片總干重實(shí)測值與模擬值的比較(虛線表示1∶1關(guān)系)

3.2.2 檉柳單株葉片總干重回歸模型

根據(jù)公式(5)計(jì)算12株檉柳的單株葉片總干重，再根據(jù)公式(6)、(7)基于冠幅周長和株高建立單株葉片總干重估算回歸模型如下：

W=3.540P2H+89.877 (R2=0.938，F(xiàn)=182.472) (10)

為了檢驗(yàn)?zāi)Ｐ偷木?，另外選取了6株檉柳用于實(shí)測值與模擬值的誤差分析，見表2。由公式(10)獲得的模擬值與實(shí)測值的絕對誤差介于29.036～67.741 g，誤差率為11.709%～52.499%，其中植株1的誤差最大，植株5的誤差最小，總體表現(xiàn)為植株葉片總干重越小，由公式(10)模擬誤差越大?？梢?，當(dāng)樣地植株葉片總干重跨度較大時(shí)，不適合用同一個(gè)公式表達(dá)單株葉片總干重與冠幅周長和株高的關(guān)系。為此，我們將植株葉片總干重劃分為兩個(gè)量級分別構(gòu)建模型。

當(dāng)W≤200 g，即P2H值小于30.3 m3時(shí)：

W=46.238 e0.049P2H(R2=0.992，F(xiàn)=500.921)

(11)

當(dāng)W>200 g，即P2H值大于30.3 m3時(shí)：

W=379.1 ln (P2H) (R2=0.986，F(xiàn)=341.226)

(12)

與公式(10)模擬結(jié)果相比，由公式(11)、(12)獲得的模擬值與實(shí)測值的絕對誤差顯著減小(表2)，平均誤差率分別為13.141%和8.008%，即模型精度達(dá)86.86%和91.99%，表明分級建立檉柳單株葉片總干重估算模型提高了模型估算精度。分級估測模型更符合植株葉片總干重與體積的實(shí)際關(guān)系：同一環(huán)境中，灌木體積越大，葉片總干重也越大，但當(dāng)灌木冠幅增加到一定程度后，葉干重、冠幅體積隨冠幅增大而增加的程度變小。

表2 檉柳單株葉片總干重實(shí)測值與模擬值比較 g

3.3 檉柳蒸騰耗水量的尺度轉(zhuǎn)換

圖2(d)的線性方程僅是由一天的枝條日液流量與葉干重?cái)M合得到的。為了提高擬合的合理性，對12株樣本中的15個(gè)枝條晴天時(shí)連續(xù)3 d的日液流量(7∶00-19∶00)平均值與葉干重的關(guān)系進(jìn)行擬合，得到如下線性回歸方程估算：

fl=14.486w+47.571 (R2=0.948，F(xiàn)=298.476)

(13)

式中：fl為檉柳單枝日液流量，即單枝日蒸騰耗水量，g；w為單枝葉片總干重，g。

用另外包裹著熱平衡液流儀的3個(gè)枝條對公式(13)進(jìn)行驗(yàn)證(見圖4)。由圖4可知莖流儀實(shí)測值與模擬值相差不大，相對誤差分別為5.358%、5.929%和4.544%，平均誤差為5.277%，說明用公式(13)估算檉柳單枝耗水量是可信的。

圖4 檉柳3個(gè)單枝莖干日液流量實(shí)測值與模擬值的比較

借助合適的尺度轉(zhuǎn)換參數(shù)可由熱平衡液流儀監(jiān)測的單枝耗水量來推導(dǎo)和獲取單株、灌叢群落的耗水量。先根據(jù)公式(9)和(13)計(jì)算12株樣本檉柳單株日液流量，再結(jié)合公式(11)和(12)計(jì)算每株的葉片總干重，建立檉柳單株液流量與葉片總干重之間的函數(shù)關(guān)系，其表達(dá)式為：

fl=15.605w+93.049 (R2=0.943，F(xiàn)=275.834)

(14)

檉柳單株耗水量線性回歸方程在P<0.01水平上的相關(guān)系數(shù)R2=0.943，F(xiàn)=275.834，達(dá)到極顯著線性相關(guān)，并且單枝耗水量估算模型[公式(13)]的估算精度在94%以上，說明利用葉干重構(gòu)建檉柳單枝、單株耗水模型是可信的。

為了進(jìn)一步驗(yàn)證檉柳單株耗水模型的可信度，隨機(jī)選擇3株小型檉柳樣本用于檢驗(yàn)。將同一時(shí)期檉柳單株日蒸騰耗水量模型計(jì)算的結(jié)果與快速稱重法擴(kuò)展獲取的同株耗水量值作比較。3株檉柳各自的葉片總干重由公式(11)或(12)估算。根據(jù)公式(8)，快速稱重法測定的3株檉柳日蒸騰量為7∶00-19∶00蒸騰速率與葉片總干重乘積的累加值。由表3可知，兩種方法的平均誤差為11.657%。單因素方差分析顯示兩種方法測定結(jié)果的Levene方差齊性檢驗(yàn)顯著水平為0.943，遠(yuǎn)遠(yuǎn)大于0.05，說明兩組測定結(jié)果在同一水平下方差一致；差異性水平為F=0.054，P=0.828>0.05，表明兩種方法測定的檉柳單株耗水量之間沒有顯著差異，故認(rèn)為用葉片總干重作為變量估算檉柳蒸騰耗水量是可行的。

表3 兩種方法測定結(jié)果的相對誤差

考慮到本文的單枝/株葉片總干重是由檉柳形態(tài)參數(shù)易測因子枝條基徑與長度，或植株冠幅與株高估測的，所以把公式(13)與(9)合并得公式(15)，公式(14)與公式(11)、(12)分別合并得公式(16)和(17)，將晴天檉柳單枝/株的耗水量回歸方程轉(zhuǎn)換為單枝/株耗水量與易測因子(枝條基徑與長度、植株冠幅和株高)之間的關(guān)系，各表達(dá)式如下：

fl=1.072 (dh)+258.168

(15)

當(dāng)PH<30.3 m3時(shí)：

fl=724.544 e0.049PH+93.049

(16)

當(dāng)PH>30.3 m3時(shí)：

fl=5 915.856 ln(PH)-18 729.702

(17)

本文的研究區(qū)為退耕還林地，檉柳的行株距均為4×2 m，所以通過估算檉柳標(biāo)準(zhǔn)木的平均日蒸騰耗水量，再乘以50 m×50 m樣地內(nèi)檉柳的株數(shù)，即可獲得樣地內(nèi)檉柳的日耗總量。50 m×50 m樣地范圍內(nèi)應(yīng)有338株檉柳，除去死亡的21株，實(shí)際共有317株。經(jīng)計(jì)算12個(gè)標(biāo)準(zhǔn)木的日平均耗水量為5.320 kg，那么50 m×50 m樣地內(nèi)檉柳日耗水總量約為1 686.44 kg，即約0.675 mm，生長季5至10月的累計(jì)耗水量約124.2 mm。另外，對灌木種植布局不統(tǒng)一，且植株個(gè)體差異較大的小尺度樣地而言，也可通過測量樣地內(nèi)每株植物的冠幅、株高，通過公式(16)和(17)計(jì)算出所有植株的日耗水量，再求和得到樣地灌木總耗水量。

4 結(jié) 語

(1)植物生態(tài)學(xué)參數(shù)是影響植物耗水特性的主要因素之一。本文以單位時(shí)間通過橫截面的液流量、單位時(shí)間通過單位橫截面的液流量和單位時(shí)間單位葉干重的液流量分析了檉柳形態(tài)學(xué)參數(shù)對不同徑級莖干液流的影響。無論用何種單位表示液流速率，3個(gè)莖干的液流變化趨勢相同，但3種單位表示中僅經(jīng)葉干重標(biāo)準(zhǔn)化處理后可以很好地消除不同枝莖干液流間的差異。同時(shí)線性回歸分析表明，檉柳莖干液流與葉干重極顯著相關(guān)，在P<0.001水平上相關(guān)度R2達(dá)0.965，而液流量與莖干邊材面積、橫截面周長和面積的相關(guān)度R2分別為0.739、0.669和0.679(P<0.001)。之所以與葉干重顯著相關(guān)，是因?yàn)槿~片是植物進(jìn)行蒸騰作用的主要器官，葉片數(shù)量直接關(guān)系到蒸騰量，且干旱區(qū)灌木分布相對稀疏葉片之間相互遮擋度不高。

(2)以枝條基徑和長度建立的檉柳單枝葉片總干重的最優(yōu)模型為：w=0.074(d2h)+14.538(R2=0.902，P<0.01)；以P2H=30.3 m3為分界點(diǎn)，利用冠幅周長與株高建立的單株檉柳葉片總干重估算模型分別為W=46.238 e0.049P2H(R2=0.992，F(xiàn)=500.921)和W=379.1 ln(P2H)-1 206.2(R2=0.986，F(xiàn)=341.226)，精度分別達(dá)86.86%和91.99%。董道瑞等[20]利用冠幅周長與株高建立的塔里木河下游檉柳生物量估算模型精度為80.09%；彭守章等[22]利用相同參數(shù)建立的黑河下游檉柳生物量估算模型精度達(dá)90%以上。本文的估算模型精度明顯高于前者，而與后者相似，主要是因?yàn)橥画h(huán)境中，當(dāng)灌木冠幅增加到一定程度后，葉干重、株高隨冠幅增大而增加的程度變小，所以在建立模型時(shí)應(yīng)依據(jù)灌木體積分級構(gòu)建，可提高模型精度。

(3)基于葉干重的單枝、單株檉柳耗水模型均達(dá)到了顯著水平，模擬值與實(shí)測值的平均誤差分為為5.277%和11.657%，說明基于葉干重的檉柳耗水估算模型是可信的。結(jié)合葉干重與冠幅、株高的函數(shù)關(guān)系，單株檉柳耗水模型可表述為當(dāng)P2H<30.3 m3時(shí)，fl=724.544 e0.049 P2H+93.049；當(dāng)P2H>30.3 m3時(shí)，fl=5 915.856 ln(P2H)-18 729.702。根據(jù)模型估算研究區(qū)檉柳生長季5-10月的累計(jì)耗水量約124.2 mm，該值處于彭守章等[22]計(jì)算的檉柳多年平均耗水量在30～386 mm研究結(jié)果的范圍內(nèi)。

本文用于建立估算模型的數(shù)據(jù)為2014年6月和8月晴天時(shí)的液流數(shù)據(jù)和同時(shí)期的植物生態(tài)學(xué)參數(shù)，因此推算的檉柳生長季耗水量與實(shí)際值可能存在誤差。借助熱平衡液流儀可實(shí)時(shí)監(jiān)測植物莖干液流變化的特點(diǎn)，建立植物不同生長階段、不同天氣條件下的檉柳單枝/株液流量與葉干重或者形態(tài)特征易測因子的關(guān)系，則有助于更加準(zhǔn)確地掌握檉柳整個(gè)生長季的蒸騰耗水量的動(dòng)態(tài)變化。