尹準(zhǔn)生， 孫長忠， 趙明揚(yáng)， 王磊

(中國林業(yè)科學(xué)研究院華北林業(yè)實(shí)驗(yàn)中心，102300，北京)

黃土高原半干旱區(qū)油松人工林冠層截留降雨模擬與分析

尹準(zhǔn)生， 孫長忠?， 趙明揚(yáng)， 王磊

(中國林業(yè)科學(xué)研究院華北林業(yè)實(shí)驗(yàn)中心，102300，北京)

基于黃土高原半干旱區(qū)2012—2013年觀測(cè)數(shù)據(jù)，采用定位研究方法對(duì)該區(qū)油松人工林林冠層降水再分配情況進(jìn)行研究；采用王彥輝模型和崔啟武模型對(duì)林冠截留量進(jìn)行模擬并對(duì)其適用性進(jìn)行對(duì)比分析。結(jié)果表明：1)油松人工林穿透雨量與林外降雨量呈顯著線性關(guān)系(R2=0.978 9)，林冠截留量與林外降雨量呈冪指數(shù)關(guān)系(R2=0.867 3)，平均林冠截留率為25.5%；2)王彥輝模型和崔啟武等模型模擬的決定系數(shù)R2分別為0.724 5和0.755 4，黃土高原半干旱區(qū)油松人工林林冠截留模擬效果崔啟武模型較王彥輝模型要好。

油松； 林冠截留模型； 模擬； 林冠截留； 黃土高原

森林水文效應(yīng)是森林生態(tài)系統(tǒng)重要的生態(tài)功能[1]。冠層截留作為森林降雨再分配的一部分，在森林水分循環(huán)過程中占據(jù)重要地位，直接影響著林分涵蓄降水能力的大小、產(chǎn)生地表徑流的能力和林地土壤有效水分的多少[2-4]。影響林冠截留量的因素很多，包括林外降雨量、降雨強(qiáng)度、植被類型和冠層結(jié)構(gòu)等[5-6]。林冠截留模型對(duì)于揭示林冠的截留作用和估算林冠截留量具有重要作用[7]；因此，根據(jù)野外試驗(yàn)樣地的實(shí)測(cè)資料，建立、檢驗(yàn)和優(yōu)化林冠截留模型，并通過模型研究林冠層截留特性仍是森林水文學(xué)的研究熱點(diǎn)[8]。

到目前為止，國內(nèi)外學(xué)者對(duì)林冠截留模型研究做了大量的工作，林冠截留模型主要分為理論模型、經(jīng)驗(yàn)?zāi)Ｐ秃桶虢?jīng)驗(yàn)半理論模型[9]。理論模型經(jīng)過嚴(yán)格的推理而來，劉家岡[10]通過描述林冠對(duì)降雨的截留過程推導(dǎo)出林冠截留降雨的理論模型，該模型機(jī)制性強(qiáng)，且模型中所有參數(shù)可測(cè)，能夠準(zhǔn)確地反映林冠截留各要素的動(dòng)態(tài)過程；但由于理論模型結(jié)構(gòu)復(fù)雜，參數(shù)的測(cè)定較為困難，所以，導(dǎo)致模型在實(shí)際應(yīng)用上受到限制。經(jīng)驗(yàn)?zāi)Ｐ褪峭ㄟ^實(shí)測(cè)數(shù)據(jù)構(gòu)建的水文模型，該類模型形式較多，結(jié)構(gòu)簡單，使用方便，通常只有林外大氣降雨量一個(gè)因素，所以不能反映林冠層截留與各生態(tài)因子之間的關(guān)系。目前，國內(nèi)外對(duì)于林冠截留的研究采用的較多的模型多為半經(jīng)驗(yàn)半理論模型，其中以J.H.C.Gash[11]模型、崔啟武等[12]模型和王彥輝[13]模型居多；但J.H.C. Gash模型需要詳細(xì)的氣象資料來確定其參數(shù)，所以導(dǎo)致模型參數(shù)精度很難保證，從而使模擬結(jié)果出現(xiàn)偏差較大[14]。而王彥輝模型和崔啟武等模型相對(duì)其他模型而言，參數(shù)較少，模型結(jié)構(gòu)簡單，使用次降雨量即可運(yùn)行，且國內(nèi)目前應(yīng)用較多，但針對(duì)這2個(gè)模型適應(yīng)性對(duì)比的報(bào)道相對(duì)較少。因此，筆者通過分析黃土高原半干旱區(qū)油松人工林降雨再分配特征，選取王彥輝模型和崔啟武等模型對(duì)該地區(qū)油松林林冠截留量進(jìn)行模擬對(duì)比，選擇最適合該地區(qū)的林冠截留模型及其參數(shù)，估算林冠截留量，旨在為黃土高原油松人工林林冠層水文效應(yīng)的深入研究提供理論參考。

1 研究區(qū)概況

試驗(yàn)地位于山西省忻州市偏關(guān)縣陳家營，E111°01′～111°22′，N39°12′～39°40′，海拔1 377 m，大陸性氣候，年平均氣溫3.7 ℃，年平均降水量425.3 mm，全年60%以上的降水都集中在7—9月份，主要植被為喬木林和灌叢，造林樹種主要為油松(PinustabulaeformisCarr.)、側(cè)柏(Platycladusorientalis(L.)Franco)和樟子松(PinussylvestrisL.var.mongolicaLitv.)等。研究區(qū)為典型的黃土高原半干旱區(qū)，所觀測(cè)油松純林林齡為30年(經(jīng)撫育的人工同齡林，林相整齊，單株差異較小)。

2 研究方法

2.1 樣地設(shè)置

設(shè)置20 m×20 m油松標(biāo)準(zhǔn)地3個(gè)，平均樹高4.8 m，平均胸徑10.6 cm，郁閉度0.8。

2.2 氣象因子測(cè)量

在林外附近空地設(shè)置小型氣象站一個(gè)進(jìn)行各氣象因子的測(cè)定，在林內(nèi)實(shí)驗(yàn)樣地進(jìn)行林內(nèi)穿透雨和樹干徑流的觀測(cè)。具體測(cè)量方法如下。

1)林外降雨的測(cè)定：林外小型氣象站內(nèi)安置一個(gè)雨量筒，測(cè)量每次降雨量即為林外降雨量。

2)林內(nèi)穿透雨的測(cè)定：在每個(gè)林內(nèi)實(shí)驗(yàn)樣地具有代表性的區(qū)域選擇一矩形(圓形)區(qū)塊，取其對(duì)角線的兩端及中點(diǎn)(圓心及任一直徑的兩端)進(jìn)行雨量筒布設(shè)，布設(shè)3個(gè)，每次降雨后觀測(cè)所得穿透雨量取其平均值即為林內(nèi)穿透雨量。

3)樹干徑流的測(cè)定：每個(gè)樣地選擇標(biāo)準(zhǔn)木3株，將直徑為10 cm的聚乙烯塑料管在樹干胸徑處纏繞2周，固定后將接縫處用密封膠封嚴(yán)，將軟管下部伸入塑料桶，收集導(dǎo)入塑料桶內(nèi)的降水，并使用量筒測(cè)量，通過公式計(jì)算即可求得樹干徑流量。計(jì)算公式為

(1)

式中：SF為樹干徑流量，mm；n為標(biāo)準(zhǔn)木株數(shù)；N為樣地內(nèi)的林木總株數(shù)；S為樣地面積，m2；Ci為第i株標(biāo)準(zhǔn)木的樹干徑流體積，mL。

4)林冠截留量的測(cè)定：林冠截留量可以通過林冠層水量平衡原理計(jì)算得到，即林外降雨量減去林內(nèi)穿透雨量和樹干徑流量計(jì)算得到。其相應(yīng)的計(jì)算公式為

Ic=P-(TF+SF)。

(2)

式中：Ic為林冠截留量，mm；P為林外降雨量，mm；TF為林內(nèi)穿透雨量，mm。

2.3 模型選取

選取王彥輝模型和崔啟武等模型對(duì)黃土高原半干旱區(qū)油松人工林進(jìn)行模擬分析。王彥輝模型和崔啟武等模型均為半經(jīng)驗(yàn)半理論模型，且只有次降雨量1個(gè)因素，在林冠截留模擬的研究中應(yīng)用較廣，故都可作為通用模型運(yùn)用于本研究中。

1)王彥輝模型。是對(duì)Horton模型的簡化，該模型將林冠截留降雨量分解為吸附截留量和附加截留量。模型方程為

(3)

2)崔啟武等模型。通過選取降雨量、相鄰2次降雨的間隔時(shí)間和林冠特性作為影響冠層截留的主要因子建立的林冠截留模型。模型方程為

(4)

式中：A為實(shí)測(cè)林分郁閉度，%；r為通過實(shí)際情況確定的經(jīng)驗(yàn)常數(shù)；I0為飽和截留量，mm，可根據(jù)實(shí)測(cè)資料或者經(jīng)驗(yàn)確定；P′為達(dá)到飽和截留量時(shí)對(duì)應(yīng)的林外降雨量，mm。

3 結(jié)果與分析

3.1 降雨再分配特征

試驗(yàn)期間(2012年4月—2013年9月)共觀測(cè)降雨事件75場(chǎng)，累計(jì)總降雨量為1 106.4 mm，單次降雨量最大值為107.6 mm，單次降雨最小值為1.3 mm，單次平均降雨量為14.8 mm。

將2012—2013年林冠層降雨再分配數(shù)據(jù)匯總，結(jié)果見表1。

根據(jù)國家氣象局規(guī)定的降雨級(jí)的分布，24 h次降雨量小于10 mm的為小雨，10～25 mm的為中雨，25～50 mm為大雨，50～100 mm為暴雨。由表1可知，試驗(yàn)觀測(cè)期間小雨事件所占的頻率為52%，中雨和大雨分別所占的頻率為33.3%和10.7%，暴雨及以上降雨事件較少，頻率僅為4.0%。觀測(cè)期間累計(jì)穿透雨量為857.7 mm，穿透雨率為77.5%，采用75場(chǎng)次降雨數(shù)據(jù)，由對(duì)穿透雨量及林外降雨量進(jìn)行回歸分析(圖1)可知，穿透雨量隨林外降雨量的增大而增大，二者之間的關(guān)系可用線性函數(shù)表示(R2=0.978 9)，由穿透雨率與林外降雨量關(guān)系圖(圖2)可知，穿透雨率隨著林外降雨量的增大呈增大的趨勢(shì)，二者之間呈冪指數(shù)函數(shù)相關(guān)關(guān)系(R2=0.471 9)，當(dāng)林外降雨量小于19.0 mm時(shí)，冠層對(duì)穿透雨的影響較大，但當(dāng)林外降雨量大于19.0 mm時(shí)，隨著林外降雨量的增大，穿透雨率逐漸趨于穩(wěn)定。樹干徑流所占比例小，試驗(yàn)區(qū)內(nèi)當(dāng)降雨量大于6.0 mm時(shí)才有樹干徑流發(fā)生，樹干徑流量所占林冠降雨再分配的比例較小，平均樹干徑流率僅為1.5%。

表1 油松林冠層降雨再分配特征

圖1 林外降雨量與穿透雨量的關(guān)系Fig.1 Relationship between precipitation and throughfall

圖2 林外降雨量與穿透雨率的關(guān)系Fig.2 Relationship between precipitation and percentage of throughfall

試驗(yàn)期間油松林累計(jì)林冠截留量為217.0 mm，占林外降雨量的19.6%，平均截留率為25.5%。由林冠截留與林外降雨量關(guān)系散點(diǎn)圖(圖3)可以看出，林冠截留量隨林外降雨量的增大而增大，當(dāng)林外降雨量增加達(dá)到一定值后林冠截留量增加幅度不明顯，二者之間的關(guān)系可采用冪指數(shù)函數(shù)擬合(R2=0.867 3)。由林冠截留率與林外降雨散點(diǎn)圖(圖4)可以看出，隨著林外降雨的增大，林冠截留率逐漸降低，二者之間呈負(fù)冪指數(shù)關(guān)系(R2=0.588 2)，當(dāng)林外降雨量級(jí)小于10 mm時(shí)，林冠層能夠吸附的雨水所占林外降雨的比例較大，平均林冠截留率為29.9%。當(dāng)林外降雨量級(jí)在10～25 mm時(shí)，林冠截留率變化在11.6%～33.0%之間，平均值為22.6%。隨著林外降雨的增大，林冠截留逐漸趨于飽和，當(dāng)降雨量大于25 mm時(shí)，林冠截留率變化幅度不大。

圖3 林外降雨量與林冠截留雨量的關(guān)系Fig.3 Relationship between precipitation and interception

3.2 林冠截留模擬

選擇崔啟武等模型進(jìn)行模擬之前必須確定經(jīng)驗(yàn)常數(shù)r和飽和截留量I0。根據(jù)本試驗(yàn)中林冠截留量與林外降雨之間的關(guān)系，林冠截留量達(dá)到6.0 mm左右則趨于穩(wěn)定，因此，選擇I0=6.0 mm；經(jīng)驗(yàn)常數(shù)r的確定一般與林冠截留率有關(guān)，根據(jù)相關(guān)文獻(xiàn)[16-18]總結(jié)的二者之間的關(guān)系，同時(shí)結(jié)合本研究的平均林冠截留率25.5%，最后確定模擬模型的r值為0.8。通過計(jì)算，得出林冠截留達(dá)到飽和時(shí)的林外降雨量P′為 16.9 mm。

利用以上確定參數(shù)后的2個(gè)模型，并結(jié)合黃土高原半干旱區(qū)油松人工林2012年4月—2013年9月的75場(chǎng)降雨資料，對(duì)林冠截留進(jìn)行模擬，2個(gè)模型的模擬結(jié)果分別見圖5和圖6。

圖5 王彥輝模型模擬油松林冠層截留量結(jié)果Fig.5 Simulation results of canopy rainfall interception for Pinus tabulaeformis stands using Wang’s model

圖6 崔啟武等模型模擬油松林冠層截留量結(jié)果Fig.6 Simulation results of canopy rainfall interception for Pinus tabulaeformis stands using Cui’s model

由圖5王彥輝模型擬合結(jié)果可知，當(dāng)截留量小于4.0 mm時(shí)，模擬值大多都大于真實(shí)值，截留量大于4.0 mm時(shí)模擬效果相對(duì)較好，模擬值與實(shí)測(cè)值基本分布在1∶1線附近。

從圖6也可以看出，當(dāng)林冠層達(dá)到飽和之前，崔啟武等模型具有較高的精度：當(dāng)林冠截留量小于6.0 mm時(shí)，模擬值與實(shí)測(cè)值基本分布在1∶1線附近；當(dāng)林冠截留大于6.0 mm時(shí)，由于林冠截留被視為一常數(shù)，但由于仍有蒸發(fā)存在，林冠截留量實(shí)際上仍處于增加狀態(tài)，所以，此時(shí)崔啟武等模型的模擬值均小于實(shí)測(cè)值。

對(duì)于所觀測(cè)到的降雨，降雨量與時(shí)間呈現(xiàn)一定的正相關(guān)趨勢(shì)，這樣就會(huì)使林冠截留在達(dá)到飽和之后，仍有降雨在發(fā)生，且降雨量越大，時(shí)間也越長，飽和之后的總蒸發(fā)量也就越大。這就造成隨著截留量的增加，偏差值也呈增大的趨勢(shì)，這與圖6顯示的結(jié)果非常一致。從總的模擬結(jié)果來看：王彥輝模型和崔啟武等模型對(duì)林冠截留模擬的總誤差分別是-17.8和11.34 mm，決定系數(shù)R2分別為0.724 5和0.755 4；王彥輝模型模擬值的相對(duì)誤差值范圍為-7.6～5.9 mm，崔啟武等模型模擬值的相對(duì)誤差值范圍為-2.13～8.05。2個(gè)模型模擬的平均絕對(duì)誤差分別為1.24和1.12 mm，但對(duì)于林冠截留量8.0 mm以上的降雨，2個(gè)模型模擬效果均一般。

王彥輝模型將林冠截留中的附加截留項(xiàng)簡化為林冠截留蒸發(fā)率和林外降雨量的乘積，并把林冠截留蒸發(fā)率假定為一個(gè)常數(shù)，這導(dǎo)致截留量小時(shí)模擬值大于實(shí)測(cè)值。崔啟武等模型則忽略了降雨期間的林冠截留蒸發(fā)，認(rèn)為林冠截留達(dá)到飽和后截留量數(shù)值不再變化，導(dǎo)致模擬達(dá)到飽和截留降雨后模擬值較實(shí)測(cè)值偏小?？傮w而言，王彥輝模型和崔啟武等模型對(duì)于黃土高原半干旱區(qū)油松人工林的模擬還是成功的，模擬精度較高。

4 討論

對(duì)2012年4月—2013年9月黃土高原半干旱區(qū)油松林冠截留進(jìn)行分析得出，平均林冠截留率為25.5%，符合國內(nèi)外學(xué)者的研究結(jié)果(平均林冠截留率變動(dòng)范圍為10%～45%[19])，略高于中國各森林生態(tài)系統(tǒng)林冠的平均截留率[20]。影響林冠截留的因素較多，主要有降雨量、降雨強(qiáng)度、空氣溫度、空氣濕度、風(fēng)速和植被特征等[21]。王彥輝模型將林冠截留降雨量分解為吸附截留量和附加截留量2部分，其中的附加截留量Ia項(xiàng)假定為降雨歷時(shí)和林冠投影面積內(nèi)的平均蒸發(fā)強(qiáng)度的乘積?？紤]到模型不宜太復(fù)雜又將附加截留項(xiàng)簡化為林冠截留蒸發(fā)率α和林外降雨量P的乘積(即Ia=αP)，并且將蒸發(fā)率α假定為一個(gè)常數(shù)；而實(shí)際降雨過程中蒸發(fā)率α隨著降雨量的增加而減小，所以必然導(dǎo)致模擬精度受到影響。崔啟武等模型考慮了降水量，植被林冠特征以及林冠干濕程度等因子，但卻忽略了林冠附加截留損失這一部分。周國逸[22]認(rèn)為，當(dāng)降雨量較大，并使林冠截留達(dá)到飽和時(shí)，林冠層的截留潛力由于水分蒸發(fā)的作用得以恢復(fù)，林冠層仍然可以繼續(xù)截留少量降雨，相應(yīng)的林冠層實(shí)際截留量將大于飽和截留量；因此，由于考慮一些實(shí)際情況，使得王彥輝模型和崔啟武等模型在模擬林冠截留量時(shí)均存在不足之處。

本研究試驗(yàn)期間，2012和2013年均為豐水年，且2012年大雨以上次數(shù)較多，而恰巧2個(gè)模型對(duì)于降雨量較大時(shí)模擬截留量均存在缺陷，故較大降雨時(shí)導(dǎo)致模擬值與實(shí)測(cè)值差別較大，模擬截留總量相應(yīng)也會(huì)產(chǎn)生誤差。這也是本研究中2個(gè)模型模擬誤差相對(duì)偏大的原因之一。對(duì)黃土高原半干旱區(qū)油松人工林，采用王彥輝模型和崔啟武等模型進(jìn)行對(duì)比發(fā)現(xiàn)，二者的決定系數(shù)R2分別為0.724 5和0.755 4，且誤差值崔啟武等模型較王彥輝模型稍好。由此可知，在該試驗(yàn)區(qū)域，考慮了降雨量和林冠特征2個(gè)因子的崔啟武等模型要略好于只考慮降雨量1個(gè)因子的王彥輝模型。崔啟武等模型的決定系數(shù)為0.755 4，說明降雨量和林冠特征只能解釋該區(qū)域油松林冠截留量的75.54%的信息，剩余的24.46%的信息需要由降雨強(qiáng)度，降雨期間的空氣溫度、空氣濕度、風(fēng)速和氣壓等因素來解釋。由于該模型的主要缺陷為林冠截留量達(dá)到飽和后忽略了林冠附加截留損失這一項(xiàng)，與現(xiàn)實(shí)情況不符；故日后的研究中可以著重對(duì)該模型這塊進(jìn)行優(yōu)化，以建立更適用、模擬更精確的林冠截留模型，為林冠層水文效應(yīng)的深入研究提供理論參考。

5 結(jié)論

1)黃土高原半干旱區(qū)油松人工林穿透雨量與林外降雨量呈顯著線性關(guān)系(R2=0.978 9)，林冠截留量與林外降雨量之間呈冪指數(shù)關(guān)系(R2=0.867 3)，平均林冠截留率為25.5%，林冠截留率隨著林外降雨的增大逐漸減少，隨著林外降雨量增加到一定值后林冠截留率值變化不大。

2)采用王彥輝模型和崔啟武等模型對(duì)油松林林冠層截留降雨量進(jìn)行模擬，決定系數(shù)R2分別為0.724 5和0.755 4，模擬效果相對(duì)較好；對(duì)比2個(gè)模型模擬結(jié)果可知，黃土高原半干旱區(qū)油松人工林林冠截留模擬效果崔啟武等模型較王彥輝模型要好。

[1] 劉世榮.中國森林生態(tài)系統(tǒng)水文生態(tài)功能規(guī)律[M].北京:中國林業(yè)出版社,1996:3-6

[2] 趙明揚(yáng),孫長忠,康磊.降雨再分配的回歸模型差異性分析[J].中南林業(yè)科技大學(xué)學(xué)報(bào),2013,33(5):79-84

[3] 王冉,王玉杰,王云琦,等.三峽庫區(qū)典型常綠闊葉林冠層降雨截留模擬研究[J].長江流域資源與環(huán)境,2014,23(2):281-286

[4] Oishi A C, Oren R, Novick K A, et al. Interannual invariability of forest evapotranspiration and its consequence to water flow downstream[J].Ecosystems,2010,13(3):421-436

[5] Dietz J, Holscher D, Leuschner C, et al. Rainfall partitioning in relation to forest structure in differently managed montane forest stands in Central Sulawesi, Indonesia[J].Forest Ecology and Management,2006,237(1):170-178

[6] Muzylo A, Llorens P, Valente F, et al. A review of rainfall interception modelling[J].Journal of Hydrology, 2009,370(1):191-206

[7] 彭煥華,趙傳燕,沈衛(wèi)華,等.祁連山北坡青海云杉林冠對(duì)降雨截留空間模擬[J].干旱區(qū)地理,2010,33(4):600-606

[8] Komatsu H ,Shinohara Y, Kume T, et a1. Relationship between annual rainfall and interception ratio for forests across Japan[J].Forest Ecology and Management,2008,256(5):1189-1197

[9] Liu Shuguang. A new model for the prediction of rain fall interception in forest canopies[J].Ecological Modelling,1997,99(9):151-159

[10] 劉家岡.林冠對(duì)降雨的截留過程[J].北京林業(yè)大學(xué)學(xué)報(bào),1987,9(2):140-144

[11] Gash J H C. An analytical model of rainfall interception by forests[J]. Quarterly Journal of Royal Meteorological Society,1979,105:43-55

[12] 崔啟武,邊履剛,史繼德,等.林冠對(duì)降雨的截留作用[J].林業(yè)科學(xué),1980,16(2):141-146

[13] 王彥輝.刺槐對(duì)降雨的截持作用[J].生態(tài)學(xué)報(bào),1987,7(1):43-49

[14] 何常清,薛建輝,吳永波,等.應(yīng)用修正的Gash解析模型對(duì)岷江上游亞高山川滇高山櫟林林冠截留的模擬[J].生態(tài)學(xué)報(bào),2010,30(5):1125-1132

[15] 王彥輝,于澎濤,徐德應(yīng),等.林冠截留降雨模型轉(zhuǎn)化和參數(shù)規(guī)律的初步研究[J].北京林業(yè)大學(xué)學(xué)報(bào),1998,20(6):25-30

[16] 衛(wèi)三平,王力,吳發(fā)啟.黃土丘陵溝壑區(qū)刺槐林冠截留模擬[J].林業(yè)科學(xué),2008,44(1):26-33

[17] 佘冬立,劉營營,邵明安,等.黃土坡面不同植被冠層降雨截留模型模擬效果及適用性評(píng)價(jià)[J].農(nóng)業(yè)工程學(xué)報(bào),2012,28(16):115-120

[18] 王佑民.中國林冠降雨再分配研究綜述[J].西北林學(xué)院學(xué)報(bào),2000,15(3):1-7

[19] 曾偉,熊彩云,肖復(fù)明,等.不同密度退耕雷竹春季林冠截留特性[J].生態(tài)學(xué)雜志,2014,33(5):1178-1182

[20] 溫遠(yuǎn)光,劉世榮.我國主要森林生態(tài)系統(tǒng)類型降水截留規(guī)律的數(shù)量分析[J].林業(yè)科學(xué),1995,31(4): 289-297

[21] 鮮靖蘋,張家洋,胡海波.森林冠層水文研究進(jìn)展[J].西北林學(xué)院學(xué)報(bào),2014,29(3):96-104

[22] 周國逸.生態(tài)系統(tǒng)水熱原理及其應(yīng)用[M].北京:氣象出版社,1997:72-73

(責(zé)任編輯：宋如華)

Simulation and analysis of canopy interception onPinustabulaeformisplantations in semi-arid areas of the Loess Plateau

Yin Zhunsheng, Sun Changzhong, Zhao Mingyang, Wang Lei

(Forestry Experiment Center of North China, Chinese Academy of Forestry, 102300, Beijing, China)

We investigated the rainfall redistribution of canopy onPinustabulaeformisplantations in the semi-arid areas of the Loess Plateau based on the data collected from 2012 to 2013 using the permanent plot method. The canopy interception was simulated separately by using Wang′s and Cui′s interception models, and the applicability of the two models were compared. The results showed that: 1) There was a linear relationship between rainfall and throughfall (R2=0.978 9), but an exponential relationship between rainfall and canopy interception (R2=0.867 3), and the average of interception rate was 25.5%. 2) The determination coefficient of Wang’s interception model was 0.724 5 and it was 0.755 4 for Cui′s. In general, Cui′s interception model was slightly better than Wang′s.

Pinustabulaeformis; canopy interception model; simulation; canopy interception; Loess Plateau

2014-09-16

2015-04-06

項(xiàng)目名稱： 國家自然科學(xué)基金“黃土高原半干旱區(qū)人工植被水分生態(tài)環(huán)境效應(yīng)規(guī)律研究”(30170769)

尹準(zhǔn)生(1988—)，男，碩士研究生。主要研究方向：森林水文。E-mail：yinzhunsheng@163.com

?通信作者簡介： 孫長忠(1957—)，男，研究員。主要研究方向：森林培育與森林生態(tài)。E-mail： sun61@163.com

S715

A

1672-3007(2015)03-0045-06