西北干旱區泡泡刺灌叢的降雨再分配特征及影響因素分析

趙文玥,吉喜斌,金博文,焦丹丹,張靖琳,郭 飛,趙麗雯

1 中國科學院西北生態環境資源研究院臨澤內陸河研究站,中國科學院生態水文與流域科學重點實驗室, 蘭州 730000 2 中國科學院大學, 北京 100049 3 北京師范大學地理科學學部,地表過程與資源生態國家重點實驗室, 北京 100875

冠層降雨再分配是指降雨進入植被冠層后被重新分配為冠層截留損失、穿透雨和樹干莖流的過程[1—2]。其中冠層截留損失是指被冠層攔截、儲存,在降雨過程中或降雨結束后以蒸發形式返回大氣的雨水[3—4],是地表蒸散發的重要組分[5],也是冠層降雨再分配過程中水分的凈損失項[6]。穿透雨是指透過冠層到達地面的降雨,包括從冠層縫隙直接滴落地面的自由穿透雨和從冠層枝葉滴落地面的釋放穿透雨[7],穿透雨約占大氣降雨的50%—80%[8—10],是冠層下表層土壤水分的主要來源[11—12]。樹干莖流是指沿樹枝、樹干向樹根流動,最終以點源的形式到達深層土壤中的雨水[13—14],雖數量較少,但集流率[15]和富集率[16]極高,是冠層下“濕島”[13, 17]和“沃島”[13, 18]效應形成發展的重要促進因子。穿透雨和樹干莖流合稱凈降雨量,是最終到達土壤可供冠層下植被利用的有效降雨量[1, 19]。因此,冠層降雨再分配對局地水文過程和水文功能具有重要影響,是局地水文循環中不可忽視的重要環節[20]。

在干旱區生態系統中,受水分和養分資源的限制,植被在空間上通常呈現斑塊或條帶狀稀疏分布[21—22]。降雨是該地區天然生態系統的主要補給水源,其數量稀少、時空分布不均且變率較大,以小降雨事件為主[2, 23]。稀疏植被冠層對于降雨的再分配作用改變了降雨特性及其時空分配格局[20],對冠層下土壤水分、養分的時空分布[12—13, 24]、生物群落的演替與分布[25]、地下水的補給[26]、區域降雨循環[27]等一系列生態水文過程產生影響。多年來,關于干旱區稀疏植被冠層降雨再分配的研究取得了重要進展,主要集中于不同種類稀疏植被冠層降雨再分配的比例構成[28—29]、變化趨勢[30]、影響因素分析[4, 29, 31]及模型構建[32—34]等。但由于干旱區植被形態結構特殊[14]、分布稀疏[2]且具有顯著的空間異質性[35],其降雨再分配過程的觀測實驗較為困難,相關長時間尺度的研究較少；在構建稀疏植被冠層截留模型時,一般基于水量平衡的原理估算冠層截留損失,假設其數值上等于林外降雨量與凈降雨量之差[7],但多數忽略了樹干莖流分量的測算；在探討植被形態結構因子對冠層降雨再分配過程的影響時,多基于葉面積指數,但對干旱區稀疏植被的觀測表明葉片的持水能力低于枝條和樹干,所以枝條和樹干在冠層降雨再分配的過程中發揮重要作用[11, 36],因此在構建干旱區稀疏植被冠層截留模型及探討相關影響因素時,應對此特征加以考慮,僅使用葉面積可能不能真實準確地反映植被冠層對降雨再分配過程的調節作用。

泡泡刺(Nitrariasphaerocarpa)為西北干旱區特有的古老荒漠殘遺種,屬超旱生灌木,是西北干旱區河西走廊中段綠洲—荒漠過渡帶的天然建群種植物,通常以灌叢沙堆的形式存在[37],在防風固沙、穩定綠洲環境等方面具有重要生態意義[38—39]。本文以泡泡刺灌叢為研究對象,通過2016年、2017年和2020年生長季對泡泡刺灌叢的冠層降雨再分配過程、植被形態特征、林外降雨特征、氣象因素進行野外定位觀測,分析了泡泡刺灌叢冠層降雨再分配的比例構成及主要影響因素(包括降雨特征、植被形態結構特征和氣象因素),構建了冠層降雨再分配要素隨降雨特征變化的經驗回歸方程,揭示了泡泡刺冠層下穿透雨的空間分布規律及樹干莖流的聚集效應。以期增進對于干旱區稀疏植被覆蓋下冠層降雨再分配過程及其影響機理的定量認識,為干旱區類似植被分布區水文循環過程中冠層截留損失分量的測算和水分收支評估提供可以借鑒的方法,以促進綠洲—荒漠過渡帶的植被恢復和生態建設實踐。

1 材料與方法

1.1 研究區概況

野外觀測實驗地位于中國科學院臨澤內陸河研究站(39°20′50″N, 100°07′48″E)綠洲-荒漠過渡帶生態水文過程試驗區,地處河西走廊中段,北鄰巴丹吉林沙漠,南鄰黑河中游臨澤綠洲,海拔1384 m。該地區屬典型的溫帶大陸性氣候,夏季高溫少雨,冬季寒冷干燥,年平均氣溫9.2℃,多年平均降雨量為124 mm,主要集中在6—9月,降雨變率大,以小降雨事件為主,年均蒸發量2390 mm[40—41]。

研究區的植被群落大多呈現斑塊狀稀疏分布,植被結構簡單,主要建群灌木有梭梭(Haloxylonammodendron)、泡泡刺(N.sphaerocarpa)、沙拐棗(Calligonummongolicum),同時零散分布著花棒(Hedysarumscoparium)、多枝怪柳(Tamarixramosissima)等,在降雨較多的年份霧冰藜(Bassiadasyphylla)、白莖鹽生草(Halogetonarachnoideus)、小畫眉草(EragrostisminorHost)等一年生草本發育。泡泡刺灌叢約占研究區總植被蓋度的45%,是該地區典型的天然稀疏植被[40]。

1.2 植被的選擇及相關植被要素的測量

首先,對研究區泡泡刺灌叢的分布位置、基徑、株高等展開植被調查[40]；隨后,在研究區隨機選擇1堆典型的、長勢良好的泡泡刺灌叢作為研究對象,并測定了影響其冠層降雨再分配過程的冠層結構參數,包括株高、冠層投影面積、樹干傾角(指枝條與水平方向的夾角)、植被面積指數(Plant area index,PAI)等(具體植被特征見表1)。其中株高、冠層厚度、基徑、樹干傾角在2016年8月、2017年7月、2020年6月觀測實驗開始前測定。冠層投影面積在2020年8月利用無人機拍攝照片,隨后運用Photoshop 2020軟件處理,再通過MATLAB R2018b計算得到。植被面積指數是指單位土地面積上冠層總表面積的一半[42],由葉面積指數、枝面積指數和莖面積指數等共同組成[43],是對植被結構和覆蓋狀況的整體描述[44—45]。研究區泡泡刺灌叢的葉面積僅占其冠層表面積的一部分,枝、莖等對雨水的截留再分配作用同樣不可忽視。因此,采用PAI能更好地反映泡泡刺冠層對降雨再分配過程的調節作用。PAI數據由2016年、2017年、2020年8月運用LAI- 2200C植物冠層分析儀(LI-COR Inc, Lincoln, Nebraska, USA)在冠層下集雨瓶所在位置貼近地表逐個測定而來。

表1 泡泡刺灌叢冠層形態結構特征的基本信息

1.3 林外降雨特征及其他氣象要素的測定

林外降雨特征由布設在泡泡刺灌叢周邊開闊地帶的翻斗式雨量筒(TE525MM, Texas Electronics Inc., Dallas, TX)測定,通過數據采集器(CR1000-XT, Campbell Scientific, UT, USA)自動記錄數據,每30 min計算平均值并儲存,精度0.1 mm。同時分別在泡泡刺灌叢0°、120°、240°方向的空地布設3個雨量筒,用于雨量訂正。雨水的收集以次降雨事件為基礎,若降雨間隔小于6 h,則視為一次降雨事件。降雨結束后立刻進行雨量的測定,若降雨發生在夜間,則第二天清晨進行收集,以減小截留雨水蒸發對觀測造成的影響。其他氣象要素包括風速、風向、溫度、濕度、凈輻射等由林外空地的自動氣象站實時監測[40]。

1.4 穿透雨的測量及計算

圖1 泡泡刺灌叢穿透雨和樹干莖流收集裝置示意圖Fig.1 The device for collecting rainfall redistribution in field, with details of homemade stemflow collectors and throughfall gauges

穿透雨通過自制的集雨瓶收集,集雨瓶由內徑6 cm、高13.5 cm的玻璃杯上部疊加內徑11 cm、外緣高2 cm的聚乙烯漏斗組成(圖1)。將集雨瓶以灌叢中心為圓心,呈輻射狀均勻放置于北、東北、東、東南、南、西南、西、西北8個方向的半徑上,相鄰兩個集雨瓶間隔20 cm,共布設93個集雨瓶(圖1)。由于泡泡刺灌叢的枝條大多呈現匍匐狀貼近沙丘表面,故將集雨瓶部分埋入土壤,頂部露出地表5—8 cm,以避免集雨瓶碰觸枝條,同時減小風對集雨瓶擺放位置的影響。穿透雨量的計算方法如下：

式中,TF是穿透雨量(mm)；VT為第i個集雨瓶中穿透雨的體積(mL)；m是冠層下集雨瓶的數量；FA為漏斗的橫截面積(cm2)。

1.5 樹干莖流的測量及計算

在灌叢8個方位上各選擇2株,共16株泡泡刺作為標準枝收集樹干莖流。具體步驟如下：首先,將小漏斗(內徑3.2 cm)側面鉆孔剖開,固定在樹干接近地表處；然后,為防止穿透雨的濺入,在漏斗敞口處覆蓋鉆孔后的鋁箔紙,鋁箔紙上圓孔的直徑稍大于樹干直徑4—5 mm確保樹干莖流可以流入；最后,用導水管連接漏斗排水口和收集瓶,并用硅膠密封接縫處(圖1)。借鑒標準枝法[46]推算整個灌叢產生的樹干莖流體積,即用標準枝的平均樹干莖流體積乘整個灌叢的枝條數。樹干莖流量的計算公式如下：

式中,SF為樹干莖流量(mm)；VS為整個泡泡刺灌叢的樹干莖流體積(mL)；SP為整個泡泡刺灌叢的冠層投影面積(m2)。此外,Herwitz等[15]提出集流率的概念,用于定量分析樹干莖流對水分的聚集作用,具體公式如下：

式中,FR為集流率；VS為樹干莖流的體積(mL)；BA為樹干基部的橫截面積(cm2)；Pg為林外降雨量(mm)。集流率即表示流經單位樹干橫截面積上的樹干莖流量與林外降雨量的比值。集流率大于1,表示除了樹干之外,冠層的其他組分對樹干莖流的產生和聚集具有促進作用[15, 17]。

1.6 冠層截留損失的計算

由于冠層截留損失難以直接測量,通常根據水量平衡的原理進行估算,其數值上等于總降雨量與凈降雨量之差：

IL=Pg-N=Pg-(TF+SF)

式中,IL為冠層截留損失量(mm)；Pg為林外降雨量(mm)；N為凈降雨量(mm)；TF為穿透雨量(mm)；SF為樹干莖流量(mm)。

1.7 數據處理

利用SPSS 22進行描述性統計和Pearson相關分析,Pearson相關分析的顯著性檢驗采用“雙尾檢驗”；使用Origin 2018進行回歸分析、曲線擬合及制圖,回歸分析和曲線擬合后,采用單因素方差分析對回歸系數進行顯著性檢驗,顯著水平設定為α=0.05；運用ArcGIS 10.2對穿透率和PAI進行克里金插值并繪制其空間分布圖。

2 結果與分析

2.1 試驗期間研究區的降雨特征

試驗期間(2016年8月、2017年7月、2020年6—9月)共觀測到可產生穿透雨的有效降雨事件23次,總降雨量為90.9 mm,單次平均降雨量為(4.0 ± 4.5)mm,最小值為0.2 mm,最大值為17.9 mm。將降雨量分為<2 mm、2—5 mm、5—10 mm、10—20 mm 4個雨量級,頻數分別為12、5、4、2,累積降雨量分別為12.3 mm、17.5 mm、31.0 mm、30.1 mm(圖2)。其中<2 mm的降雨事件發生次數最多,但其累積降雨量最小,對總降雨量的影響小。>5 mm的降雨事件僅發生6次,但其累積降雨量占總雨量的比例高達67.2%,對總降雨量的貢獻較大。實驗期間降雨強度的波動較小,平均雨強為(0.8 ± 0.5)mm/h,最小值為0.2 mm/h,最大值為2.6 mm/h。雨強<0.5 mm/h的降雨事件發生次數最多,共發生9次,但其累積降雨量僅占總降雨量的9.7%。雨強介于0.5—1.5 mm/h范圍內的降雨事件對總降雨量的貢獻最大。雨強>1.5 mm/h的降雨事件發生頻次最少,對總降雨量的貢獻也相對較小。實驗期間總降雨歷時為97 h,平均降雨歷時為(4.2 ± 3.7)h。歷時<2 h的降雨事件發生頻率最高,共發生11次,但其累積降雨量最小,約占總降雨量的15.6%。歷時>10 h的降雨事件很少發生,僅發生2次,但其累積降雨量占總降雨量的比例最大,約為33.1%。總體來看,實驗期間研究區以低雨強、短歷時的小降雨事件為主,但該類降雨事件對總降雨量的貢獻較小,主要降雨量源于發生次數較少的中雨強、長歷時的大降雨事件。

圖2 實驗期間降雨量、降雨強度、降雨歷時分布特征Fig.2 Characteristics of rainfall amount, rainfall intensity and rainfall duration distribution during the experiment

2.2 穿透雨與降雨特征和冠層結構特征的關系

試驗期間,泡泡刺灌叢的累計穿透雨量為78.89 mm,約占總降雨量的87.89%。以單次降雨事件來看,泡泡刺灌叢穿透雨量的變化范圍在0.14—16.68 mm之間,平均值約為(3.47 ± 4.16)mm。穿透雨量占次降雨量的比例(簡稱穿透率)介于73.00%—98.69%,平均值約為(84.95 ± 7.28)%。穿透雨的變異系數在14.45%—77.31%之間波動,平均值約為(29.76 ± 12.61)%。

降雨特征是影響泡泡刺灌叢穿透雨特征的重要因素,通過回歸分析和曲線擬合,建立了泡泡刺灌叢穿透雨量和穿透率與降雨特征的回歸方程(圖3)。由圖3所示,穿透雨量與降雨量、降雨歷時之間呈顯著線性正相關關系(P<0.0001),與降雨強度之間呈現冪函數關系(P<0.01)。表明穿透雨量隨降雨量、降雨強度、降雨歷時的增加而增加,當雨強>1.5 mm/h后,穿透雨量的增加速率加快,可能是由于高雨強的降雨事件更易造成雨水的飛濺,加速截留雨水向釋放穿透雨轉化導致。根據穿透雨量與降雨量的線性回歸方程推算,泡泡刺灌叢穿透雨產生的閾值約為0.2 mm,與實測結果0.2 mm基本一致。穿透率與降雨歷時之間呈線性正相關關系(P= 0.04),與降雨量、降雨強度之間無顯著相關性(圖3)。

圖3 穿透雨量、穿透率、穿透雨的空間異質性與降雨特征之間的回歸關系Fig.3 The correlation between throughfall, throughfall percentage, coefficient of throughfall spatial variability and rainfall regimesTF: 穿透雨量 Throughfall; TP: 穿透率 Throughfall percentage; CV: 穿透雨空間變異系數 Coefficient of variation; Pg: 降雨量 Rainfall amount; RI: 降雨強度 Rainfall intensity; RD: 降雨歷時 Rainfall duration; n: 樣本數量 Sample number

冠層結構特征通過改變穿透雨的傳輸路徑對冠層下穿透雨的多寡產生不可忽視的影響,相關性分析表明,平均穿透雨量與株高和PAI之間呈現極顯著的負相關關系(P<0.0001),通過對其進行曲線擬合發現線性函數的擬合效果最好(P<0.0001)(圖4)。即穿透雨量隨株高和PAI的增加而下降,可能是由于當株高和PAI較高時,雨水更易被冠層截留,進而減少了穿透雨的形成。

圖4 平均穿透雨量與株高、PAI之間的回歸關系Fig.4 The correlation between mean throughfall and shrub height, PAIMT: 單個穿透雨量筒的平均穿透雨量 Mean throughfall; SH: 株高 Shrub height; PAI: 植被面積指數 Plant area index

穿透雨的空間分布具有明顯的空間異質性,利用穿透雨的變異系數來反映穿透雨空間分布的異質性,將其與降雨特征之間進行回歸分析,結果表明,隨降雨量、降雨強度、降雨歷時的增加,泡泡刺灌叢穿透雨的空間異質性逐漸減弱,對數函數的擬合效果較好(P<0.01)(圖3)。為進一步探索穿透雨空間分布的特征,運用克里金插值繪制了不同雨量級下穿透率的空間分布圖。如圖5所示,隨降雨量的增加,穿透雨的空間異質性逐漸減弱,高強度的降雨事件也有可能減弱穿透雨的空間異質性,與回歸分析的結果一致。此外,由圖5可知穿透率的空間分布呈現多個“雨極”和“旱區”鑲嵌分布的模式,“雨極”通常出現在灌叢外緣,而“旱區”則通常出現在灌叢中部,推測可能是由于PAI的空間異質性導致,因此,運用克里金插值繪制了PAI的空間分布圖,發現PAI的高值區對應穿透率的“旱區”,而低值區則對應“雨極”。

圖5 泡泡刺灌叢穿透率和PAI的空間分布圖Fig.5 Spatial distribution of throughfall percentage and PAI for Nitraria sphaerocarpaTP: 穿透率 Throughfall percentage; PAI: 植被面積指數 Plant area index

2.3 樹干莖流與降雨特征和冠層結構的關系

試驗期間共觀測到樹干莖流13次,總樹干莖流量為1.46 mm,約占總降雨量的1.61%。從單次降雨事件來看,泡泡刺灌叢樹干莖流量的變化范圍在0.003—0.41 mm,平均值約為(0.11 ± 0.11)mm。樹干莖流量占次降雨量的比例(簡稱樹干莖流率)介于0.20%—2.41%之間,平均值約為(1.55 ± 0.78)%。

降雨特征是影響樹干莖流形成過程的重要因素,通過對泡泡刺灌叢的樹干莖流量和樹干莖流率與降雨特征進行回歸分析和曲線擬合,其變化趨勢見圖6。樹干莖流量與降雨量、降雨歷時之間呈顯著的線性正相關關系(P<0.001),與降雨強度之間呈指數函數關系(P<0.01),表明樹干莖流量先隨雨強的增加而迅速增加,當雨強>1 mm/h后,樹干莖流的增加速率逐漸減緩,可能是由于高雨強的降雨事件更易加劇截留雨水的迸濺,進而減少截留雨水向樹干莖流轉化導致。根據樹干莖流量與降雨量之間的回歸方程推算樹干莖流產生的閾值約為1.1 mm,與實際觀測值1.5 mm近似。樹干莖流率隨降雨量、降雨強度和降雨歷時的增加呈對數函數的增長趨勢(P<0.01),當降雨量<6 mm、降雨強度<1.5 mm/h、降雨歷時<6 h時,樹干莖流率隨降雨量、降雨強度和降雨歷時的增加而迅速增加,但超過這一范圍后,樹干莖流率逐漸趨于穩定。推測這一現象可能與泡泡刺灌叢樹干的最大輸水能力有關,當雨量、雨強、降雨歷時較大,超過樹干最大輸水能力時,多余的雨水可能以穿透雨的形式滴落,進而減緩了樹干莖流率的增加速率。

圖6 樹干莖流量、樹干莖流率及集流率與降雨特征之間的回歸關系Fig.6 The correlation between stemflow, stemflow percentage, funneling ratio and rainfall regimesSF: 樹干莖流量 Stemflow; SP: 樹干莖流率 Stemflow percentage; FR: 集流率 Funneling ratio

冠層結構特征同樣在樹干莖流形成的過程中發揮不可忽視的作用,單株泡泡刺產生的樹干莖流量與冠層結構特征的相關分析結果見表2。泡泡刺的樹干莖流量與樹干傾角之間呈現顯著的正相關關系(P=0.005),即樹干莖流量隨樹干傾角的增加而增加,而其與PAI、基徑、株高、冠層厚度和冠層投影面積之間沒有顯著的相關性。

表2 單株泡泡刺的樹干莖流量與冠層結構特征的相關分析

泡泡刺灌叢的集流率在3.38—306.35之間波動,平均值約為(129.66 ± 93.01),表明與空地相比,泡泡刺灌叢的樹干莖流對雨水具有明顯的聚集作用,可將129.66倍的雨水輸送到根際區。由回歸分析可知,集流率與降雨量、降雨強度、降雨歷時之間具有極顯著的相關性(P<0.001),分別呈現指數、對數和線性函數關系(P<0.001)(圖6)。當降雨量<6 mm,降雨強度<0.8 mm/h時,集流率隨雨量、雨強的增加而迅速增加,而后增加速率減緩,該現象可能與泡泡刺灌叢的冠層持水能力和樹干輸水能力有關。

2.4 冠層截留損失與降雨特征和冠層特征的關系

根據水量平衡余項法推算試驗期間共產生冠層截留損失9.54 mm,約占總降雨量的10.50%。從單次降雨事件來看,泡泡刺灌叢冠層截留損失量的變化范圍在0.01—1.53 mm,平均值約為(0.42 ± 0.39)mm。冠層截留損失占次降雨量的比例(簡稱冠層截留損失率)介于1.31%—26.93%,平均值約為(14.21 ± 7.52)%。

降雨特征在冠層截留損失形成的過程中發揮重要作用,回歸分析和曲線擬合的結果表明,泡泡刺灌叢的冠層截留損失量與降雨特征之間呈對數函數關系(P<0.01)(圖7)。當降雨量<4 mm、降雨強度<1 mm/h、降雨歷時<4 h時,冠層截留損失的增加速率較快,隨后冠層截留損失的增加速率放緩,最終趨于一定值(約為1 mm)。該現象可能與泡泡刺灌叢的冠層持水能力有關。冠層截留損失率隨降雨量和降雨歷時的增加呈現線性下降趨勢(P<0.05),與降雨強度之間無顯著相關性,表明短歷時的小降雨事件比長歷時的大降雨事件更易造成冠層截留損失。

圖7 冠層截留損失量、冠層截留損失率與降雨特征之間的回歸關系Fig.7 The correlation between canopy interception loss, canopy interception loss percentage and rainfall regimesIL: 冠層截留損失量 Interception loss; IP: 冠層截留損失率 Interception loss percentage

冠層結構特征對冠層截留損失的影響同樣不可忽視,但單株泡泡刺的冠層截留損失量與冠層結構特征的相關性分析表明,泡泡刺的冠層截留損失量與樹干傾角、PAI、基徑、株高和冠層厚度之間均無顯著的相關性(表3)。

表3 單株泡泡刺的冠層截留損失量與冠層結構特征的相關分析

2.5 冠層降雨再分配與其他氣象因子的關系

為進一步探索其他氣象因子對冠層降雨再分配過程的影響,選取平均風速、平均凈輻射、冠層表面溫度、平均大氣溫度和平均水汽壓等6個要素反映氣象特征,運用Pearson相關分析檢驗各氣象因子與冠層降雨再分配要素的相關性(表4)。結果表明,穿透雨量和冠層截留損失量與各氣象因子之間無顯著相關性；而穿透率和冠層截留損失率與冠層表面溫度、平均大氣溫度及平均水汽壓之間具有顯著的相關性(P<0.05),分別呈負相關和正相關關系,可能是由于當冠層表面溫度、平均大氣溫度、平均水汽壓較高時,冠層表面的蒸發潛力較大,促進了冠層截留雨水向冠層截留損失轉化,進而導致冠層截留損失率的增加和穿透率的減少；樹干莖流量和樹干莖流率與平均大氣濕度呈現顯著正相關關系(P<0.05),可能是由于隨大氣濕度的增加,大氣飽和差減小,在一定程度上抑制了蒸發的產生,促進了冠層截留雨水向樹干莖流轉化。

表4 冠層降雨再分配要素與氣象因子的相關性分析

3 討論

在干旱區,水分是限制植被生長發育的關鍵因素,其數量及時空分布決定著植被的發育過程和生長模式[1, 11]。降雨作為該地區的主要補給水源,植被冠層對于降雨的再分配特征將直接影響冠層下可供植被利用的有效水分[47]。因此,量化泡泡刺灌叢覆蓋下實際進入土壤的有效降雨量及其空間分布規律,對深入理解綠洲-荒漠過渡帶植被對降水再分配過程的調節及其在局地水文循環中的作用具有重要科學意義。

3.1 泡泡刺灌叢的降雨再分配特征

試驗期間泡泡刺灌叢的穿透雨量、樹干莖流量和冠層截留損失量分別為79.89 mm、1.46 mm、9.54 mm,分別占總降雨量87.89%、1.61%、10.50%(具體見圖8)。與阿根廷中部的極叉開拉瑞阿(Larreadivaricata)(78.60%、12.00%、9.40%)[29]、墨西哥南部的硬葉有刺灌叢(Matorral community)(83.30%、8.50%、8.20%)[48]、科爾沁沙地的黃柳(82.78、2.19%、15.03%)[49]及沙坡頭地區油蒿(74.83%、2.89%、22.28%)[2]的降雨再分配特征類似,但相對來說,泡泡刺灌叢的平均穿透率較高,平均樹干莖流率和冠層截留損失率較低,可能與泡泡刺灌叢特殊的植被形態特征有關。泡泡刺灌叢的樹枝較細、樹干表面粗糙且樹枝大多呈現匍匐狀貼近地表,所以其樹干莖流量較少[17,50—51]。同時,含有蠟質層的卵形葉片也不易于截留降雨[18,52],易促進雨水以穿透雨的形式滴落。此外,與其他植被相比,泡泡刺灌叢較為低矮,冠層體積和厚度較小,冠層持水能力較弱,因此,降雨在冠層內的停滯時間短,冠層截留損失和樹干莖流相對較小,大部分雨水以穿透雨的形式到達地面。

圖8 實驗期間穿透雨量、樹干莖流量、冠層截留損失量和降雨強度的變化Fig.8 Variation of throughfall, stemflow, interception loss and rainfall intensity during the experiment

3.2 泡泡刺灌叢覆蓋條件下的有效降雨量

冠層下凈降雨量的多寡及其空間分布決定著冠層下土壤水分的空間分布格局[47],進而對植被的生長發育產生重要影響[29,50—51]。因此,凈降雨量也被稱為冠層下的有效降雨量,反映了植被對雨水的利用效率[53]。實驗期間泡泡刺灌叢共產生凈降雨量81.36 mm,約占總降雨量的89.50%,其中穿透雨約占凈降雨量的98.20%,樹干莖流約占凈降雨量的1.80%。

樹干莖流雖然數量較少,卻能以點源的形式將雨水匯集并輸送至深層土壤中供給植物根部使用[11, 17, 29],通常使用集流率來反映樹干莖流的聚集效應[15]。統計發現,泡泡刺灌叢的平均集流率約為(129.66 ± 93.01),與科爾沁沙地的小葉錦雞兒(100.1)、黃柳(106.2)[50]及皋蘭地區的檸條(153.5)[17]類似,高于毛烏素沙地的籽蒿(41.5)、沙柳(64.9)[54]及沙坡頭地區的油蒿(26)[18],表明泡泡刺灌叢的樹干莖流對根際區土壤水分具有明顯的補給作用,同時,也證明除樹干外,泡泡刺灌叢的葉片及枝條有效地促進了樹干莖流的形成。此外,回歸分析發現,存在一個雨量(約為6 mm)和雨強(約為0.8 mm/h)的閾值,在此閾值之后,泡泡刺灌叢的集流率隨雨量、雨強的增加速率明顯減緩,意味著較大的雨量和雨強可能會降低樹干莖流的集流效率,這與Carlyle-Moses等[55]提出的集流率可能受降雨量閾值的限制相符。

穿透雨是有效降雨量的主要組分,對冠層下表層土壤水分的空間分布產生不可忽視的影響[11]。研究發現穿透雨的空間分布存在顯著的空間異質性,主要表現為多個“雨極”和“旱區”鑲嵌分布的模式,與毛素沙地南緣散生中間錦雞兒[56]和地中海地區松樹[12]的穿透雨空間分布格局類似。該現象在低雨強、短歷時的小降雨事件中表現最為顯著,隨雨量和雨強的增加而減弱,可能是由于高雨強、長歷時的大降雨事件更易促進冠層飽和,加速截留雨水迸濺導致。在降雨特征相同的情況下,PAI的空間異質性成為影響穿透雨空間分布格局的關鍵因素,PAI越大,雨水越易被冠層截留,不易產生穿透雨。此外,穿透雨的空間異質性也意味著我們在實驗中應盡可能多的增加雨量筒的數量以減小實驗誤差。

3.3 泡泡刺灌叢的冠層截留損失

冠層截留損失作為降雨再分配過程中的凈損失項,其大小直接影響冠層下有效降雨量的多少[12]。除了受氣象條件和降雨特征的影響外,冠層持水能力也是影響冠層截留損失的重要因素[57—58],通常將冠層持水能力定義為冠層表面完全濕潤所需要的最小降雨量[59]。實驗中發現泡泡刺灌叢的冠層持水能力可能會影響冠層截留損失隨降雨特征的變化趨勢,在低雨強、短歷時的小降雨事件中,泡泡刺灌叢具有較強的截留吸附能力,截留雨水大多以冠層截留蒸發的形式消耗,此時冠層截留損失的增加速率較快,但隨著雨量、雨強和歷時的增加,當冠層處于飽和或過飽和狀態時,超出冠層持水能力的截留雨水可能會轉化為穿透雨和樹干莖流,進而導致冠層截留損失的增加速率減緩并趨于穩定。Jian等[28, 60—61]在對我國干旱區稀疏植被的冠層截留損失進行觀測時也發現類似現象。根據Wallence等[62]提出的間接方法推算泡泡刺灌叢的冠層持水能力約為0.42 mm,顯著低于黃土高原地區的刺槐(1.34 mm)和油松(1.43 mm)[19]、青海地區的金露梅(0.55 mm)[34]、沙坡頭地區的油蒿(0.75 mm)和檸條(0.71 mm)[60]及葡萄牙中部的橄欖樹(0.97 mm)[33],可能與泡泡刺灌叢低矮的冠層、較小的冠層體積、匍匐狀的枝條及含有蠟質層的葉片有關,這也是導致泡泡刺灌叢冠層截留損失率較低的關鍵原因之一。

與其他稀疏植被相比,泡泡刺灌叢具有較高的穿透率、樹干集流率及較低的冠層截留損失率,更有利于促進大氣降雨量轉化為冠層下的有效雨量,減少冠層截留損失造成的水分凈損失,對干旱區局地水文過程和水文循環產生不可忽視的影響。除氣象因素外,植被形態結構特征也在泡泡刺灌叢降雨再分配的過程中發揮重要作用,如退化的枝條、較小且含有蠟質層的葉片等有效地促進了冠層下水分的補給、提高了對于雨水的利用效率。

4 結論

基于對河西走廊中段綠洲-荒漠過渡帶的泡泡刺灌叢開展長期冠層降雨再分配的定位觀測實驗,通過對其降雨再分配特征及影響因素進行分析,得到以下結論：

(1)實驗期間共觀測到有效降雨事件23次,總降雨量為90.9 mm,泡泡刺灌叢的平均穿透率、樹干莖流率和冠層截留損失率分別為87.89%、1.61%和10.50%。穿透雨和樹干莖流產生的閾值分別為0.2 mm和1.5 mm,冠層持水能力約為0.42 mm。

(2)降雨特征是影響泡泡刺灌叢冠層降雨再分配過程的關鍵氣象要素,回歸分析發現,穿透雨量與降雨量、降雨強度、降雨歷時之間分別呈現線性、冪函數和線性函數關系(P<0.01),穿透率與降雨歷時之間用線性函數的擬合效果較好(P=0.04)；樹干莖流量與降雨量、降雨強度、降雨歷時之間分別呈現線性、指數和線性函數關系(P<0.01),樹干莖流率、冠層截留損失量與降雨特征之間分別呈現對數函數關系(P<0.01)；冠層截留損失率與降雨量、降雨歷時之間用線性函數的擬合效果較好(P<0.05)。其他氣象要素(如平均大氣溫濕度、水汽壓、冠層表面溫度等)也在降雨再分配的過程中發揮不可忽視的作用。

(3)植被形態結構特征是影響冠層降雨再分配過程的另一重要因素,相關分析表明,穿透雨與PAI和株高之間呈現顯著的負相關關系(P<0.0001),PAI的空間異質性可能是影響穿透雨空間分布異質性的主要因素；樹干莖流與樹干傾角之間呈現顯著的正相關關系(P<0.01),隨樹干傾角的增加而增加。與其他干旱區稀疏植被相比,泡泡刺灌叢的穿透率和樹干集流率較大,冠層截留損失率較小,主要與泡泡刺灌叢特殊的植被形態特征有關。