A Simple Model for Simulating Tornado Damage in Forests

Andrew P. Holland Carolina Environmental Program, The University of North Carolina at Chapel Hill, Chapel Hill, North Carolina

Search for other papers by Andrew P. Holland in
Current site
Google Scholar
PubMed
Close
,
Allen J. Riordan Department of Marine, Earth and Atmospheric Sciences, North Carolina State University, Raleigh, North Carolina

Search for other papers by Allen J. Riordan in
Current site
Google Scholar
PubMed
Close
, and
E. C. Franklin Department of Forestry, North Carolina State University, Raleigh, North Carolina

Search for other papers by E. C. Franklin in
Current site
Google Scholar
PubMed
Close
Full access

Abstract

An analytical model is presented to describe patterns of downed trees produced by tornadic winds. The model uses a combined Rankine vortex of specified tangential and radial components to describe a simple tornado circulation. A total wind field is then computed by adding the forward motion of the vortex. The lateral and vertical forces on modeled tree stands are then computed and are compared with physical characteristics of Scots and loblolly pine. From this model, patterns of windfall are computed and are compared to reveal three basic damage patterns: cross-track symmetric, along-track asymmetric, and crisscross asymmetric. These patterns are shown to depend on forward speed, radial speed, and tree resistance. It is anticipated that this model will prove to be useful in assessing storm characteristics from damage patterns observed in forested areas.

Corresponding author address: Andrew P. Holland, The University of North Carolina at Chapel Hill, 137 East Franklin St., Chapel Hill, NC 27599-6116. Email: andy_holland@unc.edu

Abstract

An analytical model is presented to describe patterns of downed trees produced by tornadic winds. The model uses a combined Rankine vortex of specified tangential and radial components to describe a simple tornado circulation. A total wind field is then computed by adding the forward motion of the vortex. The lateral and vertical forces on modeled tree stands are then computed and are compared with physical characteristics of Scots and loblolly pine. From this model, patterns of windfall are computed and are compared to reveal three basic damage patterns: cross-track symmetric, along-track asymmetric, and crisscross asymmetric. These patterns are shown to depend on forward speed, radial speed, and tree resistance. It is anticipated that this model will prove to be useful in assessing storm characteristics from damage patterns observed in forested areas.

Corresponding author address: Andrew P. Holland, The University of North Carolina at Chapel Hill, 137 East Franklin St., Chapel Hill, NC 27599-6116. Email: andy_holland@unc.edu

1. Introduction

Official estimation of peak tornadic wind speeds is performed by the National Weather Service and relies on ground and areal surveys of damage paths. The F scale, proposed and developed by T. Fujita (e.g., Fujita 1971), is the scheme most widely used to classify tornadic wind speed based on the severity of structural damage. In this scheme, observed damage patterns are compared with a series of written descriptions and photographic examples of damage to houses, vehicles, trees, and other common objects. There are six categories of damage (F numbers), each corresponding to a range of wind speeds within the overall interval from 15 to over 130 m s−1.

The F-scale classification appears to be generally accepted as a practical means of assessing wind speeds and is routinely used by survey teams. However, the method relies primarily on damage to structures rather than to natural vegetation (Marshall 2003; Schaefer and Livingston 2003). Thus, it is difficult to classify tornadoes that occur in forested areas.

Peterson (2003) has studied tornado damaged forests to determine whether there is a way to incorporate tree damage based on species and size into an improved F-scale system. His results were encouraging when looking at forests that were similar in their composition.

Another possible way to classify tornadoes that do not pass over structures is to analyze the damage done to the trees and to try to determine the maximum wind speed of the tornado. Then the tornado could be classified using the F scale by simply putting it in the correct wind speed category.

As a first step in this endeavor, this study describes a quantitative physical model that can be readily used to assess tornadic wind speeds in forested areas. It is based on a simple Rankine vortex and a modified tree model designed by Peltola et al. (1999). The analytical model is initialized with information on tree species and stand characteristics plus estimates of tangential, radial, and forward speeds of the vortex. The model then generates a pattern of downed trees that can be compared with field observations. From this comparison, one can estimate the wind speed and, therefore, the F-number classification of a tornado.

a. An early model for the wind field

Around 1923, J. Letzmann [reviewed by Peterson (1992)] suggested that when a tornado moves through a forest it produces a damage pattern that relates to the sum of vortex rotation and forward speed. By creating a simple, analytical model of a steady-state vortex, Letzmann was able to simulate tree damage patterns that compared favorably to those observed.

In his model, Letzmann used a combined Rankine vortex to estimate the tangential flow of a tornado. Once a forward speed and radial wind component were specified, Letzmann computed the resultant wind vector for every point in and around the vortex. By experimenting with different combinations of radial, tangential, and forward speeds, Letzmann developed hand-drawn schematic illustrations of several fundamental wind field patterns. To apply the wind model to the forest, he assumed that a tree would fall at some critical velocity and that fallen stems would align with the instantaneous direction of flow.

b. Recent advances in tree response

There are, of course, limitations to Letzmann’s model. It relies on hand calculations and does not include information about the physics of tree response. Much work in the latter area has been done in the last decade. Wood (1995) used wind-tunnel experiments on model trees and pulling experiments in the field to calculate the maximum stresses that different species of trees can withstand. Mattheck and Breloer (1994) developed a simple windthrow model that calculates the additive wind and gravity forces acting upon the tree. If the sum of these forces exceeds the tree’s resistance to uprooting, the tree will fall.

Field studies have been performed to calculate the bending moments of loblolly pine trees (Pinus taeda L.) (Fredericksen et al. 1993; Hedden et al. 1995). In these studies, loblolly pines were toppled using a winch-and-cable system. The force needed to break the stem, that is, the critical turning moment, was measured for trees of different sizes and ages. Building upon earlier work, Peltola et al. (1997, 1999), Peltola and Kellomäki (1993), and Peltola (1996) developed a mechanistic tree model that predicts the critical turning moment and threshold wind speed necessary to break the stems of several tree species. The model is particularly applicable to the problem at hand because it provides a physical link between wind speed and tree breakage. As a first step, the model used in this study considers only stem breakage and does not include uprooting. The latter is a much more complex physical process controlled partly by wind stress and partly by root and soil characteristics—factors that are difficult to observe and therefore are beyond the scope of this paper.

2. Method

This work essentially combines an analytical vortex with the physical tree-response model of Peltola et al. (1999), hereinafter called the Peltola model, to produce simple visual fields of tree-fall patterns. As a starting point, we have extended the Peltola model to include information on loblolly pine, a tree species common to the southeastern United States. The resultant model produces a pattern of tree-fall somewhat analogous to the hand-drawn diagrams of Letzmann.

a. Vortex model description

The model begins with a combined Rankine vortex. Each point in a horizontal Cartesian grid is assigned x, y coordinates. These coordinates are then transformed to polar coordinates, where r is the distance from the origin to x, y:
i1558-8432-45-12-1597-e1
From (1), the tangential velocity Vtan can be calculated at each point using
i1558-8432-45-12-1597-e2
i1558-8432-45-12-1597-e3
where the maximum tangential wind velocity Vmax is specified and where Rmax is the radius of maximum tangential velocity.
Letzmann (Peterson 1992) and others (e.g., Ward 1972; Zrnić et al. 1985; Turner 1966) suggest that the radial velocity Vr (defined as positive inward) varies with radius in a similar fashion:
i1558-8432-45-12-1597-e4
i1558-8432-45-12-1597-e5
where Vrmax is the specified maximum radial wind velocity. Adding the tangential and radial velocities produces the mean wind field for a stationary tornado. To make the model simulate a moving tornado, the forward motion is then added to the wind field (Fig. 1). The wind speed at each grid point is then passed on to the tree model.

b. Tree model

The algorithm is based on the Peltola model. First, one extrapolates the two-dimensional wind field simulated by (2)(5) into a three-dimensional field in the canopy, assuming a log profile of the form
i1558-8432-45-12-1597-e6
where u(z) is the wind speed at height z, u(m) is the wind speed calculated at each grid point and is assumed to be located one tree height above the canopy, h is the tree height, and z0 is the roughness length, taken to be 0.06 times the tree height (Peltola et al. 1999). The tree height is divided into 1-m Δz segments, and, through (6), the wind speed is calculated at each segment.
Second, the model calculates the lateral and vertical forces acting on the tree. The first force considered is the force of wind drag. The mean wind force Fw (N) for each 1-m segment is given as (Monteith 1975)
i1558-8432-45-12-1597-e7
where Cd is the drag coefficient for the segment, ρ is the density of air, and A(z) is the projected cross-sectional area of the tree against the wind.

To calculate the projected area of the crown against the wind, the Peltola model assumes a crown shape simulated by triangles. For Scots and loblolly pine, two triangles share a common base, with the vertex of one triangle extending upward and the other downward. The bases of the triangles are equal to 2 times the length of the longest branch in the crown (Hakkila 1971; Peltola and Kellomäki 1993). For loblolly pine, the upper triangle is taller than the lower to differentiate between the living and dead crowns.

As the wind speed increases, the simulated canopy becomes streamlined. Peltola et al. (1999) determined that the reduction in crown area is 20% when the wind speed is less than or equal to 11 m s−1 and is 60% for wind speeds greater than 20 m s−1. Peltola and Kellomäki (1993) determine the streamlining factor St between these wind speeds as
i1558-8432-45-12-1597-e8
where u(z) is computed from (6).
Once the tree becomes streamlined and bent by the wind force, the force of gravity becomes important. This force is also calculated for each 1-m segment of the tree by
i1558-8432-45-12-1597-e9
where g is the acceleration of gravity and M(z) is the green mass of the stem and crown (Grace 1977; Jones 1983; Petty and Swain 1985; Peltola and Kellomäki 1993).

The mass of each stem segment is calculated from the green density of the stem and the volume of the segment (Clark 1991; Laasasenaho 1982; Peltola and Kellomäki 1993). The vertical distribution of crown mass is assumed to be similar to that of the crown area projected against the wind and is proportional to the stem mass (Baldwin 1987; Peltola and Kellomäki 1993).

The turning moment at the stem base can be calculated for each segment from the stem and crown mass, the wind force, and the horizontal deflection of the stem (Fig. 2), but the force must be adjusted by accounting for the gustiness of the wind and the distance between the trees. The maximum turning moment is then
i1558-8432-45-12-1597-e10
where z is the height of the segment and x(z) is the horizontal displacement of the stem, which is assumed to be directly proportional to the wind force and inversely proportional to the stem’s stiffness (Peltola and Kellomäki 1993). This displacement is given by
i1558-8432-45-12-1597-e11
and
i1558-8432-45-12-1597-e12
where l(z) is the distance from the tree top to the middle of each 1-m segment, Me is the modulus of elasticity (Pa), I is the area moment of inertia of the stem (m4), a is the distance from ground level to crown center, b is the distance from crown center to the tree top, and h is the mean tree height.
The variable Gustfactor is determined empirically from wind-tunnel experiments and is dependent on the tree height, tree spacing, and distance from the stand edge. It represents the changes in mean wind loading to the extreme normalized values (Gardiner and Stacey 1995; Gardiner et al. 1997):
i1558-8432-45-12-1597-e13
where
i1558-8432-45-12-1597-e14
and
i1558-8432-45-12-1597-e15
where s is the tree spacing and x is the distance from the forest edge in tree heights (i.e., x = 0 × h for stand edge and x = 1 × h for one tree height into the stand).
Figure 3 shows the effect that distance from stand edge has on the turning moment. This relationship assumes an infinite upwind open area. Therefore, a gap factor must account for the effect of different upwind gap sizes on the wind loading (Gardiner and Stacey 1995; Gardiner et al. 1997). Thus, we have
i1558-8432-45-12-1597-e16
where
i1558-8432-45-12-1597-e17
i1558-8432-45-12-1597-e18
and gapsize is the size of the upwind gap in tree heights. The effect of the upwind gap can be seen in Fig. 4.
The maximum turning moment at stem base can be calculated by summing the turning moments at each 1-m segment over the height of the tree:
i1558-8432-45-12-1597-e19
This total turning moment is compared with the resistance to stem breakage for each tree, with the latter based on the modulus of rupture for the tree. It is calculated at breast height (z = 1.37 m) and is given by
i1558-8432-45-12-1597-e20
where Mr is the modulus of rupture (Pa) and D is the tree diameter at breast height (DBH; m) (Jones 1983; Morgan and Cannell 1994).

The model moves the tornado vortex over the forest, calculating the turning moment for each grid point and comparing it with the resistance to stem breakage. If the turning moment exceeds the resistance, the stem is considered to be broken.

3. Results

a. Vortex wind field

As in earlier work by Letzmann (Peterson 1992), a moving Rankine vortex with variable inflow of the type described above produces an isotach field that generally consists of a family of concentric circles that enclose a calm point and an adjacent family of crescents that enclose the area of maximum wind speeds. Somewhat as defined by Letzmann, key parameters α and Gmax control the asymmetry of the wind pattern, where Gmax is the ratio of the maximum tangential wind speed to the forward speed and α is the angle between the tangential and summed component wind vectors at the point of maximum winds (Fig. 5). For the simple vortex, it can be seen in (2)(5) that α does not vary with radial distance from the vortex center.

There are several useful simple relationships that can be established for a moving counterclockwise-rotating vortex of this type for which the maxima in both the radial and tangential components occur at the same radius. One can define the axis of extremes (AB) as the straight line intersecting the center of the vortex, the point of absolute calm, and the location of the absolute maximum wind velocity Vx (Fig. 5). The orientation of AB is governed by α independent of the forward motion of the vortex. As the forward motion VS increases, the location of the absolute wind maximum remains on AB at a fixed distance Rmax from the vortex center. However, the calm point moves along AB toward the left side of the storm track and away from the center of the relative vortex. To be specific, the center of the calm circle lies at radius r0 defined by
i1558-8432-45-12-1597-e21

For a few examples, Fig. 6 combines isotach and vector fields to show the differences in the flow pattern when α and Gmax are varied. In Figs. 6a and 6b, α is held constant at 45° such that tangential wind speed is equal to radial wind speed. For this and all subsequent illustrations, Rmax is held constant at 75 m given that it has no effect on overall patterns. Figure 6a represents a slow-moving tornado, that is, one with a rotational speed that is 4.2 times the forward speed (Gmax = 4.2). In this case, the forward speed is 10 m s−1 and the maximum wind speed relative to the ground is 69 m s−1. The isotach pattern is only slightly asymmetric, but the radial inflow causes the area of maximum winds to shift 45° rearward from immediately right of center. In Fig. 6b, where Gmax = 2.1, the vortex forward speed is doubled, causing an increase in the maximum wind speed and a further leftward shift of the area of calm winds. Note that, as the forward speed increases, the location of the center of maximum winds is not displaced outward from the center of the relative vortex (located at the origin of the diagram).

In Figs. 6c and 6d the forward speed is held constant at 10 m s−1 and α is varied. When α = 36°, the tangential wind speed is greater than radial wind speed and vice versa for α = 54°. Note the rotation of the axis of extremes, AB.

b. Critical wind speed

The data used for loblolly pine in the model were taken from Wahlenberg (1960). This table includes the height, DBH, and spacing for trees aged 20–80 yr in site indices of 60–120. The site index refers to the height of 50-yr-old dominant trees of the same species in a given index area. The site index is determined for a tree species based on the soil and climatic factors characteristic of that area (Avery and Burkhart 2002, chapter 15). For every site index, the model assumes a random distribution of tree ages from 20 to 80 yr. The model calculates the critical wind speed for any height, DBH, and spacing. For example, Fig. 7 shows how the critical wind speed varies with height and form quotient, defined as the height in meters divided by the DBH in meters. For loblolly pines of given tree height, low-taper trees can withstand greater wind speeds than can high-taper trees—a result that is not surprising because a high taper value indicates a tall, thin tree and hence one that is more susceptible to wind breakage.

Spacing and crown size can also have a large effect on critical wind speed. The model assumes that spacing is equal to crown width. Therefore, tree crowns will meet but not overlap and will leave little space between trees. As the crown area increases, there is more surface area on which the wind can act. Thus, as the crown size grows with tree spacing, the wind force increases and the critical wind speed decreases. Figure 8 is a plot of critical wind speed versus tree spacing and illustrates this effect.

c. Damage patterns

The model produces a graphical plot of downed trees. The arrows are proportional to the tree height and represent trees downed in the instantaneous direction of the wind, because the model assumes no airborne stems or interaction among the trees. The gray dots represent standing trees (Fig. 9).

Suppose, for simplicity, that the radial wind component is held at zero and the maximum absolute wind speed is held constant at 60 m s−1, equivalent to an F2 tornado. For loblolly pine with a site index of 90 and an upwind gap of 10 tree heights at the forest edge, the model produces a series of damage patterns when the forward speed of the vortex is varied.

At the slowest forward speed, a broad area of intermittently downed trees lies nearly centered along the storm track and forms a semicircular pattern, as seen in Fig. 9a (cf. with Atkins and Taber 2004, their Fig. 4). A few downed trees along the left side of the path lie nearly opposite the direction of the storm motion. The general pattern is initially roughly symmetric about the axis of the track. However, as the forward speed increases, as shown in Figs. 9b–d, the damage area moves to the right side of the track and individual trees align increasingly in the direction of the storm motion (e.g., Hall and Brewer 1959, their Fig. 3; Atkins et al. 2005).

As the absolute wind speed increases to F3 strength in Fig. 10, the same general shift from symmetry to asymmetry is observed. As expected, many more trees are downed, but of interest is that progressively fewer trees are affected even though the maximum wind speed is the same for all four examples. The reason for this curious outcome is that if the maximum wind is held constant and the forward speed is increased, the tangential component must decrease. As this vortex component decreases, the initial broad ring of supercritical wind speeds shrinks and eventually shifts to the right of the track and consequently impacts fewer trees. The outcome suggests that the areal extent of damage may decrease considerably with forward speed, even for a given F number. Note, however, that the overall form of the tree damage pattern resembles that of Fig. 9.

Adding the effect of radial winds for a slow-moving storm, one again finds along the left side of the swath a tendency for tree alignment opposite that of the storm direction (Fig. 11). This tendency increases with inflow angle. As the forward motion of the vortex is increased from 3 to 18 m s−1, a new pattern develops in Fig. 12. Here the orientation of tree fall converges along a line on the right side of the track. The convergence line becomes more pronounced and moves closer to the centerline of the track as the radial component increases (cf. Figs. 12a–d). This pattern can be described as a “herringbone” (e.g., Hall and Brewer 1959, their Fig. 5).

To explain the pattern variations, one needs to consider two key controls: 1) the wind field resulting from the vortex and its forward motion and 2) the critical wind speed for the tree species in question. The first control is similar to that found by Letzmann, who illustrated the separate effects of forward speed and vortex inflow (see also Hall and Brewer 1959). However, the second effect—that determined by the tree species—can dramatically change the tree-fall pattern.

Absolute winds in a slowly moving vortex are nearly axially symmetric. If the maximum wind is slightly greater than the tree’s critical speed, the greatest damage will be located at Rmax at the upwind edge of the forest where exposure is greatest. These trees will align with the storm motion vector (see the southern boundary in Fig. 9). As the absolute wind speed increases for this slow-moving vortex, trees will fall as they encounter the approaching ring of maximum winds. The great majority of these trees will fall on their first encounter with the ring’s critical isotachs ahead of and along the sides of the vortex center. Thus, we find the semicircular pattern seen in Figs. 9a, 10a and 11.

As a radial component is added to the vortex, the situation becomes more complicated. Referring back to Fig. 6, we recall that, for a moving vortex, radial winds cause a clockwise shift in the location of the absolute wind maximum. This shift, controlled by α, causes the convergence pattern in the tree alignment observed in Fig. 12. As α increases, trees everywhere to the left of the storm track encounter maximum winds at the left-rear quadrant of the vortex (i.e., after the vortex center has passed). If wind speeds in this rear cusp of the isotach pattern exceed critical values, trees will be downed. At this location, the winds are directed from left to right: hence, the cross-track alignment of downed trees in this region. In a similar way, adjacent trees that stand between the vortex center and the track of the absolute wind maximum will experience maximum wind exposure after the vortex center has passed. These downed trees will also trend from left to right, but they will also tend to align in the forward direction. Because all trees aligned left-to-right fall after the vortex center has passed, their stems should lie atop the others. As forward speed increases, the ultimate extreme of this pattern is a pair of adjacent semicircles: the one on the right is created before the vortex center has passed, and the one on the left is created after the vortex center has passed (see Fig. 12d).

So far, the discussion has been restricted to a single species and site index and has compared damage patterns for fixed maximum absolute wind speeds. Holding the maximum winds constant requires reducing the vortex circulation as storm speed increases. Although this procedure is reasonable if one wishes to compare damage for a constant F number, it unavoidably results in a reduction in damage as the forward speed increases. However, to discern rudimentary controls of damage patterns, it might seem best to keep the vortex circulation constant and to vary the storm speed.

Figure 13 presents the results for an example of such a constant vortex. In this case, the tangential and radial wind maxima are fixed at 70 and 20 m s−1, respectively, yielding α = 16° and a peak vortex wind speed of 73 m s−1—a threshold F3 strength. The vortex is moved at varying speeds through the same stand of loblolly pine with a site index of 90 and an upwind gap of 10 tree heights. Results are generally similar to those for the trials with fixed F number. A crisscross damage pattern emerges to the right of the centerline and becomes more pronounced as the forward speed increases. This time, however, the extent of tree damage increases with forward speed because the absolute wind maximum is allowed to increase. As before, the rightward-oriented tree vectors and damage convergence increase as the radial wind component increases, as shown in Fig. 12.

The effect of critical wind speed is illustrated well by repeating the conditions shown in Fig. 13c but changing the tree species to Scots pine, which has a lower damage threshold. For this run, trees of various ages were randomly distributed through the stand. Results in Fig. 14 show not only more overall tree damage (cf. 36% for loblolly with 81% for Scots pine), but also a clear change in the character of the damage pattern. Few Scots pine survive the initial wind, and, hence, the pattern resembles that of a stronger or slower-moving vortex in a loblolly forest.

To illustrate how the tree-fall patterns vary with site index, an F4 tornado was used. The storm characteristics are identical to those of Fig. 13c. The patterns illustrated in Fig. 15 are similar: a converging pattern of damage forms to the right of the vortex track and trees nearest the upwind edge of the stand are most affected. The only striking difference is that, as the site index increases from 60 to 120, the percentage of downed trees decreases by about 10%. This reduction is simply because a site index of 120 includes stronger, bigger trees.

4. Conclusions

This study describes a physical model that produces patterns of tree fall associated with idealized tornadoes. The tornadic wind fields are specified as steady-state analytical functions that produce lateral forces on idealized trees. Empirical data from loblolly pine and Scots pine are used to simulate stem failure and tree fall for various tree spacings and site indices for a moving modified Rankin vortex with specified radial and tangential components.

It is found that characteristic tree-fall patterns for loblolly pine result from different combinations of tangential and radial components of the vortex and its forward speed. At slow forward speeds and small radial components, vortex flow that exceeds the critical value causes a semicircular damage pattern of tree fall centered along the vortex track in which most trees fall left to right. As the forward motion of the vortex increases, the damage shifts to the right and downed trees begin to align with the track. For a fixed F number, the extent of tree damage may decrease as forward speed increases.

As the radial inflow of the vortex increases, a third damage pattern emerges—one of convergence of tree alignment along a line parallel and to the right of the track. Here a unique crisscross pattern emerges with right-to-left-aligned stems superimposed upon the others. Although site index and absolute wind speed influence the amount of tree damage, they appear to have little effect on the above patterns. However, when the tree species is changed from loblolly to Scots pine, the cross-track damage pattern is found even at rapid forward speed and large radial inflow. This result illustrates the importance of the critical wind speed in determining the damage pattern. It further suggests that in a mixed forest, a combination of patterns may be found within the same damage swath.

This study is an early step in attempting to relate tornadoes to tree damage and as such has many limitations. The type of modified Rankin vortex used in this study likely represents a vast idealization of the actual near-surface flow even for a simple tornado vortex. For example, Doppler measurements by Bluestein et al. (2003) suggest that the maximum inflow in the radial wind component may occur at a radius of about 2 times that of the maximum tangential wind. Measurements by Wurman and Gill (2000) suggest that in the outer region of the vortex the tangential velocity may decrease approximately with r−0.6 rather than r−1 as assumed here. Only a limited number of actual storms have been sampled and thus it may be premature to make generalizations.

A key related question is whether a given tornado vortex flow can be realistically represented by a single set of analytical functions. Most studies and numerous cases of video footage indicate that tornadic winds are anything but steady state. Variations in tangential, radial, and forward speeds are probably large and the vortex flow may make transitions through different flow regimes (Bluestein et al. 2003; Wurman and Gill 2000; Davies-Jones 1986).

Another aspect of actual tornadoes that has been omitted here is that of airborne debris. Some trees and tree sections are picked up and deposited far from their points of origin and do not serve as markers of the instantaneous wind direction. Within the damage swath, such behavior might superimpose randomly oriented debris lying atop the trees that fall in place. This characteristic implies that sole use of areal photography would be open to misinterpretation unless accompanied by ground surveys.

The Peltola model, although based heavily on observed tree response, still represents a first-order parameterized approach. For example, it is likely that tree-stem failure most often results from a harmonic matching of wind gustiness to the frequency of tree sway. This characteristic is parameterized in the model using a gust factor. Uprooting is not included here; nor are second-order effects included. For example, partially uprooted trees often lean on their neighbors, and gaps created by uprooted or broken trees increase the wind exposure of downwind trees. It would be interesting to include some of these effects in future work.

The model uses the logarithmic profile to compute the wind speed at the tree crown from that at some distance above the treetops. Amid the sea of uncertainty, this particular assumption may seem unimportant, but it links the wind force to the tree and should be carefully considered. Studies similar to that of Oliver and Mayhead (1974), who measured wind profiles in and above forests during high winds, may prove to be useful in obtaining profile relationships.

Even in light of uncertainties, the model presented here seems to be a reasonable starting point for simulating tree damage for a simple single-vortex tornado. The model may prove useful to compare the basic simulated damage patterns with those found in the field. For example, several variables, including the tree species, the stand characteristics, and estimates of the forward speed of the vortex, may be obtainable in a given case. Once these variables are specified, observed damage patterns can be compared with those modeled for different radial and tangential components. If successful, the results might yield reliable objective estimates of the storm strength as well as vortex characteristics in forested areas.

Acknowledgments

The authors thank Dr. Jerry Davis for his comments and advice and Dr. Heli Peltola of the University of Joensuu, Finland, for her many helpful suggestions and advice in applying her tree model to this research.

REFERENCES

  • Atkins, N. T., and B. Taber, 2004: Mesoscale analysis of the Cavendish, VT bow-echo tornado on 21 July 2003. Preprints, Conf. on Severe Local Storms, Hyannis, MA, Amer. Meteor. Soc., CD-ROM, P4.2.

  • Atkins, N. T., C. S. Bouchard, R. W. Przybylinski, R. J. Trapp, and G. Schmocker, 2005: Damaging surface wind mechanisms within the 10 June 2003 Saint Louis bow echo during BAMEX. Mon. Wea. Rev., 133 , 22752296.

    • Search Google Scholar
    • Export Citation
  • Avery, T. E., and H. E. Burkhart, 2002: Forest Measurements. 5th ed. McGraw-Hill Higher Education, 456 pp.

  • Baldwin Jr., V. C., 1987: Green and dry-weight equations for above ground components of planted loblolly pine trees in the west gulf region. South. J. Appl. For., 11 , 212218.

    • Search Google Scholar
    • Export Citation
  • Bluestein, H. B., W-C. Lee, M. Bell, C. C. Weiss, and A. L. Pazmany, 2003: Mobile Doppler radar observations of a tornado in a supercell near Bassett, Nebraska, on 5 June 1999. Part II: Tornado-vortex structure. Mon. Wea. Rev., 131 , 29682984.

    • Search Google Scholar
    • Export Citation
  • Clark III, A., 1991: Weight and volume equations for loblolly pine in eastern North Carolina. USDA Forest Service Southeastern Forest Experiment Station and Federal Paper Board Co., Inc., Cooperative Study Final Rep., 57 pp.

  • Davies-Jones, R. P., 1986: Tornado dynamics. Thunderstorm Morphology and Dynamics, 2d ed. E. Kessler, Ed., University of Oklahoma Press, 197–236.

    • Search Google Scholar
    • Export Citation
  • Fredericksen, T. S., R. L. Hedden, and S. A. Williams, 1993: Testing loblolly pine wind firmness with simulated wind stress. Can. J. For. Res., 23 , 17601765.

    • Search Google Scholar
    • Export Citation
  • Fujita, T. T., 1971: Proposed characterization of tornadoes and hurricanes by area and intensity. University of Chicago Satellite and Meteorology Research Paper 91, 15 pp.

  • Gardiner, B. A., and G. R. Stacey, 1995: Designing forest edges to improve wind stability. Forestry Commission Tech. Paper 16, 8 pp.

  • Gardiner, B. A., G. R. Stacey, R. E. Belcher, and C. J. Wood, 1997: Field and wind tunnel assessments of the implications of re-spacing and thinning for tree stability. Forestry, 70 , 233252.

    • Search Google Scholar
    • Export Citation
  • Grace, J., 1977: Plant Response to Wind. Academic Press, 204 pp.

  • Hakkila, P., 1971: Coniferous branches as a raw material source. A subproject of the joint Nordic research programme for the utilization of logging residues. Communicationes Instituti Forestalis Fenniae Rep. 75, 60 pp.

  • Hall, F., and R. D. Brewer, 1959: A sequence of tornado damage patterns. Mon. Wea. Rev., 87 , 207216.

  • Hedden, R. L., T. S. Fredericksen, and S. A. Williams, 1995: Modeling the effect of crown shedding and streamlining on the survival of loblolly pine exposed to acute wind. Can. J. For. Res., 25 , 704712.

    • Search Google Scholar
    • Export Citation
  • Jones, H. G., 1983: Plants and Microclimate. A Quantitative Approach to Environmental Plant Physiology. Cambridge University Press, 323 pp.

    • Search Google Scholar
    • Export Citation
  • Laasasenaho, J., 1982: Taper curve and volume functions for pine, spruce and birch. Communicationes Instituti Forestalis Fenniae Rep. 108, 74 pp.

  • Marshall, T. P., 2003: The La Plata, MD tornado: Issues regarding the F-scale. Preprints, Symp. on the F-Scale and Severe-Weather Damage Assessment, Long Beach, CA, Amer. Meteor. Soc., CD-ROM, P1.6.

  • Mattheck, C., and H. Breloer, 1994: The Body Language of Trees: A Handbook for Failure Analysis. HMSO, 240 pp.

  • Monteith, J. H., 1975: Principles of Environmental Physics. 1st ed. Edward Arnold, 241 pp.

  • Morgan, J., and M. G. R. Cannell, 1994: Shape of tree stems: A reexamination of the uniform stress hypothesis. Tree Physiol., 5 , 6374.

    • Search Google Scholar
    • Export Citation
  • Oliver, H. R., and G. J. Mayhead, 1974: Wind measurements in a pine forest during a destructive gale. Forestry, 47 , 185194.

  • Peltola, H., 1996: Model computations on wind flow and turning moment by wind for Scots pine along the margins of clear-cut areas. For. Ecol. Manage., 83 , 203215.

    • Search Google Scholar
    • Export Citation
  • Peltola, H., and S. Kellomäki, 1993: A mechanistic model for calculating windthrow and stem breakage of Scots pine at stand edge. Silva Fennica, 27 , 99111.

    • Search Google Scholar
    • Export Citation
  • Peltola, H., M-L. Nykänen, and S. Kellomäki, 1997: Model computations on the critical combination of snow loading and windspeed for snow damage of Scots pine, Norway spruce and birch sp. at stand edge. For. Ecol. Manage., 95 , 229241.

    • Search Google Scholar
    • Export Citation
  • Peltola, H., S. Kellomäki, H. Välsänen, and V-P. Ikonen, 1999: A mechanistic model for assessing the risk of wind and snow damage to single trees and stands of Scots pine, Norway spruce and birch. Can. J. For. Res., 29 , 647661.

    • Search Google Scholar
    • Export Citation
  • Peterson, C. J., 2003: Factors influencing treefall risk in tornadoes in natural forests. Preprints, Symp. on the F-Scale and Severe-Weather Damage Assessment, Long Beach, CA, Amer. Meteor. Soc., CD-ROM, 3.1.

  • Peterson, R. E., 1992: Johannes Letzmann: A pioneer in the study of tornadoes. Wea. Forecasting, 7 , 166184.

  • Petty, J. A., and C. Swain, 1985: Factors influencing stem breakage of conifers in high winds. Forestry, 58 , 7585.

  • Schaefer, J. T., and R. L. Livingston, 2003: The consistency of F-scale ratings. Preprints, Symp. on the F-Scale and Severe-Weather Damage Assessment, Long Beach, CA, Amer. Meteor. Soc., CD-ROM, P1.1.

  • Turner, J. S., 1966: The constraints imposed on tomado-like vortices by the top and bottom boundary conditions. J. Fluid Mech., 25 , 377400.

    • Search Google Scholar
    • Export Citation
  • Wahlenberg, W. G., 1960: Loblolly Pine: Its Use, Ecology, Regeneration, Protection, Growth and Management. The School of Forestry, Duke University, 603 pp.

    • Search Google Scholar
    • Export Citation
  • Ward, N. B., 1972: The exploration of certain features of tornado dynamics using a laboratory model. J. Atmos. Sci., 29 , 11941204.

  • Wood, C. J., 1995: Understanding wind forces on trees. Wind and Trees, M. P. Coutts and J. Grace, Eds., Cambridge University Press, 133–164.

    • Search Google Scholar
    • Export Citation
  • Wurman, J., and S. Gill, 2000: Finescale radar observations of the Dimmitt, Texas (2 June 1995), tornado. Mon. Wea. Rev., 128 , 21352164.

    • Search Google Scholar
    • Export Citation
  • Zrnić, D. S., D. W. Burgess, and L. Hennington, 1985: Doppler spectra and estimated windspeed of a violent tornado. J. Climate Appl. Meteor., 24 , 10681081.

    • Search Google Scholar
    • Export Citation

Fig. 1.
Fig. 1.

Wind vectors and isotachs (shaded) representing the total wind field for a vortex with a maximum tangential velocity of 60 m s−1, a maximum radial velocity of 30 m s−1, and a forward speed of 15 m s−1. The radius of maximum vortex-relative components is 75 m. The center of the relative vortex is at the origin.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 2.
Fig. 2.

Wind force (Fwind) and gravity (Fmass) producing turning moment (T) and lateral displacement (x) on an idealized tree (after Peltola et al. 1997).

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 3.
Fig. 3.

Variation of critical turning moment with distance from stand edge for Scot pine.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 4.
Fig. 4.

Critical wind speed for stem breakage for different upwind gap sizes for Scots pine.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 5.
Fig. 5.

Rudimentary relationships controlling the isotach pattern for a moving Rankine vortex. Isotachs are curved lines (m s−1); the vortex is centered at the origin but is moving at velocity VS. Line AB intersects the calm point at radius ro, the origin, and the location X of maximum absolute wind Vx. Vectors Vt (Vtan in text) and Vr denote the tangential and radial vortex components at X, respectively.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 6.
Fig. 6.

As in Fig. 1, but for selected combinations of radial component Vr, tangential component Vt, and forward velocity VS. Note the AB axis of extremes.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 7.
Fig. 7.

Variation of critical wind speed with height and taper for loblolly pine.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 8.
Fig. 8.

Variation of critical wind speeds with tree spacing for loblolly pines with two different taper values.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 9.
Fig. 9.

Windfall patterns resulting from varying the tangential component and forward speed for an F2 tornado of fixed intensity. Values of the radial component Vr, tangential component Vt, forward motion of the vortex VS, and maximum wind Vx (m s−1) are given for each panel. The vortex center moves northward at x = 0. For all cases, the radius of maximum vortex winds is 75 m. Standing loblolly pine trees are represented by dots, and fallen trees are represented by arrows whose lengths represent the lengths of stems that were broken (see scale) and whose orientations represent stem orientations. The site index is 90, and the upwind gap (at the south edge of the stand) is 10 tree heights. The percentage for each panel refers to the percent of trees felled out of the total of 7056.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 10.
Fig. 10.

Same as Fig. 9, but for an F3 tornado.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 11.
Fig. 11.

Same as Fig. 9, but for a slow-moving F3 tornado.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 12.
Fig. 12.

Same as Fig. 9, but for a faster-moving F3 tornado.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 13.
Fig. 13.

Same as Fig. 9, but for a fixed vortex moving at increasing forward speed, resulting in increasing intensity.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 14.
Fig. 14.

Damage patterns for wind conditions identical to those of Fig. 13c, but for Scots pine with a random distribution of different aged trees. This combination resulted in 81% breakage of 12 100 trees.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Fig. 15.
Fig. 15.

Damage patterns for wind conditions identical to those of Fig. 13c, but for loblolly pine with site indices of (a) 60 and (b) 120. These site indices represent totals of 13 689 and 5184 trees, respectively, with losses of 38% and 28%.

Citation: Journal of Applied Meteorology and Climatology 45, 12; 10.1175/JAM2413.1

Save
  • Atkins, N. T., and B. Taber, 2004: Mesoscale analysis of the Cavendish, VT bow-echo tornado on 21 July 2003. Preprints, Conf. on Severe Local Storms, Hyannis, MA, Amer. Meteor. Soc., CD-ROM, P4.2.

  • Atkins, N. T., C. S. Bouchard, R. W. Przybylinski, R. J. Trapp, and G. Schmocker, 2005: Damaging surface wind mechanisms within the 10 June 2003 Saint Louis bow echo during BAMEX. Mon. Wea. Rev., 133 , 22752296.

    • Search Google Scholar
    • Export Citation
  • Avery, T. E., and H. E. Burkhart, 2002: Forest Measurements. 5th ed. McGraw-Hill Higher Education, 456 pp.

  • Baldwin Jr., V. C., 1987: Green and dry-weight equations for above ground components of planted loblolly pine trees in the west gulf region. South. J. Appl. For., 11 , 212218.

    • Search Google Scholar
    • Export Citation
  • Bluestein, H. B., W-C. Lee, M. Bell, C. C. Weiss, and A. L. Pazmany, 2003: Mobile Doppler radar observations of a tornado in a supercell near Bassett, Nebraska, on 5 June 1999. Part II: Tornado-vortex structure. Mon. Wea. Rev., 131 , 29682984.

    • Search Google Scholar
    • Export Citation
  • Clark III, A., 1991: Weight and volume equations for loblolly pine in eastern North Carolina. USDA Forest Service Southeastern Forest Experiment Station and Federal Paper Board Co., Inc., Cooperative Study Final Rep., 57 pp.

  • Davies-Jones, R. P., 1986: Tornado dynamics. Thunderstorm Morphology and Dynamics, 2d ed. E. Kessler, Ed., University of Oklahoma Press, 197–236.

    • Search Google Scholar
    • Export Citation
  • Fredericksen, T. S., R. L. Hedden, and S. A. Williams, 1993: Testing loblolly pine wind firmness with simulated wind stress. Can. J. For. Res., 23 , 17601765.

    • Search Google Scholar
    • Export Citation
  • Fujita, T. T., 1971: Proposed characterization of tornadoes and hurricanes by area and intensity. University of Chicago Satellite and Meteorology Research Paper 91, 15 pp.

  • Gardiner, B. A., and G. R. Stacey, 1995: Designing forest edges to improve wind stability. Forestry Commission Tech. Paper 16, 8 pp.

  • Gardiner, B. A., G. R. Stacey, R. E. Belcher, and C. J. Wood, 1997: Field and wind tunnel assessments of the implications of re-spacing and thinning for tree stability. Forestry, 70 , 233252.

    • Search Google Scholar
    • Export Citation
  • Grace, J., 1977: Plant Response to Wind. Academic Press, 204 pp.

  • Hakkila, P., 1971: Coniferous branches as a raw material source. A subproject of the joint Nordic research programme for the utilization of logging residues. Communicationes Instituti Forestalis Fenniae Rep. 75, 60 pp.

  • Hall, F., and R. D. Brewer, 1959: A sequence of tornado damage patterns. Mon. Wea. Rev., 87 , 207216.

  • Hedden, R. L., T. S. Fredericksen, and S. A. Williams, 1995: Modeling the effect of crown shedding and streamlining on the survival of loblolly pine exposed to acute wind. Can. J. For. Res., 25 , 704712.

    • Search Google Scholar
    • Export Citation
  • Jones, H. G., 1983: Plants and Microclimate. A Quantitative Approach to Environmental Plant Physiology. Cambridge University Press, 323 pp.

    • Search Google Scholar
    • Export Citation
  • Laasasenaho, J., 1982: Taper curve and volume functions for pine, spruce and birch. Communicationes Instituti Forestalis Fenniae Rep. 108, 74 pp.

  • Marshall, T. P., 2003: The La Plata, MD tornado: Issues regarding the F-scale. Preprints, Symp. on the F-Scale and Severe-Weather Damage Assessment, Long Beach, CA, Amer. Meteor. Soc., CD-ROM, P1.6.

  • Mattheck, C., and H. Breloer, 1994: The Body Language of Trees: A Handbook for Failure Analysis. HMSO, 240 pp.

  • Monteith, J. H., 1975: Principles of Environmental Physics. 1st ed. Edward Arnold, 241 pp.

  • Morgan, J., and M. G. R. Cannell, 1994: Shape of tree stems: A reexamination of the uniform stress hypothesis. Tree Physiol., 5 , 6374.

    • Search Google Scholar
    • Export Citation
  • Oliver, H. R., and G. J. Mayhead, 1974: Wind measurements in a pine forest during a destructive gale. Forestry, 47 , 185194.

  • Peltola, H., 1996: Model computations on wind flow and turning moment by wind for Scots pine along the margins of clear-cut areas. For. Ecol. Manage., 83 , 203215.

    • Search Google Scholar
    • Export Citation
  • Peltola, H., and S. Kellomäki, 1993: A mechanistic model for calculating windthrow and stem breakage of Scots pine at stand edge. Silva Fennica, 27 , 99111.

    • Search Google Scholar
    • Export Citation
  • Peltola, H., M-L. Nykänen, and S. Kellomäki, 1997: Model computations on the critical combination of snow loading and windspeed for snow damage of Scots pine, Norway spruce and birch sp. at stand edge. For. Ecol. Manage., 95 , 229241.

    • Search Google Scholar
    • Export Citation
  • Peltola, H., S. Kellomäki, H. Välsänen, and V-P. Ikonen, 1999: A mechanistic model for assessing the risk of wind and snow damage to single trees and stands of Scots pine, Norway spruce and birch. Can. J. For. Res., 29 , 647661.

    • Search Google Scholar
    • Export Citation
  • Peterson, C. J., 2003: Factors influencing treefall risk in tornadoes in natural forests. Preprints, Symp. on the F-Scale and Severe-Weather Damage Assessment, Long Beach, CA, Amer. Meteor. Soc., CD-ROM, 3.1.

  • Peterson, R. E., 1992: Johannes Letzmann: A pioneer in the study of tornadoes. Wea. Forecasting, 7 , 166184.

  • Petty, J. A., and C. Swain, 1985: Factors influencing stem breakage of conifers in high winds. Forestry, 58 , 7585.

  • Schaefer, J. T., and R. L. Livingston, 2003: The consistency of F-scale ratings. Preprints, Symp. on the F-Scale and Severe-Weather Damage Assessment, Long Beach, CA, Amer. Meteor. Soc., CD-ROM, P1.1.

  • Turner, J. S., 1966: The constraints imposed on tomado-like vortices by the top and bottom boundary conditions. J. Fluid Mech., 25 , 377400.

    • Search Google Scholar
    • Export Citation
  • Wahlenberg, W. G., 1960: Loblolly Pine: Its Use, Ecology, Regeneration, Protection, Growth and Management. The School of Forestry, Duke University, 603 pp.

    • Search Google Scholar
    • Export Citation
  • Ward, N. B., 1972: The exploration of certain features of tornado dynamics using a laboratory model. J. Atmos. Sci., 29 , 11941204.

  • Wood, C. J., 1995: Understanding wind forces on trees. Wind and Trees, M. P. Coutts and J. Grace, Eds., Cambridge University Press, 133–164.

    • Search Google Scholar
    • Export Citation
  • Wurman, J., and S. Gill, 2000: Finescale radar observations of the Dimmitt, Texas (2 June 1995), tornado. Mon. Wea. Rev., 128 , 21352164.

    • Search Google Scholar
    • Export Citation
  • Zrnić, D. S., D. W. Burgess, and L. Hennington, 1985: Doppler spectra and estimated windspeed of a violent tornado. J. Climate Appl. Meteor., 24 , 10681081.

    • Search Google Scholar
    • Export Citation
  • Fig. 1.

    Wind vectors and isotachs (shaded) representing the total wind field for a vortex with a maximum tangential velocity of 60 m s−1, a maximum radial velocity of 30 m s−1, and a forward speed of 15 m s−1. The radius of maximum vortex-relative components is 75 m. The center of the relative vortex is at the origin.

  • Fig. 2.

    Wind force (Fwind) and gravity (Fmass) producing turning moment (T) and lateral displacement (x) on an idealized tree (after Peltola et al. 1997).

  • Fig. 3.

    Variation of critical turning moment with distance from stand edge for Scot pine.

  • Fig. 4.

    Critical wind speed for stem breakage for different upwind gap sizes for Scots pine.

  • Fig. 5.

    Rudimentary relationships controlling the isotach pattern for a moving Rankine vortex. Isotachs are curved lines (m s−1); the vortex is centered at the origin but is moving at velocity VS. Line AB intersects the calm point at radius ro, the origin, and the location X of maximum absolute wind Vx. Vectors Vt (Vtan in text) and Vr denote the tangential and radial vortex components at X, respectively.

  • Fig. 6.

    As in Fig. 1, but for selected combinations of radial component Vr, tangential component Vt, and forward velocity VS. Note the AB axis of extremes.

  • Fig. 7.

    Variation of critical wind speed with height and taper for loblolly pine.

  • Fig. 8.

    Variation of critical wind speeds with tree spacing for loblolly pines with two different taper values.

  • Fig. 9.

    Windfall patterns resulting from varying the tangential component and forward speed for an F2 tornado of fixed intensity. Values of the radial component Vr, tangential component Vt, forward motion of the vortex VS, and maximum wind Vx (m s−1) are given for each panel. The vortex center moves northward at x = 0. For all cases, the radius of maximum vortex winds is 75 m. Standing loblolly pine trees are represented by dots, and fallen trees are represented by arrows whose lengths represent the lengths of stems that were broken (see scale) and whose orientations represent stem orientations. The site index is 90, and the upwind gap (at the south edge of the stand) is 10 tree heights. The percentage for each panel refers to the percent of trees felled out of the total of 7056.

  • Fig. 10.

    Same as Fig. 9, but for an F3 tornado.

  • Fig. 11.

    Same as Fig. 9, but for a slow-moving F3 tornado.

  • Fig. 12.

    Same as Fig. 9, but for a faster-moving F3 tornado.

  • Fig. 13.

    Same as Fig. 9, but for a fixed vortex moving at increasing forward speed, resulting in increasing intensity.

  • Fig. 14.

    Damage patterns for wind conditions identical to those of Fig. 13c, but for Scots pine with a random distribution of different aged trees. This combination resulted in 81% breakage of 12 100 trees.

  • Fig. 15.

    Damage patterns for wind conditions identical to those of Fig. 13c, but for loblolly pine with site indices of (a) 60 and (b) 120. These site indices represent totals of 13 689 and 5184 trees, respectively, with losses of 38% and 28%.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 623 207 71
PDF Downloads 465 137 18