|
|
||||||||
Structure and Development |
Department of Plant Biology, Cornell University, Ithaca, New York 14853-5908 USA
Received for publication August 27, 2002. Accepted for publication October 25, 2002.
| ABSTRACT |
|---|
|
|
|---|
L4/3 and EI
L11/4
Ms5/3, whereas Mf /Ms
L1.15, i.e., peduncles disproportionately increase in their biomass as they increase in length, but mechanically support a disproportionately smaller floral biomass relative to their biomass. Calculations show that the tall peduncles from wind-sheltered sites have a larger fruit dispersal range and a lower factor of safety than the shorter peduncles produced in open sites. These and other observations are interpreted to indicate that tubular peduncles enhance relative fitness in terms of propagule dispersal (but not propagule number per stem) while maintaining a sufficient factor of safety against mechanical failure.
Key Words: allometry Asteraceae Brazier buckling flexural stiffness fruit dispersal Hieracium pilosella plant biomechanics winding loading
| INTRODUCTION |
|---|
|
|
|---|
These trade-offs may be evaluated theoretically in terms of the allometric (scaling) relationships among stem diameter D, biomass (mass) M, and flexural stiffness EI with respect to length L, because, for any stem biomass investment, mechanical theory predicts very specific scaling exponents (Table 1). For example, three models for stem taper (geometric, elastic, and stress self-similarity; see Appendix), one of which is equivalent to the Euler buckling formula (elastic self-similarity; see McMahon, 1973
), predict very different exponents for the relationship between stem diameter and length (1, 3/2, and 2, respectively), stem mass (proportional to D2L) and length (3, 4, and 5, respectively), and stem axial second moment of area (proportional to D4) and length (4, 6, and 8, respectively).
|
DLu2). Assuming that stems are "designed" to maintain a constant stress at their base
b such that
b
F/D2
constant, it follows that D
Lu2. If ambient wind speeds scale as the 1/2- or 1-power of length, stem diameter is predicted to scale as 2- or 3-power of length, respectively, whereas much larger exponents are predicted for the scaling of stem biomass or second moment of area with respect to length (Table 1).
Hollow beams with very thin walls are also susceptible to Brazier buckling, i.e., localized ovalization and crimping (Niklas, 1989
, 1990
, 1992
; Spatz et al., 1997
). Engineering theory shows that the critical Brazier bending moment MB is proportional to D3 for thin-walled tubes. For wind-loaded stems, this bending moment is proportional to the product of drag and stem length, i.e., FL
D3. Thus, if wind speeds scale as the 1/2- or 1-power of stem length, then stem diameter is predicted to scale as the 3/2- or 2-power of length (Table 1).
These theoretical expectations are based on a number of assumptions predicated on the supposition that self-loading and wind-induced drag forces govern the allometry of peduncles. Nevertheless, empirically determined scaling exponents, significantly lower or higher than those predicted for the relationships among stem diameter, biomass, and second moment of area, would provide circumstantial evidence that trade-offs among economy, dispersal, and safety have occurred.
This study was designed to examine these trade-offs empirically by comparing the morphometry and mechanical allometry of the vertical flower stalks (peduncles) of Hieracium pilosella drawn from populations found in windy and sheltered sites against the scaling relationships predicted by mechanical allometric theory. Hieracium pilosella was selected because its peduncles are hollow and non-septate, its aggregated (composite) flowers are borne distally, and its dry fruits (achenes) are wind dispersed. To develop a sense of the interspecific allometric trends, additional species in the Asteraceae with similar peduncle morphology were studied (i.e., H. aurantiacum, Taraxicum officinale, Tragopogon pratensis).
| MATERIALS AND METHODS |
|---|
|
|
|---|
Reciprocal transplants were made between pop. 1 and pop. 3 (in May of the year preceding this study to maximize transplant acclimation to new environmental conditions) to determine whether genotypic or environmental differences affected peduncle morphometry or mechanical properties. Twenty plants of comparable size (rosette diameters) from pop. 1 and pop. 3 were randomly selected along a 6-m linear transect, removed with as much of their roots as possible, and randomly planted among members of the other population. Transplants were watered for 2 wk and then left undisturbed until the following year during which the herbaceous vegetation surrounding pops. 12 and pop. 3 was regularly mowed or hand clipped to an average height of 10 cm and 2 cm, respectively. All transplants from pop. 1 to pop. 3 failed to flower or died for unknown reasons. Therefore, data are reported for plants transplanted from pop. 3 to pop. 1 (designated as Tr. 3-1; see Table 3).
|
A total of 251 peduncles was collected. Each peduncle was cut at its base at ground level and placed in water to reduce wilting. All morphometric measurements and mechanical tests were performed within 24 h of collection. Hieracium pilosella peduncles (n = 195) were assigned an identification number that did not identify the population from which it was collected to provide double-blind comparisons of peduncle morphometry and mechanical properties among pops. 13.
The flowers of each peduncle were removed and weighed. The remaining stem was segmented to determine its outer and inner diameters along its length using a handheld microcaliper. Segments were also weighed to determine fresh mass per unit length. The inner and outer diameters of segments were used to calculate peduncle wall thickness. The peduncles of some species were branched. The dimensions and masses of branches were also measured and recorded (n = 100). The mass of lateral branches was added to that of the main stems to calculate loading conditions at the base of each stem.
Mechanical tests
The composite-tissue elastic modulus E of the basal 8 cm long segment of representative H. pilosella peduncles was determined by means of three-point-bending tests using an Instron 4502 testing machine. Before each bending test, the outer and inner diameters of each specimen were measured at both ends using a handheld electronic microcaliper to determine an average axial second moment of area I using the formula I = (Ib + Ia)/ 2 where I =
(do4 di4)/64, subscripts b and a refer to the basal and apical portions of the peduncle, and do and di are averaged outer and inner diameters (measured at two locations), respectively. Alternative methods for calculating I were explored; the method used was least sensitive to error.
Each basal segment was then suspended horizontally between two vertical supports 5 cm apart and bent by applying a blunt-edge loading point at its mid-length at a strain rate of 25 mm/min. No transverse compression at the supported ends of peduncles was observed during loading experiments. In each case, the length to outer average diameter of each specimen
20. Ten data points per second were recorded continuously during each test to construct a force vs. displacement graph. A bending test was terminated once the applied (maximum) force Fmax exceeded the proportional limit of the peduncle; Fmax was determined by visual inspection of the force vs. displacement graph. The yield stress
y (the onset of plastic deformation) at this stage of testing was calculated using the formula
y = 16 LdoFmax/[
(do4 di4)]. The composite-tissue elastic modulus of each basal segment was calculated using the formula E = 4mL3/[3
(do4 di4)], where m is the slope of the initial (linear) portion of the force vs. displacement graph and L is the distance between the two vertical supports (=5 cm).
Simulated wind speed profiles and estimated fruit dispersal distances
Fruit dispersal distance was estimated based on the average wind speed utop recorded at 0.5 m above ground within pop. 1 and pop. 3 when flower mass was maximum (June). Wind speeds were simultaneously measured in the field at 0.5-h intervals for 48 h using Sierra-Misco model WSD335 anemometers and wind vanes and Omnidata model EL824-MS data loggers. Each wind vane was located 0.5 m above ground within each population.
A logarithmic vertical wind speed profile was simulated using the formula ui = [utop/ln(h/zo)]ln[(h xi)/zo], where ui is the estimated wind speed at any distance xi from where utop was measured, h equals 0.5 m, and zo is the roughness parameter (which was taken as 0.05 times the height of surrounding vegetation; see Niklas and Spatz, 2000
). The surrounding vegetation of pop. 1 and pop. 3 was 0.20 m and 0.10 m, respectively; therefore, zo was taken as 0.01 and 0.005, respectively.
Fruit dispersal distance y was estimated using the formula y = Luia/uT, where L is the average peduncle length for each population, uia is the averaged wind speed between L and ground level (where ui = 0), and uT is the terminal settling velocity of fruits (see Okubo and Levin, 1989
; Niklas, 1994
). Values for uia were determined from simulated wind speed profiles for pop. 1 and pop. 3 (see Fig. 6). Values for uT were determined by dropping 10 representative fruits from pop. 1 and 3 into a cylinder measuring 1.5 m in length and 0.20 m in diameter fitted with two orthogonally positioned holes located 0.20 m above the cylinder base to provide access for a stroboscopic light and a camera lens. Photographic prints of descending stroboscopically illuminated fruits (500 flashes/s) were measured to calculate rates of descent (see Niklas, 1992
, pp. 459461).
|
y/
c, and (3) the yield stress of stem tissues divided by the bending stress generated by wind drag measured at the peduncle base,
y/
W. All of these parameters were directly measured based on morphometric and mechanical testing data with the exception of
W. The wind-induced bending stress was estimated based on the wind speed profiles measured for the sheltered and open populations in the following way. The profile for each population was used to calculate the mean wind speeds along the length of each peduncle. The drag force per unit length was then calculated along the length of each peduncle based on stem and flower projected areas. This force was then integrated along stem length to calculate the bending moment at the base of each peduncle. Finally, this bending moment MW was then used to calculate
W using the formula
W = MWDb/2IB , where Db is the external diameter and Ib is the second moment of area measured at the base of the peduncle.
Statistical analyses
Model Type II (reduced major axis) regression analysis was used to determine the scaling exponents and allometric constants (regression slopes and Y-intercepts,
RMA and ßRMA, respectively) of pair-wise comparisons of log10-transformed data. This protocol is recommended when functional rather than predictive relationships are sought among variables that are biologically interdependent and subject to unknown measurement error (Sokal and Rohlf, 1981
; Niklas, 1994
).
All statistical analyses used the formulas log Y2 = log ßRMA +
RMAlogY1, where Y2 and Y1 are interdependent variables (e.g., reproductive and standing leaf biomass per plant),
RMA =
OLS/r, where
OLS and r are the slope and correlation coefficient determined from ordinary least squares (Model Type I) regression analysis, and log ßRMA = log(cap Y2)
RMAlog(cap Y1), where (cap Y) denotes the mean of variable Y. The 95% confidence intervals for ßRMA values were computed based on the corresponding 95% CI values of
RMA (Sokal and Rohlf, 1981
).
| RESULTS |
|---|
|
|
|---|
|
0.10 even for beams constructed of anisotropic materials (Ades, 1957
Mechanical properties
The axial second moments of area I at the base of H. pilosella peduncles were statistically significantly correlated with peduncle length L, whereas the composite-tissue elastic modulus E was not (Fig. 2AB). Nevertheless, the flexural rigidity EI and thus the ability to resist bending forces at the base of peduncles correlated significantly with L (Fig. 2C). Across all peduncles, I scaled roughly as the 2.23-power of L, whereas EI scaled as the 2.75-power of L and as the 1.66-power of peduncle biomass Ms (Table 2). No statistically significant differences were evident among H. pilosella growing in different populations.
|
|
y measured at the base of peduncles was noticeably different between transplants and non-transplants (Fig. 2D). A strong correlation was observed between E and
y for the peduncles of non-transplants; these two mechanical properties were poorly correlated among the peduncles of transplants (Table 2). We have no explanation for this difference.
Allometric relations
Hieracium pilosella peduncle biomass MS scaled as the 1.36-power of peduncle length L. Based on comparisons of the 95% confidence intervals of the scaling exponents (
RMA), this relationship held for transplants and non-transplants (Fig. 3A; Table 3). Across all species, including H. pilosella, MS scaled as the 1.56-power of L (r2 = 0.893, n = 335, F = 2794, P < 0.0001) (Fig. 3B).
|
Hieracium pilosella inflorescence biomass relative to peduncle biomass, i.e., Mf /MS, decreased with increasing peduncle length (Fig. 4A). Mf /MS scaled as the 1.15-power of L (Table 3). The scaling exponent and Y-intercept for this relationship were statistically invariant across pops. 13, and, based on the allometry of transplants, not affected by wind conditions (Table 3). Because H. pilosella inflorescence biomass was poorly correlated with peduncle length (0.14
r2
0.22), we interpret these data to indicate that taller and more massive peduncles do not provide an advantage in terms of supporting a larger number of progeny per stem.
|
Wind speeds and estimated fruit dispersal distances
The mean ± SE wind speed measured at 0.5 above ground for pop. 1 and pop. 3 was 0.288 ± 0.043 m/s and 0.518 ± 0.023 m/s, respectively; maximum (minimum) wind speeds measured were 1.1 (0.01) m/s and 2.2 (0.01) m/s, respectively. Based on a logarithmic vertical wind speed profile and the average peduncle lengths, wind speeds at the top of peduncles in pop. 1 and pop. 3 were 0.26 m and 0.24 m/s, and 0.12 m and 0. 36 m/s, respectively (Fig. 5). The average wind speed from the top of peduncles to ground level in pop. 1 and pop. 3 (i.e., uia) was 0.18 m/s ad 0.27 m/s, respectively. The average terminal settling velocity of mature fruits did not differ between the two populations (uT = 0.013 ± 0.006 m/s).
|
Factors of safety
Three criteria for stem mechanical failure were selected to calculate factors of safety for stem mechanical failure: the quotient of the Euler critical buckling biomass and floral biomass, Pcr /Mf , the quotient of the stem tissue yield stress and the stem compressive stress,
y /
c, and the quotient of the stem tissue yield stress and the wind-induced bending stress at the base of the stem,
y /
W. For each of these criteria, all peduncles had factors of safety greater than unity (Fig. 6). On average, the peduncles of plants growing in windy sites had significantly higher factors of safety than those produced by wind-sheltered plants (P < 0.001). The lowest factors of safety were those sheltered peduncles computed for
y/
W.
| DISCUSSION |
|---|
|
|
|---|
In terms of elevating reproductive organs, it is easy to show that tapered stems can elevate flowers and fruits higher than their untapered counterparts with an equivalent biomass, basal diameter, and bulk material (tissue) density. For example, if Lstc denotes the length of the solid truncated conical stem with apical and basal diameters Da and Db and if Lcy is the length of its cylindrical counterpart with an equivalent Db and bulk density, then Lstc will be 2.4 times longer than Lcy when Da /Db = 0.2 (see Appendix). Hollow, tapered peduncles are even more cost-effective. For example, if Lhtc denotes the length of a hollow truncated conical stem with wall thickness t, then Lhtc will be 17 times longer than Lstc when t/Db = 0.01 (see Appendix).
However, Euler's formula for elastic buckling and the formula for the critical Brazier bending moment show that trade-offs exist between the potential for mechanical failure and the extent to which a distal load can be elevated (see Wainwright et al., 1976
; Niklas, 1992
; Speck et al., 1990a
, b
). Euler's formula states that the critical buckling load Pcr of a cylinder with flexural rigidity EI, diameter D, and length L is proportional to EI/L2. Because I is proportional to D4 and E is constant for any tissue construction, we see that Pcr
L2/r4, where r = L/D, i.e., for any biomass investment, the critical buckling load increases rapidly as a stem elevates a distal load. Likewise, the critical Brazier buckling moment Mbr for a hollow tube is proportional to EDt2, where t is wall thickness (see Ades, 1957
). Because Mbr is proportional to the product of a bending force F and some function of tube length L and because the mass of a tube is proportional to DLt, we see that F
Et/L2, i.e., for any mass investment, the force required to locally crimp a thin-walled, tubular stem decreases dramatically with increasing length and decreasing wall thickness.
The susceptibility of tubular peduncles to mechanical failure is exacerbated under windy conditions because plant aerial organs experience wind-induced drag forces in addition to the bending forces incurred by their own mass and the mass of the flowers they support (see Vogel, 1981
; Niklas and Spatz, 2000
; Spatz and Brüchert, 2000
). Because mechanical forces are additive, either stem length must be reduced or wall thickness must be increased to maintain a sufficient factor of safety against failure. The allometry of H. pilosella peduncles complies with this expectation. Plants growing in open sites produce shorter peduncles with significantly higher factors of safety than those produced in wind-sheltered habitats.
Nonetheless, the factor of safety steadily erodes with increasing peduncle length and increasing fruit dispersal range, which is mechanically tenable for two reasons. First, taller peduncles do not support more massive inflorescences than their shorter counterparts, and, second, calculations indicate that the taller peduncles produced in wind-sheltered sites experience lower average wind speeds and thus smaller drag forces. This phenomenology may not hold true for all populations. Our measurements of ambient wind conditions are limited and thus idiosyncratic. However, for the habitats studied here, our data indicate that the potential for wind-induced mechanical stem failure is balanced against the benefits of propagule dispersal. Thus, the allometry and mechanics of hollow non-septate peduncles reflect a number of trade-offs that are reconciled in ways that are likely to be reproductively beneficial.
Finally, we draw attention to the radical departures between the scaling exponents predicted by simple mechanical theory for stems (Table 1) and those empirically observed for peduncles (Tables 23). For example, both elastic self-similarity (which follows directly from Euler's formula) and Brazier buckling theory (for wind-loaded peduncles) predict that stem biomass and second moment of area should scale as the 4- and 6-power of length, respectively. Likewise, a constant stress model for self- or wind-loaded peduncles predicts that stem biomass and second moment of area should scale as the 5- and 8-power of stem length, respectively. In contrast, we observe that M
L1.36 and I
L2.75. These and other significant differences between observed and predicted allometric scaling exponents indicate that the morphometry of peduncles is far more efficient with respect to elevating reproductive organs than would be expected based on naive expectations about self-loading or wind-induced drag. Although mechanical theory can be used effectively to postulate form-function relationships, it must always be challenged empirically. Plants cannot obviate physical principles, but they have evolved adaptively in often unexpected ways.
| APPENDIX |
|---|
|
|
|---|
for these three models (see Table 1). The derivation for the
governing geometric self-similarity is mathematically trivial, because the model stipulates that radius R measured anywhere along length L is proportional to L, i.e., R
L1.
To determine the stress self-similarity
, we must find the numerical value of
for R = kL
such that stress
measured anywhere along L is constant. Noting that the maximum stress
max in any cross section is given by the formula
max = EKRmax, where E is the elastic modulus, K is curvature, and Rmax is the maximum distance from the centroid axis, the bending moment at any distant x from the tip of a leaning stem Mx is given by the formula Mx = EKIx. Because
max = MxRx/Ix, where Ix is the second moment of area measured at x, it follows that Mx = FxLx = (
/3)
VxLx =
Rx2Lx2 = k1Lx2
+2, where Fx is the self-generated bending force,
is bulk tissue density, Vx is stem volume measured at x, and k is a constant of proportionality. Because Ix = (
/4)Rx4 = k2Rx4, we see that
max = MxRx/Ix = (k1Lx2
+2kL
/k2Rx4) = k3Lx3
+2/kL4
= k5 Lx3
+24
. If
max is constant anywhere along Lx, then 3
+ 2 4
= 0. Therefore,
= 2, i.e., R = kL2. This derivation does not hold true for perfectly vertical stems, because
x is a compressive stress and thus proportional only to the mass of tissues distal to x. Assuming that R does not vary as a function of x, it follows that
x
L.
The scaling exponent for the elastic self-similarity model follows directly from Euler's formula for the critical buckling load Pcr of a vertical columnar support member with length L and uniform radius R, i.e., Pcr = kEI/L2. Noting that Pcr
V
R2L, we see that R2L
I/L2. Because I
R4, it follows that R
L3/2.
Volume and length relations
Peduncle geometry conforms to a hollow, double-tapered, truncated cone (see Fig. 1). However, for simplicity, peduncle geometry is here approximated as that of a hollow single tapered truncated cone to draw heuristic comparisons with solid cones and cylinders. The objective is to estimate the ability of hollow peduncles to extend in length beyond that of their solid counterparts with equivalent volume and bulk material (tissue) densities, i.e., equivalent mass. The volume of a solid cylinder Vsc with diameter D and length L equals
D2L/4. The volume of a solid truncated cone Vstc with length Lstc and basal and apical diameters Da and Db equals
Lstc(Db2 + DbDa + Da2)/12, or
Lstc(1 + k + k2) Db2/12, where k = Da /Db. Assuming both geometries are composed of the same materials with equivalent bulk densities and that both have equivalent basal diameters (i.e., D = Db), it follows that Lstc/Lsc = 3/(1 + k + k2) such that, when k = 0.2, Lstc = 2.419Lsc.
The volume of a hollow truncated cone with length Lhtc, outer and inner basal diameters Db and db, and outer and inner apical diameters Da and da is given by the formula Vhtc =
Lhtc [(1 + k + k2)Db2 (db2 + dbda + db2)]/12, where k = Da/Db. Assuming that the structure has a uniform wall thickness t, we see that db = Db 2t and da = Da 2t = kDb 2t. Therefore, Vhtc =
Lhtc {(1 + k + k2)Db2 [(Db 2t)2 + (Db 2t)(kDb 2t) + (kDb 2t)2]}/12 =
Lhtc(6Dbt + 6 kDbt 12t2)/12.
Provided that a hollow and a solid truncated cone have the same volume and material bulk density, Lhtc(6Dbt + 6kDbt 12 t2) = Lstc(1 + k + k2) Db2. Denoting t/Db as
(where
0.5 and k
2
), we see that Lhtc[6
(1 + k 2
)] = Lstc(1 + k + k2) such that Lhtc/Lstc = (1 + k + k2)/[6
(1 + k 2
)]. If k = 2
(no inner apical diameter, i.e., Da = 2 t), then Lhtc/Lstc = (1 + 2
+ 4
2)/6
. Therefore, when
= 0.01, Lhtc = 17Lstc.
| FOOTNOTES |
|---|
2 Author for reprint requests (kjn2{at}cornell.edu
) ![]()
| LITERATURE CITED |
|---|
|
|
|---|
McMahon T. A. 1973 Size and shape in biology. Science 179: 1201-1204
Monteith J. L. 1973 Principles of environmental physics. Elsevier, New York, New York, USA
Niklas K. J. 1989 Nodal septa and the rigidity of aerial shoots of Equisetum hyemale. American Journal of Botany 76: 521-531[CrossRef][ISI]
Niklas K. J. 1990 Safety factors in vertical stems: evidence from Equisetum hyemale. Evolution 43: 1625-1636[ISI]
Niklas K. J. 1992 Plant biomechanics. University of Chicago Press, Chicago, Illinois, USA
Niklas K. J. 1994 Plant allometry. University of Chicago Press, Chicago, Illinois, USA
Niklas K. J. H.-C. Spatz 2000 Wind-induced stresses in cherry trees: evidence against the hypothesis of constant stress levels. Trees 14: 230-237[CrossRef]
Okubo A. S. A. Levin 1989 A theoretical framework for data analysis of wind dispersal of seeds and pollen. Ecology 70: 329-338[CrossRef][ISI]
Sokal S. S. F. J. Rohlf 1981 Biometry. Freeman, San Francisco, California, USA
Spatz H.-C. H. Beismann F. Brüchert A. Emanns T. Speck 1997 Biomechanics of the giant reed Arundo donax. Philosophical Transactions of the Royal Society (London) B 352: 1-10[CrossRef]
Spatz H.-C. F. Brüchert 2000 Basic biomechanics of self-supporting plants: wind loads and gravitational loads on a Norway spruce tree. Forest Ecology and Management 135: 3-44[CrossRef][ISI]
Spatz H.-C. T. Speck 1994 Local buckling and other modes of failure in hollow plant stems. Biomimetics 2: 149-173
Speck T. H.-C. Spatz D. Vogellehner 1990a Contributions to the biomechanics of plants. I. Stabilities of plant stems with strengthening elements of different cross-sections against weight and wind forces. Botanica Acta 103: 111-122[ISI]
Speck T. H.-C. Spatz D. Vogellehner 1990b Contributions to the biomechanics of plants. II. Stability against local buckling in hollow plant stems. Botanica Acta 103: 123-130[ISI]
Vogel S. 1981 Life in moving fluids. Willard Grant, Boston, Massachusetts, USA
Wainwright S. A. W. B. Biggs J. D. Currey J. M. Gosline 1976 Mechanical design in organisms. John Wiley and Sons, New York, New York, USA
This article has been cited by other articles:
![]() |
A. M. Cooley, A. Reich, and P. Rundel Leaf support biomechanics of neotropical understory herbs Am. J. Botany, April 1, 2004; 91(4): 573 - 581. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |