213 Appendix A. ongwave radiation in a canopy stand The longwave radiation transfer equations can be written as (Ross, 1981): di ld (,r) d cosθ = G(,r)i ld (,r)+η,r ( ) (A.1) di lu ( ),r d cosφ = G,r ( )i lu (,r) η(,r) (A.2) where i ld (,r), i lu (,r) are the longwave downward and upward radiation intensities inside a plant stand, r is the direction in space, θ is the inclination angle and φ is equal to (π θ ), G(,r) is the foliage area orientation function, η(,r) is the leaf emission coefficient. The boundary conditions are given as: i ld ( ) = ε A i lu ( )= ε B + A BJ π Ω 1 i d (,r')cosθ'dω' (A.3) where ε A is the intensity of longwave radiation of the atmosphere, ε B is that emitted from the ground surface below the stand, A BJ is the albedo of the ground surface for longwave radiation. By assuming the leaf temperature is the same in a horizontal layer, one can obtain the following expressions: R ld π /2 ( )= 2π ε A ()exp θ G() θ cosθ sinθ dθ + cosθ R lc π /2 ( ' ) 2 G()exp θ ( ' ) G() θ sinθ dθ cosθ d' (A.4)
214 R lu π /2 ( )= 2R lg exp ( φ cosφsinφdφ + cosφ π /2 R ( ' ) 2 G()exp φ ( φ cosφ sinφdφ lc cosφ d' + π /2 ( ) exp ( φ 2π 1 δ B cosφ sinφdφ cosφ π /2 G() θ 2π ε A ()exp θ cosθ sinθ dθ + cosθ π /2 2π ( 1 δ B ) exp ( φ cosφsinφ dφ cosφ E π /2 ( ) 2 G()exp θ ( θ sinθdθ cosθ d' (A.5) where R ld () and R lu () are the longwave downward and upward irradiances. The first term in equation (A.4) defines that part of longave radiation of the atmosphere which penetrates the layer (, ) of the stand and reaches the level, the second term defines the downward flux of the longwave radiation of leaves in the layer (, ), which reaches the level. In equation (A.5), the first term defines that part of longwave radiation of the ground which penetrates the layer (, ) and reaches the level, the second term describes the upward flux of longwave radiation of leaves in the layer (, ) which reaches the level, and the third term defines that part of longwave radiation of the atmosphere which penetrates the whole layer of the stand, is reflected from the ground and, penetrating through the layer (, ), reaches the level. The last term represents that part of the downward flux of radiation of the leaves off the whole stand which reaches the ground, is reflected from it, penetrates the layer (, ) and reaches the level. Since ( 1 δ B ).5, the last two terms in expression (A.5) can be neglected as a first approximation and in the case of isotropic longwave radiation of the atmosphere, we can obtain the following expressions from equations (A.4) and (A.5): R ld ( )= R lc ( R lc R la )a D ( ) (A.6) R lu ( )= R lc ( R lc R lg )a D ( ) (A.7) where R lc is the longwave radiation emitted from the stand, R la is the longwave radiation from the atmosphere, R lc is the longwave radiation from ground surface, a D () is the penetration function π /2 a D () = 2 exp( G(θ) cosθ )cosθ sinθ dθ (A.8)
The net longwave radiation can be written as: R ln ( )= ( R lc R la )a D 215 ( )a D () R lg R lc For relatively dense stand, R R R = R and one obtains: lg lc, lu lc ( ) (A.9) R ln = ( R la R lc )a D ( ) (A.1) The penetration function for longwave radiation can be defined as: π /2 a D () = 2 exp( G(θ) cosθ )cosθ sinθ dθ (A.11) where is the leaf area index, G(θ) is the so-called G-function which is a complicated double integral of inclination (θ) and azimuth orientation. This equation can be equally applied to diffuse radiation in a plant canopy. The penetration function for shortwave radiation in a canopy with uniform orientation (spherical orientation) can be expressed as: a s ()= exp( 2sinβ ) (A.12) where β is the solar elevation. It is clear that both penetration functions are exponential. Data from Ross (1981, Table II.5.1) presented in Figure A1 shows that both shortwave and longwave radiation can be approximated by the same function (Beer s law) with the same parameter value (bulk extinction coefficient). 1 Penetration function.8.6.4 isotropic overcast 3deg 6deg 9deg maize sunflower k=.66 k=.7.2 1 2 3 4 5 AI Figure A1. Penetration function for shortwave and longwave radiation in a spherical canopy under various sky conditions.
216 Appendix B: Details of Analytical Solutions We here formulate u and u' (i.e. u/ ζ) for the analytical solutions presented in section 4.9.4, equations (4.13) and (4.14). The most general formulation (Broadbridge, 199) has not been given succinctly, or in a form that permits exact evaluation. We express u as: u = u +u s1 + u s2 (B.1) where u is required regardless of initial and lower boundary conditions, u s1 and u s2 are non-zero summations when there is a finite-depth profile, and u s2 is non-zero only for a finite-depth profile with non-zero initial soil-water content. To represent the effect of non-zero initial soil-water content Θ i, we define: Q = ρ +.5 Θ i 1 Θ (B.2) i With a completely dry initial condition Θ i =, equation (B.2) yields Q = ρ. Expressions for u and u' /u u =.5exp( ζ 2 / τ)a ( + b c + d e) (B.3) where: a = 2exp( ζ Qτ )/ τ ) 2 (B.4) b = f c = f ( ζ gτ)/ τ ) (B.5) ( ζ Qτ )/ τ) (B.6) d = f( ζ + gτ)/ τ ) (B.7) e = f ( ζ +Qτ)/ τ) (B.8) g = ρ 1+ 1 ρ (B.9) f(x)= exp( x 2 )erfc(x) (B.1)
217 Differentiating (B.3) with respect to ζ yields: u' = exp( ζ 2 /τ)qa+ ( c e) gb ( d) ( ) (B.11) Equation (4.13) requires u' /u. Because b, in (B.5), is some orders of magnitude larger than the other terms in (B.3) and (B.11), Θ-values computed directly from these equations may slightly exceed the theoretical maximum at the trailing edge of the wetting front, using 64-bit arithmetic. By taking the quotient of (B.3) and (B.11) and rearranging the algebra, b may be eliminated from the numerator, and the numerical error is minimised. The resulting expression is: u' /u = 2 Qa+ ( c d )+ gd 2g 1 a + b c + d e a c + d e a + b c + d e (B.12) Equations (B.3) and (B.12) can be used to solve (4.13) and (4.14). Expressions for summation terms of u and u' 4 ( ) u s1 = 1 S j exp( c 1 )fc ( 2 ) exp( c 3 )f( c 4 ) (B.13) 2 n= 1 j = 1 3 ( ) u s2 = ω j exp( c 5 )fc ( 6 ) exp( c 7 )f( c 8 ) (B.14) n= j = 1 where: c 1 = 4n 2 ( 1 Θ i )l ρ 2n( 1 Θ i )l + ζ ( ) 2 /τ (B.15) c 2 = ( 2n( 1 Θ i )l +ζ)/ τ ( 2nρ + b j ) τ (B.16) c 3 = 4n 2 ( 1 Θ i )l ρ 2n( 1 Θ i )l ζ ( ) 2 /τ (B.17) c 4 = ( 2n( 1 Θ i )l ζ)/ τ ( 2nρ b j ) τ (B.18) c 5 = 4n 2 ( + 4n)1 ( Θ i )l ρ ( 2n + 1) ( 1 Θ i )l + ζ c 6 = ( 2n + 1)1 ( Θ i )l +ζ)/ τ 2nρ + a j ( ) 2 / τ (B.19) ( ) τ (B.2) c 7 = 4n 2 ( + 4n)1 ( Θ i )l ρ ( 2n + 1) ( 1 Θ i )l ζ ( ) 2 /τ (B.21)
218 c 8 = ( 2n + 1)1 ( Θ i )l ζ)/ τ 2nρ + a j ( ) τ (B.22) [ ] (B.23) ω j = exp( 2( 1 Θ i )l Q), ω 1 /2, ω 1 /2 S j = [ +1,+1, 1, 1] (B.24) a j = [ ρ,2ρ, Q,Q] (B.25) [ ] (B.26) b j = ρρ+ ( 1), ρρ+ ( 1),Q, Q The derivatives of the summation terms with respect to ζ are: 4 u' s1 = S j ( 2nρ + b j )exp ( ()f c 1 ( c 2 )+ exp( c 3 )fc ( 4 ) (B.27) n= 1 j =1 3 u' s2 = 2ω j ( 2nρ + a j )exp ( ( c 5 )fc ( 6 )+ exp( c 7 )fc ( 8 ) (B.28) n = j = 1 For a finite-depth profile, a test of computational roundoff error at the critical lower boundary may be made by checking the deviation of z (ζ max ) from depth l, where ζ max = (1 Θ i )l ρτ represents the lower boundary at any point in time. Analytic solutions are valid in principle for ζ max >. However this test is required because summation terms in the expressions for u/ ζ and u may be too large for accurate computation, even with 64-bit arithmetic, because the greatest number that can be represented is 1 ±38. This problem occurs if ρ is too large in relation to l and soil-water content is large at the lower boundary. It cannot be solved simply by rearranging the algebra (Broadbridge, 199). The test provides a good direct check on the accuracy of z as computed by (4.14), and an indirect check on the accuracy of Θ as computed by (4.13). For a semi-infinite profile, no fixed value of ζ max will correspond to an arbitrary scaled depth l. For computation to that depth, ζ max lies in the range of (1 Θ i )l ρτ ζ max l. In general, the ζ-range must be searched to obtain a solution for a given value of z.