Abstract
Tree biomechanics studies using dynamic methods of analysis are reviewed. The emphasis in this review is on the biomechanics of open-grown trees typically found in urban areas, rather than trees in forests or plantations. The distinction is not based on species but on their form, because open-grown trees usually grow with considerable branch mass and the dynamic response in winds may be different to other tree forms. Methods of dynamic analysis applied to trees are reviewed. Simple tree models have been developed to understand tree dynamic responses, but these largely ignore the dynamics of branches. More complex models and finite element analyses are developing a multimodal approach to represent the dynamics of branches on trees. Results indicate that material properties play only a limited role in tree dynamics and it is the form and morphology of the tree and branches that can influence the dynamics of trees.
The biomechanical studies on trees that have taken a dynamic approach to their analysis are reviewed in this paper. Studies published in the last 20 years are mainly considered, with older, seminal studies included where appropriate. The field of biomechanics is often broken into two complementary approaches: statics and dynamics. Both methods are useful in studying the structure of trees under mechanical loading. This review of tree dynamics is part of a project on tree biomechanics that includes a review of tree statics presented as a separate paper (Dahle et al. 2013, in review).
This paper includes an introduction to biomechanical studies of trees; a review of the main literature on tree dynamics; the different dynamic methods of analysis that have been used on trees, including approaches used in forestry, urban forestry, and wind tunnels; and a summary of the complex multimodal studies that consider both tree and branch dynamics. The final section presents research on tree dynamics that is indicating how the form and morphology of trees influences the dynamic response in winds because the branch dynamics can be important.
The emphasis in this review is on the biomechanics of open-grown trees rather than trees in forests or plantations. The distinction is not based on species, but on their form, because open-grown trees, both excurrent and decurrent, usually have considerable branch mass. The term open-grown trees is used rather than urban trees, because the biomechanical principles are not unique to urban trees but rather to all trees of the open-grown form.
Despite the abundance of literature describing the interaction of wind and trees, particularly as it relates to tree dynamics (Moore and Maguire 2004; de Langre 2008; Gardiner et al. 2008; Sellier et al. 2008; de Langre 2012), the literature is almost nonexistent regarding recommendations for pruning open-grown trees to reduce wind damage (Smiley and Kane 2006; Gilman et al. 2008a; Gilman et al. 2008b; Pavlis et al. 2008). Therefore, this review concludes with a list of perceived knowledge gaps in the field of tree biomechanics.
Biomechanical Studies of Trees
Biomechanics applies the basic principles of structural engineering theory to the study of plant forms, including trees. A fundamental premise is that plants cannot violate the laws of physics (Niklas 1992). Biomechanics studies trees as mechanical objects (de Langre 2008), using engineering and physical principles in an attempt to understand the structural properties of trees and how they interact with the environment. The growth rate of trees is largely determined by physiological constraints, particularly those affecting photosynthesis and water transport. But regardless if these are optimal, tree size and shape are still limited by biomechanical constraints (Spatz and Bruechert 2000). Wood in trees is flexible and behaves as neither an ideal solid nor an ideal fluid (Vogel 1996). Wood and most plant materials are described as viscoelastic because their mechanical properties are both elastic and viscous (fluid like). These properties result in non-linear behavior (Miller 2005), and under mechanical loading, plant material does not act like steel or concrete and may not conform exactly to current mechanical models. For this reason it is important to be aware of the limitations of trying to get an exact value for a plant parameter, and to recognize when theory and reality fail to coincide (Niklas 1992). Furthermore, biological materials acclimate and can change their material properties as they age and grow (Lindström et al. 1998; Lichtenegger et al. 1999; Reiterer et al. 1999; Brüchert et al. 2000; Spatz and Brüchert 2000; Lundström et al. 2008; Dahle and Grabosky 2010b; Speck and Burgert 2011), so the dynamic responses can be difficult to predict.
Dynamics and Trees
Wind exerts the largest dynamic forces on trees and is the most important factor for dynamic loading on plants in the terrestrial environment (Niklas 1992). Tree response to wind is ultimately a dynamic process (de Langre 2008), and although understanding the static behavior of trees provides a good basis for understanding their overall behavior, it is a simplification of reality (Moore and Maguire 2004). One of the main reasons for studying trees in winds is to assess their stability and some of the earliest studies recognized that windthrow is also a dynamic process (Coutts 1986).
A dynamic analysis is more complicated than a static analysis because it includes all the static forces and additional components of inertial forces due to the motion, the damping forces and the dissipation of energy, the displacement and phase differences, the natural frequencies, and the consequent changes in motion (Den Hartog 1956). A force applied in a static manner will result in a deflection of a certain magnitude. The same force applied in a dynamic or cyclic manner, at a certain frequency, may increase or amplify the motion and produce a larger effect than the same force applied statically. This effect is called the dynamic amplification factor (DAF) or the dynamic response factor (DRF) and has been applied to trees in only a few studies (Sellier and Fourcaud 2009; James 2010; Ciftci 2012; Ciftci et al. 2013).
Open-grown Trees
The shape or morphology of the tree and the distribution of oscillating branch masses becomes important during dynamic studies (Rodriguez et al. 2008; Sellier and Fourcaud 2009; Ciftci et al. 2013). Slender forest conifers sway in a relatively simple manner, whereas open-grown trees, with many independent and larger branch masses, sway in a complex manner that is different from forest conifers and not yet fully understood (James et al. 2006). The dynamic interaction of branches in winds can significantly modify the frequency and damping of a tree (Moore and Maguire 2005).
Dynamic studies of trees aims to understand how individual trees respond in winds (Baker and Bell 1992; Roodbaraky et al. 1994; James et al. 2006; Kane and Smiley 2006; Baker 1997; Castro-Garcia et al. 2008; Kane et al. 2008; Kane and James 2011; Ciftci 2012) with a few studies investigating the effect of pruning on wind loads (Smiley and Kane 2006; Pavlis et al. 2008; Gilman et al. 2008a; Gilman et al. 2008b; James 2010; Ciftci et al. 2013). Studies examining tree failure due to winds in urban areas have been undertaken after wind storms (Duryea et al. 2007; Lopes et al. 2007; Kane 2008; Matheny and Clark 2009) but with only limited correlation to actual wind velocity and gustiness.
MAIN LITERATURE AND CONFERENCE PROCEEDINGS
Tree Dynamics
There have been several reviews on trees and wind, but they are mainly focused on forest trees (Moore and Maguire 2004; de Langre 2008; Gardiner et al. 2008). A bibliography for tree care professionals was published by Cullen (2002a). Moore and Maguire (2004) reviewed the concepts and dynamic studies by examining the natural frequencies and damping ratios of trees in winds. Gardiner et al. (2008) reviewed the mechanistic modeling of forest trees and the risk of damage in plantations. De Langre (2008) reviewed the literature on plants and examined the more complex fluid mechanics and multimodal models that are being developed to describe the complex dynamic responses of plants and trees. While not exclusively on dynamics of trees, there have been several major conferences on wind and trees that have published proceedings or books with contributions from many authors (Coutts and Grace 1995; Ruck et al. 2003; Mitchell 2007; Mitchell 2008; Schindler et al. 2012) and also conferences on plant biomechanics (Telewski et al. 2003; Salmen 2006; Speck and Burgert 2011; Moulia and Fournier 2012; Thibaut 2012). The first urban tree biomechanics conference was held in Savannah, Georgia, U.S. (Smiley and Coder 2002).
Winds and Tree Damage
Wind data can be expressed in a number of ways, including scales, such as the Beaufort Scale, and more commonly as wind speeds that use a variety of units (e.g., miles per hour, kilometers per hour, knots, m s−1). This can be an obstacle to disseminating knowledge and for practical tree risk management (Cullen 2002b). Instantaneous wind speed is usually not available, and it is customary to quote an average wind speed (either 10-minute or one-hour average) and a gust wind speed taken as a three-second average (Holmes 2007). The wind speed at which tree failure begins to occur is defined in forest studies of trees as the critical wind speed (Oliver and Mayhead 1974; Petty and Swain 1985; Coutts 1986; Blackburn et al. 1988; Peltola and Kellomaki 1993; Hedden et al. 1995; Peltola 1996b; England et al. 2000; Gardiner et al. 2000; Zhu et al. 2000; Cullen 2002b; Zeng et al. 2007; Gardiner et al. 2008; Schelhaas 2008; Wood et al. 2008). Mechanistic models in forestry research use the critical wind speed value to calculate the percentage of failure likely to occur in a forest stand and the approach does not focus on individual tree failure (Gardiner et al. 2008).
Mayer (1987) cautioned that no tree can survive violent storms and posed the question of how the results from investigations on tree sways can be used in practice. How trees fail under dynamic wind loading is not known because the actual dynamic process has never been verified in field experiments due to lack of measurements (Hale et al. 2010). The assumption that the extreme (maximum) wind loading in any particular storm is the key factor in determining whether damage occurs has never been verified in field experiments and it is possible that root fatigue (Rodgers et al. 1995) from a number of storms could actually be more important (Hale et al. 2010).
Studies of the impact of hurricane force winds on urban trees in Florida, U.S., found failure by trunk breakage exceeded overturning in some species [Pinus elliottii (slash pine), 64% broke during Hurricane Jeanne], but other species [Pinus clausa (sand pine)] had 71% breakage during Hurricane Jeanne (Duryea et al. 2007). Following another hurricane (Ivan), uprooting was the main mechanism of failure. Other than post-storm surveys that relate estimated wind speed to tree failures (Kane 2008), there are at present no definitive methods that predict tree failure at a specific wind speed.
Dynamic Analysis Methods
The principles of the dynamic behavior of structures was first published by Den Hartog (1956), and most subsequent texts (e.g., Clough and Penzien 1993; Chopra 1995; Balachandran and Magrab 2004) still use the fundamental equations described in this book. Dynamic analysis examines the forces and displacements of moving structures and considers the inertial forces of mass (m), the elastic forces (k) as in a spring, and the damping forces (c) that dissipate energy. A static analysis considers only the spring forces (k).
When studying the dynamic behavior of a structure, three different approaches are commonly used (Clough and Penzien 1993):
lumped-mass procedure, mass concentrated at a discrete point
generalized displacements for uniformly distributed mass, where a trunk is treated as a beam
the finite element method (FEM)
The Lumped-mass Procedure
The lumped-mass procedure assumes the mass is concentrated at a discrete point as it oscillates dynamically. This greatly simplifies the analysis because inertial forces develop only at these mass points. This method has been used to develop spring-mass-damper models for trees as a single mass (e.g. Milne 1991; Miller 2005), as seen in Figure 1, or as a complex system of coupled masses that represent the trunk and branches (James et al. 2006; Theckes et al. 2011; Murphy and Rudnicki 2012).
Dynamic models using a spring-mass-damper system representing: (a) a tree as a single mass (Miller 2005), and (b) as multiple masses with a trunk and branches (James et al. 2006).
A simple spring-mass-damper system (Figure 1a) is described by a second-order differential equation:
1
where c, m, and k are damping coefficient, mass, and stiffness, respectively; x, ẋ, and ẍ are the displacement, velocity, and acceleration, respectively; and f(t) is the wind-induced time varying (dynamic) force. Equation 1 describes the motion of a single degree of freedom system (SDOF), and if used for analysis of a tree (Figure 1a), approximates the tree to a single oscillating mass (m) with a stiffness (k) and a damping (c) (Miller 2005). A more complex mass model representing branches as oscillating masses attached to a main trunk (Figure 1b) extends this concept to consider branches as oscillating masses attached to the main trunk (James et al 2006).
The oscillating lumped-mass model has been used for trees (Milne 1991; Baker and Bell 1992; Peltola and Kellomaki 1993; Guitard and Castera 1995; Peltola 1996a; Baker 1997; Kerzenmacher and Gardiner 1998; Saunderson et al. 1999; Flesch and Wilson 1999b; England et al. 2000; Miller 2005; James et al. 2006; Jonsson et al. 2007; James 2010; Thekes et al. 2011; Murphy and Rudnicki 2012). Analyses of the mass-spring-damper model of a tree may include a spectral analysis approach using Fourier transformations and transfer functions based on a SDOF model that is often not explicitly stated (Peltola 1996b; Rudnicki et al. 2008; Schindler 2008).
A simple model of a tree (Figure 1a) has a dynamic response whose response amplitude is frequency dependent. Depending on the frequency of sway, the dynamic response is dominated by stiffness, damping, or inertia (Balachandran and Magrab 2004). At low frequencies, the response is dominated by stiffness. As the frequency of the applied force increases, the dynamic response increases until it equals the natural frequency of the system. At this point resonance occurs and there is an amplification of the sway, which depends on the damping, and is known as the damping-dominated region. As frequencies increase further, the rapid force impulses do not cause the mass to move because of its inertia; this is known as the inertial region.
In the damping-dominated region, at frequencies close to the natural frequency, the amplification of sway response has been described for trees as a DAF (Sellier and Fourcaud 2009; Ciftci et al. 2013) (James 2010).
The DAF applied to trees by Sellier and Fourcaud (2009) was defined as the ratio of the maximum displacement under turbulent wind to the displacement caused by the static, instantaneous wind force. DAF was calculated at breast height and at the base of the live crown of a 35-year-old maritime pine (Pinus pinaster Ait.), with values between 0.98 and 1.19. These values seem low due to the DAF being based on displacements, which at breast height would always be small. Ciftci (2012) used FEM to investigate the effect of branches on DAF of a large sugar maple (Acer saccharum L.), also finding that changes to tree geometry induced greater changes in DAF. However, recent studies have indicated that different growth forms in woody plants show distinct ontogenetic trends in mechanical properties (Dahle and Grabosky 2010b; Speck and Burgert 2011), so material properties cannot be ignored in dynamic analyses (Moore and Maguire 2008; Ciftci 2012).
The DRF (James 2010) was defined as the ratio of maximum base moment to mean base moment. It varied among species; more flexible trees (Cupressus sempervirens L., Washingtonia robusta H. Wendl.) exhibited higher values than stiffer trees (Agathis australis D. Don).
Damping has the effect of reducing the amplitude of oscillation and is most effective around the natural frequency region. Damping has little effect at lower frequencies, shown as the static region, and also has little effect at the higher frequencies where the inertia of the mass is the dominant effect on the response. Damping is usually not well understood in vibrating structures (Clough and Penzien 1993) and may be more complex in nature as it may have a non-linear response to produce soft and hard spring mass systems (Miller 2005). In trees, damping forces are considered velocity dependent (Kollmann and Krech 1960; Moore and Maguire 2004; Jonsson et al. 2007) and include frictional forces, aerodynamic drag, collisions, and internal (viscoelastic) forces (Milne 1991).
Because the amplitude response of a dynamic structure is frequency dependent, the natural frequencies of trees have been investigated by either (a) inducing sway in still air conditions, usually with an attached rope (Sugden 1962; Mayhead 1973a; Mayhead et al. 1975; Milne 1991; Gardiner 1992; Roodbaraky et al. 1994; Guitard and Castera 1995; Baker 1997; Flesch and Wilson 1999; Moore and Maguire 2004; Jonsson et al. 2007; Kane and James 2011) or (b) by measuring the tree response in wind conditions and using a power spectrum approach (Holbo et al. 1980; Peltola et al. 1993; Gardiner 1995; Hassinen et al. 1998; James et al. 2006; Moore 2008; Rudnicki et al. 2008).
Complex mass models of trees (e.g., Figure 1b) produce more complex dynamic responses known as multimodal responses, which are discussed later in this review, and also develop mass damping when two or more coupled masses oscillate.
Dynamics of Beams with Distributed Mass
Another method of dynamic analysis considers the structure of a beam or column with the mass distributed along its length. The dynamic equation for a uniform vibrating beam is a fourth-order partial differential equation that is accurate for small deflections, and has been used to study the oscillations and damping of the stems of woody and non-woody plants (Finnigan and Mulhearn 1978; Mayer 1987; Spatz and Speck 2002; Brüchert et al. 2003; Speck and Spatz 2004). Using this equation for tree analysis assumes the tree is like a beam with mass distributed along its length. The first structural model of a tree (Greenhill’s model, Spatz 2000) considered the tree as a pole (Figure 2), and used a static analysis to calculate how tall a tree could grow before it buckled under its own weight (Spatz 2000). There was no consideration of dynamic loads from winds.
Plant stems considered as a beam with distributed mass (Brüchert et al. 2003).
Greenhill’s (1881) simple pole model for trees has been the conceptual basis for both static and dynamic analyses and has been used to analyze dynamics of trees growing in closely spaced plantations or forests (Papesch 1974; Finnigan and Mulhearn 1978; Mayer, 1987; Wood 1995; Peltola 1996b; Flesch and Wilson 1999; Gardiner et al. 2000; Spatz 2000; Novak et al. 2001; Spatz and Speck 2002; Bruchert et al. 2003; Gardiner et al. 2005; Jonsson et al. 2007; Spatz et al. 2007; Moore and Maguire 2008; Rudnicki et al. 2008; Schindler 2008).
The dynamic response of an oscillating beam is more complex than for a single mass because the beam can vibrate in many modes. The first mode is a simple back and forth sway of the whole beam at a frequency known as the natural or fundamental frequency. Other sway responses are possible and the beam may deflect in different shapes (known as mode shapes) that occur at different frequencies (Figure 3a). In theory, a beam considered a uniform continuous structure has an infinite number of vibrating modes, but in practice, most of the energy of vibration occurs in the first few modes. The first or fundamental mode occurs at the lowest frequency, and has the most energy and amplitude.
Dynamic modes applied to trees: (a) modes of a beam (Schindler et al. 2010) and (b) modes of branched structures (Rodriguez et al. 2008).
Finite Element Method
In dynamic analysis, FEM combines features of both the lumped mass and uniformly distributed mass procedures. It is applicable to all structures and requires computer analysis due to the complex calculations (Sellier et al. 2006, Dupuy et al. 2007; Rodriguez et al. 2008; Moore and Maguire 2008; Sellier and Fourcaud 2009; Theckes et al. 2011; Ciftci 2012; Ciftci et al. 2013).
FEM divides a structure or beam into an appropriate number of elements whose sizes may vary, and the ends of each element (nodes) become the generalized coordinates. The deflection of the complete structure can then be expressed in terms of generalized coordinates. This method is good for one- and three-dimensional structures and has the advantage of being able to select the desired number of generalized coordinates by dividing the structure into the appropriate number of segments. For uniform materials, such as steel and concrete, interpolation functions of each segment may be identical and computations are simplified (Figure 4).
Finite element models showing the crown structure of three trees (Moore and Maguire 2008).
An advantage of FEM is that complex wind-loading scenarios can be modeled. The dynamic response of the structure (i.e., the tree) is important, but an equally important factor is the wind loading, which can be quite complex (Finnigan and Brunet 1995; Belcher et al. 2012). Recent FEM studies have investigated tree response to different wind-loading scenarios (Sellier et al. 2008; Sellier and Fourcaud 2009). Use of FEM to explore the complex structural dynamics of decurrent trees holds great promise, but it requires accurate empirical measurements of many parameters peculiar to the tree and loading conditions to produce a reliable result.
STRATEGIES USED IN DYNAMICS RESEARCH ON TREES
Different strategies have been used by various researchers to study tree biomechanics, and tree dynamics in winds. In this review the strategies are broadly grouped as:
Forestry – (trees in groups) economic damage on plantation grown trees,
Open-grown trees – (individual trees with branches, both excurrent and decurrent) tree stability and risk assessment, predominantly in urban areas,
Wind tunnels – small trees to measure drag coefficients in constant velocity winds, and
Modeling – (dynamic models of trees) computer studies, finite element methods and mathematical modeling.
Forestry
Forestry studies examine plantation trees (mainly conifers) and the economic losses caused by damaging winds (Moore and Maguire 2005; Peltola 2006), usually with the aim to determine threshold values of storm damage. Threshold values include wind speed, gustiness, duration of storm, terrain, soil type, soil moisture, stand characteristics (e.g., height, density, diameter at breast height, crown length), and the physical condition of a tree (Mayer 1987). The threshold value of wind speed at which damage to trees occurs, termed the critical wind speed, is an important variable for forest managers (Peltola 2006) and in forest modeling (Gardiner 1995; Moore and Maguire 2004; Frank and Ruck 2008). The factors of site, tree species, soil, wind climate, critical wind speed, and silvicultural treatments, such as thinning, are considered together in order to calculate the risk of damage to a naturally regenerated forest or plantation. The forestry studies estimate the percent damage to a group average of trees, rather than explicitly predicting failure of any individual tree.
To predict the percentage of trees in a forest stand likely to fail during a storm, mechanistic models have been developed (Gardiner et al. 2008; Frank and Ruck 2008; Schelhaas 2008; Wood et al. 2008). The models calculate the critical wind speed required to break or overturn trees, and then determine the probability of damage at the geographic location of the trees, based on some assessment of local wind climatology and empirical relationships. Models have been shown to be valid in certain circumstances (Gardiner et al. 2000), but their deterministic nature is sometimes at odds with field observations of wind throw. By definition, the models are restricted to excurrent trees in plantations and are really an application of statics to a dynamic phenomenon, and so are not yet applicable to open-grown trees in urban areas.
There has been considerable work using static pulling tests, mainly on forest conifers (Nicoll et al. 2006) and on small and young trees (Lundström et al. 2007) with a high slenderness ratio usually above 50 and often over 100 (e.g., slenderness values [46-136] Hale et al. 2010; [58-94] Jonsson et al. 2006). Tree-pulling tests have had an important role in providing valuable information on mechanical stability of trees of varying size and tree species, and the information is useful in mechanistic modeling, but the simulation of static loading by tree pulling alone is not enough to explain the mechanical stability of trees (Peltola 2006).
The critical wind speeds that cause failure depend on tree species, growth pattern, and location, and estimates vary. However, ultimately few tree species can survive violent storms with mean wind speeds over a period of 10 minutes, exceeding 30 m s−1 near the top of the canopy without damage (Peltola 1996a).
Open-grown Trees
The distinction between open-grown trees and forest trees is made in this review because of differences in the growth and form of the trees, particularly with respect to their canopy architecture. Research on dynamic response of forest trees may not be applicable to open-grown trees that develop a complex distribution of branch masses. When applying dynamic methods to tree sway in winds, recent research has indicated that the branches and the form of the tree are important in understanding how trees respond in winds (James et al. 2006; Spatz et al. 2007; Rodriguez et al. 2008; Sellier and Fourcaud 2009; Theckes et al. 2011; Ciftci et al. 2013).
Research on open-grown trees in winds aims to understand how individual trees respond in winds, and investigates aerodynamic properties (Baker and Bell 1992; Roodbaraky et al. 1994; Baker 1997; Ennos 1999), dynamic properties of frequency and drag (Kane and Smiley 2006; Kane and James 2011), effect of pruning dose on wind response (Gilman et al. 2008a; Gilman et al. 2008b; Pavlis et al. 2008), and wind loads (James 2006; James 2010).
There is very little data on the wind loading of open-grown trees during storms, and much of what we know about how trees fail comes from post-storm tree damage surveys (Duryea et al. 2007; Lopes et al. 2007; Kane 2008; Matheny and Clark 2009). Extreme European wind storms on December 26–28, 1999, were directly responsible for killing 95 people in France; 15 in Germany; 11 in Switzerland; 11 in the United Kingdom; and 5 in Spain. Damage was estimated at more than USD $10 billion where wind speeds exceeding 160 km/h were recorded along the French coast (Lopes et al. 2007).
There is currently no definitive method to predict failure of an individual tree. Arboricultural assessments of trees include visual tree assessment (Mattheck and Breloer 1994), tree risk assessment methodology (Smiley et al. 2011), quantified tree risk assessment (Ellison 2005), and statics integrated methods that combine static pulling with dynamic wind load assessment (Wessolly 1991; Brudi and van Wassenaer 2002; Detter and Rust 2013). The effect of pruning dose and trunk movement in tropical storm winds has been investigated (Gilman et al. 2008a), using artificially generated winds at speeds up to 26.8 ms−1 on trees of 6.1 m average height. Trees generally moved similarly in wind regardless of ANSI pruning type applied, although crown/branch thinning may be more effective in reducing motion than other pruning types. Gilman et al. (2008a) suggested that it may not be wise to extrapolate these results to larger trees, and that further testing is required to examine the pruning effect of individual branches when they are coupled as a continuous dynamic structure. Conflicting results were obtained in further studies using similar trees and methods, where crown thinning was less effective at reducing trunk movement (Gilman et al. 2008b). In this study it was noted that branches on thinned trees appeared to move more than branches on other treatments but not in the same direction. This complex branch movement indicates that the dynamic effects of branches may play an important role in acting as a buffer to dampen and reduce motion (Moore and Maguire 2005; James et al. 2006).
Wind Tunnels
Wind tunnel tests have been used to study wind effects on trees, but there are serious limitations due to the size of the wind tunnel and the trees that can fit into them. In general, derived results are only strictly applicable under similar conditions, but the forces are of the right magnitude for mechanistic models, similar to static pulling experiments (Peltola 2006). Scale models of trees have been used in wind tunnels to represent forest trees, to study the dynamics of wind turbulence on forest canopies and to examine the effects of commercial practices such as thinning and spacing (Stacey et al. 1994; Gardiner and Stacey 1996; Gardiner et al. 1997; Gilman et al. 2008a).
By necessity, trees in wind tunnels are small and the wind flow conditions are steady state or quasi-static (Holmes 2007) and drag dominated. Because conditions in wind tunnels are quasi-static, tests on trees have previously been reported in static papers on trees (Peltola 2006) rather than in dynamic reviews. One of the problems with results from wind tunnel tests is the question of scale, and how to select the appropriate wind speed in relation to the scale of the model and of full-sized trees (Peltola 2006).
Wind tunnel tests on trees have been performed on individual scale models (Tevar Sanz et al. 2003; Gromke and Ruck 2008), on model canopies (Finnigan and Mulhearn 1978; Wood 1995; Gardiner et al. 1997; Novak et al. 2001; Gardiner et al. 2005), and on small trees (Fraser 1967; Mayhead 1973b; Rudnicki et al. 2004; Vollsinger et al. 2005; Cao et al. 2012) and individual leaves (Vogel 1989). Some studies conducted in wind tunnels have investigated the effect of pruning on drag of conifers (Fraser 1967; Mayhead et al. 1975; Rudnicki et al. 2004) and deciduous trees (Vollsinger et al. 2005), but interpretation of results is limited because few replications were used (Fraser 1967; Mayhead et al. 1975) and the trees were small (less than 2 m tall) (Rudnicki et al. 2004; Vollsinger et al. 2005). Smiley and Kane (2006) and Pavlis et al. (2008) simulated wind tunnel conditions by placing small trees on the back of a truck and driving at high speed. They examined the effect of pruning on drag reduction by applying different pruning methods. Reduction in drag induced bending moment differed by pruning type, mainly due to the mass of foliage removed, but predicting the reduction in drag was not reliable based on area of crown removed. Tree mass was the best predictor of drag for red maple (Acer rubrum), but these results were on small trees, and the authors recommended caution when extrapolating drag values to larger red maples.
A frequently cited study of drag on British forest trees (Mayhead 1973b) used a wind tunnel to determine drag coefficients, but Mayhead commented that it is probably unsound to test trees less than 3–4.5 m high because larger trees have a different morphology (Niklas 1994a; Niklas 1995; Osunkoya et al. 2007; Dahle and Grabosky 2010a). Mayhead (1973b) found large variations in results both between and within species and suggested that the range of variation was either natural or a result of poor technique.
Wind tunnels are used by civil engineers to determine drag on solid objects known as bluff bodies (Holmes 2007), and constant wind velocity is used to create steady state or quasi-static conditions. Bluff bodies have a fixed frontal area exposed to the wind, a fixed shape that has a set value of streamlining, and a constant drag coefficient that is proportional to the square of velocity. Results from small bluff body models may be scaled up for large structures such, as tall buildings, where additional factors, such as the aerodynamic admittance function, may need to be considered (Holmes 2007). These methods may not be suitable for flexible objects, such as trees, because the response of small-scale models under constant wind speed conditions may not represent the response of large trees under actual wind conditions (Mayhead 1973b). The drag coefficient for trees may not be a constant value and could be proportional to wind speed (v) or the square of wind speed (v2) (Cullen 2002b). Average values of drag for trees are often quoted, but the large range and variability is often overlooked. Mayhead (1973b) reported drag coefficients for several conifer species of importance to British forestry and the results have become standard values used in windthrow risk modeling, despite very small sample sizes (e.g., Gardiner et al. 2000). Gardiner et al. (2005) cautioned that these wind tunnel tests are a simplification of a real forest, and in some instances can only provide a rough approximation to reality.
COMPLEX TREE MODELS
All models used for dynamic analysis of trees make assumptions, but some assumptions (e.g., ignoring branches or treating them as rigid, lumped masses) may not adequately represent the complex dynamic response of trees (Moore and Maguire 2004). More complex models are needed to account for different tree shapes and species, and in particular account for the dynamic influence of branches (Kerzenmacher and Gardiner 1998; England et al. 2000; Sellier and Fourcaud 2009; Theckes et al. 2011; Ciftci 2012; Ciftci et al. 2013). The dynamic effect of branches on frequency and damping became increasingly important as crown architecture deviated from a slender, cantilevered beam (Sellier and Fourcaud 2009).
Models of trees must account for the dynamic contribution of branches, particularly in trees where the mass of branches is significant. For complex botanical structures, such as trees, a multi-degree of freedom system or a multimodal analysis is required to account for complex dynamic interaction of the branches and trunk (de Langre 2008; Rodriguez et al. 2008). The dynamics of trees with many large branches is complex because the swaying branches are attached to other swaying branches and then to the trunk. Each of the swaying masses influences the other swaying masses to create different modes of sway and also has an effect on the frequencies and damping of the overall structure. How the branched architecture and tree geometry influences the dynamics of the tree is therefore a central question to be investigated (Rodriguez et al. 2008). Multimodal analysis has only been used in a few studies to analyze the dynamic characteristics of trees (Fournier et al. 1993; Moore and Maguire 2005; Sellier et al. 2006; de Langre 2008; Rodriguez et al. 2008; Ciftci 2012; Murphy and Rudnicki 2012).
Multimodal response can occur in two different ways: (a) in beams (Figure 3a) and (b) in branched structures (Figure 3b). Multimodal response in beams occurs where the distributed mass along a single beam flexes in a number of modal shapes (Figure 3a) and is described previously in the section on beams. Although it is a multimodal dynamic analysis, the beam model does not account for oscillating branches.
Multimodal response in branched structures (Figure 3b) occurs when several coupled masses (branches) oscillate in a complex manner, often with an in-phase and out-of-phase response so that several modal sway responses are possible. The coupled masses, with their individual oscillation response, are connected to another oscillating mass, resulting in a coupled response of the combined masses. The branched multimodal method has been applied to trees where the branches are considered as coupled masses that oscillate on the trunk, which itself is an oscillating mass (James et al. 2006; Spatz, 2007; Rodriguez et al. 2008; Thekes et al. 2011; Ciftci 2012; Murphy and Rudnicki 2012; Ciftci et al. 2013). It is possible to extend this branched concept to second- and third-order branches where the complexity could be expected to increase further.
Complex models of trees that represent the dynamic oscillations of branches have used either (a) a multiple spring-mass-damper model (James 2003; James et al. 2006; Spatz 2007; Thekes et al. 2011; Murphy and Rudnicki 2012) or (b) a FEM approach (Rodriguez et al. 2008; Ciftci 2012; Ciftci et al. 2013) where modes are generated from branches moving together or apart in a complex manner (as in a fractal tree) (Rodriguez et al. 2008).
Where multimodal response occurs due to the swaying branches oscillating with each other, a damping effect known as mass damping may occur. A mass damping system described by Den Hartog (1956) has been defined for trees (James et al. 2006), and occurs when the branches sway together (in phase) or against each other (out of phase) in a complex manner. Damping from branches has been identified for a two degree-of-freedom system in a T- or Y-shaped branched structure (Spatz et al. 2007; de Langre 2008; James 2010; Thekes et al. 2011; Murphy and Rudnicki 2012; Spatz and Thekes 2013) based on a tuned mass damper system and in trees creates a modal energy transfer (de Langre 2008; Thekes et al. 2011; Spatz and Thekes 2013) as a protective mechanism against large sways. Complex dynamics that include branches could be beneficial to the tree by enhancing wind energy dissipation through a mechanism called multiple resonance damping (Spatz et al. 2007), multiple mass damping (James et al. 2006), or branch damping (Spatz and Thekes 2013). A prerequisite for this mechanism to occur is a multimodal behavior of the tree, with high modal density in the frequency range and significant branch deformations.
This dynamic response was found for trees with contrasting architectures in a three-dimensional modal analysis and FEM modeling (Rodriguez et al. 2008). Branch oscillations influence the dynamic behavior of trees to a greater extent than can be explained simply by their additional mass (Moore and Maguire 2008; Ciftci 2012).
TREE MORPHOLOGY AND MATERIAL PROPERTIES
Size and morphology of trees need to be considered when using complex dynamic analyses because the dynamics of branches affects the oscillating frequency and damping of the whole tree (Rodriguez et al. 2008; Speck and Burgert 2011). In a study of tree aerodynamic behavior it was found that material properties play only a limited role in tree dynamics (Sellier and Fourcaud 2009). In contrast, small morphological variations can produce extreme behaviors, such as either very little or nearly critical dissipation of stem oscillations. Effects of branch geometry on dynamic amplification are substantial yet not linear (Sellier and Fourcaud 2009). Recent studies in the biomechanics of plant stems indicate that different growth forms in woody plants show distinct ontogenetic trends in mechanical properties (Dahle and Grabosky 2010b; Speck and Burgert 2011).
Open-grown trees have diverse branch morphology, as shown in a survey of 40 woody tree and shrub species in New York (Evans et al. 2008). The size of a tree is also an important parameter because large trees have a different morphology to small trees, and it is probably unsound to test trees less than 3–4.5 m tall (Mayhead 1973b). The morphology of branches also changes with size (Bertram 1989; Dahle and Grabosky 2010a), which must be taken into account. Natural morphological variation within and across species of open-grown trees need to be considered, and care taken when attempting to scale up results and extrapolate data. (Mayhead 1973b; Gilman et al. 2008a)
Dynamics studies of olive trees (Castro-Garcia et al. 2008) and walnut trees (Rodriguez et al. 2012) swaying under forced vibration during harvesting also show the multimodal response, similar to the wind excitation results, which is due to the dynamic interaction of branches on the tree.
CONCLUDING REMARKS
The dynamic response of open-grown trees in winds is greatly influenced by the size and form of the tree, and at least partly due to the dynamics of branches. Simple models have been used for forest and plantation trees and are useful for dynamic analysis of slender trees with few branches. However, more complexity, such as through a multimodal approach, is needed for a dynamic analysis of open-grown trees, because the dynamic coupling of branches has an influence on the response of the tree (de Langre 2008; Rodriguez et al. 2008).
There appear to be gaps in the literature on several topics relating to dynamic analysis of open-grown trees and their response in winds, including:
Recommendations for pruning open-grown trees to reduce wind damage (Gilman et al. 2008b).
The dynamic contribution and the damping effects of branches. Studies indicating that the form of the tree has a greater influence than the material properties (Sellier and Fourcaud 2009) have implications for branch removal and future pruning practices.
Modeling of open-grown trees should account for the multimodal branch response. The complexity of dynamic analysis is likely to increase in the near future but will need to be condensed into simpler methods for practical use.
Tree failure under actual wind conditions has not yet been measured (Hale et al. 2010), and so extending the results from current research is difficult, especially when trying to determine tree failure and stability in winds.
Associated with tree failure is the understanding of the energy transfer from the wind to the tree. There is little published data on actual wind loads on trees and understanding the energy transfer process may assist in understanding how trees and branches dissipate energy and dampen the wind energy. This may be an important factor in understanding how trees survive otherwise damaging winds.
Finally, the topic of torsional forces and loads on trunks and branches has not yet been investigated, yet may be critical in understanding the total loads on trees (Niklas 1992). Torsional forces that twist trunks and branches are observed in trees during winds, but no method has yet been developed to measure the dynamic torsional loads experienced by trees during winds.
Acknowledgments
We would like to thank the International Society of Arboriculture and the ISA Science and Research Committee for funding this literature review. Thanks to Aaron Carpenter who assisted in formatting of the manuscript and to two anonymous reviewers for valuable comments.
- © 2014, International Society of Arboriculture. All rights reserved.