Abstract

This paper presents the Growth-Flow method for animating the effect of wind on the motion and growth of virtual plant branches and leaves. The method incorporates changes to the growth rate when a plant is exposed to winds with speeds higher than a threshold. In particular, growth rate is reduced in branch elongation, increased in the branch radius, reduced in leaf length, and increased in leaf thickness. In addition, when a plant is exposed to wind for long time periods, the branch growth angle is changed to align more closely with the wind vector. The Growth-Flow method incorporates all these effects on growth and motion due to wind in one algorithm.

1. Introduction

To enhance the realism of virtual plants, the effects of wind on the motion and growth of plants should be modeled. The effects of wind blowing through trees and especially its effects on branches have been extensively studied in computer animation. To date, only a few computer animation researchers have studied the effect of wind on leaf motion, with no consideration of the effects of wind on the growth of virtual plants.

Plants are affected by wind immediately through physical motion and over time through changes in growth. Wind contact with a plant immediately causes physical motion by bending, swaying, or twisting the stem or leaves. Studies on branch movement in trees have approached the issue using physically based methods [14], procedurally based methods [58], motion capture [9], simulation [10], and shell dynamics [11]. Studies have examined leaf motion with respect to flow [12], biomechanics [13], spring and rotational movement [14], wilting [15], response to wind by broadleaf types [16], response to wind by both broadleaf and needle leaves [14], as a component of overall plant motion [10], and as part of shell space [11]. In addition, animation with leaf growth effects has been performed using a physically based approach without any consideration of wind [17]. Over time, plant growth is affected by wind, according to its strength and duration. Plants respond in three major ways: the girth of the stem is increased, the length of the stem is decreased, and the direction of growth is changed. The leaves of the plant also change in two ways: their thickness is increased and their shape becomes closer to round and stunted [18, 19].

This paper proposes the Growth-Flow method for animating the effect of wind on the motion and growth of stems and leaves of virtual plants. Incorporating the effects of wind on growth into a virtual plant that supports motion is an original contribution of this paper. One might assume that plant growth is too slow a process to be considered jointly with plant motion. However, in many uses of virtual plants, the growth is shown at such an accelerated rate that it is reasonable to combine the two. For example, a fast-growing plant may provide a challenge in a computer game, and it will appear more authentic if it also moves with the wind and its growth is affected by the wind. Moreover, a lecturer explaining the effect of wind on long-term growth may want to use an accelerated model to illustrate the interrelationships.

When animating plants, a technique illustrating the combined effect of wind on motion and growth is generally not needed, since only motion is typically of interest. However, when depicting the effect of wind on the growth of plants in an animated scenario, it is likely that one will also want to depict motion. Our experience is that when the wind is acknowledged to have an effect on growth, users want to see evidence of the wind, that is, simultaneous motion due to the wind. The Growth. Flow technique described in this paper satisfies this need.

The remainder of this paper is organized as follows. Background information on wind, plant motion, plant growth, and the responses of plants to wind is given in Section 2. Formulae and algorithms for the motion and growth of both stems and leaves in response to wind are given in Section 3. Experimental results are presented in Section 4, a discussion related to current literature is given in Section 5, and conclusions are drawn in Section 6.

2. Background

This section describes the motion and growth effects caused by wind. Section 2.1 provides background information concerning wind and the basis of its animation. Section 2.2 provides the background for the effect of wind on plants. Section 2.2.1 describes branch movement, Section 2.2.2 describes branch growth, Section 2.2.3 describes leaf movement, and Section 2.2.4 describes leaf growth.

2.1. Wind Motion

Wind is air in motion. It is produced when atmospheric (air) pressure is higher in one area and lower in another. The velocity of the wind results from the pressure difference between an area of high pressure and an area of low pressure, with higher velocities from greater differences. Wind velocity is characterized by its speed and direction [20]. Existing methods simulate wind by distinguishing between gusts and stable winds [14, 21] and by treating it as a fluid flow modeled by the Navier-Stokes equations [12, 2224]. We assume that , the velocity of the wind at time everywhere in the simulation, is set according to a wind model. Throughout this paper, all variables shown in bold face represent vectors. Gusts could be modeled in the wind model by rapidly increasing and decreasing the wind speed in particular areas from time to time. Instead, for simplicity, we give leaf movements a random component that can be interpreted by the viewer as wind gusts. This random component is used to decide which of the bending, twisting, and vibrating motions to apply to each leaf.

2.2. Plant Response to Wind

As mentioned in Section 1, plants respond to wind immediately with motion and over time with changes in growth patterns. Plant motion in response to wind has been studied most thoroughly with respect to trees. The response of a tree branch or leaf to wind depends upon its tensile strength, the size of its surface area facing the wind, and its mass [25, 26]. If a branch has a large surface area facing the wind, a low-tensile strength, and a large mass, it can fracture when the force of the wind exceeds the strength of the branch. Leaves respond in a similar fashion to branches, although with less influence from the leaf’s mass, which in negligible, and more from the surface area.

Branch movement in response to wind is observed as a swaying motion, which is also evident in leaf movement. In addition, leaves display a twisting motion in response to the wind and a vibration motion when the wind is tangent to the blade [8].

Plant growth is affected by wind, among many other environmental factors. Several studies have been conducted to discover the response of various plant species to wind. One study observed that elongation was significantly retarded in young plants of Hordeum vulgare, Bryonia dioica, Cucumis sativus, Phaseolus vulgaris, Mimosa pudica, and Ricinus communis [27]. Another observed that leaf growth was inhibited by 30% in cocklebur plants and grape vines [28]. In addition, when trees of Liquidambar Styraciflua were shaken daily, elongation was reduced by 70%–80% [29]. Recently, a study regarding airflow and plant growth discovered that increased airflow enhanced stem elongation when the plant did not flex its stem and restricted elongation otherwise [30]. Another recent study on wind and mechanical stress on leaves observed that wind restricted growth and also made leaves thicker and closer to being round [18]. As mentioned in Section 1, the overall effect of wind on stems (or branches in trees) is to thicken them, restrict their elongation, and change the direction of their growth and the overall effect on leaves is to thicken them and change their shape closer to being round [18, 19].

The remainder of this section describes previous research on plant motion in more detail by considering branch movement, branch growth, leaf movement, and leaf growth.

2.2.1. Branch Movement

Branches are typically modeled using a segment-based approach, exemplified by the work of Sakaguchi and Ohya [1], in which each segment represents an approximately linear part of a branch between nodes. To avoid confusion, in this paper, we reserve the term “node” for its sense in botany, where it refers to a point on a plant where new growth is possible. Typically, a node corresponds to a joint in an animated model. A branch is represented by a computer science tree data structure where the vertices represent segments and the arcs represent the existence of (botany) nodes. The motion of a particular segment depends on the motions of its child segments (i.e., connected segments closer to the external end of the branch) and the effect of the wind on the particular segment. The effect of the wind on a particular segment is calculated based on the combination of wind velocity, axial damping, and restoration force. The axial damping is influenced by the tensile strength of the branch, which depends on its thickness. The restoration force (resistivity) is the force influencing a branch that is not in its rest position to return to that position. Equation (1) can be used to calculate motion of a segment [1]: where is the moment of force (torque), is the external force (i.e., wind), is the resistivity, is the axial damping force at a node, is the force that propagates back from the child branches (if any), indicates the cross product, and is a vector from the node at the beginning of the segment (where it attaches to its parent segment) to the center of gravity of the segment [1]. The branch movement method described by Sakaguchi and Ohya [1] is illustrated in Figure 1. On the left, three segments joined at two nodes are shown, and on the right, the same segments are shown after the wind has moved the top segments to the right. The segment of interest is the middle one. The forces listed in (1) act in the directions shown.

2.2.2. Branch Growth

Growth is an irreversible increase in size, mass, or number of cells [31]. Ignoring other environmental factors, growth for a living thing, such as a plant, shows no logarithmic, linear, or exponential increase with respect to time, depending on the current time during the year. Three existing ways of measuring growth are (1) absolute growth rate (AGR), which is the total gain in size (height, weight, and radius) or mass (dry or fresh weight) over time, (2) relative growth rate (RGR), which is the change in the logarithm of mass over time, and (3) percentage growth rate (PGR), which is the percentage change in length or mass over time. AGR and RGR are used for measuring growth in plant physiology, whereas PGR is used in horticulture. For illustration AGR, RGR, and PGR can be characterized as shown in (2), (3), and (4), respectively, where and are the relevant size or mass measurements and and , where , are the respective times at which these measurements were taken.

Branch growth is seasonal. When it occurs, a sigmoidal growth curve is observed such that the growth rate slowly increases during the early part of the growing season, remains constant (linear growth) during the main part of the growing season, and then declines so that it reaches zero when the segment achieves its maximum size [31]. Stem elongation proceeds differently among plants: some grow at the apical meristem (the branch tip) resulting in long older growth, some exhibit varying rates of elongation over several segments, and some grow evenly along all stem segments. For simplicity, we assume the later. Given this assumption, the principal differentiating characteristic is the growth rate during the linear-growth period. Absolute and percentage growth rates for various plants are provided in Table 1 in the exact form given in the listed references, along with wind inhibition percentages, which will be discussed shortly. We are presenting a scattering of examples because apparently no large-scale-comparative study has been reported.

Branch growth is affected by wind through a thickening of the stem [3236] and reduced elongation [3739]. The direction of the branch is altered to align more closely with the wind vector if the wind is constant and strong over the lifetime of the plant. In general, the effect of wind on growth is not significant unless the wind persists over long periods and has sufficient speed, for example, over 4 meters per second for aspen [19]. The amount of effect on growth is correlated with wind speed [19]. Table 1 also includes the observed effect of wind on growth rate for the sample plants where available.

Another common growth formula is based on the Richards growth model [4749] which considered exponential growth as described in where is time, is the size of the growing part (e.g., amount of dry matter) at time , is the asymptotic value of the size, and , , and are constant coefficients, where determines the curve type, describes the linear component of growth, and determines the rate at which growth changes over time [47]. According to Richards, “ is the “rate constant” which determines the spread of the curve along the time axis, is usually unimportant biologically, since it reflects only the choice of the zero of the time” ([47], p. 290). (This should not to be confused with the vector used elsewhere in this paper.) The constant influences the shape of the curve. The formula used within the Growth-Flow method resembles this growth formula.

2.2.3. Leaf Movement

A leaf is attached to a petiole, which is in turn attached to a branch, as illustrated in Figure 2. A petiole is a smaller branch structure that ensures proper attachment of the leaf to the branch. The movement of the petiole influences the movement of the leaf and vice versa. The tensile strength of the petiole influences how far a leaf bends and whether the petiole breaks off its branch. The surface area of the leaf blade is also a contributing factor to the motion of the leaf in response to wind, with a larger surface area resulting in a stronger response. The tensile strength of the blade, which depends on the strength of the vein, also influences the movement of the leaf. If a leaf has been exposed to drought and its tensile strength is low, its blade can bend more in response to the wind than it would otherwise. All the previously discussed factors, that is, wind speed, tensile strength of petiole, tensile strength of blade, and leaf blade surface area, contribute to the overall leaf motion. The first three of these factors are considered in more details in Section 3.

2.2.4. Leaf Growth

Steady wind of sufficient speed has been observed to affect the growth of a leaf through an increase in the thickness of the leaf blade. For aspen, this effect was observed for wind with a speed of more than 4 meters per second [19]. In addition, leaves tend to become smaller in surface area [19] and closer to being round [18] due to wind stresses on their surfaces.

3. The Growth-Flow Method

The Growth-Flow method has two components resulting from wind influence: physical motion and growth effects on branches and leaves. The methods used to simulate the physical motion and growth for a branch are described in Section 3.1, and those for a leaf are described in Section 3.2.

3.1. Branch Response to Wind

To simulate the physical motion and growth for a branch, we consider the movement in Section 3.1.1 and the growth in Section 3.1.2.

3.1.1. Branch Movement

Branch movement in response to wind was modeled by swaying the branch toward the wind vector and back over time. To explain in a simple fashion, consider a 2D example, as illustrated in Figure 3, which uses a right-handed coordinate system: part (a) shows the centers of rotation of a plant (which correspond to the nodes); part (b) shows that as the wind pushes left on a branch segment, the angle of rotation increases by around the center of rotation for that branch segment, and part (c) shows that the angle of rotation is decreased as the branch returns to its original position due to the restorative force and tensile strength of the branch. To produce the rotation, the value for is increased over half a time period until a certain angle is reached and then it is decreased back to its original value over the other half of the time period. To achieve a natural slowing down and speeding up at the ends of the arcs, a sinusoidal ease-in/ease-out function is used [50].

Although the example was shown in 2D, the wind vector is actually 3D. Given the branch direction and wind direction defined in the same coordinate system, the axis of rotation for the wind effect is obtained by taking their cross product . The amount to rotate around this axis is determined based on the angular acceleration, calculated by dividing the torque from (1) by the mass of the branch. The maximum rotation, which would cause the branch to align with the wind vector, is . The mass of a branch is arbitrarily set to 10 times the mass of a leaf. A leaf is thus affected 10 times more than a branch by the wind.

3.1.2. Branch Growth

To represent the effect of wind on growth, we assume that there is a species-specific minimum speed , and maximum speed , such that growth is not affected if the wind is less than or equal to and that the full effect occurs if the wind is greater than or equal to . We assume that the effect is linearly interpolated between these speeds. To quantify the effect of wind on growth, we define the wind growth factor , which yields a value in the range , as follows:

The three effects on branch growth are simulated: girth increase, elongation reduction, and cumulative change in growth angle (henceforth called the growth angle). For a wind with uniform velocity (and constant speed ), we use (7), (8), and (9) for these three factors: where is the elapsed time; , , and are the radius (girth), length, and growth angle at time ; , , and are the original radius, length, and growth angle at time 0; and are the normal rates of radius growth and branch elongation; , , and are the proportional changes in radius growth, branch elongation and growth angle due to the wind, respectively; and is a conversion constant from meters per second to degrees per second. The species-specific constants , , and , are assumed to be in the range . Equations (7) and (8) specify that the changes are affected by both the speed and the duration of the wind but the effect of the wind direction is negligible. If the length reaches a maximum determined by the species, no further increase in length occurs. If the wind speed is less than a species-specific minimum speed , then growth is not affected by the wind.

3.1.3. Branch Growth and Motion

See Algorithm 1.

Input: time , branch vector b, and wind velocity w.
Output: branch girth , branch length , and branch growth angle , all at time .
  Set the wind growth factor (see (6)).
  Set the branch radius (see (7)).
  Set the branch length (see (8)).
  If , then .
  Set the growth angle (see (9)).
  Determine the axis of rotation for the branch as .
  Rotate the branch around axis according to the growth angle .
  Additionally rotate the branch around axis according to wind velocity w (1).

3.2. Leaf Response to Wind

The leaf response to the wind is modeled by considering movement and growth, as discussed in Sections 3.2.1 and 3.2.2, respectively.

3.2.1. Leaf Movement

The leaf blade moves in response to wind in one of three ways: bending, twisting, and vibrating. Bending is the motion when the leaf blade is facing the wind (perpendicular to the wind). Twisting is the motion where the leaf blade turns in response to the wind, which occurs when the blade is not perpendicular to the wind. Vibration is tiny movements along the petiole/leaf axis when the leaf blade is parallel to the wind. The petiole response to wind is limited to bending, because vibration and twisting are negligible due to the small surface area. Another consideration is the effect of neighbouring leaves on the motion of the leaf where there is collision and bounce. For the Growth-Flow approach, we ignore these effects based on our assumption that they are relatively small.

As mentioned in Section 2, to model the effects of gusts of wind, the Growth-Flow approach includes a stochastic parameter that influences leaf motion. Let represent the fraction of the wind’s velocity that affects a leaf at the current time. This factor is incorporated in the equations for leaf motion. The remainder of this subsection describes the Growth-Flow approach to leaf motion by considering each of bending, twisting, and vibrating. Throughout, we use time to represent the time interval, to represent the tensile strength of the petiole, and to represent the tensile strength of the blade.

(1) Bending. In the Growth-Flow approach, if a leaf blade is facing (perpendicular to) the wind, then bending occurs either at the petiole center of rotation, as illustrated in Figure 4, or at the blade center of rotation, as illustrated in Figure 4. If the tensile strength of the petiole is less than that of the blade, that is, , then any bending will be petiole bending; otherwise it will be blade bending. We denote the amount of petiole bending by and the amount of blade bending by . Blade bending occurs if the tensile strength of the blade is less than the force of the wind, as illustrated in Figure 4. In either case, we model the bending as a rotation at a single pointas follows: The direction of bending, , stays constant until the sum of the previous angle of rotation, , plus an update of in the same direction would be higher than the maximum angle or lower than the minimum angle , in which case switches, as described in (10). The update to the bending, , is computed based on the sign, , the wind speed, , a conversion constant, , over a change in time . The complete motion, shown in Figures 5 and 6, is accomplished by incrementing the current angle with the change in angle at the center of rotation, as described by (11). If the leaf is oriented on the plane, bending occurs in the dimension.

(2) Twisting. Another motion leaves exhibit is twisting around the petiole. We modeled twisting as a rotation of the petiole at the point where it joins the branch. We assume that in the local frame of reference the leaf is initially lying on the plane and the branch is aligned roughly with the -axis. The amount of twisting, , is incremented/decremented as specified by the sign, , over a change in time, , influenced by the wind speed, , subject to a conversion constant , in the range from to , as described in (12). The complete motion, which is shown in Figure 7, is accomplished by incrementing the current angle with the change in angle, as described by

(3) Vibrating. The third motion a leaf exhibits is vibration, as illustrated in Figure 8. If the leaf is oriented flat on the plane, vibration is in the dimension. The amount of vibration, , is incremented/decremented as specified by the sign, , over a change in time, , influenced by the wind speed, , subject to a conversion constant , in the range from to , as described in (14). The complete motion, which is shown in Figure 8, is accomplished by incrementing the current angle with the change in angle, as described by (15)

The three equations are the same except for varying angles of rotation and varying values for the maximum and minimum angles of rotation. They can be implemented with a single function with different values for its parameters.

3.2.2. Leaf Growth

Changes in leaf thickness and leaf radius are modeled in (16) and (17), respectively, as follows: where is the elapsed time, is the resulting leaf thickness at time , is the original leaf thickness, is the normal leaf thickness growth rate, is a species specific constant, is the resulting leaf radius at time , is the original leaf radius, is the normal leaf growth rate, is the proportional change in the radius, and is a species specific constant.

3.2.3. Leaf Growth and Motion

See Algorithm 2.

Global Initialize:
      Choose a random type of leaf motion
Input: Time , old leaf parameters ( ), wind speed
Output: New leaf parameters ( ), thickness ( ), and radius ( )
 If or , that is, if the wind exceeds the maximum tensile strengths of the
 petiole ( ) or blade ( ), then exit because the leaf is no longer attached.
, , , and
 If
   If ,
     Bend the petiole at the petiole center of rotation
     (See (10) and (11))
   else
     Bend the leaf blade at the leaf center of rotation
      ((11) with β instead of θ)
 else if
   Twist the leaf (see (12) and (13))
 else if
   Vibrate the leaf at the leaf center of rotation (see (14) and (15))
 Calculate the leaf thickness (see (16))
 Calculate the leaf radius (see (17))

4. Results

The tree branch, petiole, and vein structure were generated using a nested L-System [51], and the leaves were generated using the Radius Variation method [52]. Figure 11 demonstrates the effect of a 10 km/h, left-to-right wind on a tiny branch with four leaves attached by petioles over one full cycle of swaying motion. To simplify viewing, in this figure, the branch is shown at one-third scale and the leaves and petioles have their normal size. During the first 16 frames, the branch is moving from right to left. For example, the leftmost petiole visible in the first frame can be seen (where visible) to gradually move to the left in the first 16 frames. Then it moves back to the center over the next 16 frames, it continues to the right over 16 more frames, and finally back to the left over the remaining 16 frames to complete the cycle. The leaves also exhibit a variety of bending, twisting, and vibrating behaviors. For each frame, one of the three behaviors is randomly selected for updating (1/3 probability each). We set the maximum and minimum angles for these behaviors as , , , , , , , and . Experiments with other ranges for the three behaviors also resulted in apparently natural plant behavior, although the ranges chosen seemed the most realistic visually. Overall, as Figure 9 illustrates, the responses of the leaves and the branches to wind do mimic nature.

Further illustration of the response to wind is provided in Figure 10 for a seven-year-old tree. As in Figure 9, the tree trunk and branches are shown at a smaller scale.

The Growth-Flow method, in addition to supporting full motion response of branches and leaves to the wind, also limits growth as illustrated in Figure 11. Parts (a) and (b) of this figure demonstrate growth over a two-second sample without wind and with wind, respectively. There is a marked difference in the plant’s growth with the branch length and leaf size both reduced when wind is present, which matches observations about natural trees.

As shown by the above results, the Growth-Flow method is effective at calculating leaf behavior in response to wind, both with respect to motion (Figures 9 and 10) and growth (Figure 11). The wind effect on plant motion results in the bending, swaying, and twisting of the plant. The wind effect on plant growth results in thickening of leaf blade, reduction of leaf radius, reduction in stem elongation, and increase in stem radius and influences the direction of growth to align more closely with the wind vector, when its direction is consistent for a length of time.

5. Discussion

The Growth-Flow method for animating plants while considering the influence of both motion and growth, can be contrasted with previous work. Crucially, none of the previous work considers the effect of growth on plant animation. In particular, the most similar approaches to plant motion to the one adopted here, including the segment-based [1], stochastic-based [6], and physics-based [13, 24] methods, did not consider growth factors. Feng et al. [14] proposed a method based on only two rotations local to the leaf instead of the three rotations we consider. Although growth is not considered, Ota et al. [53] use noise functions with three movements (horizontal, vertical, and rotational) to provide a more sophisticated stochastic model of plant motions. Our method focuses on the local behavior of the leaf in response to wind. In contrast, the method proposed by Zioma [8] considers the lifting effects of wind on groups of leaves, but does not consider the response of individual leaves. Akagi and Kitajima [2] devised a boundary condition map to locate tree parts in a simplified grid, which reduces the number of calculations needed for motion simulation, but ignored the effect of wind on growth.

The current paper considers branch elongation, branch radius increase, leaf thickening, leaf radius decrease, and change in branch growth angle. Other diverse approaches have been devised for animating plant motion, including simulating motion with filtered uncorrelated noise in a frequency domain [7, 53, 54], simulating overall motion by using rotations based on torque and oscillation [55, 56], using mass-spring networks [57, 58], displaying texels with billboarding techniques [59], and using a finite element model (FEM) [6062]. None of these diverse methods included growth.

6. Conclusion

Both motion and growth of branches and leaves are influenced by the wind. The Growth-Flow method proposed in this paper is the only method known for plant animation that considers both. This method, as illustrated in Figures 9 to 11, includes the following changes to growth: decrease in branch elongation, increase in branch radius, change in branch growth angle, increase in leaf thickness, and decrease in leaf size. All of these effects were simulated and match the observed response of plants to wind duration and velocity over extended periods of time.

Further study is required to detect collisions between plant parts and to add the influence of petiole and leaf motions to those of the branches. In addition utilization of a global approach to incorporate leaf-petiole-branch interactions requires investigation. Growth is also influenced by other environmental factors, such as temperature and nutrient, carbon dioxide, and water availability, which could be incorporated in to the algorithm as well.