北京山區(qū)側(cè)柏林冠層對降雨動力學特征的影響

史 宇，余新曉，張建輝，羅海江，張佳音

(1. 中國環(huán)境監(jiān)測總站, 北京 100012; 2. 北京林業(yè)大學水土保持學院, 水土保持與荒漠化防治教育部重點實驗室, 北京 100083; 3. 中國環(huán)境科學研究院, 北京 100012)

北京山區(qū)側(cè)柏林冠層對降雨動力學特征的影響

史 宇1，余新曉2,*，張建輝1，羅海江1，張佳音3

(1. 中國環(huán)境監(jiān)測總站, 北京 100012; 2. 北京林業(yè)大學水土保持學院, 水土保持與荒漠化防治教育部重點實驗室, 北京 100083; 3. 中國環(huán)境科學研究院, 北京 100012)

林冠層枝葉對雨滴的能量要進行再次分配，從而在降雨雨滴大小分布、降落速度及動能等性質(zhì)上都發(fā)生了變化，是森林水文生態(tài)功能的重要環(huán)節(jié)。然而目前國內(nèi)外在森林影響降雨動力學特征方面的研究較少且研究方法不夠完善。研究地點設(shè)在北京妙峰山林場，利用激光雨滴譜儀對側(cè)柏人工林林內(nèi)和林外進行了長期同步觀測，對典型場降雨中林內(nèi)與林外雨滴的數(shù)量、大小、速度、動能等動力學特征值進行了對比分析。結(jié)果表明：林內(nèi)的雨滴總數(shù)高于林外37.8%，林內(nèi)小徑級和大徑級的雨滴數(shù)目都有所增加，而中間徑級雨滴減少。林外降雨體積貢獻率隨雨滴直徑增加呈單峰曲線變化，而林內(nèi)降雨呈多峰山狀變化；在試驗側(cè)柏林下雨滴的終點速度有所減小，特別是對直徑0.75—4.5mm范圍的雨滴更為明顯。林內(nèi)降雨相比于林外其雨滴直徑與雨滴終點速度的相關(guān)性降低，雨滴譜分布更加分散；林外降雨動能貢獻值隨雨滴直徑增加呈單峰曲線變化趨勢，而林內(nèi)降雨呈多峰山狀變化。在雨滴終點速度為4—9m/s的區(qū)間內(nèi)林外雨滴貢獻的動能明顯高于林內(nèi)降雨。側(cè)柏林地雨滴降落到地面的動能比林外空曠地減少了27.0%，林內(nèi)雨滴動能減小的最主要原因是林冠層削減了林內(nèi)降雨量。

雨滴動能；雨滴大??；終點速度；雨強；北京山區(qū)；側(cè)柏林

土壤侵蝕指地表土壤受外營力作用而被分散、搬運、堆積的物理過程，水力侵蝕的動力主要來自降雨動能，特別是較大降雨強度的大雨或暴雨更是造成水土流失的關(guān)鍵因素[1- 2]。雨滴擊濺引起土壤結(jié)構(gòu)破壞、顆粒分散及土粒濺蝕，形成結(jié)皮層和導致入滲下降及地表徑流增加，還會引起地表徑流紊動和土體沖刷。林冠是森林對雨滴動能產(chǎn)生影響的第一個作用層。降雨在進入林地前具有一定的動能與勢能，降雨通過林冠層后，由于受到冠層的再分配作用而以多種形式下落，出現(xiàn)林冠截持和樹干流等現(xiàn)象，林冠層枝葉對雨滴的能量要進行再次分配，從而在降雨雨滴大小分布及降落速度、動能等性質(zhì)上都發(fā)生了變化[3- 6]。森林植被對降雨動能的削減是其水土保持功能的重要體現(xiàn)[7]。因此，定量評價植被對降雨動能的影響和研究植被調(diào)節(jié)降雨動能的機制是植被的生態(tài)水文功能的研究中非常有意義的工作。

雨滴大小組成和雨滴速度是研究森林對雨滴動能的影響的兩個要素。降雨的雨滴大小組成(即雨滴譜)直接決定了雨滴的速度及動能，與降雨的濺蝕力密切相關(guān)[8]。關(guān)于雨滴大小的測量方法主要有色斑法、面粉球法[9]、油滴法[10]、高速攝影或照像法、雨滴測量儀法等，國內(nèi)之前的相關(guān)研究以濾紙色斑法為主。國內(nèi)外的學者對降雨的地面終點速度提出了不同的計算方法[12- 14]，但這些方法有的計算結(jié)果相差比較大，對不同直徑的雨滴計算方法也存在一定的局限性。

國內(nèi)的對林冠層調(diào)控雨滴動能的研究從20世紀80年代陸續(xù)開展研究對象以人工林居多[5,15- 19]，天然林較少[6,8]，在儀器設(shè)備和觀測計算方法上都遠落后于國際水準。本研究利用了激光雨滴譜儀對北京山區(qū)側(cè)柏人工林林內(nèi)和林外同步進行了連續(xù)觀測。此種方法的優(yōu)勢在于通過激光原理可直接測定某段時間內(nèi)通過儀器受雨面積的雨滴大小、數(shù)量和速度，并直接利用動能公式計算雨滴動能，避免了利用經(jīng)驗公式推算帶來的不確定性，測定結(jié)果較為精準。本文利用激光雨滴譜儀長時間定點觀測得到的數(shù)據(jù)對2008—2010年多場典型降雨中林內(nèi)與林外雨滴的數(shù)量、大小、速度、動能等動力學特征值進行了對比分析，初步探討了林冠層對雨滴結(jié)構(gòu)和能量調(diào)控的特征和機理。

1 研究區(qū)概況

本文研究數(shù)據(jù)出自位于北京妙峰山林場的國家林業(yè)局首都圈森林生態(tài)系統(tǒng)定位觀測研究站，該研究站地處北京市西北郊太行山北部，燕山東端，北緯39°54′，東經(jīng)116°28′。妙峰山林場屬華北大陸型季風氣候，春季干旱多風，夏季涼爽多雨，冬季干燥寒冷。年平均氣溫12.2℃，最高氣溫39.7℃，最低氣溫19.6℃，無霜期180天，年降水量近700mm，多集中在7、8月份。

監(jiān)測樣地海拔170m，坡向為南偏東68°，坡度15°，土層瘠薄，A+B層厚度為20 cm左右。樣地為側(cè)柏林，林分密度為1531株/hm2，平均年齡53 a，郁閉度0.7，林分樹高和胸徑特征如表1所示。林下灌木層蓋度40%，主要有構(gòu)樹(Broussonetiapapyrifera)、酸棗(ZizyphusjujubaMill)、小葉鼠李(Rhamnusparvifolia)、荊條(Vitexnegundovar.heterophylla)、孩兒拳頭(Grewiabiloba)、胡枝子(Lespedezabicolor)等。草本層蓋度20%，主要有鴨跖草(Commelinacommunis)、羊胡子草(Eriophorumvaginatum)、菅草(Themedajaponica)、中華卷柏(Selaginellasinensis)等。監(jiān)測樣地林分為研究地區(qū)典型的側(cè)柏人工中齡林的結(jié)構(gòu)特征。

表1 側(cè)柏樣地喬木胸徑、樹高特征

2 研究方法

2.1 儀器布設(shè)

2008年6月至2010年10月期間，在北京妙峰山側(cè)柏林監(jiān)測樣地內(nèi)和該樣地50m外的空曠地處布設(shè)兩臺LPM激光雨譜儀(Laser Precipitation Monitor，德國，thiesclima公司)同時進行林內(nèi)和林外降雨監(jiān)測。該儀器應用激光原理對高速運動物體進行測定，儀器觀測高度為91.0cm，探測面積為45.6cm2，測定對象最小直徑達到0.16mm，測定時間間隔為1s，可實時測定雨滴的總量，大小，強度，和運動速度，可體現(xiàn)場降雨的整個變化過程。同時利用激光雨譜儀配套的自動氣象站監(jiān)測空氣溫度(℃)、空氣相對濕度(%)，風速(m/s)和風向等。用CI- 110植物冠層圖像分析儀測定雨滴譜儀上方75°天頂角的葉面積指數(shù)為1.65。

2.2 數(shù)據(jù)選取

研究區(qū)降雨有較強的地帶性特征，主要雨型可分為兩種：一是由局部地形和氣候影響產(chǎn)生的來勢猛、歷時短(1 h左右)的小面積降雨，稱為短陣性雨型；二是主要是季風(鋒面)影響的大面積普通降雨雨型。前者出現(xiàn)次數(shù)較多但降雨量的貢獻率不大，而后者出現(xiàn)次數(shù)偏少但對年降水量的貢獻率可占80%以上，是地表徑流和土壤侵蝕的主要驅(qū)動因子，因此將后一種雨型作為主要的研究對象。為保證數(shù)據(jù)的典型性，分別選取了2008年1場，2009年3場和2010年2場較長時間降雨，從每場降雨經(jīng)過的首個整點時間開始每隔30min取樣1次，共獲得100組1min內(nèi)雨滴特征數(shù)據(jù)值，具體取樣信息見表2。

表2 雨滴數(shù)據(jù)選取情況

2.3 降雨動能的計算

根據(jù)通用動能計算公式，單個雨滴的動能Ed(J)主要取決于雨滴的直徑D(mm)和下落時的終點速度μT(D)(m/s)，Ed可以根據(jù)下式來計算：

式中,ρ為水的密度，標準狀態(tài)下為1.0×10-6kg/mm3。

雨滴終點速度與雨滴形狀、大小和空氣狀況等因素緊密相關(guān)，但由于這些因素都存在很大的變異性，并且很難通過測量的手段獲得，從而給雨滴動能的計算帶來了困難。實際計算中通常使用經(jīng)驗的統(tǒng)計關(guān)系來代替復雜的基于物理過程的計算。本次研究采用LPM激光雨譜儀探測可以直接測定出林內(nèi)和林外的雨滴終點速度，不需要通過經(jīng)驗公式來進行換算。將雨譜儀測得的100組雨滴速度和雨滴直徑來代入能量計算公式，可測得每個徑級的雨滴在各次降雨中的動能情況，將同一次降雨中所有徑級的雨滴相加可得出每個1min取樣時間內(nèi)的總降雨動能，再用總動能除以儀器采樣面積0.00456m2就能得出單位面積每分鐘降雨動能值Es(J m-2min-1)。用Es除以每分鐘降雨量得出了單位降雨量、單位時間、單位面積的降雨動能，用Ep(J m-2min-1mm-1)來表示。

3 結(jié)果與分析

3.1 林冠對雨滴大小分布的影響

3.1.1 林內(nèi)外雨滴大小分布差異

表1反映了林內(nèi)和林外所有樣本中各直徑級別的雨滴數(shù)目分布情況，可以看出林內(nèi)外雨滴直徑大都分布在0.125—0.5mm之間，但林內(nèi)與林外的各徑級雨滴分布規(guī)律有著一定的差異。林內(nèi)的雨滴總數(shù)為149053個而林外雨滴總數(shù)為108130個，林內(nèi)的雨滴總數(shù)高于林外37.8%，林內(nèi)≤0.25mm的雨滴無論從數(shù)量還是百分比上都要高于林外，而林外雨滴減少幅度最大的是直徑處于0.75—2.00mm之間的雨滴，可以確定林冠能將體積較大的雨滴分散形成了更多小雨滴，這與羅德和王琛的研究結(jié)果近似[18- 19]。

表3 雨滴大小分布

雨滴個數(shù)為雨譜儀的有效測定面積上(0.0046m2)測得的數(shù)值

林冠層除了分散雨滴作用外還有著匯集雨滴的作用。從表3中可以看出，林內(nèi)的雨滴數(shù)量增加幅度最大的是0.375mm徑級，對比林外此徑級的雨滴數(shù)量增加了114%，占雨滴總數(shù)的百分比也從11.25%增加到了17.51%。據(jù)雷瑞德測定[16]，一枚針葉聚集自然下落的雨滴直徑為2.8—3.4mm，而兩枚針葉為3.4—3.8mm，從本研究結(jié)果上看側(cè)柏葉聚集雨滴的特征與兩枚針葉的更為近似。同時也可以看出，大于3mm的每個徑級林內(nèi)的雨滴數(shù)目都多于林外。林外降雨幾乎沒有直徑大于4mm的雨滴，而林內(nèi)則出現(xiàn)了一些7mm以上的大雨滴。

3.1.2 不同徑級雨滴對降雨的貢獻率

林冠層對雨滴大小分布的改變對降雨的影響需要分析各徑級雨滴對總降雨的體積貢獻率才能有更全面的認識。圖1反映了全部樣本中林內(nèi)和林外不同徑級雨滴對總降雨的體積貢獻率變化，圖2則反映了累積體積貢獻率的變化情況。從圖1可以看出，林外降雨體積貢獻率最大的雨滴直徑范圍是0.5—2mm，整體上看呈單峰曲線變化，其中貢獻率最大的為1mm徑級雨滴，達到了15.1%；而林內(nèi)降雨體積貢獻率隨雨滴直徑增加呈多峰山狀變化，在0.5、4、5mm處為體積貢獻率的峰值，貢獻率最大的為徑級在4mm的雨滴，達到了10.9%。圖2反映了林內(nèi)外不同徑級雨滴對總降雨累積體積貢獻率，從圖中可以看出林內(nèi)外雨滴體積累計分布曲線有一定差別，林內(nèi)曲線較為平緩，雨滴體積積累隨雨滴直徑呈緩慢均勻上升狀態(tài)，而林外曲線則上升急促，到雨滴直徑為2mm處時的體積累積已超過了90%。林內(nèi)外雨滴各徑級累積體積的分布規(guī)律表明，由于林冠對降雨的攔截再分配功能，林下降雨雨滴在各個徑級間更為分散，中間直徑的雨滴數(shù)量減少，體積貢獻率變低，而小雨滴和大雨滴的數(shù)量都有所增加，其中大雨滴的體積貢獻率有了明顯增加，這與張穎和王琛等人的研究結(jié)論相同[5,19]。

圖1 林內(nèi)外不同徑級雨滴對總降雨體積貢獻率Fig.1 Total rainfall volume contribution of different diameter level raindrops inside and outside forest

圖2 林內(nèi)外不同徑級雨滴對總降雨累積體積貢獻率Fig.2 Total rainfall volume accumulation contribution of different diameter level raindrops inside and outside forest

圖3 林內(nèi)外雨滴終點速度隨雨滴直徑變化情況 Fig.3 Change of raindrops terminal speed along with raindrops diameter outside and inside forest

3.2 林冠對雨滴終點速度的影響

3.2.1 林內(nèi)外雨滴終點速度差異

將選取的100個1min內(nèi)雨滴數(shù)據(jù)樣本中林內(nèi)外各徑級的雨滴平均速度隨雨滴直徑繪制成折線圖(圖3)。從圖中可以看出，林內(nèi)外雨滴終點速度都隨著雨滴直徑的增加而增加。當雨滴直徑小于0.75mm時林內(nèi)和林外的相同直徑的雨滴終點速度差異很小，當雨滴直徑在大于0.75mm時，同一徑級林內(nèi)的雨滴速度開始低于林外并且隨著雨滴直徑的增加林內(nèi)外雨滴終點速度差異逐漸加大。林外雨滴在直徑達到3mm后速度穩(wěn)定在7m/s左右，而直徑大于4mm的林內(nèi)雨滴隨著直徑增加其雨滴終點速度一直在不斷增加，當雨滴直徑大于6mm后還出現(xiàn)了一段較快的增加趨勢，林內(nèi)大雨滴的最高速度能達到7.8m/s左右。

圖4 林外雨滴直徑-速度分布Fig.4 Distribution of raindrops diameter-velocity outside forest圖中實線為Gunn-Kinzer曲線

圖5 林內(nèi)雨滴直徑-速度分布Fig.5 Distribution of raindrops diameter-velocity inside forest圖中實線為Gunn-Kinzer曲線

3.2.1 雨滴終點速度與雨滴大小分布的相關(guān)性

為研究林內(nèi)與林外雨滴大小和雨滴終點速度分布的差異性，選取了2010年8月21日降雨中時間接近、雨強相似的兩組雨滴數(shù)據(jù)樣本，儀器測定并由自帶軟件生成的林內(nèi)與林外雨滴分布譜圖見圖4和圖5。Atlas利用Gunn和Kinzer在海平面標準大氣的實驗觀測結(jié)果[14,20]，提出了雨滴直徑與雨滴終點速度的經(jīng)驗公式V=9.65-10.3e-0.6D[20,14]，式中V為雨滴終點速度(m/s)，D為雨滴直徑(mm)。從圖4中可以看出，林外的雨滴直徑和終點速度分布基本上符合Gunn-Kinzer曲線的規(guī)律，可以認為Gunn-Kinzer公式對于研究區(qū)的降雨有著較好的適應性。而林內(nèi)雨滴直徑和終點速度的分布較為較分散，出現(xiàn)了低速(4.2—5.8m/s)的較大雨滴(4—4.5mm)，和少量速度大于5.8m/s的較小雨滴(0.25—0.375m/s)。

3.3 林冠對雨滴動能的影響

3.3.1 林內(nèi)外不同雨滴直徑和終點速度對降雨動能的貢獻

圖6 林內(nèi)外不同直徑雨滴的動能貢獻 Fig.6 Contribution for rainfall kinetic energy by different diameter level raindrops outside and inside forest

圖6表示所有樣本降雨數(shù)據(jù)中不同徑級雨滴貢獻的雨滴動能Es情況。從圖中可以看出林外降雨總動能值Es隨雨滴直徑增加呈單峰曲線變化趨勢，0.5—3mm徑級的雨滴貢獻了95%以上的降雨動能，小于0.5mm的雨滴數(shù)量眾多但由于其質(zhì)量小而且終點速度較低，其對動能的貢獻率不大，而大于3mm的雨滴數(shù)量過少對總動能的貢獻也很小。而林內(nèi)降雨各徑級對雨滴總動能的貢獻隨雨滴徑級的增加呈多峰山狀變化，各徑級之間雨滴提供的動能較林外降雨來說差異較小，直徑在2—7.5mm范圍的雨滴貢獻的動能占總動能的絕大部分。可見林冠層的存在削減了林外雨滴直徑0.5—3mm的降雨動能，并將削減動能的大部分轉(zhuǎn)化給聚合生成的大徑級雨滴。

圖7 林內(nèi)外不同終點速度雨滴的動能貢獻 Fig.7 Contribution for rainfall kinetic energy by different raindrops terminal speed outside and inside forest

圖7表示了所有樣本降雨數(shù)據(jù)中不同終點速度的雨滴貢獻動能Es情況。從圖中可以看出林內(nèi)和林外終點速度小于2.6m/s的雨滴貢獻動能的情況差異不大，這個速度范圍內(nèi)的雨滴動能對總降雨動能的貢獻程度都很微小。在雨滴終點速度處于2.6—11m/s的區(qū)間內(nèi)林外雨滴動能隨終點速度增加總體上呈單峰曲線變化趨勢，而林內(nèi)降雨的動能呈多峰山狀變化。在雨滴終點速度為4—9m/s的區(qū)間內(nèi)林外雨滴的動能明顯高于林內(nèi)降雨，林內(nèi)降雨只在速度為10mm處雨滴動能高于林外降雨?？偨Y(jié)林內(nèi)外的變化規(guī)律，可以認為林冠層主要削減了雨滴終點速度為4—9m/s的林外降雨的雨滴動能值。

3.3.2 林內(nèi)外降雨動能對比

通過對雨譜儀實測的100組數(shù)據(jù)得出的林內(nèi)外雨滴動能各個相關(guān)參數(shù)的對比情況如表4所示。從表中可以看出，100組數(shù)據(jù)的林內(nèi)平均雨強為4.35mm/h，而林外為5.68mm/h，這種差異反映了林冠層對降雨的截持作用，相當于占總量23.4%的降雨無法降落到地面。林內(nèi)平均每分鐘每平方米的雨滴個數(shù)為324028.26個，而林外此數(shù)值為235065.21個，林內(nèi)比林外高38%。林外單個雨滴的動能平均值Ed為29.02μJ，而林內(nèi)Ed值為24.56μJ。從雨滴總數(shù)量上看林內(nèi)要高于林外，但由于林內(nèi)大量的小雨滴其動能值非常微小，導致降低了林內(nèi)雨滴的平均動能值有所降低。平均Es值是林內(nèi)外雨滴動能情況的綜合反映，林外平均Es值為0.859 J m-2min-1，而林內(nèi)為0.627 J m-2min-1，相當于側(cè)柏林地雨滴降落到地面的動能比林外空曠地減少了27%，這個結(jié)果說明林冠層對降雨動能的削減作用較為明顯。而平均Ep值反映了假定林內(nèi)外降雨量完全相同情況下的動能差異情況，可以看出林外平均Ep值為6.821 J m-2min-1，而林內(nèi)為7.957 J m-2min-1，林內(nèi)要高于林外。這反映了林冠層枝葉積聚雨滴以及其他的相關(guān)物理過程能夠使等降雨量的降雨動能有所增加。

表4 林內(nèi)外降雨動能對比

4 討論

(1)降雨通過林冠層時被分解為許多分量，一部分是未經(jīng)樹體接觸而從林冠空隙中穿過的自由穿透降雨，其降雨特征同林外降雨相比沒有發(fā)生變化；一部分是林冠截持降雨，它被林冠枝葉截持并最終直接蒸發(fā)到空氣中，不會對地面發(fā)生作用，可認為這部分降雨的動能被林冠層完全消減了；另一部分是沿樹干緩緩流到林地的干流，樹干流的不存在對地面的擊濺作用，其降雨動能可忽略不計；還有一部分降雨，雖被林冠攔截過，但最終還是降落到林地，稱為林冠降雨。實際觀測中不能把直接穿透降雨和林冠降雨分開，常合稱為穿透降雨。本研究所測定的林內(nèi)降雨動能即為廣義上的穿透降雨，其與林外降雨的差值能夠比較準確的反映林冠層對降雨動能削減的程度。

本研究發(fā)現(xiàn)林下降雨雨滴在各個徑級間更為分散，中間直徑的雨滴數(shù)量減少，體積貢獻率變低，而小雨滴和大雨滴的數(shù)量都有所增加，其中大雨滴的體積貢獻率有了明顯增加。通過對林內(nèi)降雨中雨滴大小的分布規(guī)律研究以及與林外降雨的對比分析，可認為林冠降雨也可分為兩個部分：一部分是林冠層通過枝葉和雨滴的碰撞將自然降雨中大量的雨滴分散成了更小的雨滴，可稱作飛濺液滴；另一部分是由于葉片的匯聚作用將一部分的小雨滴聚合成直徑較大的雨滴降落在林下地面，可稱為葉片滴落雨。而研究得出的林內(nèi)降雨雨滴分布規(guī)律就是自由穿透降雨、飛濺液滴和葉片滴落雨的共同反映，這與Nanko等得出的研究結(jié)論基本一致[21]。

(2)雨滴終點速度是降雨參數(shù)中一個重要的微物理量，決定了降水動能的大小。雨滴下落的過程中同時受到3種力共同作用，即重力，空氣阻力和空氣浮力。在這3種力的共同作用下，雨滴的下落運動實際上是一個加速度逐漸減小直至為零的變加速運動，當加速度為零時，雨滴作勻速運動。雨滴速度從0到勻速的過程所經(jīng)歷的時間稱為弛豫時間，所走過的路程稱為慣性路程，最后所達到的均勻速度稱為雨滴終點速度[22]。而林冠層通過截留、分散和匯聚作用改變了雨滴的原有運動規(guī)律，會使林內(nèi)與林外的雨滴終點速度分布規(guī)律出現(xiàn)一定的差異性。

研究結(jié)果表明，當雨滴直徑小于0.75mm時林內(nèi)和林外的相同直徑的雨滴終點速度差異很小，直徑較小的雨滴能在較快時間和較短的運行距離內(nèi)達到勻速運動。當雨滴直徑在大于0.75mm時，相同徑級的雨滴林內(nèi)速度開始低于林外并且隨著雨滴直徑的增加林內(nèi)外雨滴終點速度差異逐漸加大，這說明在此直徑范圍內(nèi)林內(nèi)雨滴無法在8m左右的降落高度內(nèi)達到終點速度，并且林內(nèi)雨滴直徑越大其終點速度與勻速運動速度的差異值越大。而在林內(nèi)通過林冠匯聚作用重新生成的葉片滴落雨隨著直徑增加其雨滴終點速度一直在不斷增加，說明葉片滴落雨可通過自身質(zhì)量的增加不斷提高下落時的加速度從而在相同的運動距離內(nèi)達到更大的終點速度。林內(nèi)雨滴由于林冠的截留、分散和匯聚作用，相比于林外其雨滴直徑與雨滴終點速度的相關(guān)性降低，雨滴譜分布更加分散，這與Frasson和Krajewski在玉米冠層下得出的研究結(jié)論相似[23]。

(3)從能量守恒的角度來講，雨滴從天空形成直至降落地面的過程中林冠層對降雨的阻攔相當于增加了雨滴下降的阻力，林冠層對雨滴做負功，則雨滴的總動能必然減小。通過研究得出相同受雨面積下林內(nèi)降雨動能比林外減小了27.0%，證明了側(cè)柏林對降雨的動能有著有效的削減效果。同時，研究得出林內(nèi)降水量比林外減小了23.4%，而減少的這部分降雨被截留在枝葉表面或沿樹干緩緩流下，它們的動能基本上完全損失，說明了林冠層削減降雨動能主要還是通過林冠截留作用實現(xiàn)的。

研究中發(fā)現(xiàn)林內(nèi)降雨單位雨滴的平均動能Ed和單位面積單位降水量的動能Ep平均值均要高于林外降雨。經(jīng)分析可知，造成這種現(xiàn)象的原因在于林冠枝葉能將降雨中小雨滴匯集成了大雨滴最終降落至地面，這部分經(jīng)林冠匯聚作用新生成的葉片滴落雨同時具有較大的質(zhì)量和速度，是林內(nèi)降雨動能的主要組成部分。之前的一些研究也認為由于林冠枝葉積聚雨滴作用，林下雨滴動能或濺蝕力在某些情況下有所增加[16, 24- 25]。因此，植被覆蓋可以在削減林內(nèi)降雨總量方面降低林內(nèi)土壤的水力侵蝕，但在局部可能會加重雨滴對土壤表面的擊濺侵蝕作用。

林冠層對降雨的能量調(diào)節(jié)是一個非常復雜的過程，這與植被的枝葉特性、降落高度、降雨類型、降雨強度、氣象因素等都有著密切關(guān)系。本研究由于條件有限，僅限于研究單一的樹種和林分結(jié)構(gòu)下林內(nèi)與林外降雨動能的變化情況，并不能反映林冠結(jié)構(gòu)以及降雨期間氣象條件對動能的影響。希望今后此類研究能夠繼續(xù)深入開展，并與先進的傳感技術(shù)相結(jié)合，進一步揭示林冠層對降雨能量調(diào)節(jié)作用的規(guī)律。

[1] Zhang H X, Wang W Z. The characteristics of hard rain and its distribution over the Loess Plateau. Bulletin of Soil and Water Conservation, 1983, (1): 35- 44.

[2] Wang W Z. Study on the relations between rain fall characteristics and loss of soil in loess region. Bulletin of Soil and Water Conservation, 1983, 3(4): 7- 13, 65- 65.

[3] Lei A L, Zhang X D, Tang K L. Contrast of a few calculating methods on fal1 velocity of water drops. Bulletin of Soil and Water Conservation, 1995, 15(4): 43- 47.

[4] Xu R. Calculation method of natural rainfall and artificial rainfall kinetic energy. Soil and Water Conservation in China, 1983, (3): 37- 39.

[5] Zhang Y, Xie B Y, Yu X X, Niu J Z, Li W, Yue Y J. Influence of young growth canopy of several typical tree species on rainfall characteristics in the Loess Plateau. Journal of Beijing Forestry University, 2009, 31(4): 70- 76.

[6] Yin H, Guan W B, Xue X X, Xie C H. Redistribution of rainfall energy by canopy of dark coniferous forests in Gongga Mountain, southwestern China. Journal of Beijing Forestry University, 2010, 32(2): 1- 5.

[7] Luo W X, Bai L Q, Song X D, Liu T Y, Ma Z M. Runoff and scouring amount in forest and grass land with different cover rate. Journal of Soil and Water Conservation, 1990, 4(1): 30- 35.

[8] Li Z X, Ouyang Z Y, Zheng H, Jiang L J, Liu X L, Su Y M. Changes of raindrop composition and rainfall energy under coniferous forest and sub- alpine shrubbery in Wolong nature reserve region. Journal of Soil and Water Conservation, 2004, 18(4): 125- 133.

[9] Zhou P H, Wu C L. The research method of soil anti-scourability experiment in Loess Plateau. Journal of Soil and Water Conservation, 1993, 7(1): 29- 34.

[10] Zhou G Y. Hydrothermal Principle and Application of Ecosystem. Beijing: China Meteorological Press, 1997: 84- 85.

[11] Luo D, Yu X X, Dong L. Changes of nature rainfall energy through forest canopy in Miyun mountain region. Journal of Soil and Water Conservation, 2008, 22(3): 60- 63.

[12] Lei A L, Zhang X D, Tang K L. Contrast of a few calculating methods on fal1 velocity of water drops. Bulletin of Soil and Water Conservation, 1995, 15(4): 43- 47.

[13] Beard K V, Pruppacher H R. A determination of the terminal velocity and drag of small water drops by means of a wind tunnel. Journal of Atmospheric Sciences, 1969, 26(5): 1066- 1072.

[14] Gunn R, Kinzer G D. The terminal velocity of fall for water droplets in stagnant air. Journal of Meteorology, 1949, 6(4): 243- 248.

[15] Zhou G Y. Distribution of rainfall kinetic energy by canopies of artificial forest tree species, and its ecological effects. Acta Phytoecologica Sinica, 1997, 21(3): 250- 259.

[16] Lei R D. Influences of canopy of pine upon kinetic energy of precipitation. Journal of Soil and Water Conservation, 1988, 2(2): 31- 39.

[17] Wu G Y, Cheng J, Zhu Z H, Wu F Q, Yu X L. Rainfall characteristics under canopy of plantation forest in the southern Loess Plateau. Research of Soil and Water Conservation, 2011, 18(4): 32- 35.

[18] Luo D. Studies on the Dynamic Characteristics of Rainfall under the Forest in Beijing Mountainous Area [D]. Beijing: Beijing Forestry University, 2008.

[19] Wang C. Research of Microclimatic Characteristics of Forest in Beijing [D]. Beijing: Beijing Forestry University, 2010.

[20] Atlas D, Srivastava R C, Sekhon R S. Doppler radar characteristics of precipitation at vertical incidence. Reviews of Geophysics, 1973, 11(1): 1- 35.

[21] Nanko K, Hotta N, Suzuki M. Evaluating the influence of canopy species and meteorological factors on throughfall drop size distribution. Journal of Hydrology, 2006, 329(3/4): 422- 431.

[22] Jia X C. Observational Studies on Summer Rain Clouds from Different Area Raindrop Spectrum Characteristics [D]. Nanjing: Nanjing University of Information and Science Technology, 2009.

[23] de Moraes Frasson R P, Krajewski W F. Characterization of the drop-size distribution and velocity-diameter relation of the throughfall under the maize canopy. Agricultural and Forest Meteorology, 2011, 151(9): 1244- 1251.

[24] Zhou Y, Wats D. Hydrological effect of Pinus Yunnanensis Forest on soil erosion control in the alpine gorge region. Journal of Soil and Water Conservation, 1998, 4(3): 31- 38.

[25] Wang L X, Xie M Z. Mountain Shelter Forest Ecological Benefits of Soil and Water Conservation in Hydrology and Information System. Beijing: China Forestry Press, 1997.

參考文獻：

[1] 張漢雄, 王萬忠. 黃土高原的暴雨特性及分布規(guī)律. 水土保持通報, 1983, (1): 35- 44.

[2] 王萬忠. 黃土地區(qū)降雨特性與土壤流失關(guān)系的研究. 水土保持通報, 1983, 3(4): 7- 13, 65- 65.

[3] 雷阿林, 張學棟, 唐克麗. 幾種計算水滴降落速度方法的比較. 水土保持通報, 1995, 15(4): 43- 47.

[4] 徐銳. 關(guān)于天然降雨和人工降雨的動能計算方法. 中國水土保持, 1983, (3): 37- 39.

[5] 張穎, 謝寶元, 余新曉, 牛健植, 李維, 岳永杰. 黃土高原典型樹種幼樹冠層對降雨雨滴特性的影響. 北京林業(yè)大學學報, 2009, 31(4): 70- 76.

[6] 殷暉, 關(guān)文彬, 薛肖肖, 謝春華. 貢嘎山暗針葉林林冠對降雨能量再分配的影響研究. 北京林業(yè)大學學報, 2010, 32(2): 1- 5.

[7] 羅偉祥, 白立強, 宋西德, 劉天義, 馬志明. 不同覆蓋度林地和草地的徑流量與沖刷量. 水土保持學報, 1990, 4(1): 30- 35.

[8] 李振新, 歐陽志云, 鄭華, 姜立軍, 劉興良, 宿以明. 臥龍地區(qū)針葉林及亞高山灌叢對降雨的雨滴譜及能量的影響. 水土保持學報, 2004, 18(4): 125- 133.

[9] 周佩華, 武春龍. 黃土高原土壤抗沖性的試驗研究方法探討. 水土保持學報, 1993, 7(1): 29- 34.

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

[11] 羅德, 余新曉, 董磊. 密云山區(qū)林冠層對天然降雨能量影響的初步研究. 水土保持學報, 2008, 22(3): 60- 63.

[12] 雷阿林, 張學棟, 唐克麗. 幾種計算水滴降落速度方法的比較. 水土保持通報, 1995, 15(4): 43- 47.

[15] 周國逸. 幾種常用造林樹種冠層對降水動能分配及其生態(tài)效應分析. 植物生態(tài)學報, 1997, 21(3): 250- 259.

[16] 雷瑞德. 華山松林冠層對降雨動能的影響. 水土保持學報, 1988, 2(2): 31- 39.

[17] 吳光艷, 成婧, 祝振華, 吳發(fā)啟, 于小玲. 黃土高原南部人工林林冠對降雨特征的影響分析. 水土保持研究, 2011, 18(4): 32- 35.

[18] 羅德. 北京山區(qū)森林植被影響下的降雨動力學特性研究 [D]. 北京: 北京林業(yè)大學, 2008.

[19] 王琛. 北京地區(qū)森林小氣候特征研究 [D]. 北京: 北京林業(yè)大學, 2010.

[22] 賈星燦. 夏季不同地區(qū)降水云系雨滴譜特征的觀測研究 [D]. 南京: 南京信息工程大學, 2009.

[24] 周躍, Wats D. 高山峽谷區(qū)云南松林土壤侵蝕控制的水文效應. 水土保持學報, 1998, 4(3): 31- 38.

[25] 王禮先, 解明曙. 山地防護林水土保持水文生態(tài)效益及其信息系統(tǒng). 北京: 中國林業(yè)出版社, 1997.

ResearchonchangesofdynamiccharacteristicsofrainfallthoughPlatycladusOrientalisplantationcanopyinBeijingMountainArea

SHI Yu1，YU Xinxiao2,*，ZHANG Jianhui1, LUO Haijiang1, ZHANG Jiayin3

1ChinaNationalEnvironmentalMonitoringCenter,Beijing100012,China2CollegeofSoilandWaterConservationBeijingForestryUniversity&KeyLaboratoryofSoilandWaterConservationDeserticationCombatingofMinistryofEducation,Beijing100083,China3ChineseResearchAcademyofEnvironmentalSciences,Beijing100012,China

Due to forest canopy layer redistributes the energy of rainfall, the size distribution, terminal velocity and kinetic energy of raindrops change through the canopy. These changes are important part of the ecological hydrology functions of forest. However, at home and abroad, there is less studies on forest affecting the dynamic characteristic of rainfall is less, and the research methods are always imperfect. This study, performed at the Miaofeng Mountain forest farm in Beijing, used two laser raindrop spectrometers to synchronously monitor precipitation both inside and outside aPlatycladusorientalisplantation, and analyzed the size, amount, velocity, kinetic energy and other dynamic characteristics of raindrops outside and inside forest for the typical rainfall events. The results show the amount of raindrops outside forest was 37.8% above inside forest, the amount of small and big diameter raindrops increased inside forest, but the middle diameter raindrops decreased. The rainfall volume contribution outside forest showed unimodal curve with the increase of raindrop diameter change, but that showed polymodal curve inside forest. The terminal velocity reduced inside forest, especially for the raindrops with diameter in the range of 0.75—4.5 mm. Compared with outside forest, the correlation between the diameter and terminal velocity of raindrops inside forest was reduced, and the raindrop spectrum distribution was more scattered. The rainfall kinetic energy contribution outside forest showed unimodal curve with the increase of raindrop diameter change, but that showed polymodal curve inside forest. When the raindrops terminal velocity was in the range of 4—9m/s, the kinetic energy contribution of the raindrops outside forest was obviously higher than inside forest. The total kinetic energy of the raindrop fallen on the ground inside thePlatycladusorientalisplantation was 27.0% less than outside forest, and the main reason of that is the forest canopy cut down the precipitation inside forest.

raindrop kinetic energy; raindrop size; terminal velocity; Beijing mountain area;Platycladusorientalis

林業(yè)公益性行業(yè)科研資助項目(201104005)

2013- 07- 08;

2013- 10- 10

*通訊作者Corresponding author.E-mail: yuxinixao111@126.com

10.5846/stxb201307081861

史宇，余新曉，張建輝，羅海江，張佳音北京山區(qū)側(cè)柏林冠層對降雨動力學特征的影響.生態(tài)學報,2013,33(24):7898- 7907.

Shi Y，Yu X X，Zhang J H, Luo H J, Zhang J Y.Research on changes of dynamic characteristics of rainfall thoughPlatycladusOrientalisplantation canopy in Beijing Mountain Area.Acta Ecologica Sinica,2013,33(24):7898- 7907.