北京山區側柏林冠層對降雨動力學特征的影響

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

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

北京山區側柏林冠層對降雨動力學特征的影響

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

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

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

雨滴動能；雨滴大小；終點速度；雨強；北京山區；側柏林

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

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

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

1 研究區概況

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

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

表1 側柏樣地喬木胸徑、樹高特征

2 研究方法

2.1 儀器布設

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

2.2 數據選取

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

表2 雨滴數據選取情況

2.3 降雨動能的計算

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

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

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

3 結果與分析

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

3.1.1 林內外雨滴大小分布差異

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

表3 雨滴大小分布

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

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

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

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

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

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

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

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

3.2.1 林內外雨滴終點速度差異

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

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

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

3.2.1 雨滴終點速度與雨滴大小分布的相關性

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

3.3 林冠對雨滴動能的影響

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

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

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

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

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

3.3.2 林內外降雨動能對比

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

表4 林內外降雨動能對比

4 討論

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

[25] 王禮先, 解明曙. 山地防護林水土保持水文生態效益及其信息系統. 北京: 中國林業出版社, 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

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

2013- 07- 08;

2013- 10- 10

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

10.5846/stxb201307081861

史宇，余新曉，張建輝，羅海江，張佳音北京山區側柏林冠層對降雨動力學特征的影響.生態學報,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.