Effect of tree size and competition on tension wood production over time in beech plantations and assessing relative gravitropic response with a biomechanical model

? Premise of the study: Gravitropic movements are unexpected mechanical processes that could disturb tree design allometries derived from the physics of nonliving bodies. We investigated whether the scaling law of gravitropic performance (power of -2 of stem diameter) derived from integrative biomechanical modeling is disturbed by ontogeny or environment, then discuss the silvicultural and dendroecological consequences. ? Methods: In a beech (Fagus sylvatica) plantation, four plots with different initial planting densities evolved without any intervention for 26 yr. Regular tree inventories and a silvicultural model were used to monitor competition over time in each plot. The radial production of tension wood was quantified using a cross-section of the stems at 1.30-m height, and an integrative biomechanical model computed the tree gravitropic performance over time. ? Key results: All trees developed tension wood over the whole period, with higher amounts at the youngest age, resulting in theoretical lean corrections of ca. 20-30° on the first 4 m of the stem over the whole period. The scaling law of gravitropic performance is slightly larger than the power of -2 of stem diameter. ? Conclusions: Gravitropic performance in forest ecosystems is mainly limited by size (diameter). Ontogenic acclimation of tension wood formation allows the youngest trees to be more reactive. No additional effect of spacing was found. However, silviculture influences size and, therefore, tree reactivity at a given age. Such results will be helpful for dendroecological approaches that use wood as a marker of environmental disturbances or a trait linked to plant strategies.     

Biomechanical modeling of gravitropic response of branches: roles of asymmetric periphery growth strain versus self-weight bending effect

Bending movement of a branch depends on the mutual interaction of gravitational disturbance, phototropic response, and gravitropic correction. Four factors are involved in gravitropic correction: asymmetric growth strain, eccentric growth increment, heterogeneous longitudinal elasticity (MOE), and initial radius which are associated with reaction wood production. In this context, we have developed a simplified model to calculate the rate of curvature change by combination of these factors. Experimental data from Taiwan red cypress were used to test the validity of the model. Our results show clearly that asymmetric growth strain is the main factor involved in correction. Eccentric growth increment has positive efficiency and increases correction in addition to growth strain, while the total effect of longitudinal MOE variation has negative efficiency and decreases correction. Spring-back strain measurement is found to be useful for the measurement of self-weight bending moment of branch. The branches studied are essentially close to a biomechanical equilibrium which maintains branches in horizontal positions. In the case of a deciduous dicotyledonous tree, flamegold, the effect of defoliating behavior on measured growth strain and curvature change was formulated by a modified model. Growth strain is the sum of measured growth strain after defoliation and spring-back strain during defoliation. The curvature change can be calculated by using measured growth strain and spring-back strain after defoliation. These results show that full-leaf branches of flamegold have a tendency to bend downward, but defoliated branches have a tendency for upward bending. The efficiency of correction increased after defoliation due to weight loss.     

Identification of biomechanical factors involved in stem shape variability between apricot tree varieties

BACKGROUND AND AIMS: Stem shape in angiosperms depends on several growth traits such as elongation direction, amount and position of axillary loads, stem dimensions, wood elasticity, radial growth dynamics and active re-orientation due to tension wood. This paper analyses the relationship between these biomechanical factors and stem shape variability. METHODS: Three apricot tree varieties with contrasting stem shape were studied. Growth and bending dynamics, mechanical properties and amount of tension wood were measured on 40 1-year-old stems of each variety during one growth season. Formulae derived from simple biomechanical models are proposed to quantify the relationship between biomechanical factors and re-orientation of the stems. The effect of biomechanical factors is quantified combining their mechanical sensitivity and their actual variability. RESULTS: Re-orientations happened in three main periods, involving distinct biomechanical phenomena: (a) passive bending due to the increase of shoot and fruit load at the start of the season; (b) passive uprighting at the fall of fruits; (c) active uprighting due tension wood production at the end of the season. Differences between varieties mainly happened during periods (a) and (b). CONCLUSIONS: The main factors causing differences between varieties are the length/diameter and the load/cross-sectional area ratios during period (a). Wood elasticity does not play an important role because of its low inter-variety variability. Differences during period (b) are related to the dynamics of radial growth: varieties with early radial growth bend weakly upward because the new wood layers tend to set them in a bent position. The action of tension wood during period (c) is low when compared with passive phenomena involved in periods (a) and (b).     

Effects of bending stress on taper and growth of stems of youngEucalyptus regnans trees

A glasshouse experiment investigated the effect of bending stress on stem radial and height growth and stem taper ofEucalyptus regnans seedlings. Eighteen-week-old, potted seedlings were bent continuously for 8 weeks with a static bending stress. The bending treatment was then removed and the seedlings grown for another 12 weeks. Other seedlings were stayed vertically throughout the experiment whilst control seedlings were neither bent nor stayed. Seedlings were rotated every 2 days to prevent reaction wood developing asymmetrically in the stems of bent trees. Bent trees had higher radial growth rates, developed more tapered stems and had higher safety factors (the ratio of stem radius to the minimum radius required to prevent the tree toppling over) than unbent seedlings. They produced a band of tension wood in their stems and ceased height growth whilst bent. When bending ceased, they resumed normal radial and height growth. Unbent trees developed more cylindrical stems. There were no differences in growth behaviour between stayed and control trees. Bent and unbent trees all developed a butt swell, the taper of which was not affected by treatment. It was concluded that bending stress has substantial effects on both the size and taper of tree stems. However, the development of butt swell is independent of the bending stress applied. The results were considered in relation to biomechanical theories of tree stem development.     

Stem biomechanics of three columnar cacti from the Sonoran Desert

The allometric relationship of stem length L with respect to mean stem diameter D was determined for 80 shoots of each of three columnar cactus species (Stenocereus thurberi, Lophocereus schottii, and S. gummosus) to determine whether this relationship accords with that predicted by each of three contending models purporting to describe the mechanical architecture of vertical shoots (i.e., geometric, stress, and elastic similitude, which predict L proportional to D(alpha), with alpha = 1/1, 1/2, and 2/3, respectively). In addition, anatomical, physical, and biomechanical stem properties were measured to determine how the stems of these three species maintain their elastic stability as they increase in size. Reduced major axis regression of L with respect to D showed that alpha = 2.82 ± 0.14 for S. thurberi, 2.32 ± 0.19 for L. schottii, and 4.21 ± 0.31 for S. gummosus. Thus, the scaling exponents for the allometry of L differed significantly from that predicted by each of the three biomechanical models. In contrast, these exponents were similar to that for the allometry previously reported for saguaro. Analyses of biomechanical data derived from bending tests performed on 30 stems selected from each of the three species indicated that the bulk stem tissue stiffness was roughly proportional to L2, while stem flexural rigidity (i.e., the ability to resist a bending force) scaled roughly as L3. Stem length was significantly and positively correlated with the volume fraction of wood, while regression analysis of the pooled data from the three species (i.e., 90 stems) indicated that bulk tissue stiffness scaled roughly as the 5/3-power of the volume fraction of wood in stems. These data were interpreted to indicate that wood served as the major stiffening agent in stems and that this tissue accumulates at a sufficient rate to afford unusually high scaling exponents tot stem length with respect to stem diameter (i.e., disproportionately large increments of stem length with respect to increments in stem diameter). Nevertheless, the safety factor against the elastic failure of stems (computed on the basis of the critical buckling height divided by actual stem length) decreased with increasing stem size tot each species, even though each species maintained an average safety factor equal to two. We speculate that the apparent upper limit to plant height calculated for each species may serve as a biomechanical mechanism for vegetative propagation and the establishment of dense plant colonies by means of extreme stem flexure and ultimate breakage, especially for S. gummosus.     

Wood mechanics, allometry, and life-history variation in a tropical rain forest tree community

Wood density plays a central role in the life-history variation of trees, and has important consequences for mechanical properties of wood, stem and branches, and tree architecture. Wood density, modulus of rupture, modulus of elasticity, and safety factors for buckling and bending were determined for saplings of 30 Bolivian rain forest tree species, and related to two important life-history axes: juvenile light demand and maximum adult stature. Wood density was strongly positively related to wood strength and stiffness. Species safety factor for buckling was positively related to wood density and stiffness, but tree architecture (height : diameter ratio) was the strongest determinant of mechanical safety. Shade-tolerant species had dense and tough wood to enhance survival in the understorey, whereas pioneer species had low-density wood and low safety margins to enhance growth in gaps. Pioneer and shade-tolerant species showed opposite relationships between species traits and adult stature. Light demand and adult stature affect wood properties, tree architecture and plant performance in different ways, contributing to the coexistence of rain forest species.     

Study of tension wood in the artificially inclined seedlings of <Emphasis Type="Italic">Koelreuteria henryi</Emphasis> Dummer and its biomechanical function of negative gravitropism

Key message

Stem reorientation is critical to tree survival. With anatomical observation and strain measurement, the tension wood formation and biomechanical behavior were studied to gain insights into tree uprighting process.

Abstract

Tension wood plays a role in maintaining the mechanical stability of angiosperm trees. Both biological and physical aspects of tension wood are essential in understanding the mechanism of trunk or branch reorientation. In this study, we worked on both tension wood formation and its biomechanical function in artificially inclined 2-year-old Koelreuteria henryi seedlings. The tension wood formation and reorientation process of the trunk last for about 3 months. With pinning method, we confirmed that at the beginning of inclination the cambial zone including the vascular cambium and the developing normal wood fibers on the upper side of the inclined trunk perceives the onset of mechanical change and starts to produce G-fibers that generate a strong contractile released growth strain (RGS) for gravitropic correction. Stronger contractile RGS and more tension wood were found at the trunk base than at the half-height, suggesting that the trunk base plays a key role in trunk uprighting of K. henryi seedlings. The eccentric cambial growth in the tension wood side increases the efficiency of gravitropic correction and the compressive strains measured in the opposite wood of some inclined seedlings also help the upright movement.

     

TWIG: a model to simulate the gravitropic response of a tree axis in the frame of elasticity and viscoelasticity, at intra-annual time scale

Trees are able to maintain or modify the orientation of their axes (trunks or branches) by tropic movements. For axes in which elongation is achieved but cambial growth active, the tropic movements are due to the production of a particular wood, called reaction wood which is prestressed within the growing tree. Several models have been developed to simulate the gravitropic response of axes in trees due to the formation of reaction wood, all within the frame of linear elasticity and considering the wood maturation as instantaneous. The effect viscoelasticity of wood has, to our knowledge, never been considered. The TWIG model presented in this paper aims at simulating the gravitropic movement of a tree axis at the intra-annual scale. In this work we studied both the effect of a non-instantaneous maturation process and of viscoelasticity. For this purpose, we considered the elastic case with maturation considered as an instantaneous process as the reference. The introduction of viscoelasticity in TWIG has been done by coupling TWIG to a model developed for bridges. Indeed from a purely mechanical point of view, bridges and trees are very similar: they are structures which are built in stages, they are made of several materials (composite structures), their materials are prestressed (wood is prestressed during the maturation process as a result of polymerisation of lignin and cellulose to form the secondary cell wall and concrete is prestressed during drying). Simulations gave evidence that the reorientation process of axes can be significantly influenced by the kinetics of maturation. Moreover the model has now to be tested with more experimental data of wood viscoelasticity but it appears that in the range of a relaxation time from 0 to 50 days, viscoelasticity has an important effect on the evolution of tree shape as well as on the values of prestresses.     

Growth stress controls negative gravitropism in woody plant stems

In the shoots of woody plant species, reaction-wood fibers are formed on the upper or lower side of the secondary xylem of a leaning trunk or branch wherever large, internal growth stress is generated. Negative gravitropic movement in woody plant stems is proposed to be the result of growth stress generated in the reaction-wood tissue. This study examines the interaction between bending moment due to increasing self-weight and recovery moment resulting from asymmetric growth stress, and tests a hypothesis that describes the relationship based on the structural mechanics "beam theory". Simulations of observed tree branch morphology of Magnolia kobus DC., Juniperus chinensis L., Abies saccharinensis Fr. Schum., and Prunus spachiana Kitamura f. spachiana cv. Plenarosea showed that (i) the growth stress generated in the reaction wood is sufficient to counteract the gravitropic response to increasing self-weight, and (ii) the specific directional angle of the shoot apex or preferred angle of the elongation zone plays an important role in controlling the spatial shape of the branch stem that is peculiar to plant species with large growth stress generated in the reaction-wood tissue.     

A model to simulate the gravitropic response and internal stresses in trees,considering the progressive maturation of wood

Key message

The developed model of gravitropism takes non-instantaneous maturation of wood into account which enabled to correctly simulate different gravitropic phases and realistic internal stress profiles.

Abstract

A new biomechanical model of tree movement in relation to gravity (gravitropism) is proposed in this study. The modelling of the progressive maturation of wood cells is taken into account, as well as spatio-temporal variations in maturation strains (MS) and mechanical properties. MS were identified using an inverse method that allows the model to fit the gravitropic reaction observed experimentally. For this purpose, the curvature during righting movement, the geometry and the mass distribution of a two-year-old poplar tree was measured. The identified MS are higher than expected, which shows the underestimation of MS by usual measurements. By using the same mechanical parameters and MS as an input, the model gives satisfying results in terms of shape modelling for different trees up to 32 days after tree tilting. The model is able to simulate the latency phase observed in the tree righting movement, and the internal stress profile in the trunk is realistic (low compressive value in the central part of the trunk and zero stress in newly formed cells). The next development of the model will aim to simulate the end of the gravitropic phase in relation with the regulation of MS by the tree.     

Safety and streamlining of woody shoots in wind: an empirical study across 39 species in tropical Australia

? Wind is a key mechanical stress for woody plants, so how do shoot traits affect performance in wind? ? We used a vehicle mounted apparatus to measure drag, streamlining and mechanical safety in 127 vertical lead-shoots, 1.2 m long, across 39 species in tropical Australia. ? Shoot dimensions and stem tissue properties were closely coupled so that shoots with low stem specific gravity or larger projected area had thicker stems. Thicker stems provide larger second moment of area (I), which increased shoot safety and bending stiffness but impeded shoot reconfiguration in strong winds, including frontal area reduction. Nonetheless, increasing I also improved streamlining. Streamlining was unrelated to traits except I. Stem tissue material properties only had small effects. Higher modulus of rupture increased shoot safety and higher Young's modulus impeded shoot reconfiguration. ? We found no conflict between bending stiffness and streamlining for woody shoots. Stiffness might help streamlining by increasing damping and stability, thereby reducing flagging in wind. Tissue-level traits did influence shoot-level mechanical safety and behaviour, but shoot geometry was much more important. Variable shoot and stem traits, which all influenced shoot biomechanics, were integrated in shoots to yield a relatively narrow range of outcomes in wind.     

Numerical modelling of shape regulation and growth stresses in trees

The main objective of this paper is to present the results of a study of the interactions between the growth and design of a tree with regards to biomechanical factors at the plant level. A numerical incremental model dedicated to the calculation of tree mechanical behaviour has been integrated in the plant architecture simulation software AMAPpara. At any stage of tree growth, a new equilibrium was calculated considering the weight increment applied on the structure, i.e. the mass of new wood layers and vegetative elements, as well as the biomechanical reaction caused by cell maturation strains in both normal and reaction wood. The resulting incremental displacements allowed the tree shape to be modified. The field of growth stresses was calculated within the stem, using a cumulative process taking into consideration the past history of each growth ring. The simulation results of trunk and branch shape, as well as internal stresses, were examined after consideration of different growth strategies. A block of trees was also simulated in order to show the influence of spatial competition on stem curvature and the variability in growth stress.     

The gravitropic response of poplar trunks: key roles of prestressed wood regulation and the relative kinetics of cambial growth versus wood maturation

In tree trunks, the motor of gravitropism involves radial growth and differentiation of reaction wood (Archer, 1986). The first aim of this study was to quantify the kinematics of gravitropic response in young poplar (Populus nigra x Populus deltoides, 'I4551') by measuring the kinematics of curvature fields along trunks. Three phases were identified, including latency, upward curving, and an anticipative autotropic decurving, which has been overlooked in research on trees. The biological and mechanical bases of these processes were investigated by assessing the biomechanical model of Fournier et al. (1994). Its application at two different time spans of integration made it possible to test hypotheses on maturation, separating the effects of radial growth and cross section size from those of wood prestressing. A significant correlation between trunk curvature and Fournier's model integrated over the growing season was found, but only explained 32% of the total variance. Moreover, over a week's time period, the model failed due to a clear out phasing of the kinetics of radial growth and curvature that the model does not take into account. This demonstrates a key role of the relative kinetics of radial growth and the maturation process during gravitropism. Moreover, the degree of maturation strains appears to differ in the tension woods produced during the upward curving and decurving phases. Cell wall maturation seems to be regulated to achieve control over the degree of prestressing of tension wood, providing effective control of trunk shape.     

The mechanical role of bark

The ability of stem bark to resist bending forces was examined by testing in bending segments of Acer saccharum, Fraxinus americana, and Quercus robur branches with and without their bark. For each species, the bark contributed significantly to the ability of stem segments differing in age to resist bending forces, but its contribution was age-dependent and differed among the three species. The importance of the mechanical role of the bark decreased basipetally with increasing age of F. americana and Q. robur stem segments and was superceded by that of the wood for segments ≥ 6 yr old. A. saccharum bark was as mechanically important as the wood for stem segments 7 yr old but was not a significant stiffening agent for younger or older portions of stems. On average, the stiffness of the bark from all three species was 50% that of the wood. However, the geometric contribution to the flexural rigidity of stems made by the bark (i.e., the bark's second moment of area) was sufficiently large to offset its lower stiffness (Young's modulus) relative to that of the wood. A simple model is presented that shows that the bark must be as mechanically important as the wood when its radial thickness equals 32% that of the wood and its stiffness is 50% that of the wood. Based on this model, which is shown to comply with the data from three species purported to have stiff woods, it is evident that the role of the bark cannot be neglected when considering the mechanical behavior of juvenile woody stems subjected to externally applied bending forces.     

Biomechanics of the columnar cactus Pachycereus pringlei

We report the longitudinal variations in stiffness and bulk density of tissue samples drawn from along the length of two Pachycereus pringlei plants measuring 3.69 and 5.9 m in height to determine how different tissues contribute to the mechanical stability of these massive vertical organs. Each of the two stems was cut into segments of uniform length and subsequently dissected to obtain and mechanically test portions of xylem strands, stem ribs, and a limited number of pith and cortex samples. In each case, morphometric measurements were taken to determine the geometric contribution each tissue likely made to the ability of whole stems to resist bending forces. The stiffness of each xylem strand increased basipetally toward the base of each plant where stiffness sharply decreased, reaching a magnitude comparable to that of strands 1 m beneath the stem apex. The xylem was anisotropic in behavior, i.e., its stiffness measured in the radial and in the tangential directions differed significantly. Despite the abrupt decrease in xylem strand stiffness at the stem base, the contribution made by this tissue to resist bending forces increased exponentially from the tip to the base of each plant due to the accumulation of wood. A basipetal increase in the stiffness of the pith (and, to limited extent, that of the cortex) was also observed. In contrast, the stiffness of stem rib tissues varied little as a function of stem length. These tissues were stiffer than the xylem in the corresponding portions of the stem along the upper two-fifths of the length of either plant. Tissue stiffness and bulk density were not significantly correlated within or across tissue types. However, a weak inverse relationship was observed for these properties in the case of the xylem and stem rib tissues. We present a simple formula that predicts when stem ribs rather than the xylem strands serve as the principal stiffening agents in stems. This formula successfully predicted the observed aspect ratio of the stem ribs (the average quotient of the radial and tangential dimensions of rib transections), and thus provided circumstantial evidence that the ribs are important for mechanical stability for the distal and younger regions of the stems examined.     

Development and Growth Form of the Neotropical Liana <Emphasis Type="Italic">Croton nuntians</Emphasis>: The Effect of Light and Mode of Attachment on the Biomechanics of the Stem

The neotropical liana Croton nuntians (Euphorbiaceae) can occur in a variety of different growth habits. Juvenile freestanding plants are mechanically stable without support and resemble morphologically young trees or shrubs, whereas adult plants are climbers. Ontogenetic variation of bending and torsion properties of different growth phases are analyzed by measurements of flexural stiffness, structural bending modulus, torsional stiffness and structural torsional modulus. Mechanical and anatomical data show two fundamentally different patterns for juvenile freestanding and adult climbing plants. In freestanding plants, mechanical properties and the contribution of cortex, wood, and pith to the stem cross-section vary only little during ontogeny as is typical for semi-self-supporting plants. In contrast, climbing plants become significantly more flexible during ontogeny, as is characteristic for lianas. This is accompanied by a transition to the formation of a less dense wood type with large diameter vessels and an increasing contribution of flexible tissues (less dense wood and cortex) to the cross-sectional area and the axial second moment of area of the stems. Depending on the environmental conditions, freestanding plants can differ considerably in their appearance due to differences in branching system or stem taper. Therefore the influence of light quantity, measured as percentage of canopy opening, on the mechanical properties and the stem anatomy was tested. Freestanding plants grown with strong shade are significantly more stiff in bending compared with plants grown with a moderate light environment.     

Effect of circumferential heterogeneity of wood maturation strain, modulus of elasticity and radial growth on the regulation of stem orientation in trees

Active mechanisms of re-orientation are necessary to maintain the verticality of tree stems. They are achieved through the production of reaction wood, associated with circumferential variations of three factors related to cambial activity: maturation strain, longitudinal modulus of elasticity (MOE) and eccentric growth. These factors were measured on 17 mature trees from different botanical families and geographical locations. Various patterns of circumferential variation of these factors were identified. A biomechanical analysis based on beam theory was performed to quantify the individual impact of each factor. The main factor of re-orientation is the circumferential variation of maturation strains. However, this factor alone explains only 57% of the re-orientations. Other factors also have an effect through their interaction with maturation strains. Eccentric growth is generally associated with heterogeneity of maturation strains, and has an important complementary role, by increasing the width of wood with high maturation strain. Without this factor, the efficiency of re-orientations would be reduced by 31% for angiosperms and 26% for gymnosperms. In the case of angiosperms, MOE is often larger in tension wood than in normal wood. Without these variations, the efficiency of re-orientations would be reduced by 13%. In the case of gymnosperm trees, MOE of compression wood is lower than that of normal wood, so that re-orientation efficiency would be increased by 24% without this factor of variations.     

Patterns of auxin distribution during gravitational induction of reaction wood in poplar and pine

Gravistimulation of tree stems affects wood development by unilaterally inducing wood with modified properties, called reaction wood. Commonly, it also stimulates cambial growth on the reaction wood side. Numerous experiments involving applications of indole-3-acetic acid (IAA) or IAA-transport inhibitors have suggested that reaction wood is induced by a redistribution of IAA around the stem. However, in planta proof for this model is lacking. Therefore, we have mapped endogenous IAA distribution across the cambial region tissues in both aspen (Populus tremula, denoted poplar) and Scots pine (Pinus sylvestris) trees forming reaction wood, using tangential cryosectioning combined with sensitive gas chromatography-mass spectrometry analysis. Moreover, we have documented the kinetics of IAA during reaction wood induction in these species. Our analysis of endogenous IAA demonstrates that reaction wood is formed without any obvious alterations in IAA balance. This is in contrast to gravitropic responses in roots and shoots where a redistribution of IAA has been documented. It is also of interest that cambial growth on the tension wood side was stimulated without an increase in IAA. Taken together, our results suggest a role for signals other than IAA in the reaction wood response, or that the gravitational stimulus interacts with the IAA signal transduction pathway.     

Mathematical modelling of morphogenesis in fungi: a key role for curvature compensation ('autotropism') in the local curvature distribution model

The assumption that the mushroom stem has the ability to undergo autonomic straightening enables a mathematical model to be written that accurately mimics the gravitropic reaction of the stems of Coprinus cinereus . The straightening mechanism is called curvature compensation here, but is equivalent to the 'autotropism' that often accompanies the gravitropic reactions of axial organs in plants. In the consequently revised local curvature distribution model, local bending rate is determined by the difference between the 'bending signal' (generated by gravitropic signal perception systems) and the 'straightening signal' (proportional to the local curvature at the given point). The model describes gravitropic stem bending in the standard assay with great accuracy but has the virtue of operating well outside the experimental data set used in its derivation. It is shown, for example, that the mathematical model can be fitted to the gravitropic reactions of stems treated with metabolic inhibitors by a change of parameters that parallel the independently derived physiological interpretation of inhibitor action. The revised local curvature distribution model promises to be a predictive tool in the further analysis of gravitropism in mushrooms.     

The role of the epidermis as a stiffening agent in Tulipa (Liliaceae) stems

We investigated the hypothesis that the epidermis is a tension-stressed "skin' whose contribution to stem stiffness depends on the turgor pressure exerted on it by an hydrostatically inflated inner "core' of tissues. This hypothesis was tested by relying on the intensities of bending stresses due to stem flexure, which must reach their maximum levels at the outer surface of epidermis such that damage to the surface of the stem should produce the most significant decrease in overall flexural stiffness. We discerned whether the principal tension supporting members at the stem surface (cellulosic microfibrils) were oriented parallel or normal to stem length by comparing the bending stiffness of stems before and after their surface cells first received three parallel longitudinal incisions followed by one helical incision, and by comparing the bending stiffness of stems for which the sequence of cuts was reversed. The same protocol was also applied to stems with various water potentials to determine the effect of hydrostatic pressure on stem stiffness contributed by the surface. Based on the behavior of 82 turgid Tulipa stems, parallel cuts reduced, on average, stem stiffness by 8%, whereas a subsequent helical incision further reduced stiffness by 42%. In contrast, an initial helical incision reduced stem stiffness by 50%, while three subsequent parallel cuts through the same stems did not significantly further reduce stiffness. These results suggested that the net orientation of cellulose microfibrils in the outer epidermal walls was parallel to stem length. This was confirmed by microscopic observations of cells with dichroic staining and polarized light. The responses to surgical damage were directly proportional to stem water potential. We thus conclude that the epidermis, probably in conjunction with a single layer of subepidermal collenchyma cells, acts as a tension-stiffening agent that can contribute as much as 50% to overall stem stiffness We present a simple mechanical model that can account for all our observations.