免费文献传递   相关文献

Theoretical research on a model for predicting the shadow boundary of and individual conical crown on a slope

坡面单木圆锥形树冠投影边界的模型理论



全 文 :第 26 卷第 10 期
2006 年 10 月
生   态   学   报
ACTA ECOLOGICA SINICA
Vol. 26 ,No. 10
Oct. ,2006
坡面单木圆锥形树冠投影边界的模型理论
张 蕾 ,于海业 3 ,杨昊谕 ,张志平
(吉林大学生物与农业工程学院 ,地面机械仿生技术教育部重点实验室 ,长春 130022)
基金项目 :吉林省科技厅重点资助项目 (2004020721)
收稿日期 :2006201211 ;修订日期 :2006207229
作者简介 :张蕾 (1979~) ,女 ,吉林人 ,博士生 ,主要从事农业生物环境研究. E2mail : zhanglei20902 @sina. com3 通讯作者 Corresponding author. E-mail :haiye @jlu. edu. cn
Foundation item :The project was supported by Emphases projects of Department of Science & Technology of Jilin Province (No. 2004020721)
Received date :2006201211 ;Accepted date :2006207229
Biography :ZHANGLei , Ph. D. , mainly engaged in agriculture biological environment . E2mail : zhanglei20902 @sina. com
摘要 :遮荫作为树冠对辐射截留的集中表现 ,与林下环境 ,特别是光环境特征具有密切关系。为满足坡地人工林下作物、中药栽
培环境模拟、预测和控制的理论需要 ,以圆锥形树冠为研究对象 ,借助几何光学中光的直线传播原理 ,根据水平投影与坡面投影
的亲似对应关系 ,在坡平面直角坐标系内讨论树冠的坡面投影形状、确定投影关键点坐标 ,在测树因子、经纬度、坡度、坡向、时
间 (包括时刻和日期)等变化要素基础上建立坡面投影理论边界模型 Z ( x′, y′) = 0。经模拟验证 ,各时刻实测与计算阴影边界
曲线相关系数大于等于 01928、显著水平 ( P)均达到 0101 ,表明该模型可较精确地描述正圆锥的坡面投影。在此基础上 ,进行模
型的实用检验和分析 ,选择北半球中纬度单株樟子松 ,计算得到 5 月 12 日 8 :00~14 :00 树木遮荫影响最远到达东 0169 倍树高、
西 1118 倍树高、南 012 倍树高、北 014 倍树高 ,绝大多数阴影位于近似扇环的区域内 ,计算结果符合北半球太阳视运动规律且实
际观测结果与计算值之间相关系数为 01983 ( P = 01008) ,证明该模型适用于树冠的遮荫分析 ,并可能在邻体遮光干扰的进一步
研究中发挥作用。最后讨论该模型的优点及局限性 ,明确模型在实际研究中合理运用的条件。
关键词 :坡面 ;圆锥形树冠 ;亲似对应 ;投影边界 ;模型
文章编号 :100020933(2006) 1023317207  中图分类号 :O185 ,Q141 ,Q948. 1 ,S181  文献标识码 :A
Theoretical research on a model for predicting the shadow boundary of an individual
conical crown on a slope
ZHANGLei , YU Hai2Ye 3 , YANG Hao2Yu , ZHANG Zhi2Ping  ( Key Laboratory for Terrain2Machine Bionics Engineering of MOE , School
of Biological and Agricultural Engineering , Jilin University , Changchun 130022 , China) . Acta Ecologica Sinica ,2006 ,26( 10) :3317~3323.
Abstract :Tree2shading could represent solar radiation intercepted by the tree crown. Understory light environment also had a close
relation with tree2shading. The conical crown was studied because it was familiar and its shadow , whether on a flat plane or slope ,
was more complex than other crown profiles. The shadow profile of an individual conical crown on a slope was an ellipse or a
closed curve composed of an ellipse and its two tangents , which depended on the relationship between the solar altitude ( hs ) and
the base angle ( A) of the cone. In the rectangular plane coordinate system ,a theoretical model was developed , which took tree
height ( H) , height below branch ( h) , crown radius( r) , longitude (λ) , latitude (φ) , height ( ht) , slope (α) , aspect (β) ,
date ( d) and time ( t) into account to describe the shadow boundary of an individual conical crown on the slope. The model was
based on geometrical optics theory and the corresponding and imitating principle between shadows on a flat plane and a slope. In
the simulation test , the results of Spearman and Kendall’s tau2b test showed all correlations between measured and calculated
values of boundaries were not less than 01928 and were significant at the 0. 01 level (22tailed) from 8 :00 to 14 :00. The model ,
therefore , was suitable to describe the shadow boundary of an individual cone , and the more the actual profile of the crown
resembled a cone , the more accurate the simulated results. Furthermore , this model was practically used to forecast the variable
scale of tree2shading of an individual tree species , Pinus sylvestris , in which the shadow of the crown always lied , from 8 :00 to
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
14 :00 on May 12th . Theoretically , this scale reached 0169H on the east side of the tree , 1118H on the west , 012H on the south
and 014H on the north. The result was in accordance with the apparent motion rules of the sun and the correlation was 01983
between observations and calculations , though it was different from existing research. There were two possible reasons for this
difference. One was that two scales were in different time spans respectively. The other reason was that this paper was based on
the apparent motion rules of the sun and did not consider any aspects of plant ecology. Advantages and limitations of the model are
listed which define how to use this model well with interrelated research.
Key words :slope ; conical crown ; the corresponding and imitating principle ; shadow boundary ; model
我国的人工林面积居世界首位 ,并以前所未有的速度递增 ;同时 ,我国是多山国家 ,多数人工林营建于坡
地[1 ] ,这些为开展人工林下作物、中药栽培提供了有利条件。
森林群落中 ,由于上层林冠的存在 ,直接导致林下的辐射量降低[2 ] ,而光生态因子又是林下植被发育、生
态系统作用过程关键的决定性因素 ,其空间变化远大于其他任何植物可利用资源[3 ] 。树冠遮荫作为树冠辐射
截留的集中表现 ,与林下包括太阳辐射在内的小气候特征具有密切关系[4 ] 。
随着生态学的发展 ,人们越来越倾向于用数学模型或生态模型来描述各种生态过程 ,其特点是关系简明 ,
精确性较高[5 ] 。20 世纪 80 年代以来已有一些关于树冠遮荫相关模型的报道[6~11 ] ,这些研究主要集中于单株
树冠投影长度的定量计算、投影形状和变化的定性讨论、单株树冠投影面积的近似几何运算及采用简单累加
或矩形法近似计算非单株树冠投影面积。所有这些研究均未考虑地形因素 ,也未考虑地面树木分布及阴影分
布的位置要素 ,缺乏对植物、空间和时间尺度多个影响因子的综合考虑 ,造成其在分析坡面非单株树冠投影连
续变化、邻体影响及环境模拟等方面的应用受到限制。
针对目前研究现状 ,本文着重在二维直角坐标系内 ,建立基于测树因子、地形、时间 (包括时刻和日期) 等
变化要素的单木树冠坡面投影边界模型 ,并对其模拟效果和应用进行了验证与说明 ,为坡面人工林下环境 ,尤
其是光环境动态变化的计算机模拟、预测和精确控制提供理论依据 ,同时也为邻体遮光干扰等研究提供一定
的理论支持。
1  模型的建立
111  理论基础
常见树冠基本形状包括[12 ] : ①圆锥形 ,多为壮年针叶树 ,如落叶松、黑松等。②圆球形 ,如槐树等。③圆
柱形 ,如黑杨等。3 种冠形中 ,以圆锥形树冠投影情况最为复杂 ,以下讨论以此为例。
11111  树冠的水平面投影  对任意圆锥形树冠 (如图 1 (a) ) ,树高 H ,枝下高 h ,半冠幅 r ,树冠锥顶 S ,圆锥底
面圆心 P ,圆锥底角 A ( A = arctg[ ( H - h)Πr ]) ; Ib 、Ic 分别为 A ≤hs < 90°和 0°< hs < A 时的入射光线 , hs 为水
平面太阳高度角 ;树冠剖面ΔSAD ( A 点位于阴面)与入射光线平行且过圆锥轴线 ,树冠剖面ΔSCE 与ΔSAD 垂
直相交于圆锥轴线 ; S F、SB 是光线照射下树冠阴、阳面交线 (即阴线) , G 为 B 、F 中点。当 A ≤hs < 90°时 , G、
B 、F 分别与 P、C、E 重合。
A ≤hs < 90°时 ,水平面投影是 8 点为圆心 ,半径等于 r 的圆 ;0°< hs < A 时 ,水平面投影是以 8 点为圆心 , r
为半径的扇形和过圆外点 1 与圆相切的两直线围成区域[13 ] 。图 1 (b) 、(c) 中 1、A′及 2~8 点分别是圆锥顶点
S 、圆锥底面 A~ G 及 P 点投影。1、4、7、8、A′点共线 ,且位于树木主干投影上 ,2、6 点是直线 l12 、l16与圆 8 切
点。其中 ,8 点为 4、A′中点 ,7 点为 2、6 中点。
11112  坡面投影与水平面投影关系  从光的直线传播原理出发 ,考虑到太阳光线入射的平行性 ,树木位于坡
面上时 ,假设光线射向坡面形成阴影 P′r′后 ,还可穿过坡面在水平地面上形成阴影 Pr ,则 Pr 与 P′r′应为一对
透视仿射对应 (亲似对应) ,即图 1 (d) 、图 1 (e)所示投影分别为图 1 (b) 、(c)的亲似对应。平面圆心 8 点投影为
椭圆圆心 8′点 ,圆的共轭直径投影后 ,共轭关系仍存在但垂直关系被破坏 ,平面投影上直线分点间线段比例、
线与线的平行及点与线的从属关系在坡面上仍然成立[14 ] 。
8133  生  态  学  报 26 卷
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
图 1  圆锥形树冠平面、坡面投影示意图
Fig. 1  Shadow of an individual conical crown on a flat plane and a slope respectively
(a)树冠结构参数 structural parameters of a crown ; (b) A ≤hs < 90°时的水平投影 shadow on a flat plane when A ≤hs < 90°; (c) 0°< hs < A 时的水平
投影 shadow on a flat plane when 0°< hs < A ; (d) (b) 对应的坡面投影 shadow corresponding to ( b) on a slope ; (e) (c) 对应的坡面投影 shadow
corresponding to (c) on a slope
图 2  坐标系的建立
Fig. 2  Rectangular coordinate system
112  直角坐标系的建立 (图 2)
对任意坡度为α的坡面 ,取树木 OS 所在 O 点为原
点 ,过 O 作坡面垂线OP ,令 OP 在水平面w 上落影所在
直线为 x 轴 ,落影反方向为正方向。x 轴逆时针旋转
90°得到 y 轴。xOy 为 w 上的直角坐标系。将水平面 w
绕 y 轴逆时针旋转α,得到与坡面重合的平面 w′(即坡
面 w′) 及该平面上的直角坐标系 x′Oy′。这里 , y 与 y′
轴重合 ,且为 w 与 w′上投影对应的亲似轴。
113  坡面 w′上任意时刻投影边界函数
11311  确定平面 w 上关键点坐标
图 3  树冠水平面投影示意图
Fig. 3  Shadow of a crown on a flat plane
(1) 参数介绍 (图 3)
①水平面太阳高度角 hs ∈[0°,90°] ;太阳方位角 As ∈[0°,360°] ,南为 0°,顺时针方向至东为 270°。hs 、As
均为地理纬度 (φ) 、经度 (λ) 、和时刻 ( t)的函数[14 ] 。这里 , t 取值范围的计算与海拔 ( ht)有关[14 ] 。
②坡向角β∈[0°,360°] ,南为 0°,顺时针方向至东
为 270°。需实测。
③A0 ∈[0°,360°] ,为 x 轴逆时针旋转至树木主干
投影的角度 ,可借助β、As 计算。
④ A′及 1~8 点是求投影边界的关键点。
(2) 0°< hs < A 时 ,点坐标求算 (图 3)
①1 点 x1 = Hctg hs cosA0 = u1 , y1 = Hctg hs sinA0 =
v1 ;
②4 点 x4 = ( hctg hs - r) cosA0 = u4 , y4 = ( hctg hs -
r) sinA0 = v4 ;
③8 点 x8 = hctg hs cosA0 = u8 , y8 = hctg hs sinA0 =
v8 ;
913310 期 张蕾  等 :坡面单木圆锥形树冠投影边界的模型理论  
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
④ A′点 xA′= ( hctg hs + r) cosA0 = uA′, yA′= ( hctg hs + r) sinA0 = vA′;
连接 O5 与 O3 , ∠1 O5 = ∠1 O3 = arctg[ r/ | O8| ] ∈[0°,90°] ,则
⑤3 点 x3 = ( r2 + h2 ctg2 hs ) 1Π2 cos ( A0 + ∠1 O3) = u3 ,
y3 = ( r2 + h2 ctg2 hs ) 1Π2 sin ( A0 + ∠1 O3) = v3 ;
⑥5 点 x5 = ( r2 + h2 ctg2 hs ) 1Π2 cos ( A0 - ∠1 O5) = u5 ,
y5 = ( r2 + h2 ctg2 hs ) 1Π2 sin ( A0 - ∠1 O5) = v5 ;
根据过圆外 1 点与圆 8 相切的条件 ,切点 2、6 坐标可通过 x1 、x8 、y1 、y8 、r 的函数计算 :
⑦2 点 x2 = f 2 ( u1 , u8 , v1 , v8 , r) = u2 , y2 = g2 ( u1 , u8 , v1 , v8 , r) = v2 ;
⑧6 点 x6 = f 6 ( u1 , u8 , v1 , v8 , r) = u6 , y6 = g6 ( u1 , u8 , v1 , v8 , r) = v6 ;
⑨7 点为 2、6 中点 x7 = ( u2 + u6 )Π2 = u7 , y7 = ( v2 + v6 )Π2 = v7 ;
(3) A ≤hs < 90°时 ,投影边界满足 ( x - u8 ) 2 + ( y - v8 ) 2 = r2 , x ∈[ u8 - r , u8 + r ]。
11312  确定坡面 w′上投影关键点坐标  坡面 w′上 1′~8′、A″点分别为 w 平面上对应点的投影。
一般情况下 ,基株 OS 的主干投影不会与 y ( y′)轴重合 ,即 ui ≠u′i′≠0 (如图 4) 。
图 4  u′i′≠0 时的坡面投影点与平面投影点关系示意图
Fig. 4  Relationship between the shadow boundary points on a plat plane and a slope when u′i′≠0
(1) i′( i = 1、4、7、8、A′)点坐标
①图 4 (a)中坡面太阳高度角 h′s ∈[0°,90°] ,是地理纬度 (φ) 、坡度 (α) 、坡向 (β)和时刻 ( t)的函数[15 ] 。
② i′点 : Di = | Oi | = ( u2i + v2i ) 1Π2 ,那么 , D′i′= | Oi′| = Di sin hsΠsin h′s
x′i′ = [ D′i′cos ( h′s - hs ) cosA0 ]Πcosα = u′i′, y′i′ = D′i′cos( h′s - hs ) sinA0 = v′i′
(2) j′( j = 2、6、3、5)点坐标
2′与 6′、3′与 5′点分别对称 ,对称点 7′、8′已知 ,则 :
直线 l16与 y 轴交于点 O6 (0 , vO6 ) ,坡面上 1′、6′亦与 O6 共线 ,6′点满足于直线 l1′O6 ;同理 ,2′点满足于直线
l1′O2 ;再由对称条件 ,可建立 4 个方程构成的线性方程组。据此 , j′( u′j′、v′j′) ,均可求得。
对于主干投影在 y ( y′) 轴上 ,即 u′i′= 0 的特殊情况 , i′点与 i 点重合 ;同时 ,由几何关系不难推出 , u′j′=
ujΠcosα, v′j′= vj [ ( h - uj tgα)Πh ] , j = 2、3、5、6。推导过程略。
由于圆 8 上所有点均可求 ,而椭圆 8′上点与其一一对应 ,利用 11312 (2)方法可求 j′以外的任意其他点 ,那
么在任何情况下均足以满足求椭圆一般方程的需要。
11313  边界模型 Z ( x′, y′)  在坡面直角坐标系 x′Oy′内 ,借助椭圆上 5 个互不重合的点 ,解线性方程组 ,得
圆心为 8′的椭圆 F ( x′, y′) = a1 x′2 + a2 x′y′+ a3 y′2 + a4 x′+ a5 y′+ 1 = 0 ;借助直线上已求得的两点 ,确定直
线 y′= l1′2′( x′)与 y′= l1′6′( x′) 相交函数为 G( x′, y′) = [ l1′2′( x′) - y′]·[ l1′6′( x′) - y′] = 0 ,那么阴影边界满
足 :
0233  生  态  学  报 26 卷
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
(1) A ≤hs < 90°时 , Z ( x′, y′) = F ( x′, y′) = 0 ,定义域 DZ = 椭圆定义域 DF 。
(2) 0 < hs < A 时 , Z ( x′, y′) 是由椭圆 F ( x′, y′) 与直线相交函数 G( x′, y′) 构成的封闭曲线 ,通过直线 y′=
l2′6′( x′)进行如下分割 :
① u′2′≠u′6′且 v′1′< 0 ,椭圆 F ( x′, y′)在直线 y′= l2′6′( x′)上方
Z ( x′, y′) =
F( x′, y′) = 0 , y′> l2′6′( x′)
G( x′, y′) = 0 , y′≤ l2′6′( x′)
② u′2′≠u′6′且 v′1′> 0 ,椭圆 F ( x′, y′)在直线 y′= l2′6′( x′)下方
Z ( x′, y′) =
F( x′, y′) = 0 , y′> l2′6′( x′)
G( x′, y′) = 0 , y′≤ l2′6′( x′)
③ u′2′= u′6′且 u′1′> 0 ,椭圆 F ( x′, y′)在直线 y′= l2′6′( x′)左侧
Z ( x′, y′) =
F ( x′, y′) = 0 , x′< u′2′
G( x′, y′) = 0 , x′≥ u′2′
④ u′2′= u′6′且 u′1′< 0 ,椭圆 F ( x′, y′)在直线 y′= l2′6′( x′)右侧
Z ( x′, y′) =
F ( x′, y′) = 0 , x′> u′2′
G( x′, y′) = 0 , x′≤ u′2′
  定义域 DZ 由 DF 和 1′点决定。
图 5  模拟圆锥树冠实测与计算投影比较
Fig. 5  Comparison of measured and calculated values of boundaries
2  模型的验证与分析
211  边界模型的模拟验证与分析
由于实际树木阴影面积较大 ,特别是形体高大的树木 ,而目前尚无准确测量树木阴影面积的合适方法 ,手
工描绘投影边界、测量投影均会产生较大的测量误差 ,用这些实测值来验证理论计算结果 ,很难区分模拟效果
的好坏是模型表达的问题还是测量误差造成的。所以本文以小型正圆锥体模拟树冠进行模型的验证。
任意选定 :2005 年 12 月 15 日 ,以母线长 518cm 的正圆锥体模拟 r = 219cm 的树冠 , H = 9cm , h = 4cm。实
验地点位于东经 125°21′,北纬 43°52′,坡向东偏南 22°
(即β= 338°) 、α= 18°, ht = 240m。当日日照时间约为
7 :00~16 : 00 ,实验时间选择日照强度较大的 8 : 00~
14 :00 ,每隔 1h 在坐标纸上绘出阴影边界 ,与计算机数
值计算的理论值比较 (如图 5) 。
边界按步长 015 对实测与计算曲线插值 ,应用
SPSS1110 进行 Spearman 及 Kendall’s tau2b 相关性检验 ,
结果一致 :各时段两曲线插值结果相关系数均大于等于
01928 ( P ≤0101) ,表明该模型可较精确的模拟正圆锥体
的坡面投影。由于实验以正圆锥模拟树冠 ,边界差异主
要来自计算误差及太阳辐射半影效应对实测边界描绘
的影响。这一理由可解释图 5 中椭圆部分模拟效果较
直线部分差且 8 : 00 与 14 : 00 的模拟差异相对较大的
问题。
由于模拟实验的坡度、坡向为任意选定 ,虽无法反
映所有可能情况 ,但具备一定的随机性。对于实际地
形、实际树冠 ,投影规律与模拟圆锥一致 ,只是树冠锥体
体积更大。所以 ,该模型对实际树冠坡面投影模拟效果
123310 期 张蕾  等 :坡面单木圆锥形树冠投影边界的模型理论  
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
的精确与否主要取决于树冠与正圆锥体的接近程度。
212  模型的实用验证
Z ( x′, y′) = 0 最直接的应用是通过积分进行遮荫面积计算 ,但除此之外还可根据其在坐标系内位置的变
化 ,确定树木遮荫影响范围等。
图 6  8 :00~14 :00 坡面树冠投影边界散点图
Fig. 6  Scatter plot of shadow boundaries from 8 :00 to 14 :00 on a slope
长春南湖公园 (东经 125°18′,北纬 43°51′) 内 ,坡度
α= 10°的向南坡面 (β = 0°) 上 12a 生樟子松 ( Pinus
sylvestris) , H = 3155m、h = 1129m、r = 1m。按本文方法计
算 ,单株树木在 2006 年 5 月 12 日 8 :00~14 :00 间树冠
投影最远到达东 0169 H (约 2145m) 、西 1118 H (约
412m) 、南 012 H (约 0171m) 、北 0145 H (约 116m) 。从这
些数值上看 ,绝大多数阴影位于 8 :00~14 :00 近似扇环
的理论投影边界区域内 (如图 6) 。这一范围虽与李树
人等[7 ]以基株为中心 ,东西 2 H 为长轴、南北 115 H 为短
轴的近椭圆范围存在差异 ,但本文结果与北半球太阳视
运动规律相符 ,且当日实际树冠投影观测结果为阴影最
东到达 2111m、西 4m、南 0152m、北 1174m ,与理论值之
间相关系数为 01983 ( P = 01008) ,即计算与实际基本符
合。本结果与已有研究存在差异的可能原因有两个 : ①已有椭圆遮光范围是针对一年中每天 7 :00~17 :00 均
适用而提出的 ,时间跨度大于本文 ,所以南、北、东、西投影均较长 ; ②已有研究从生态场角度出发 ,考虑植物生
态场作用特点和特征 ,而本文仅从太阳视运动规律出发 ,单纯反映树冠随太阳运动产生坡面投影实际位置的
变化。
边界模型对遮荫影响范围的预测可在邻体遮光干扰、林下辐射分布等研究中发挥作用。
3  讨论与结论
综上 ,与其他现有遮荫模型相比 , Z ( x′, y′) = 0 具有以下特点 :
(1) 灵活性强  该模型考虑测树因子 ( H、h、r) 、地理位置特征 (φ、λ、ht、α、β) 及时间 (包括时刻 t 及日期
d)等几乎所有对树冠投影有重要影响的变化因素 ,每输入一组变量值 ,都能求出与之对应的确定的边界表达
式。虽然多变量也同时带来计算复杂、工作量大的问题 ,不过 ,在实际工作中 ,基于研究目的的不同 ,函数中包
含的变量总有一些可以取为定值 ,甚至只剩一个变量。如研究选定地点阴影日变化 ,则 H、h、r、φ、λ、ht、α、β
及 d 均为定值 ,仅 t 是变量。确定的研究内容使变量数目减少、计算简单、模型表达具体化。
(2) 实用性强  Z ( x′, y′) = 0 弥补了已有相关模型未考虑地形因素和树木位置因素的不足 ,无论是从我
国森林多分布于山地的实际出发 ,还是从环境动态模拟、量化的角度来讲都更具实用性 :除了可通过积分求得
阴影面积外 ,还可通过坐标平移得出林内任意与基株相对位置已知的树木投影边界 ,进而对林下阴影区、透光
区、阴影交叠区域进行划分 ;可对林内任意选定区域的未来遮荫面积、邻体遮荫影响等进行理论计算 ;有了确
定的坐标范围 ,还可将树冠遮荫的动态与林下辐射、温度、湿度变化乃至光、热传输联系起来 ,对确定位置的太
阳辐射等环境因素进行定量 ,最终实现林下环境的模拟和预测。
  (3) 局限性
①该模型对阴影模拟效果的好坏主要取决于树冠与正圆锥体的相似程度。
②模型在同龄人工纯林内应用更广泛。天然林和异龄林林分组成多样化 ,冠层结构存在较大的异质性 ,
遮荫情况复杂 ,造成该模型的适用有限。
③由于日出后、日落前半影效应明显 ,模型在这段时期对实际阴影的模拟不可避免的存在较大误差 ,但
这段时间内太阳辐射强度较小 ,植物的光合作用较弱 ,对于林下植物生长环境研究的意义不大 ,所以并不影响
2233  生  态  学  报 26 卷
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net
模型在林下植物生长环境研究中的实际应用。
了解模型的优点及局限性 ,灵活运用该模型进行合理的科学计算 ,可为坡面人工林下作物、中药栽培环境
特别是光环境动态变化研究提供重要的理论依据。同时 ,本文对单株树冠投影边界函数的推导及结果也可作
为进一步开展天然林、异龄林遮荫动态变化研究的基础。
References :
[ 1 ]  Zhang X Q , Xu D Y. A canopy radiation model for even2aged plantations  Ⅰ. theoretical calculations. J . of Acta Phytoecologica Sinica , 2001 , 25 (5) :
609~615.
[ 2 ]  Bellow J G, Nair P K R. Comparing common methods for assessing understory light availability in shaded2parennial agroforetry systems. Agriculture and
Forest Meteorology , 2003 , 114 (324) : 197~211.
[ 3 ]  Bazzaz F A. Plants in changing environments : linking physiological , population , and community ecology. Cambrige : Cambridge University Press , 1996.
[ 4 ]  Richard Lee. Forest microclimatology. New York : Columbia University Press , 1978.
[ 5 ]  Wang D L. Some comments on ecological field. Chinese Journal of Applied Ecology , 2000 ,11 (3) : 472~476.
[ 6 ]  Wu L L , Wang Z C. A study on tree shadow movements. Journal of Nanjing Forestry University , 1991 , 15 (2) : 61~66.
[ 7 ]  Li S R , Zhao Y. Study on the shading mathematic model of tree crown. Acta Agriculture University Henanensis , 1994 , 28 (4) : 361~366.
[ 8 ]  Wu L L , Wang Z C. Tree screen effects on space. Journal of Nanjing Forestry University , 2000 , 24(supplement) : 47~49.
[ 9 ]  Li S R , ZhaoY, Yan Z P. Light2ecological field of Larix leoptolois canopy . Chinese Journal of Applied Ecology , 1997 , 8 (2) : 123~126.
[10 ]  Shen Q , Tang R N. A theoretical study on the number of trees in an interplanted tea stand. J . Biomath. , 1994 , 9 (5) : 66~70.
[11 ]  Lu G Z , Bu Z D , Tian F J , et al . Study on shading of apricot intercropping with agricultural crops on the terrace land in the semiarid area of the western
part of Liao ning. Protection Forest Science and Technology , 2005 , 1 (64) : 18~20.
[12 ]  Zheng W J . Arbores sinicae. Beijing : China Forestry Press , 1985.
[13 ]  Xu Z Q , Wu H X. Foundations of projection geometry. Beijing : Higher Education Press , 1987. 1~15.
[14 ]  Yao L H. Agriculture meteorology. Beijing : China Forestry Press , 1992. 198~199.
[15 ]  Bao F P. Mountain climates. Beijing : Science Press , 1983. 4.
[16 ]  Ma Q Y, Liu Z G, Pan XL , et al . A study on the canopy structure and light distribution of Larix principis2ruppurechtii during the growing season. Journal
of Beijing Forestry University , 2000 , 22 (4) : 18~21.
[17 ]  Eshetu , Yirdaw , Olavi Luukkanen. Photosynthetically active radiation transmittance of forest plantation canopies in the Ethiopian highlands . Forest Ecology
and Management , 2004 , 188 (123) : 17~24.
参考文献 :
[ 1 ]  张小全 , 徐德应. 人工同龄林冠层辐射场模拟模型 Ⅰ. 理论计算. 植物生态学报 , 2001 , 25 (5) : 609~615.
[ 5 ]  王德利. 关于生态场的几点评述. 应用生态学报 , 2000 , 11 (3) : 472~476.
[ 6 ]  吴力立 , 王宗淳. 树冠遮阴动态研究. 南京林业大学学报 , 1991 , 15 (2) : 61~66.
[ 7 ]  李树人 , 赵勇. 树冠遮光数学模型的研究. 河南农业大学学报 , 1994 , 28 (4) : 361~366.
[ 8 ]  吴力立 , 王宗淳. 林木对空间的遮蔽效应. 南京林业大学学报 , 2000 , 24 , 增刊 : 47~49.
[ 9 ]  李树人 , 赵勇 , 阎志平. 日本落叶松林冠层光生态场研究. 应用生态学报 , 1997 , 8 (2) : 123~126.
[10 ]  沈泉 , 唐荣南. 林茶间作的理论基础研究. 生物数学学报 , 1994 , 9 (5) : 66~70.
[11 ]  卢国珍 , 步兆东 , 田福军 ,等. 辽西半干旱区梯田地杏农复合系统树木遮荫作用规律研究. 防护林科技 , 2005 , 1 (64) : 18~20.
[12 ]  郑万钧主编. 中国树木志. 北京 : 中国林业出版社 , 1983.
[13 ]  许志群 , 吴海霞编. 射影几何基础. 北京 : 高等教育出版社 , 1987. 1~15.
[14 ]  姚丽华主编. 农业气象学. 北京 : 中国林业出版社 , 1992. 198~199.
[15 ]  傅抱璞著. 山地气候学. 北京 : 科学出版社 , 1983. 4.
[16 ]  马钦彦 , 刘志刚 , 潘向丽 ,等. 华北落叶人工林生长季内的林冠结构和光分布. 北京林业大学学报 , 2000 , 22 (4) : 18~21.
323310 期 张蕾  等 :坡面单木圆锥形树冠投影边界的模型理论  
© 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. http://www.cnki.net