ABSTRACT
We have focused on the giant kelp Nereocystis luetkeana to examine the mechanical scaling of benthic marine organisms loaded in tension by hydrodynamic forces. If we consider simply the allometry of the kelp’s morphological characters, we conclude that their stipes are underscaled relative to the blade area they support (i.e. that the kelp do not maintain stress or elastic similarity as they grow). However, a closer look at the characteristics of these kelp in the field reveals (1) that they have different blade shapes (and hence drag coefficients) and stipe material properties in different hydrodynamic environments, and (2) that they show a decrease in drag coefficient as they become larger. One consequence of these adjustments of blade and stipe morphology is that the maximum stresses in N. luetkeana stipes, when the kelp are pulled by peak tidal currents in their respective habitats, are similar for kelp of different sizes and for kelp from different sites. Hence, sessile organisms such as these kelp can, via their growth responses in different mechanical environments, show a phenomenon analogous to dynamic strain similarity. In addition, N. luetkeana also maintain a constant environmental stress factor, the ratio of the stress required to break a component of an organism (in this case the stipe) at some stage in its life to the maximum stress normally encountered in the habitat by that component during that stage (in this case, stress due to the drag on a kelp exposed to the peak tidal currents typical of the site at which it lives), both between habitats and as they grow.
INTRODUCTION
The consequences of body size for how organisms function has long intrigued biologists (reviewed in Alexander, 1971; Pedley, 1977; Banse and Mosher, 1980; McMahon and Bonner, 1983; Peters, 1983; Calder, 1984; Schmidt-Nielsen, 1984; Denny et al. 1985; Reiss, 1989; Patterson, 1992). Changes in body size are often accompanied by changes in shape, both during the ontogeny of an individual (e.g. Calder, 1984) and during the evolution of a lineage (e.g. Gould, 1966; Bonner, 1968; Stanley, 1973). Most of the attention on how such shape changes affect the mechanical functioning of organisms has focused on the allometry of terrestrial creatures, whose support structures are loaded as buckling columns or bending cantilevers by their weight (proportional to volume) (e.g. McMahon, 1973, 1975). In contrast, little is known about the scaling of benthic aquatic organisms, for whom hydrodynamic forces are more important than gravity. The few studies of scaling that have examined benthic organisms have focused on structures loaded in bending or shear (e.g. Koehl, 1977; Denny et al. 1985); however, many benthic organisms are loaded in tension. Although some general rules have been proposed about how tensile structures such as kelp stipes (stems) should scale (Peterson et al. 1982), field measurements of how forces on these structures vary with size have not been made. Furthermore, although some investigators have recognized that the gaits used by locomoting animals can affect their mechanical scaling (e.g. Rubin and Lanyon, 1984), researchers have not yet explored the ways in which the responses of sessile organisms to environmentally induced loads might affect their scaling.
The overall goal of this study was to analyze the mechanical scaling of aquatic tensile biological structures. We focused our analysis on the stipes (stems) of the giant bull kelp Nereocystis luetkeana (Mertens) Postels and Ruprecht (see Fig. 1A).
Objectives
Scaling of tensile structures
A variety of biological structures (such as spider silk, fruit stems or the stipes of kelp with gas-filled floats) are loaded in simple tension. One objective of the present study was to measure how the load (F) on the stipe of a tensile aquatic macrophyte scales with the size of the organism and, in turn, how the dimensions of the stipe vary with respect to F.
If two structures are scaled to show ‘stress similarity’, then the stress (force per unit cross-sectional area of material bearing that force) at some position in the large structure is the same as the stress at a comparable position in the smaller structure when both perform the same function, such as supporting their own weight (e.g. McMahon, 1975). Structures of different sizes loaded in tension show stress similarity when d ∝F0.5, where F is the force on the structure and d is its diameter (Peterson et al. 1982). Since stress in a tensile structure is simply force per unit cross-sectional area, diameter can be independent of length (Wainwright, 1970) as long as the weight of the tensile element is negligible relative to F.
If two structures are scaled to show ‘elastic similarity’ in bending, then the deflections of the large and small structures (normalized to their respective lengths) are the same when they perform comparable functions (McMahon, 1973). We define elastic similarity for a tensile structure as maintenance of ΔL/L, where ΔL is the linear displacement of the end of the tensile structure as it is stretched by load F, and L is the length of the structure before F is applied. For a structure loaded in tension, extension per unit length is a function of stress (e.g. Wainwright et al. 1976); therefore, a tensile structure that maintains stress similarity also maintains elastic similarity.
While the force on a tensile fruit stem is proportional to the volume of the fruit, Peterson et al. (1982) suggested that the force on a tensile kelp stipe is due to hydrodynamic drag on the blades (they assumed drag to be proportional to blade area, AB). Therefore, they predicted that d∝AB0.5. Although Peterson et al. (1982) did find that stipe diameter was independent of stipe length for the elk kelp Pelagophycus porra, they also found that stipe diameter appeared to be underscaled relative to blade area (d∝AB0.268) for stress (and therefore, elastic) similarity.
We tested the assumption that the forces (F) on macrophytes are simply proportional to blade area (AB), and we measured how the diameters (d) of the stipes scale with F. An organism in a unidirectional current is exposed to drag force (proportional to projected area), whereas an organism in ocean waves is also subjected to acceleration reaction force (proportional to volume) (Koehl, 1977; Denny et al. 1985). These hydrodynamic forces depend not only on the area or volume of the organism but also on the velocities and accelerations of the water flow in its habitat. Therefore, we measured the forces on kelp at similar velocities to those measured at the sites from which the kelp were collected.
Size-dependent material properties
A second objective of the present study was to incorporate measurements of size-dependent material properties into an analysis of mechanical scaling. The scaling models mentioned above have been based on the assumption that tissue material properties are independent of the size of an organism. This assumption, however, is invalid for many terrestrial plants (Niklas, 1992) and vertebrates (Currey, 1984) and some species of seaweed (Delf, 1932; Dudgeon and Johnson, 1992), all of which show changes in the mechanical properties of their support tissues as they grow and age.
Dynamic strain similarity
A third objective of the present study was to explore how the responses of kelp to different water flow environments might affect how the load (F) on the stipe scales with thallus (body) size. Rubin and Lanyon (1984) refined classical scaling ideas by considering the dynamic load imposed upon a structure during use. They found that during vigorous locomotion, the strains (for small deformations, strain is the change in length divided by the undeformed length) in the bones of vertebrates representing a range of body sizes were similar because large animals use gentler gaits than do small ones. They called this maintenance of peak functional strain independent of body size ‘dynamic strain similarity’. The concept of dynamic strain similarity should also apply to non-locomoting creatures subjected to environmental forces (e.g. kelp in flowing water) if the organisms respond to the environment in ways that affect the magnitudes of such forces. For example, passive reconfiguration of sessile organisms in moving fluids can reduce drag (e.g. Koehl, 1977; Vogel, 1984). Furthermore, some sessile organisms can change the drag they experience by changing their shapes (via muscle contraction or growth) in response to their water-flow habitats (e.g. Koehl, 1977).
Environmental stress factor
A fourth objective of this study was to investigate whether aquatic macrophytes maintain their environmental stress factor as they grow in different water-flow habitats.
By ‘environmental stress factor’ we mean the ratio of the strength (σbrk, the stress required to break the tissue) of a component of an organism at some stage in its life (in this case, the stipe of a non-senescent kelp sporophyte) to the typical maximum stress (σmax) experienced by that component as it is used day-to-day in the habitat of the organism during that life stage (in this case, the stress in the stipe due to the maximum hydrodynamic force on the kelp during times of peak tidal current). In contrast, the safety factor of a component of an organism is the ratio of the strength of the tissue from which the component is made to the peak stress that it experiences during its lifetime (Alexander, 1981). For example, the environmental stress factor for a leg bone would be the ratio of the strength of the bone material in a healthy adult to the stress in the bone during running at peak velocity or jumping to maximum height. In contrast, the safety factor for that bone might be the ratio of the strength of the portion of the bone tissue most weakened by osteoporosis in an old animal to the stress experienced when that animal fell off a cliff.
There is variability within a population in the strength (σbrk) of the support tissues of different individuals, and there is also variability in the peak stress (σpeak) encountered by different individuals in their lifetimes. Therefore, the safety factor is often considered statistically, both for organisms and for man-made structures (see Alexander, 1981); such a statistical safety factor is the mean of the σbrk values of the individuals in a population divided by the mean of the σpeak values they encounter in their lifetimes. In practice, the safety factor is very difficult to determine because long-term field studies of the populations are required to determine how σbrk values change with age. Similarly, measuring the lifetime σpeak values experienced by members of a population can be a formidable undertaking, although Denny and Gaines (1990) and Denny (1991) have described how the statistics of extremes can be used to estimate peak forces on wave-swept organisms. In contrast, environmental stress factor is easier to determine because it involves measurements of the strength of an organism’s tissues at a particular stage in its life and of the peak stresses experienced by those tissues as used by the organism at that stage in its habitat. Note that environmental stress factor can change with season and with the age of the organism, whereas safety factor represents the entire lifetime of the organism.
Several lines of evidence suggest that sessile organisms might be able to maintain a constant environmental stress factor or safety factor as they grow in different habitats. Some sessile marine organisms have been shown to change their breaking force (Etter, 1988; Palumbi, 1984) in response to their hydrodynamic environment, and some terrestrial plants have been found to maintain a constant safety factor as they increase in size (Tateno and Bae, 1990).
We have focused on the giant kelp Nereocystis luetkeana to determine how the thallus allometry, stipe material properties and blade responses to hydrodynamic environment affect environmental stress factors as the kelp grow in different water-flow habitats.
Nereocystis luetkeana
Nereocystis luetkeana sporophytes form extensive beds along the Pacific coast of North America, from California to Alaska. N. leutkeana blades, which can be up to 4 m long, are held near the water surface by a gas-filled pneumatocyst atop a long, slim stipe attached to the substratum by a holdfast (Abbott and Hollenberg, 1976). These kelp, which occur in water about 3–17 m deep, are exposed to tidal currents and, in some habitats, to non-breaking waves (Koehl and Alberte, 1988).
N. luetkeana are essentially annual kelp: young sporophytes appear in the early spring, grow to the water surface by mid-summer, and wash away during winter storms (Rigg, 1917; Markham, 1969; Foreman, 1970; Nicholson, 1970; Vadas, 1972; Abbott and Hollenberg, 1976; Hawkes, 1981; Rosell and Srivastava, 1984; Wheeler et al. 1984; Miller and Estes, 1989), although a few individuals may persist through the winter (Rigg, 1917; Foreman, 1970; Rosell and Srivastava, 1984; Wheeler et al. 1984). N. leutkeana blades grow like conveyor belts: new tissue is added in the meristematic region near the blade base, sori (patches on the blades in which spores are produced) develop distally in older tissue, and sori are shed and old tissue is eroded at the blade tip. The main spore-producing season appears to be mid to late summer (Rigg, 1917): although spore production begins after the kelp reach the surface (Foreman, 1970; Duncan, 1973) and continues through autumn (Foreman, 1970; Abbott and Hollenberg, 1976; Amsler and Neushul, 1991), blade growth rate begins to slow in late summer (Wheeler et al. 1984), photosynthetic rate decreases (Rosell and Srivastava, 1984) and blades become badly frayed by winter (Foreman, 1970). Furthermore, epiphyte load on the blades increases from late summer into the winter (Peters, 1913; Rigg, 1917; Markham, 1969; Foreman, 1970; Nicholson, 1970; Hawkes, 1981; Rosell and Srivastava, 1984).
MATERIALS AND METHODS
Field sites
Nereocystis luetkeana were studied from three sites near San Juan Island, Washington, USA. The ‘protected’ site was a stretch of shore between Shady Cove and Cantilever Point near Friday Harbor Laboratories that experienced peak velocities of 0.5 m s−1 (maximum tidal current plus maximum velocity fluctuation measured during a storm, Koehl and Alberte, 1988). The ‘current-swept’ site was Turn Rock in San Juan Channel that experienced peak velocities of 1.5 m s−1 (the maximum tidal current at this site during a year was 1.1 m s−1 (National Oceanic and Atmospheric Administration, 1981); the maximum velocity fluctuation measured during a storm at this site was 0.4 m s−1 (Koehl and Alberte, 1988)). The ‘wavy’ site was Cattle Point, which was subjected to small waves superimposed on tidal currents; daily peak velocities ranged between 0.6 and 1.7 m s−1 on non-stormy days at a site about 100 m from ours (D. Pentcheff, unpublished data; recorded over nine separate days, April–July).
N. luetkeana were collected from the kelp bed at each of the sites described above. A boat was anchored at two haphazardly selected positions at Turn Rock (one near the edge and another near the middle of the bed), and all the kelp within reach of the boat at each position were collected. This technique was also used to sample kelp from Shady Cove, but because the kelp bed was narrow at this site, there was no distinction made between the edge and the middle of the bed. Kelp at Cattle Point were collected by wading from the shore to a haphazardly selected point in the middle of the kelp bed; all kelp encountered within 1 m of a transect run from that point parallel to the shore were taken until the collecting bags were full. Collections were made at Shady Cove on August 23 and September 29, at Turn Rock on July 19, and at Cattle Point on July 6, 1982.
Drag measurements
The drag (D) on individual N. luetkeana was determined (Koehl and Wainwright, 1985) to the nearest 0.1 N using an Ametek model LKG5 force gauge. Each kelp was towed just below the water surface at a number of different velocities (measured with a Marsh–McBirney model 511 electromagnetic flow meter) outside the wake of a boat. We focused on velocities at the high end of the range encountered in situ to examine the susceptibility of kelp to breakage by currents: N. luetkeana from the protected and the wavy sites were towed at 0.25, 0.50 and 0.75 m s−1, and kelp from the current-swept site at 0.50, 1.00 and 1.50 m s−1.
Morphological measurements
Each morphological variable quantified in this study (shown in Fig. 1 and described below) was measured at least three times to determine the precision of our techniques. Symbols used to characterize morphology throughout the manuscript are listed in Table 1.
Blade characteristics
All blades on a kelp were removed, spread out along with grids marked in cm2, and photographed from an overhanging balcony one storey (3.5 m) above them. The resulting slides were projected onto paper and the perimeter of each blade was traced. Vernier calipers were used to measure blade length (LB) and blade width (WB) at the widest part of the blade on these tracings (Fig. 1B); these measurements were made to the nearest millimeter after conversion. Since there were several tiny or broken blades on each kelp, the mean length and width of only the longest 50% of the blades on each kelp was used to characterize the LB and WB for each kelp. The projected areas of the blades were measured to the nearest 0.001 m2 by weighing paper cut-outs of the tracings. The total blade projected area (AB) for a kelp was calculated by summing the areas of all of the blades of the kelp. To obtain an estimate of the error in our first method (see Allometric analysis below), projected area measurements were repeated on seven of the kelp (ranging in size from 0.1 to 2.3 m2) using another technique: the slides of blades were projected onto a digitizing tablet (Jandel model 220/2210) interfaced with a Toptec XT computer, and their areas were calculated using Sigma-Scan software (Jandel Scientific).
Stipe and holdfast characteristics
Stipe diameter was measured to the nearest millimeter using vernier calipers at the holdfast (d1) (see Fig. 1D), and at 15 cm intervals (lA, measured with a ruler) along the entire length of the stipe; stipe diameter (dN) just proximal to the pneumatocyst was also measured. Stipe length (LS) was determined to the nearest centimeter by summing the values of lA for the entire stipe plus lN and the length of the pneumatocyst (see Fig. 1D,F). The narrowest stipe diameter (d) measured for each kelp was used in subsequent statistical comparisons, for calculation of stress due to drag (see below) and for scaling arguments.
The volume of the gas-filled cavity (VG) in each N. luetkeana from Turn Rock and Shady Cove was determined by cutting the kelp transversely at the position of dN (Fig. 1D), filling the pneumatocyst and the stipe cavities with water, and then measuring the volume of this water to the nearest milliliter using a graduated cylinder.
The areas of intact, undamaged holdfasts (H) were determined for kelp from Cattle Point and Shady Cove (no undamaged holdfasts were available for kelp from Turn Rock). The longest dimension (Ha) of each holdfast, as well as the width (Hb) of the holdfast perpendicular to Ha were measured to the nearest 0.5 cm with a ruler (Fig. 1E). Holdfasts were assumed to be elliptical and their areas were calculated using Ha and Hb as the major and minor axes.
Principal components analysis
To compare the morphologies of the kelp from different sites, we used principal components analysis (Chatfield and Collins, 1980) to summarize a number of the variables we measured to characterize stipe morphology (stipe length and stipe smallest diameter) and blade morphology (total projected blade area, blade length, blade width and drag, which depends on blade ruffliness as well as size; Koehl and Alberte, 1988). The principal components analysis was performed using SPSSX software (version 9; SPSS, Inc) on a VACS mainframe computer. Each principal component represents a linear summary (similar to a linear multiple regression) of the data that combines several variables into one descriptive variable (the principal component). Several principal components can be constructed from each set of variables. The first principal component is the linear equation for the line that best fits the data; that is, that accounts for most of the variation. The second principal component is the line constructed from the data, perpendicular to the first principal component, that accounts for most of the remaining residual variation.
Principal components analysis not only allows us to summarize many variables with only one or two new variables (principal components), but it also allows us to look for patterns in the clumping of the data when the principal components that are derived from the data are plotted against one another.
Mechanical properties of stipe tissues
Mechanical properties of stipe tissues (material properties; symbols listed in Table 2) from each kelp were measured by conducting tensile tests using techniques described by Koehl and Wainwright (1977, 1985). These tests were conducted on a Houndsfield Tensometer, model W, arranged to record force and extension electronically as described by LaBarbera (1985). Tests were carried out on at least three replicate pieces of stipe taken from the end nearest the holdfast, where kelp in situ usually break (Foreman, 1970; M. A. R. Koehl, personal observation). Material properties were only determined for specimens that did not slip or break at the tensometer grips. When more than one test was successful for a kelp, we used the mean value for that kelp. The pre-test length between the grips of each specimen (L0), which was measured to the nearest 0.1 mm with vernier calipers, was between 11 and 14 cm. Stipe specimens were cylindrical, and the diameter of the narrowest portion of each was measured to the nearest 0.1 mm using vernier calipers; cross-sectional areas (A0) were calculated using this diameter. Strain rates (strain rate=(ΔL/L0)/t, where ΔL is the increase in length of the specimen as it is pulled, and t is the time interval over which ΔL occurred) ranged from 0.025 to 0.084 s−1; strain rate within this range had no significant effect (i.e. P>0.05; linear regression analysis) on any of the material properties used in our analysis. [Note that a test for zero slope is unbiased in the presence of error in the independent variable (Fuller, 1987). Therefore, no correction of slope was necessary for this result.]
The material properties we calculated are illustrated in Fig. 2. The extension ratio (λ) of the specimen is given by:
We also conducted a series of stress–extension tests in which we pulled N. luetkeana stipes to various extensions and then returned them to their initial lengths (L0). The area under a plot of stress versus extension ratio for a specimen being pulled is a measure of the work (Wpull) per volume required to stretch the specimen to a given extension ratio, while the area under the curve as the specimen is being brought back to L0 is a measure of the strain energy (Wrecoil) stored in the specimen available for elastic recoil (Fig. 3B). A measure of the resilience of stipe tissue is Wrecoil/Wpull. Such resilience measurements were conducted for N. luetkeana collected from the protected site during the summer of 1975.
Site-specific mechanical features
We used estimates of typical peak velocities for the two sites subjected to unidirectional currents (0.5 m s−1 for the protected site and 1.5 m s−1 for the current-swept site, as described above in the section on field sites) to determine the ‘peak drag force’ on each kelp in situ. The peak drag on each kelp was divided by the narrowest cross-sectional area of the stipe of that kelp to yield an estimate of the peak stress due to drag in the stipe of each kelp (σdrag, Fig. 2). The extension ratio corresponding to that stress was determined from each stress/extension ratio curve for that kelp, and the mean of those extension ratios for each kelp was used as an estimate of the peak extension ratio due to drag experienced by that kelp in situ (λdrag, Fig. 2). For each curve, we noted whether the point on the curve defined by (σdrag, λdrag) fell in the region of E1 or E2 (defined above; Fig. 2).
The environmental stress factor (ESF) for each kelp was estimated as σbrk/σdrag, where σbrk is the mean breaking stress of the stipe tissue of that kelp.
To compare how kelp from different sites might fare under the same flow conditions, we also calculated the σdrag and the λdrag for a velocity of 0.5 m s−1 for each kelp from each of the three sites.
Allometric analysis
Standard least-squares linear regression will generally underestimate the slope (b) for allometric relationships with a low r2 because of error in the measurement of the independent variable (LaBarbera, 1989). Reduced major axis estimates are often used to avoid such underestimation of the slope; however, the use of reduced major axis is not recommended because it can yield meaningless results, such as a slope between two uncorrelated variables (Harvey and Pagel, 1991). Alternatively, the effect on the slope of the error in the independent variable can be estimated if the measurement error of the independent variable is known. We used the reliability ratio, κ (where κ=r, and r is the correlation between repeated measurements of a variable; Fuller, 1987) to estimate the measurement error in our independent variables, and then used this κ to correct the slopes of our regression lines, as described below. To determine κ for drag measurements, we made three separate measurements of drag on each kelp, and we compared regressions of the first set with those for the second and the first set with the third set. The value κ was calculated as the mean r for these two regressions. To determine κ for blade area, we measured blade area using two different techniques (described above) and calculated the regression of the results of technique one with those of technique two. The value κ was determined from the r for this regression. We calculated the corrected slope (β) by multiplying the uncorrected slope (b) of the double logarithmic plot by 1/κ for plots of stipe diameter as a function of drag and of stipe diameter as a function of blade area.
Statistical analysis
Analyses of variance (ANOVAs), t-tests, Kendal’s T and linear regression analyses were performed using StatView software (version 1.03; Abacus Concepts, Inc.) on a MacIntosh IIci. Homogeneity of variances was tested using the Fmax-test. Analysis of covariance (ANCOVA) comparisons were performed using ANCOVA software created by R. Etter (version 1.0).
RESULTS AND DISCUSSION
In this section we will discuss the allometry of N. leutkeana in different habitats, but will then move beyond this standard morphological scaling analysis to consider how the kelp actually perform in the field. To this end, we will first incorporate measurements of the hydrodynamic forces on the kelp into our scaling analysis, and we will then determine whether stipe material properties change with size or habitat. All these data will then be used to consider whether N. leutkeana in diverse habitats show dynamic strain similarity and whether they are scaled in such a way as to maintain a constant environmental stress factor. First, however, we must consider the seasons for which our data for N. luetkeana are relevant.
Deterioration of kelp in the autumn
As described in the Introduction, Nereocystis luetkeana are generally considered to be annual kelp, the main source of mortality being breakage during autumn and winter storms. To examine whether changes in morphology or material properties occurred as kelp were damaged or removed from a kelp bed in the autumn, we measured autumn kelp as well as summer kelp at the protected site (the site that was most easily accessible during bad weather). We found no significant difference in the material properties of stipe tissues between kelp collected in May at this site (data from Koehl and Wainwright, 1977) and those of kelp collected in August (t-tests, all P>0.10). However, by the end of September, the population of kelp remaining at this site were shorter (Table 1) and their stipe tissues had lower breaking stresses, elastic moduli and W/V ratios (Table 2) than did summer kelp. The autumn kelp were also more abraded and bore a heavier load of epiphytes. Therefore, except for those variables whose values did not change between August and September, all between-site comparisons were made using only data for healthy summer kelp from each site.
Morphological comparisons of kelp from different sites Blade morphology
N. luetkeana at the current-swept site had flat, narrow, very long strap-like blades, whereas those at the protected site had significantly wider (Table 1), undulate (i.e. ruffled) blades (Fig. 1B,C). This pattern of strap-like blades in current-swept habitats versus ruffled blades in protected habitats has already been reported for N. luetkeana (Koehl and Alberte, 1988), as well as for a number of other species of macroalgae (reviewed in Koehl and Alberte, 1988; McEachreon and Thomas, 1987; Armstrong, 1989; Jackelman and Bolton, 1990; Gutierrez and Fernandez, 1992). These blade shape differences in N. luetkeana (M. A. R. Koehl, unpublished data) and Laminaria saccharina (Gerard, 1987) are due to plastic growth responses to tensile loads on the blades.
Kelp at the wavy site showed an intermediate morphology: they were undulate, but narrow (Table 1). Since long, flexible kelp in the oscillatory flow associated with waves move with the flow and experience little force before the flow direction reverses (as explained in Koehl, 1984, 1986; and measured by Koehl et al. 1991), it is not surprising that these long kelp moving back and forth in small waves showed the ‘protected’ type of undulate blade morphology.
Stipe morphology
The morphology of stipes from N. luetkeana at all three sites was quite similar. Although other brown macrophytes have been found to have wider stipes (Rice et al. 1985) or larger holdfasts (Norton, 1986) in exposed habitats, there was no significant difference between the narrowest stipe diameters or between the holdfast areas (Table 1) of N. luetkeana from different flow habitats. Furthermore, even though the current-swept kelp had larger blade areas than did the protected kelp, there was no significant difference between these sites in the volume of the gas-filled cavity (Fig. 1D) that provides the buoyancy to hold the blades near the water surface. In contrast, stipe length did vary between sites, being significantly lower at the wavy site. N. luetkeana stipes grow in length until the blades approach the water surface (Hurd, 1916; Foreman, 1970; Nicholson, 1970; Duncan, 1973; Duncan and Foreman, 1980); hence, the shorter stipes at the wavy site simply reflected the fact that the water was shallower at this site.
Principal components analysis
The morphological features of kelp from the different sites were summarized using principal components analysis, as described in Materials and methods section. We included drag force (D) in the principal components analysis because drag reflects the morphology of the blades (large size and high degree of ‘ruffliness’ result in high drag; Koehl and Alberte, 1988). The first principal component (PC1 = 0.89LB + 0.87WB + 0.77D + 0.75d + 0.75LS + 0.65AB) accounted for 61% of the variation in these variables, while the second principal component (PC2 = 0.55WB + 0.50D + 0.14d - 0.30AB - 0.34LB - 0.37LS) accounted for 15% of the variation in these variables (abbreviations defined in Table 1). The two variables that carried the most weight in the first principal component (LB, WB) are both related to the size of the blades, so we call PC1 the ‘size variable’. The two variables that carried the most weight in the second principal component (WB, D) are both related to blade shape (wide blades are undulate in shape, and blades with high drag are undulate; Koehl and Alberte, 1988), so we call PC2 the ‘shape variable’.
A plot of PC2 as a function of PC1 (Fig. 4) shows that kelp from different sites cluster in different regions of this morphospace. Kelp from the shallow wavy site were smaller (had significantly lower PC1 values) than were those from the other two deeper sites, and autumn kelp were significantly smaller than summer kelp from the protected site. Blades at the current-swept site were flatter and narrower (had significantly lower PC2 values) than were those at the other two sites (Table 1).
The relationship of PC2 to PC1 gives a rough summary of how shape changes with size for N. luetkeana from the different sites. Shape was independent of size for kelp from the protected site, whereas PC2 increased with size for kelp from the wavy site (for interpretation, see discussion of blade morphology above). In contrast, at the current-swept site, larger kelp had lower PC2 values (i.e. had flatter, narrower blades) than did smaller kelp (statistics reported in Fig. 4). Hence, at the current-swept site, where hydrodynamic forces on N. luetkeana are likely to be the largest, the kelp took on a more streamlined blade morphology as they grew (see discussion of drag below).
Allometry
Since a N. luetkeana stipe is a tensile structure, stipe diameter would not be expected to scale with stipe length (Wainwright, 1970; Peterson et al. 1982). Indeed, we found that stipe diameter was independent of stipe length at all three sites (linear regression analysis: current-swept, P=0.59, N=9; wavy, P=0.24, N=21; summer protected, P=0.35, N=7; autumn protected, P=0.49, N=11).
For both elastic similarity and stress similarity between small and large tensile structures, the diameter of the structure should be proportional to the force on the structure raised to the power 0.5 (Peterson et al. 1982). Since the mechanical load on the stipe of a N. luetkeana is due to the drag force on the blades of the kelp (Koehl and Wainwright, 1977), stipe diameter might be expected to scale with AB0.5 (Peterson et al. 1982). The results of our analysis, however, indicate that, according to this criterion, stipe diameter (d) is underscaled with respect to blade area (AB) at all sites: at the protected site, d∝AB0.1; at the wavy site, d∝AB0.2; and at the current-swept site, d∝AB0.3 (Fig. 5). All these exponents are significantly less than 0.5 (statistics reported in Fig. 5). Similarly, Peterson et al. (1982) found for the stipes of the elk kelp Pelagophycus porra that d∝AB0.268, and several species of red algae also show such underscaling of stipe diameter with blade area (Carrington, 1990; Dudgeon and Johnson, 1992). Thus, it appears that N. luetkeana stipe diameters do not increase enough as the kelp grow to maintain elastic or stress similarity and that this underscaling is most pronounced at the more protected site.
Although the preceding analysis suggests that N. luetkeana are not scaled according to elastic or stress similarity, this analysis was based on the commonly used assumptions that (1) the mechanical load on the structure was proportional to some morphological variable (in our case, blade area), and (2) the structure’s material properties were constant. In the following sections, we will explore the validity of each of these assumptions.
Drag
Drag on kelp of different morphologies and sizes
At a given velocity, N. luetkeana from the current-swept site experienced lower drag than did kelp of comparable blade area from the protected or the wavy habitats (Fig. 6A).
The CD values of N. luetkeana (Fig. 7) are similar to those reported for other species of deformable seaweeds (e.g. Denny, 1988; Carrington, 1990; Dudgeon and Johnson, 1992). The current-swept kelp were more streamlined (i.e. had lower drag coefficients, CD; Fig. 7) than those from the other sites because the narrow, flat blades on a current-swept kelp flapped with lower amplitude and collapsed together into a narrower bundle in flowing water than did the undulate blades of kelp from the other sites. (This phenomenon was quantified by Koehl and Alberte, 1988.) While several other studies of macroalgae have also shown that thallus shape affects drag (Armstrong, 1987; Koehl and Alberte, 1988; Sheath and Hambrook, 1988; Dudgeon and Johnson, 1992), Carrington (1990) found that shape did not affect CD values of very collapsible seaweeds at the high velocities characteristic of shores exposed to large breaking waves.
Vogel (1984) has proposed a number called the ‘figure of merit’ (B) to describe the relative reduction in drag experienced by flexible structures as they reconfigure as flow velocity increases, where B is the slope of a double logarithmic plot of speed-specific drag (D/U2, where D is drag and U is velocity) as a function of velocity; the greater the absolute value of the negative slope, the greater the relative drag reduction experienced with an increase in velocity. We determined B for each individual kelp by calculating the slope of a linear regression of a plot of log(D/U2) versus logU. N. luetkeana from the protected site had significantly steeper values of B (mean B=-1.2, S.D.=0.28, N=16) than did those from the current-swept site (mean B=-0.85, S.D.=0.16, N=8) or the wavy site (mean B=-0.75, S.D.=0.25, N=19) (ANOVA and Fisher PLSD, P<0.05). These results are similar to those obtained by Armstrong (1989), who found for the kelp Hedophyllum sessile that ruffled individuals from a protected site had steeper values of B than did strap-like individuals from a more exposed site. In the case of N. luetkeana, the strap-like kelp from the current-swept site collapsed into streamlined bundles at lower velocities than did the ruffled kelp; therefore, the kelp with strap-like blades did not show as much drag reduction with further increases in velocity as did the ruffled kelp (which had much greater initial drag and which required higher water speeds to achieve passive streamlining). Our values for B are comparable with those reported for other macroalgae: H. sessile (-0.57 to -1.2; Armstrong, 1989), Sargassum filipendula (-1.06 to -1.47; Pentcheff, cited in Vogel, 1984); various species of small intertidal brown and red algae (-0.28 to -0.76; Carrington 1990); and freshwater red algae (-0.33 to -1.27; Sheath and Hambrook, 1988).
Kelp from all three sites showed a decrease in CD as blade area increased (i.e. the bundles of blades on large kelp were more streamlined than those on small kelp) (Fig. 7). Since drag is proportional to the product of CD and blade area, our data indicate that the consequences for drag of growth in blade area are greater for small kelp than for large ones.
The total drag per kelp at a given velocity (0.5 m s−1) increased with blade area for kelp from all three sites, with the rate of increase being greatest for kelp from the wavy site and lowest for kelp from the current-swept site (Fig. 6A). Therefore, although blade area is a good predictor of the force that a N. luetkeana stipe must bear, the relationship of force to blade area also depends on the morphology of the blades, which is site-dependent.
Drag forces borne by stipes in different habitats
Not only do blade morphologies differ between sites, but ambient current velocities differ as well. In this study, we could only assess the forces on stipes of kelp from the sites subjected to unidirectional flow, since kelp from the wavy site can move with the flow and hence transmit little force to their stipes (see discussion of blade morphology above). When we considered drag at the relevant ‘typical’ peak velocity for each of the sites subjected to unidirectional currents (0.5 m s−1 at Shady Cove and 1.5 m s−1 at Turn Rock; determination of typical peak velocities explained above in Materials and methods section on field sites), we found: (1) that the drag forces on kelp bearing similar blade areas were comparable at the two sites, and (2) that the rate of increase of drag with an increase in blade area did not differ between the sites (Fig. 6B). Therefore, the load that a N. luetkeana stipe must bear due to the drag on the blades is similar between different flow habitats because the kelp adjust their blade morphology.
The relationship of stipe diameter (d) to site-relevant drag (drag at the water velocity relevant to that site, Ds) is plotted in Fig. 8 for kelp from the sites exposed to unidirectional flow. For elastic similarity and for stress similarity, we would expect d∝ Ds0.5 if material properties were constant (Peterson et al. 1982). Indeed, for kelp from the current-swept site, the relationship of d to Ds is not significantly different from d∝Ds0.5 (statistics reported in Fig. 8). However, for kelp from the protected site, d∝ Ds0.2, which is significantly different from the expectation for elastic or stress similarity if we assume that stipe material properties are constant. In the next section, we shall examine this assumption.
Material properties
A typical stress/extension ratio curve (σ/λ curve) for N. luetkeana stipe tissue is given in Fig. 2, and the material properties for stipe tissues from the three sites are summarized in Table 2. Our measurements of stipe material properties are consistent with those reported by Koehl and Wainwright (1977) for N. luetkeana. This kelp’s stipe tissue is stiffer at low extension ratios (λ) than at high values of λ, with the transition in modulus (from E1 to E2) beginning at λ values in the range 1.1–1.2. The stipe tissues of another kelp, Postelsia palmaeformis, also have σ/λ curves of this shape (Holbrook et al. 1991).
When N. luetkeana stipe tissue was pulled to λ values within the E1 range, it was very resilient (Wrecoil/Wpull=0.93, S.D.=0.06, N=6) and maintained stress until it was returned to its resting length (Fig. 3A). The high resilience of N. luetkeana stipes at λ values in the E1 region permits them to snap back to resting length between bursts of high force (see discussion below of N. luetkeana as shock absorbers).
In contrast, stipe tissue pulled to λ values in the E2 range was significantly less resilient (Wrecoil/Wpull=0.52, S.D.=0.05, N=19; t-test, P=0.0001), and stress in the tissue dropped to zero before the specimens were fully returned to resting length (Fig. 3B). Vincent and Gravell (1986) also found hysteresis when they pulled Laminaria blade tissues to high extensions. Although highly resilient materials (such as N. luetkeana tissue in the E1 region of its σ/λ curve) are brittle because they store strain energy that can be readily transferred to an advancing crack (see Biedka et al. 1987; Denny et al. 1989; Vincent, 1990; Holbrook et al. 1991), the reduced resilience of N. luetkeana stipe in the E2 range of its σ/λ curve may lessen such susceptibility to brittle failure. Furthermore, the plastic deformation and great extensibility of N. luetkeana stipe tissue in the E2 region of its σ/λ curve may help to blunt the tips of cracks, which should reduce local stresses at the crack tips (Gordon, 1968). Indeed, crack-tip blunting has been shown to reduce susceptibility to brittle failure in the blades of the red macroalga Iridaea flaccida (Denny et al. 1989).
The low elastic moduli (E values) and breaking stresses (σbrk values) we measured for N. luetkeana stipes (Table 2) are comparable with those of stipe and blade tissues of other brown and red macroalgae (Delf, 1932; Charters et al. 1969; Koehl, 1982, 1986; Vincent and Gravell, 1986; Armstrong, 1987; Denny et al. 1989; Carrington, 1990; Holbrook et al. 1991; Lowell et al. 1991; Dudgeon and Johnson, 1992). Although the breaking extensions (λbrk values) of N. leutkeana stipes are at the high end of the range reported for other species of macroalgae (Koehl, 1986; Armstrong, 1987; Sheath and Hambrook, 1988; Denny et al. 1989; Holbrook et al. 1991), the work per volume (W/V) required to break N. luetkeana stipes is slightly lower than that measured for other kelp (Holbrook et al. 1991; Lowell et al. 1991). When compared with other biomaterials, the stipe tissues of N. luetkeana and other species of macroalgae are neither stiff nor strong (i.e. they have low E values and σbrk values; Wainwright et al. 1976; Koehl, 1986; Denny et al. 1989), nor are algal tissues tough (i.e. they have low tearing energy; Vincent, 1990; Denny et al. 1989). However, the high extensibility (i.e. high λbrk) of algal stipes renders these structures able to absorb a lot of strain energy before failing (i.e. they have relatively high values of W/V when compared with other biomaterials; Wainwright et al. 1976; Koehl, 1986).
Comparisons of Nereocystis luetkeana from different sites
The results of an ANOVA comparing material properties of stipes from N. luetkeana from the different sites are summarized in Table 2. Stipe tissues of kelp collected during the summer from the protected site were significantly stronger (σbrk) than were those from the current-swept and wavy sites. Although there was no significant difference in breaking extension ratio (λbrk) between sites, the work per volume required to break a stipe was higher for the protected kelp than for those from the other two sites. Stipes from the protected habitat were stiffer at low extensions (E1) than were those from the other sites, whereas at higher extensions (E2), stipes from both the protected and current-swept habitats were stiffer than stipes from the wavy habitat. Because our results indicate that the material properties of N. luetkeana can differ between sites, we cannot assume a single value for such variables as modulus or strength when conducting allometric analyses for kelp from different habitats. Other species of kelp also show differences in material properties in different water-flow habitats (Charters et al. 1969; Armstrong 1987; McEacheron and Thomas, 1987). Furthermore, several studies have shown that such changes in kelp material properties can be environmentally induced (Kraemer and Chapman, 1991; Lowell et al. 1991).
Relationship of material properties to thallus size
There was no correlation between the various material properties of stipe tissues we measured and the blade areas (a mechanically relevant measure of thallus size) of the N. luetkeana kelp from which they were taken (Table 3). The only exception to this was the positive association between breaking extension ratio and blade area for the kelp from the wavy site. However, this result could have been due to chance, since the P<0.05 criterion for significance implies that one out of twenty comparisons may yield a correlation simply due to chance alone (we ran twenty statistical tests and found only one significant association). Therefore, the assumption that N. luetkeana stipe tissue properties remain constant as the kelp grow appears to be valid (although as kelp age and become damaged in autumn their material properties deteriorate, as described above).
Scaling of Nereocystis leutkeana in different water-flow habitats
Since stipe material properties, blade morphology and ambient flow conditions differ between habitats for N. luetkeana, we must consider the stresses and extensions experienced by stipes in nature to evaluate the mechanical consequences of size. At site-relevant peak velocities, σdrag and λdrag did not change with size (Fig. 9A,B). However, although stress due to site-relevant drag (σdrag) did not differ between sites, extension ratio due to site-relevant drag (λdrag) was greater at the current-swept site than at the protected site (Table 2).
Peak stresses in the stipes of N. luetkeana are maintained, not by allometric scaling of stipe dimensions, but rather by adjustments in blade morphology in different flow habitats. Therefore, we see within this species of kelp a phenomenon analogous to the ‘dynamic strain similarity’ (sensuRubin and Lanyon, 1984) between species of vertebrate animals: peak stresses in the bones of vertebrates representing a range of body sizes are maintained, not entirely by allometric scaling of bone dimensions, but rather by the use of more gentle gaits by large animals. Vertebrates (for which the loads on limb bones are mainly due to locomotion) can adjust their behavior, whereas kelp (for which the loads on stipes are mainly due to environmental hydrodynamic forces) can adjust their morphology. Similarly, sessile animals subjected to hydrodynamic forces can alter the magnitude of those forces by shape changes due to growth (e.g. Hunter, 1988) or behavior (Koehl, 1977). Spiders provide another example of how an organism’s behavior can compensate for an insufficiently robust support structure: although the drag line is not strong enough to bear the force of a falling spider, the ability of the animal to draw drag-line silk from its spinnerets during a fall absorbs enough energy to prevent breakage of the drag line (Brandwood, 1985). All these examples of dynamic strain similarity illustrate the importance of determining the mechanical loads on structures as they are used by organisms in nature.
Nereocystis luetkeana kelp as shock absorbers
N. luetkeana can be subjected to short (<1 s) pulses of rapid water flow in the turbulent currents they encounter in both protected and exposed habitats (Koehl and Wainwright, 1977; Koehl and Alberte, 1988). The low modulus of the stipes of these kelp permits them to ‘give’ when subjected to a burst of high drag (Koehl and Wainwright, 1977; Koehl, 1984, 1986). We found that N. luetkeana stipes in the field operate in the E1 range of their stress/extension ratio curves (all 37 tests for nine kelp from the current-swept site and all 72 tests for eight summer kelp from the protected site), where the tissues have very high resilience (Fig. 3A). This high resilience permits them to snap back to resting length before the next pulse of high force. In contrast, some deteriorating autumn kelp operated in the E2 region of their stress/extension ratio curves (four out of 11 kelp tested, eight out 37 mechanical tests), where their resilience is lower (Fig. 3B).
The ability of a N. luetkeana to withstand transient high forces without breaking depends on the mechanical work required to break the stipe. There are two mechanisms by which the performance of a stipe as a ‘shock-absorber’ can be improved: the tissues can require a high work per volume (W/V) for breakage or the stipe can be long (i.e. the total volume of tissue absorbing the strain energy is large). There was no significant difference in the length (which depends on water depth; Duncan, 1973) or total tissue volume of the stipes from the protected and the current-swept sites (Table 1) but, because the protected stipes had a higher W/V, the total work required to break healthy summer kelp was higher at this site than at the current-swept site (Table 2). Although N. luetkeana from the shallow wavy site had the same W/V as the current-swept kelp (Table 2), we would expect these shorter kelp from shallow water to be less effective shock-absorbers than the longer N. luetkeana from the other two deeper sites. It is interesting to note that at the current-swept site there was a positive association between the total work required to break the stipes (W) and the blade area (AB) that the stipes supported (Table 3).
Environmental stress factors
The environmental stress factor of a component of an individual organism at a particular stage in its life is the ratio of the breaking stress of that component to the peak stress that it experiences at that stage in its habitat. In the case of N. luetkeana, we used stipe breaking stress (σbrk) and the stress due to drag (σdrag) at the peak site-relevant velocity (defined above) to calculate environmental stress factor for each kelp from the current-swept and the protected sites (the sites exposed to unidirectional currents). Environmental stress factor was independent of site (Table 2) and of blade area (Fig. 9C), suggesting that the blade morphology, stipe diameter and tissue properties of N. luetkeana were scaled in such a way that environmental stress factor was maintained between habitats and as the kelp grew.
Cost of maintaining the environmental stress factor (Cs). N. luetkeana maintain the environmental stress factor, in part, by the streamlining of blades. Such streamlining reduces photosynthetic rate due to self-shading as moving water pushes the strap-like blades on top of each other into a bundle (Koehl and Alberte, 1988). Reduced photosynthesis should translate into reduced blade growth rate and, hence, into a lowering of spore production rate. Furthermore, trade-offs between tissue strength and photosynthesis have been noted for other species of macroalgae (reviewed by Koehl, 1986); however, such costs have not been assessed yet for N. luetkeana.
In contrast to N. leutkeana, which maintain their environmental stress factor as they increase in size during the summer, several species of intertidal red algae show a decrease in this factor with increasing size (Carrington, 1990; Dudgeon and Johnson, 1992). Unlike N. luetkeana, these red algae can regenerate from their holdfasts after breakage; hence, their cost of failure (Cf) might be expected to be lower. Moreover, among these red algae, a species that regenerates slowly, Mastocarpus stellatus, has a higher environmental stress factor than does a co-occurring species, Chondrus crispus, that regenerates more quickly (i.e. that may have a lower Cf) but that is broken more frequently (and therefore has a higher P) in winter storms (Dudgeon and Johnson, 1992). Koehl (1984, 1986) has noted a similar pattern among kelp species, where those with weak tissues tend to have life histories characterized by rapid growth, early reproduction and susceptibility to seasonal dislodgment. Although these observations are consistent with the predictions of Alexander’s (1981, 1988) theory, we cannot evaluate whether they represent optimization of the environmental stress factor until the reproductive output or success of these macrophytes has been quantified.
Conclusions
We have investigated the mechanical scaling of Nereocystis luetkeana, a giant kelp whose stipes are loaded in tension due to hydrodynamic forces on the blades of the kelp. If we simply consider the allometry of the morphological characters of these kelp, we conclude that the stipes are underscaled relative to the blade area they support (i.e. that the kelp do not maintain stress or elastic similarity as they grow). However, we must examine the assumptions underlying such an allometric analysis: (1) that load on a support structure (the stipe) is simply proportional to a morphological feature (blade area), and (2) that material properties of the stipe are constant. We found that kelp alter their blade morphologies (hence their drag coefficients) and their stipe material properties in different hydrodynamic environments such that they maintain elastic similarity and environmental stress factor as they grow. Thus, sessile organisms such as these kelp can, via their growth responses to environmental stresses, show dynamic strain similarity. Our study illustrates the importance of moving beyond simple allometric analyses of the mechanical support structures of organisms: if we are to understand the functional consequences of changes in organism size, we must also measure how organisms actually perform in the field in different habitats.
Acknowledgements
This research was supported by N.S.F. Grants OCE-9217338 and OCE-8352459 to M.A.R.K. and by a NATO Postdoctoral Fellowship and a Northeastern University Postdoctoral Fellowship to A.S.J. We are grateful to A. O. D. Willows for the use of facilities at Friday Harbor Laboratories, to J. Ahouse for help with data analysis, to D. Pentcheff for the use of his unpublished flow data, to P. Spowart for her sketches and to S. Wainwright for collaboration on the resilience measurements. We thank the members of the Berkeley biomechanics group for helpful discussions and advice, O. Ellers and S. Worcester for their critical reading of this manuscript and M. Denny and an anonymous reviewer for helpful suggestions.