周 威，康永祥，劉建軍，夏國(guó)威，劉小林

（1.西北農(nóng)林科技大學(xué) 林學(xué)院，陜西 楊凌712100；2.賀蘭山森林生態(tài)系統(tǒng)定位研究站，寧夏 銀川750021；3.甘肅小隴山森林生態(tài)系統(tǒng)定位研究站，甘肅 天水741022）

水在森林生態(tài)系統(tǒng)中的循環(huán)與分配整合了能量流動(dòng)和養(yǎng)分循環(huán)等生態(tài)過(guò)程［1］，因而水文功能在森林生態(tài)系統(tǒng)的作用中是非常重要的。張焜［2］、王波等［3］、孔維健等［4］、劉小林等［5］分別在不同地域?qū)μ烊会橀熁旖涣值乃墓δ苓M(jìn)行了研究，為探求天然針闊混交林的水文作用提供了有價(jià)值的參考。但由于研究對(duì)象的林分特征、氣候條件、研究側(cè)重點(diǎn)等存在差異，要全面認(rèn)識(shí)天然針闊混交林的水文效應(yīng)還任重而道遠(yuǎn)。

森林引起生態(tài)系統(tǒng)內(nèi)降雨分配的時(shí)空差異，對(duì)林地水文過(guò)程產(chǎn)生顯著影響［6］。國(guó)內(nèi)外針對(duì)林冠層分配降水研究較多，但在秦嶺西段小隴山的研究鮮見(jiàn)報(bào)道［7］。小隴山林區(qū)是兼有中國(guó)南北氣候特點(diǎn)的典型天然次生林區(qū)，天然次生針闊混交林是一個(gè)重要森林類(lèi)型。研究小隴山林區(qū)森林冠層下穿透降雨規(guī)律及與大氣降雨之間的相關(guān)性，對(duì)揭示區(qū)域森林生態(tài)水文作用具有重要意義。

在森林生態(tài)系統(tǒng)中，不同冠層對(duì)降雨的截留效應(yīng)響應(yīng)不同，導(dǎo)致了林內(nèi)穿透降雨的空間分布異質(zhì)性［8］。地統(tǒng)計(jì)學(xué)是一種空間分析方法，用于描述區(qū)域化變量的空間分布特征，并提供了一種無(wú)偏最優(yōu)的空間插值方法，可用來(lái)充分認(rèn)識(shí)空間結(jié)構(gòu)的特征，已廣泛應(yīng)用于礦產(chǎn)、土壤、生態(tài)、氣象等因子分布圖的繪制［9］。本文試圖將地統(tǒng)計(jì)學(xué)的理論方法應(yīng)用到穿透降雨的研究中，探討穿透降雨空間分布的特征。

本文通過(guò)對(duì)小隴山天然次生針闊混交林下穿透降雨特征及其空間分布格局與降雨量關(guān)系的分析研究，旨在進(jìn)一步明確針闊混交天然林生態(tài)系統(tǒng)的水分運(yùn)轉(zhuǎn)過(guò)程及功能的生態(tài)學(xué)機(jī)制，以求更深入地探究森林水文的過(guò)程及作用機(jī)理，以期為正確認(rèn)識(shí)小隴山林區(qū)的森林水文效應(yīng)及流域森林生態(tài)水文功能提供一定的觀測(cè)數(shù)據(jù)和經(jīng)驗(yàn)參數(shù)，對(duì)制定森林管理措施提供科學(xué)依據(jù)。

1 研究區(qū)概況

研究區(qū)位于甘肅省小隴山林業(yè)科學(xué)研究所沙壩實(shí)驗(yàn)基地，地處小隴山林區(qū)中心地帶，行政區(qū)劃屬天水市秦州區(qū)娘娘壩鎮(zhèn)，東接觀音林場(chǎng)，南接高橋林場(chǎng)，西鄰麻沿林場(chǎng)，北連李子園林場(chǎng)。屬秦嶺南坡石質(zhì)山地，海拔在1 550～2 100m，年平均氣溫1.2～8.8℃，最高氣溫30.3℃，最低氣溫－22.4℃，≥10℃的年積溫2 480.4℃，初霜期10月16日，終霜期5月4日，無(wú)霜期154d，年均降雨量757mm，平均年蒸發(fā)量1 012.2mm，平均相對(duì)濕度78%，屬大陸性季風(fēng)氣候，為暖溫帶濕潤(rùn)區(qū)。土壤以山地褐土和山地棕壤為主，成土母質(zhì)為巖石，土層厚度30—55cm，土壤質(zhì)地為砂壤，表土層腐殖質(zhì)含量豐富，礦物質(zhì)養(yǎng)分含量中等，呈微酸性到中性。

樣地設(shè)在沙壩實(shí)驗(yàn)基地亂石窖溝，地理坐標(biāo)為105°53′72″—105°53′73″E，34°8′56″—34°8′57″N，海拔1 630～1 660m，樣地坡度15°～25°，南坡中坡位。樣地為以銳齒櫟（Quercus aliena var.acuteserrata）為主的天然次生針闊混交林，幾乎無(wú)人為干擾，植被繁密。樣地內(nèi)胸徑大于5cm的喬木總計(jì)120株，平均胸徑16.8cm，平均樹(shù)高19.9m，平均郁閉度0.93，喬木林群落的樹(shù)種組成為4銳齒櫟，2華山松（Pinus armandii），1四照花（Cronus japonica var.chinen－sis），1燈臺(tái)樹(shù)（Bothrocaryum controversum），2軟闊葉樹(shù)種。樣地內(nèi)下木主要有白檀（Symplocos paniculata）、千金榆（Carpinus cordata）、鵝耳櫪（Carpinus turczaninowii）、葛棗獼猴桃（Actinidia polygama）、灰栒子（Cotoneaster acutifolius）等，平均高度2.5m，蓋度為60%。草本主要有菝葜（SmilaxL.）、羊胡子草（EriophorumL.）、懸鉤子（Rubus L.）等，蓋度15%。

2 研究方法

2.1 試驗(yàn)布設(shè)與數(shù)據(jù)收集

林外大氣降雨量資料來(lái)源于小隴山森林生態(tài)系統(tǒng)定位研究站，位于甘肅省小隴山林業(yè)實(shí)驗(yàn)局林業(yè)科學(xué)研究所沙壩實(shí)驗(yàn)基地，地理坐標(biāo)為105°54′E，34°07′N(xiāo)。在沙壩實(shí)驗(yàn)基地天然次生林中選取一塊投影28.5m×27m的樣地，借助測(cè)繩和羅盤(pán)儀沿等高線(xiàn)精確圈出樣地。在樣地內(nèi)網(wǎng)格狀機(jī)械布設(shè)100個(gè)雨量桶（直徑20cm），保持雨量筒高度距地面60cm，按品字形水平布設(shè)，行與行、列與列之間投影距離均為3m，用于收集樣地的穿透降雨量。每次降雨后及時(shí)測(cè)定各個(gè)雨量桶內(nèi)的水量，計(jì)算雨量時(shí)根據(jù)收集面積將雨量桶內(nèi)每次降雨量換算為mm單位。

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

采用Microsoft Excel 2007和SPSS 18.0對(duì)數(shù)據(jù)進(jìn)行一般處理分析，利用穿透降雨的測(cè)定數(shù)據(jù)，計(jì)算得到每次降雨事件每個(gè)雨量筒的穿透率（相同時(shí)間下林內(nèi)穿透降雨量與林外大氣降雨量的比值）。完成基本統(tǒng)計(jì)分析后，用SPSS 18.0對(duì)數(shù)據(jù)通過(guò)直方圖和正態(tài)QQ圖進(jìn)行正態(tài)分布檢驗(yàn)，確定各穿透率數(shù)據(jù)均符合正態(tài)分布，不需要進(jìn)行數(shù)據(jù)轉(zhuǎn)換。

然后利用地統(tǒng)計(jì)軟件GS＋9.0對(duì)穿透率數(shù)據(jù)進(jìn)行半方差函數(shù)分析和參數(shù)計(jì)算，所謂半方差函數(shù)，又叫半變異函數(shù)，就是兩點(diǎn)間插值的方差的1／2，其表達(dá)式為：

式中：γ（h）——間距h的半方差函數(shù)值；N（h）——以h為間距的穿透雨觀測(cè)點(diǎn)對(duì)的數(shù)目；Z（xi）——樣點(diǎn)在xi處的穿透率；Z（xi＋h）——與xi相隔間距h處樣點(diǎn)的穿透率。下同。半變異函數(shù)值隨樣點(diǎn)間距的增加而增大，并在一定的間距（稱(chēng)為變程，range）達(dá)到一個(gè)基本穩(wěn)定的常數(shù)（稱(chēng)為基臺(tái)，sill）［10］。

采用GS＋9.0中不同類(lèi)型半變異模型進(jìn)行擬合，因相比于圓形模型、指數(shù)模型、高斯模型和線(xiàn)性模型而言，球面模型擬合效果最好，故本文所用半變異函數(shù)模型均為球面模型［11］。

式中：C0——塊金，即間距為0時(shí)的半方差函數(shù)值，是由實(shí)驗(yàn)誤差和隨機(jī)因素引起的變異；C——結(jié)構(gòu)方差，指空間自相關(guān)部分引起的變異；（C0＋C）——基臺(tái)，表示系統(tǒng)內(nèi)的總變異，一般情況下，基臺(tái)值越大表示總的空間異質(zhì)性程度越高，反之越??；a——變程，指空間自相關(guān)距，表示觀測(cè)值之間的距離大于該值時(shí)是相互獨(dú)立的，小于該值時(shí)存在空間自相關(guān)性［12］。

采用ArcGIS 10.0這一GIS軟件平臺(tái)，利用地統(tǒng)計(jì)學(xué)中常用的數(shù)據(jù)插值空間變異分析方法——Kriging插值法，進(jìn)行穿透降雨空間分布的預(yù)測(cè)模擬分析。

Kriging法是建立在半變異函數(shù)理論和結(jié)構(gòu)分析的基礎(chǔ)上，在有限的區(qū)域內(nèi)對(duì)區(qū)域化變量的取值進(jìn)行無(wú)偏最優(yōu)估計(jì)的一種方法，可對(duì)周?chē)膶?shí)測(cè)值進(jìn)行加權(quán)以得出未測(cè)位置的預(yù)測(cè)，公式由數(shù)據(jù)的加權(quán)總和組成：

式中：Z（si）——第i個(gè)位置處的測(cè)量值；λi——第i個(gè)位置處的測(cè)量值的未知權(quán)重；s0——預(yù)測(cè)位置；N——測(cè)量值數(shù)量。權(quán)重λi取決于測(cè)量點(diǎn)、預(yù)測(cè)位置的距離和預(yù)測(cè)位置周?chē)臏y(cè)量值之間空間關(guān)系的擬合模型。

使用ArcGIS 10.0的地統(tǒng)計(jì)分析模塊中普通Kriging法的球面模型對(duì)穿透率數(shù)據(jù)進(jìn)行插值，預(yù)測(cè)出樣地中300多萬(wàn)個(gè)點(diǎn)的穿透率值，繪制成柵格圖像數(shù)據(jù)，直觀地模擬樣地的降雨分布。

通過(guò)計(jì)算得到各降雨事件的穿透率平均值、標(biāo)準(zhǔn)差（S）和變異系數(shù)（Cv），作為評(píng)價(jià)不同降雨量時(shí)穿透率的平均狀況和總體不均勻度的指標(biāo)。

3 結(jié)果與分析

3.1 各降雨事件穿透降雨特征

在本實(shí)驗(yàn)觀測(cè)期2011年8月至2012年7月共收集了29次降雨事件數(shù)據(jù)，總降雨量為947.3mm，總的穿透雨量和穿透率分別為742.6mm和0.783 9。通過(guò)對(duì)這29次降雨事件實(shí)測(cè)數(shù)據(jù)的分析，發(fā)現(xiàn)當(dāng)大氣降雨量小于1.9mm時(shí)，林冠截留了幾乎全部降水，沒(méi)有穿透雨產(chǎn)生，當(dāng)大氣降雨量大于1.9mm時(shí)，樣地中100個(gè)樣點(diǎn)平均穿透降雨量與大氣降雨量緊密相關(guān)，穿透降雨量隨大氣降雨量的增加而增大，呈顯著的線(xiàn)性正相關(guān)關(guān)系（圖1a）。擬合的方程如下：

式中：T——穿透降雨量（mm）；P——大氣降雨量（mm）。

隨著降雨量的增加，穿透率也呈增加的趨勢(shì)。對(duì)穿透率與降雨量的關(guān)系進(jìn)行回歸分析，比較后得出對(duì)數(shù)函數(shù)具有較好的擬合性（圖1b），其關(guān)系式為：

式中：Tr——穿透率；P——大氣降雨量（mm）。

圖1 穿透雨量、穿透率及穿透雨量變異系數(shù)與大氣降雨量的關(guān)系

可見(jiàn)，穿透率隨降雨量的變化基本分為3個(gè)階段：在小雨雨量級(jí)（單場(chǎng)降雨量10mm以下），穿透率隨著降雨量的增加而迅速增大，為快變期；在中雨雨量級(jí)（單場(chǎng)降雨量10～25mm），穿透率增加速率明顯變緩，為漸變期；在大雨雨量級(jí)（單場(chǎng)降雨量25～60mm）和暴雨雨量級(jí)（單場(chǎng)降雨量60mm以上），穿透率逐漸趨于平穩(wěn)，波動(dòng)范圍在0.793 9～0.907 5，為穩(wěn)定期。表1的穿透率平均數(shù)也可以印證這一結(jié)論。

穿透降雨變異系數(shù)（降雨事件的穿透降雨標(biāo)準(zhǔn)差與穿透降雨量的比值）揭示穿透降雨的空間分布不均勻程度，其隨著降雨量的增加而呈下降的趨勢(shì)。對(duì)穿透降雨變異系數(shù)與降雨量的關(guān)系進(jìn)行回歸分析，經(jīng)過(guò)比較得出Mitscherlich模型具有較好的擬合性（圖1c），Mitscherlich模型是Richards模型的一個(gè)特例，最初是描述生物生長(zhǎng)量的模型。其關(guān)系式為：

式中：Tv——穿透降雨變異系數(shù)；P——大氣降雨量（mm）。

穿透降雨變異系數(shù)隨降雨量的變化也可分為快變期、漸變期和穩(wěn)定期3個(gè)階段，但變化規(guī)律與穿透率隨降雨量的變化相反。在小雨雨量級(jí)內(nèi)，穿透降雨變異系數(shù)隨著降雨量的增加而迅速減小，在中雨雨量級(jí)之后逐漸趨于穩(wěn)定，波動(dòng)范圍在0.193 0～0.193 2。

表1 不同雨量級(jí)穿透率的不同指標(biāo)特征

3.2 不同雨量級(jí)下穿透降雨分布格局

穿透率觀測(cè)值的半方差函數(shù)分析（表2）表明：基臺(tái)值（C0＋C）隨雨量級(jí)的增大而減小，說(shuō)明降雨量增大使得降雨穿透率在空間上趨向于均勻化；各穿透率半方差函數(shù)模型空間自相關(guān)范圍（變程a）大致介于5.5～6.2m；模型殘差值RSS是選擇模型的主要依據(jù)，所有擬合的半變異函數(shù)模型殘差值均極小，RSS越小則擬合效果越好，說(shuō)明各模型擬合效果都很好；決定系數(shù)R2表示穿透率半方差模型的解釋效率，即模型擬合精確度的數(shù)字表示，可見(jiàn)暴雨的模擬精度最高，小雨的模擬精度最低，模擬精度隨雨量級(jí)的增大而升高，這可能與穿透降雨不確定度隨降雨量的增加而減小有關(guān)。

表2 不同雨量級(jí)下穿透率半變異函數(shù)理論模型及其參數(shù)

表2中總體指各樣點(diǎn)29次穿透降雨率平均值的集合，小雨、中雨、大雨、暴雨分別指各樣點(diǎn)在不同雨量級(jí)下的穿透降雨率平均值的集合。

塊金與基臺(tái)的比值C0／（C0＋C）表示實(shí)驗(yàn)誤差和隨機(jī)因素引起的空間變異占系統(tǒng)總變異的比例。從表2可見(jiàn)，每組穿透率數(shù)據(jù)C0／（C0＋C）值都很小，說(shuō)明實(shí)驗(yàn)誤差和隨機(jī)因素引起的空間變異很小，主要是空間自相關(guān)部分引起的變異。但C0／（C0＋C）值隨雨量級(jí)的增大而增大，說(shuō)明隨著降雨量的增加，隨機(jī)因素引起的空間異質(zhì)性逐漸增加，由結(jié)構(gòu)性的空間因素引起的空間異質(zhì)性逐漸下降，即與樹(shù)冠結(jié)構(gòu)的相關(guān)性逐漸降低。

從空間自相關(guān)的角度來(lái)看，C0／（C0＋C）表示變量的空間相關(guān)性程度，如果C0／（C0＋C）＜25%，說(shuō)明系統(tǒng)具有強(qiáng)烈的空間相關(guān)性，如果25%≤C0／（C0＋C）≤75%，說(shuō)明系統(tǒng)具有中等的空間相關(guān)性，如果C0／（C0＋C）＞75%，則說(shuō)明系統(tǒng)空間相關(guān)性很弱［13］。而表2中每組數(shù)據(jù)C0／（C0＋C）都顯著小于25%，因此均表現(xiàn)出強(qiáng)烈的空間相關(guān)性。其中，小雨的空間自相關(guān)最顯著。

用Kriging插值法繪制穿透降雨率在各雨量級(jí)下的空間分布圖（圖2），由圖2中可直觀地看出穿透降雨在不同雨量級(jí)的空間分布格局及其差異，穿透率的空間分布受雨量級(jí)的影響呈現(xiàn)出不同的特點(diǎn)。在小雨條件下，穿透率普遍偏小，其分布東北部及個(gè)別區(qū)域穿透率相對(duì)略高，其他各處較為平均，穿透率處于0.75以下（0.6～0.75）的占到樣地總面積的80%以上。中雨條件下的穿透率略有增大，穿透率在整個(gè)樣地中的分布比較平滑勻稱(chēng)，其分布更具有代表性。大雨條件下，穿透率整體較高，但從表2來(lái)看，其分布的不均勻性并不很強(qiáng)，變化相對(duì)比較平緩，85%以上的區(qū)域穿透率在0.8～0.95。相對(duì)來(lái)說(shuō)，在暴雨條件下的穿透降雨分布呈現(xiàn)出一些差異，其分布的不均勻性更明顯更突出。但在不同雨量級(jí)下，穿透降雨率分布形態(tài)也有不同程度的相似，如東北部均存在一個(gè)相對(duì)高穿透率區(qū)域，中部偏西及西南角的同一位置均存在相對(duì)低穿透率區(qū)域。因?yàn)槌涤炅看笮⊥?，林冠下穿透降雨的空間異質(zhì)性同時(shí)受多個(gè)其他因素綜合影響，包括葉面積指數(shù)、冠層厚度、離樹(shù)干距離、冠層郁閉度、冠形、雨強(qiáng)、風(fēng)速、風(fēng)向、坡度等。

圖2 基于100個(gè)樣點(diǎn)在不同雨量級(jí)下的穿透降雨率平均值預(yù)測(cè)的穿透降雨率在樣地內(nèi)分布

3.3 穿透率實(shí)測(cè)數(shù)據(jù)與預(yù)測(cè)數(shù)據(jù)的比較

由表1可知，實(shí)測(cè)與Kriging預(yù)測(cè)的穿透率的平均數(shù)差別很小，可以看成是一致的，但是實(shí)測(cè)數(shù)據(jù)與預(yù)測(cè)數(shù)據(jù)的標(biāo)準(zhǔn)差和變異系數(shù)差別很大，預(yù)測(cè)數(shù)據(jù)的標(biāo)準(zhǔn)差和變異系數(shù)明顯小于實(shí)測(cè)數(shù)據(jù)。且根據(jù)表1實(shí)測(cè)數(shù)據(jù)可知，雨量級(jí)從小雨增大時(shí)，穿透降雨率的變異系數(shù)減小。由此可知，穿透降雨分布的不均勻性隨雨量級(jí)的增加而下降。但預(yù)測(cè)數(shù)據(jù)不符合上述規(guī)律，在預(yù)測(cè)數(shù)據(jù)中，大雨的標(biāo)準(zhǔn)差和變異系數(shù)最小，其次是小雨以及中雨，暴雨的標(biāo)準(zhǔn)差和變異系數(shù)最大，且其遠(yuǎn)高于小雨、中雨和大雨，小雨、中雨和大雨的標(biāo)準(zhǔn)差和變異系數(shù)差別不大。由此推斷，除大雨外，穿透降雨分布的不均勻性隨雨量級(jí)的增加而增大。實(shí)測(cè)數(shù)據(jù)與預(yù)測(cè)數(shù)據(jù)產(chǎn)生矛盾。

雖然圖3中預(yù)測(cè)值和實(shí)測(cè)值頻數(shù)分布的基本趨勢(shì)是一致的，但預(yù)測(cè)值的最大值比實(shí)測(cè)最大值小，而預(yù)測(cè)值的最小值比實(shí)測(cè)最小值大，這意味著預(yù)測(cè)值比實(shí)測(cè)值的極差縮小。因此，實(shí)測(cè)數(shù)據(jù)頻數(shù)分布圖相對(duì)矮寬，預(yù)測(cè)數(shù)據(jù)的頻數(shù)分布圖相對(duì)高窄，Kriging預(yù)測(cè)數(shù)據(jù)的離散度下降，分布更為集中，于是預(yù)測(cè)值的分布不能完全再現(xiàn)樣本分布，即存在“平滑效應(yīng)”。李超等［14］對(duì)平滑效應(yīng)的產(chǎn)生原因做了詳盡闡釋?zhuān)交?yīng)的存在會(huì)導(dǎo)致預(yù)測(cè)值的空間連續(xù)程度增強(qiáng)，變異程度下降。

盡管Kriging預(yù)測(cè)值存在條件偏差，但由于在Kriging方程組中對(duì)無(wú)偏性進(jìn)行了強(qiáng)制約束，使其仍然是全局最優(yōu)的。對(duì)于平滑效應(yīng)的認(rèn)識(shí)也要全面、客觀，要針對(duì)具體的研究目的而定。

雖然平滑效應(yīng)導(dǎo)致了以上的問(wèn)題，但正是Kriging預(yù)測(cè)值存在平滑效應(yīng)才使得Kriging法能夠用來(lái)繪制穿透率的等值線(xiàn)和趨勢(shì)面，對(duì)預(yù)測(cè)結(jié)果的離散程度沒(méi)有任何約束的插值方法是無(wú)法繪制等值線(xiàn)或等值面的。

因此，Kriging預(yù)測(cè)數(shù)據(jù)在反映總體分布規(guī)律時(shí)數(shù)據(jù)不一定精確，但可以表示相對(duì)大小的趨勢(shì)，統(tǒng)計(jì)分析應(yīng)以實(shí)測(cè)數(shù)據(jù)為主。同時(shí)，對(duì)插值結(jié)果的可靠性可以用Kriging方差來(lái)指示，Kriging方差越大說(shuō)明結(jié)果的可靠性越差，否則結(jié)果就越可靠［15］?？梢?jiàn)大雨條件下的穿透率插值結(jié)果最可靠，暴雨條件下的最不可靠。

圖3 不同雨量級(jí)穿透率頻數(shù)分布

4 討論

本試驗(yàn)在觀測(cè)期內(nèi)，降水總量為947.3mm，總穿透雨量和穿透率分別為742.6mm和0.7839。穿透降雨量隨大氣降雨量的增加而增大，呈顯著線(xiàn)性正相關(guān)關(guān)系。很多相似研究從不同尺度、不同林分、不同氣候等角度均認(rèn)為線(xiàn)形方程能較好地?cái)M合穿透雨量和降雨量之間的關(guān)系［16－18］。

穿透率隨著大氣降雨量的增加而逐漸增大，最后趨于穩(wěn)定值，對(duì)數(shù)函數(shù)可以較好地描述穿透率隨降雨量的變化。這一結(jié)果與翟杰等［19］對(duì)紫金山麻櫟的研究，鮑文等［20］對(duì)岷江上游油松的研究結(jié)果相一致，而劉章文等［21］對(duì)祁連山灌叢的研究認(rèn)為指數(shù)函數(shù)擬合效果更好。實(shí)際上劉章文等所采用的是Mitscherlich模型，采用Mitscherlich模型擬合本實(shí)驗(yàn)數(shù)據(jù)的穿透率與降雨量的方程為：

可見(jiàn)Mitscherlich模型與對(duì)數(shù)函數(shù)擬合精度完全相同（R2＝0.439）。將 Mitscherlich模型式（7）與對(duì)數(shù)函數(shù)公式（5）中的Tr用T／P代替，將其分別經(jīng)過(guò)轉(zhuǎn)換得到穿透雨量與降雨量的關(guān)系公式：

以上兩式與線(xiàn)性方程（4）擬合精度完全相同（R2＝0.994），表明其與線(xiàn)性方程均可描述穿透雨量與降雨量的關(guān)系。因而可知，本實(shí)驗(yàn)采用對(duì)數(shù)函數(shù)與Mitscherlich模型也可較好擬合穿透率與降雨量之間的關(guān)系。李振新等［22］認(rèn)為穿透雨率同林外降雨量之間的關(guān)系可以用logistic曲線(xiàn)方程模擬的結(jié)論不適用于本研究，模擬精度相對(duì)稍差（R2＝0.428），但其穿透率存在一個(gè)最大值的特點(diǎn)與Mitscherlich模型相似。

穿透降雨變異系數(shù)與降雨量呈負(fù)相關(guān)關(guān)系，Mitscherlich模型可以較好地描述穿透降雨變異系數(shù)隨降雨量的變化。說(shuō)明穿透降雨空間分布的不均勻度隨大氣降雨量的增加有下降的趨勢(shì)，最終逐漸趨于穩(wěn)定值。這與戰(zhàn) 偉 慶 等［23］、Rodrigo A 等［24］的 研 究觀點(diǎn)一致，但所用擬合模型不同，他們分別認(rèn)為指數(shù)模型和對(duì)數(shù)模型較為合適。戰(zhàn)偉慶［23］認(rèn)為造成這一現(xiàn)象的原因是，當(dāng)林冠達(dá)到飽和持水量時(shí)，林外降雨幾乎全部轉(zhuǎn)化為穿透雨。而我們認(rèn)為此時(shí)林外降雨并未全部轉(zhuǎn)化為穿透雨，仍然有樹(shù)體吸收和蒸發(fā)等過(guò)程發(fā)生，只是林外降雨轉(zhuǎn)化為穿透雨的比例趨于穩(wěn)定。

5 結(jié)論

穿透降雨率的半變異函數(shù)分析表明，天然針闊混交林冠下穿透降雨率的不均勻度隨雨量級(jí)的增大而減弱，空間分布趨向于均勻化。且隨著雨量級(jí)的增大，隨機(jī)因素對(duì)穿透率的影響相對(duì)增加，而林冠結(jié)構(gòu)等結(jié)構(gòu)性因素對(duì)穿透降雨率產(chǎn)生的影響有逐漸變小的趨勢(shì)，但結(jié)構(gòu)性因素始終是最主要的影響因素。從穿透率的插值分布圖可以直觀地看到穿透率的空間分布在不同雨量級(jí)存在著明顯的差別。

由于Kriging預(yù)測(cè)數(shù)據(jù)存在“平滑效應(yīng)”，使得預(yù)測(cè)值的分布不能完全體現(xiàn)實(shí)際分布，故Kriging預(yù)測(cè)數(shù)據(jù)只用來(lái)反映總體分布規(guī)律，數(shù)據(jù)不精確，統(tǒng)計(jì)分析應(yīng)以實(shí)測(cè)數(shù)據(jù)為主。

天然次生針闊混交林下穿透雨空間分布的差異對(duì)林地的土壤水分分布和養(yǎng)分循環(huán)利用及生物多樣性必然會(huì)產(chǎn)生重要影響，有待于進(jìn)一步研究。

［1］ Anna A，Anselm R.Trace metal fluxes in bulk deposition，through－fall and stem－flow at two evergreen oak stands in NE Spain subject to different exposure to the industrial environment［J］.Atmospheric Environment，2004（38）：171－180.

［2］ 張焜.重慶四面山4種類(lèi)型天然林水文功能研究［D］.北京：北京林業(yè)大學(xué)，2012.

［3］ 王波，張洪江，杜士才，等.三峽庫(kù)區(qū)天然次生林凋落物森林水文效應(yīng)研究［J］.水土保持通報(bào)，2009，29（3）：83－87.

［4］ 孔維健，周本智，安艷飛，等.天然次生林和人工毛竹林水文生態(tài)特征比較［J］.水土保持研究，2010，17（1）：113－116.

［5］ 劉小林，鄭子龍，藺巖雄，等.甘肅小隴山林區(qū)主要林分類(lèi)型土壤水分物理性質(zhì)研究［J］.西北林學(xué)院學(xué)報(bào)，2013，28（1）：7－11.

［6］ Sun Ge，Mcnulty S G，Amaty A D M，et al.A comparison of the watershed hydrology of coastal forested wetlands and the mountainous uplands in the Southern US［J］.Journal of Hydrology，2002，263（1／4）：92－104.

［7］ 馬熙淵，馬立琦.小隴山林區(qū)主要林分對(duì)降水的再分配［J］.甘肅林業(yè)科技，2003，28（2）：1－4.

［8］ 譚俊磊，馬明國(guó)，車(chē)濤，等.基于不同郁閉度的青海云杉冠層截留特征研究［J］.地球科學(xué)進(jìn)展，2009，24（7）：825－833.

［9］ 時(shí)忠杰，王彥輝，熊偉，等.單株華北落葉松樹(shù)冠穿透降雨的空間異質(zhì)性［J］.生態(tài)學(xué)報(bào)，2006，26（9）：2877－2886.

［10］ 史利江.基于GIS和地統(tǒng)計(jì)學(xué)的土壤養(yǎng)分空間變異特征研究［D］.上海：上海師范大學(xué)，2006.

［11］ Esr.ArcGIS Resources［OL］.［2013－03－10］.http：∥resources.arcgis.com／en／help／previous－h(huán)elp／index.html.

［12］ 姜城.不同經(jīng)營(yíng)體制下土壤養(yǎng)分空間變異規(guī)律及管理技術(shù)的研究［D］.北京：中國(guó)農(nóng)業(yè)科學(xué)院研究生院，2000.

［13］ 王麗霞，段文標(biāo)，陳立新，等.紅松闊葉混交林林隙大小對(duì)土壤水分空間異質(zhì)性的影響［J］.應(yīng)用生態(tài)學(xué)報(bào)，2013，24（1）：17－24.

［14］ 李超.土壤水分的空間變異特性分析和農(nóng)田干旱評(píng)估［D］.北京：清華大學(xué)，2009.

［15］ Yamamoto J K.An alternative measure of the reliability of ordinary Kriging estimates［J］.Mathematical Geology，2000，32（4）：489－509.

［16］ 胡珊珊，于靜潔，胡堃，等.華北石質(zhì)山區(qū)油松林對(duì)降水再分配過(guò)程的影響［J］.生態(tài)學(xué)報(bào)，2010，30（7）：1751－1757.

［17］ 劉世海，余新曉.北京密云水庫(kù)庫(kù)區(qū)水源涵養(yǎng)林冠層水文特征研究［J］.林業(yè)科學(xué)，2005，41（1）：194－197.

［18］ 薛建輝，郝奇林，吳永波，等.3種亞高山森林群落林冠截留量及穿透雨量與降雨量的關(guān)系［J］.南京林業(yè)大學(xué)學(xué)報(bào)：自然科學(xué)版，2008，32（3）：9－13.

［19］ 翟杰，張志民，胡海波，等.紫金山麻櫟林降水分配格局研究［J］.林業(yè)科技開(kāi)發(fā)，2011，25（4）：63－67.

［20］ 鮑文，包維楷，何丙輝，等.岷江上游油松人工林對(duì)降水的截留分配效應(yīng)［J］.北京林業(yè)大學(xué)學(xué)報(bào)，2004，26（5）：10－16.

［21］ 劉章文，陳仁升，宋耀選，等.祁連山典型灌叢降雨截留特征［J］.生態(tài)學(xué)報(bào)，2012，32（4）：1337－1346.

［22］ 李振新，鄭華，歐陽(yáng)志云，等.岷江冷杉針葉林下穿透雨空間分布特征［J］.生態(tài)學(xué)報(bào)，2004，24（5）：1015－1021.

［23］ 戰(zhàn)偉慶，張志強(qiáng)，武軍，等.華北油松人工林冠層穿透雨空間變異性研究［J］.中國(guó)水土保持科學(xué)，2006，4（3）：26－30.

［24］ Rodrigo A，Avila A.Influence of sampling size in the estimation of mean throughfall in two Mediterranean holm oak forests［J］.Journal of Hydrology，2001，243（3／4）：216－227.