Wood density as a proxy for vulnerability to cavitation : Size matters

In this study, vulnerability to cavitation, P50 (i.e. the water potential causing 50 % loss of hydraulic conductivity), of Norway spruce trunkwood at different cambial age was related to wood density. Wood density was calculated from mass in the oven dry state related either to volume at the oven-dry state (dry wood density) or to volume at full saturation (basic wood density). Dry wood density and basic wood density were strongly linearly related (r2 = 0.99); there was however a shift from the 1:1 reference line with increasing dry wood density. Dry wood density as well as basic wood density had similar high predictive qualities for P50 (r2 = 0.85). The quadratic regression lines took however a quite different course below -4 MPa because volume shrinkage increased with increasing dry wood density. For Norway spruce sapwood with high dry wood density, mixing up different wood density traits would thus result in a predicted overestimation or underestimation of vulnerability to cavitation. Gravimetrically wood density measurements at different moisture contents (starting at full saturation) are easily to achieve on standard size specimens and the conversion curves obtained will be of high value for future ecological studies on other species and across species.


Introduction
Wood density is regarded as a functional trait in ecophysiological research and has been used so far in many metaanalyses in order to learn more about the different strategies of tree species to survive drought (e.g.Chave et al., 2009;Schuldt et al., 2013;Anderegg et al., 2016;Gleason et al., 2016;Jupa et al., 2016).In conifers, strong relationships between vulnerability to cavitation and wood density were found across species (Hacke et al., 2001;Hacke and Jansen, 2009;Lachenbruch and McCulloh, 2014) and within a given species (Domec et al., 2009;Rosner, 2013, Rosner et al., 2014).Several studies report however either weak or no such within-species- (Martínez-Vilalta et al., 2009;Corcuera et al., 2011;Lamy et al., 2012) or interspecies correlations (Gleason et al., 2016).Hacke et al. (2001) hypothesized that greater resistance to cavitation requires a safer design for resisting implosion since cell walls have to withstand higher tensile strain before the water column breaks.Based on the assumption that mechanical strength and stiffness increase with increasing wood density, tensile stress in a water filled conduit at a given water potential should increase with decreasing squared double wall/lumen ratio (Hacke et al., 2001;Pittermann et al., 2006;Domec et al., 2009).Wood density is affected by anatomical and chemical traits that determine how much mass can be found in a given volume.These relationships are much simpler in conifers than in angiosperms.Conifer wood consists mainly of tracheids and has less than 10 % parenchyma tissue (Morris et al., 2016) and its wood density is strongly related to latewood percentage (Domec et al., 2009).Angiosperm wood can have -even within a narrow density range -variable amounts of fibre wall-, fiber lumen-and parenchyma fractions (Ziemińska et al., 2015).Moreover, in some species, wood density is mainly driven by the density of tissues outside the vessels rather than by the vessel lumen fraction (Ziemińska et al., 2013).The characteristics of the vessels (diameter, wall thickness) should however determine hydraulic functioning (Awad et al., 2010).Density has no direct effect on hydraulic functioning; empirical relationships exist because some of these traits that affect density also affect hydraulics (Lachenbruch and McCulloh, 2014).Thus, the statistical relationship between vulnerability to cavitation and wood density is no direct relationship.In fact, in conifers, vulnerability to cavitation is related to the characteristics of the membranes of bordered pits (Domec andGartner, 2002b, Delzon et al., 2010;Pittermann et al., 2010;Jansen et al., 2012;Bouche et al., 2014).Seal capillary-seeding, i.e. when the torus is not tightly sealed against the pit border, either due to weak aspiration or poor air tightness of the torus/pit aperture interface, is regarded as the © The Authors most likely mode of air-seeding in bordered pits of conifers (Delzon et al., 2010).Earlywood tracheids with thicker cell walls may have bordered pits that are less prone to seal capillary-seeding.In that regard, " It is entirely appropriate to take advantage of these indirect correlations as long as the proxy is presented as what it is, and the correlation is expected to be valid over the conditions used."(Lachenbruch and McCulloh, 2014).
In ecological studies, the functional trait "wood density" is referred to "basic wood density", that is the wood mass in the oven dry state divided by the volume in the green state.Basic wood density datasets for many species have been made available (e.g.Zanne et al., 2009).On the other hand, there is an increasing trend to analyze within-ring differences of wood density by means of x-ray transmission and relate that to hydraulic performance of a given species (e.g.Dalla-Salda et al., 2014;Martinez-Meier et al., 2015;Wilkinson et al., 2015).Wood density assessed by X-ray transmission is based on the wood mass in the oven dry state divided by the volume in the oven dry state (or at equilibrium moisture content below fiber saturation) and is thus a different trait than "basic wood density".In wood technology research, density where mass and dimensions are measured either in the oven dry state or at equilibrium moisture content (in general 12%) is traditionally termed "(dry) wood density".The incompatibility of "basic wood density" and "dry wood density" is due to dimensional changes during moisture loss that become more relevant with increasing density of a given wood specimen (Simpson, 1993) because "moisture-induced dimensional change increases as wood density increases."(Forest Products Laboratory, 2010).Moreover, dimensional changes in detached sapwood can be observed above fiber saturation (i.e. about 30 % moisture content), thus when there is still free water in the conduit lumen (Rosner et al., 2009).
Aim of this study is to emphasize the importance of a clear terminology when "wood density" is used as a functional trait or as a proxy for hydraulic vulnerability and to create awareness for its high variability within a conifer species.P50 (i.e. the water potential causing 50 % loss of hydraulic conductivity) of Norway spruce (Picea abies [L.] Karsten) trunkwood of different cambial age from mature and young trees was related to wood density calculated either from volume at the oven dry state (dry wood density) or from volume at full saturation (basic wood density) and the results are compared.

Plant material
The dataset for P50, dry wood density and basic wood density of Norway spruce (Picea abies (L.) Karst.)comprised existing datasets as well as new measurements.Information on the origin of the plant material, age of trees and sampling location within the trees can be found in Table 1.The parameter set is described in detail in Table 2.

Vulnerability to cavitation (P50)
P50 (air pressure that results in 50 % conductivity loss) was assessed by the air injection method.The air injection technique enables the construction of a whole vulnerability curve on a given sample within a few hours (Cochard et al., 2013) up to one (working) day.Sapwood specimens of mature and young trunkwood had a size of 6 mm (radial), 6 mm (tangential) and 110 mm (longitudinal).Specimen from small trees had a length of 110 mm and a variable diameter that was smaller than 7 mm.Specimen for hydraulic testing comprised 1-3 complete annual rings, those were rings 1-3 in specimens from Austria and Germany, rings 8-10 in specimens from Norway and rings 17-19 in specimens from Sweden.
For hydraulic testing, specimens from small trees were debarked (more details in Rosner et al., 2007); whereas small wood beams were produced form the trunkwood samples of older trees with a chisel and by planing on a sliding microtome (Rosner et al., 2008).Specimens were kept wet during all preparation steps.Wood specimens were soaked in distilled water under partial vacuum at 4°C for at least 48 h to refill embolized tracheids.Hydraulic flow measurements were carried out under a pressure head of 8.0 kPa (specimen from small trees) or 5.4 kPa (wood beams) with distilled, filtered (0.22 µm), and degassed water containing 0.005 vol.% Micropur (Katadyn Products Inc., Wallisellen, Switzerland).After the first conductivity measurement at full saturation, air overpressure was applied by means of a double-ended pressure chamber (PMS Instruments Co., Corvallis, Oregon).After a relaxation period in distilled water of 30 minutes, hydraulic conductivity was measured again.The first pressure application was 1.0 MPa, thereafter the pressure was subsequently increased in steps of 0.5 -1.0 MPa until more than 90 % loss of conductivity was reached.P50 values for each sapwood beam were calculated as described in Pammenter and Vander Willigen (1998).In specimen from Norway, conductivity measurements were performed only until 70% loss of hydraulic conductivity, since that proved sufficient to calculate reliable P50 values for each tree segment (Rosner et al., 2014).

Basis wood density and dry wood density
Basic wood density was defined as the mass (weight) in the oven dry state divided by the volume in the fully saturated green state [kg m -³].Dimensions for calculating the volume at the fully saturated green (never dried) state were measured with a caliper (resolution of 0.1 mm) after saturation at partial vacuum under water (Hietz et al., 2008).The volume was determined before the hydraulic flow measurements.Dry wood density was defined as the volume in the oven dry state divided by the mass in the oven dry state [kg m -³].Specific gravity is often mixed up with density.Specific gravity of wood can be calculated by dividing dry wood density by the density of water at 4.4 °C (1000 kg m -³).From the standpoint of physics, specific gravity is defined as the density of wood where mass and volume are determined in the same (oven dry) state relative to the density of water (Williamson and Wiemann, 2010).After flow measurements, specimens were dried for 48 hours at 103°C in order to determine weight and volume at the oven dry state.Lower temperatures can result in remainders of bound water in the cell walls (Williamson and Wiemann, 2010).After drying, specimens were immediately transferred from the oven into an exsiccator containing dried silica gel.As soon as the wood specimens had reached room temperature, they were weighed (one by one, the rest remained in the exsiccator) and their dimensions were determined by means of a caliper.Volume shrinkage was calculated as the volume loss from the fully saturated state to the oven dry state related to the volume in the fully saturated state (in %).

Results
Stem wood from small trees had much lower P50 values (-4.88 ± 0.62 MPa) than sapwood from the crown (-2.55 ± 0.45 MPa) or from the lower trunk (-2.44 ± 0.21 MPa).Dry wood density and basic wood density were strongly linearly related (Fig. 1).With increasing wood density there was however a shift from the 1:1 reference line.The equation for calculating basic wood density from dry wood density of Norway spruce sapwood across cambial age was: Basic wood density = 37.39 + 0.815 * dry wood density [kg m -³] [Equation 1] The relationship between P50 and dry wood density as well as basic wood density is shown in Fig. 2. Basic wood density data for young trees from Germany were calculated from dry wood density data with Equation 1. Due to the shift from the 1:1 reference that increased with increasing wood density (Fig. 1), the quadratic regression lines took a quite different course below -4 MPa.For sapwood with high wood density, mixing up different wood density traits would result in an overestimation (basic wood density used in an empirical model for P50 and dry wood density) or underestimation (dry wood density used in an empirical model for P50 and basic wood density) of vulnerability to cavitation.Overestimation is meant as the prediction of higher vulnerability to cavitation, thus a lower hydraulic safety of up to 2 MPa (Fig. 3), underestimation the prediction of a much higher hydraulic safety.The relationship between dry wood density and volume shrinkage of Norway spruce sapwood is presented in Fig. 4. Volume shrinkage and thus the shift from the 1:1 relationship between dry wood density and basic wood density (Fig. 1) increased with increasing dry wood density.P50 predicted for wood density data using either empirical models for basic wood density or dry wood density.Regression lines for the models are given in Fig. 2. The quadratic regression is significant at the 0.01 % level (P < 0.0001, n = 94).Light green dots represent young trees (Austria), green dots young trunkwood (Norway) and dark green dots mature trunkwood (Sweden).

Discussion
Basic wood density was strongly linearly related to dry wood density in Norway spruce sapwood (Fig. 1); shifts from a 1:1 relationship could be explained by differences in volume shrinkage (Fig. 4).Volume shrinkage of Norway spruce sapwood presented in this study (Fig. 4) corresponded to that reviewed in literature (Kärkkäinen and Marcus, 1985).Wood structure is variable within a tree and between different trees and so is volumetric shrinkage.A general concept is that volume shrinkage below fiber saturation increases with dry wood density (Simpson, 1993;Forest Products Laboratory, 2010).Fiber saturation is defined as a state during a dehydration cycle when cell lumina contain no longer free water but when cell walls are fully saturated with liquid (Stamm, 1971;Skaar, 1988;Berry and Roderick, 2005).
When conduits reach moisture contents below fiber saturation, the cell wall shrinkage process is initiated.Radial and tangential wall shrinkage, summarized as transverse shrinkage, account for most of the volumetric shrinkage, whereas shrinkage in the longitudinal direction is often much below 1% (Kärkkäinen and Marcus, 1985).Transverse cell wall shrinkage decreases with decreasing dry wood density (Stamm, 1971;Skaar 1988;Chafe and Ilic, 1992;Pang and Herritsch, 2005).Cellulose microfibrils of the tracheid walls are embedded in an amorphous matrix.This matrix tends to contract isotropically when it loses water.The microfibrils deform the matrix when resisting this contraction (Barber and Meylan, 1964).In this study, highest wood density values were found for stem wood of very young Norway spruce trees.Strictly spoken, wood density of these specimens is "stem density without bark" because such specimens (also branches) include pith tissue.However, the linear relationship between dry wood density and basic wood density was not influenced by the presence of pith tissue.The regression line drawn across cambial age had exactly the same course as those for wood of young trees and older trees (Fig. 1).Care was taken that only specimens with a small amount of pith were used for the study (Rosner et al., 2014).In addition, the presence of the pith in specimens from small trees led to no exponential increase in volumetric shrinkage (Fig. 4).Rather the opposite was the case, as indicated by the depletion of the quadratic regression curve when volumetric shrinkage was plotted against dry wood density.For young branch wood or trunkwood of species with a very high amount of pith tissue results might differ.
P50 showed a quite wide range in Norway spruce sapwood specimens.This is not surprising, because sapwood with extreme differences in wood structure originating from small trees (Fig. 5a) and the mature trunk (Fig. 5b) were investigated.When conifer branches (e.g.Dalla-Salda et al., 2009;Corcuera et al., 2011;Lamy et al., 2014), young stems (e.g.Rosner et al., 2007;Sáenz-Romero et al., 2013) or mature trunkwood of similar cambial age (e.g.Rosner et al., 2008;Dalla-Salda et al., 2009) is tested much narrower P50 ranges are observed.Within-species variation in P50 of sapwood of a given organ at similar cambial age can be a result of genetic determination and site conditions (Rosner et al., 2007;Rosner et al., 2008;Corcuera et al., 2011;Sáenz-Romero et al., 2013) or driven by sudden drought events (Rosner et al., 2014).Differences in P50 of a given organ observed so far made up however not much more than 1.5 MPa in conifers (e.g.Dalla-Salda et al., 2009;Lamy et al., 2011;Corcuera et al., 2011;Rosner et al., 2016) and in angiosperms (e.g.Awad et al., 2010).Some conifer species, especially pines, exhibit a much lower genotypic variability and phenotypic plasticity (Sáenz-Romero et al., 2013;Lamy et al., 2014) than Norway spruce (Rosner et al., 2007;2008).As mentioned above, the sample set comprised trunkwood from young small trees with a very high hydraulic safety and sapwood from the mature trunk that had much lower hydraulic safety.High intraspecific variability in P50 has been reported for other conifer species too; Domec et al. (2009) found P50 values in Douglas fir ranging in roots from -1.4 MPa to -1.0MPa, in the main trunk from -4.1 to -3.2 MPa and in branches from -6.1 to -5.9 MPa.In Ponderosa pine, P50 ranged from -1.5 to -0.7 MPa in roots, from -4.1 to -3.2 MPa in the main trunk and from -4.7 to -4.5 MPa in branches.Recently, Bouche et al. (2016) found however low inter-organ variability in P50 in four different Pinaceae species, including Douglas fir.Stem wood of young small spruce trees (seedling size) investigated in the present study had a mean P50 of -4.88 MPa.This is slightly lower than P50 of -4.38 MPa found by Mayr et al. (2003) for leader shoots of older Norway spruce trees.The high hydraulic safety demand in the sapwood of small young trees might be explained by the shallow root system of seedlings that implies less access to ground water (Rosner, 2013).Accordingly, in Douglas fir seedlings (2nd internode from top) P50 can be much lower (-5.6 MPa) than at the base of the trunk of old trees (-3.3 MPa) of the same species (Domec and Gartner, 2002a).  1 and 2.
Whereas strong relationships were found between P50 and basic wood density within conifer species across cambial age (Domec et al., 2009;Rosner et al., 2014) no such relationships have been reported for conifer branch wood (Martínez-Vilalta et al., 2009;Corcuera et al., 2011;Lamy et al., 2012).Branch wood contains however variable amounts of compression wood that has higher wood density but also higher hydraulic vulnerability than opposite wood (Mayr and Cochard, 2003).In branches, the mechanical demands might mask the empirical relationship between wood density and P50 (Rosner et al., 2007;Rosner, 2013).In the dataset presented in this study, dry wood density (Rosner et al., 2014) and basic wood density were strongly related to P50 across cambial age in Norway spruce trunkwood (Fig. 2).The regression lines shifted however from each other with increasing wood density.Mixing up different wood density parameters can thus lead to wrong prediction of drought sensitivity of a given species (Fig. 3).Basic wood density is oven-dry mass related to the volume in the green wood state and should be thus a more reliable ecological trait than dry wood density or density at defined equilibrium moisture content (below fiber saturation).However, determination of the volume in the green state can bear sources of errors, because dimensional changes in detached sapwood can be observed already above fiber saturation (i.e. about 30 % moisture content), thus when there is still free water in the conduits (Rosner et al., 2009).In that regard, volume at the "green state" is a rather vague definition and shall be replaced by volume at the fully saturated state.
Across species, conversion formulas between wood densities at variable moisture contents are available (e.g.Forest Products Laboratory, 2010, Chapter 4, Tables 4-6a and 4-6b); the best conversion formula is however the one that is species-specific and that covers the whole wood density range of between-tree and within-tree differences.Such relationships are unavailable yet for most tropical wood species (Williamson and Wiemann, 2010).In species with high wood density, volume shrinkage of more than 20 % can be expected (Forest Products Laboratory, 2010).Conversion formulas have been developed for lumber (mainly mature trunkwood) in order to calculate e.g.shipping weights and they might not take into account within-tree variability.In species where branches or very young trunks have quite high densities (compared to the main mature trunk) these formulas are not applicable.Moreover, the formulas include heartwood and sapwood and do not cover the whole moisture content range that is possible in sapwood of living trees at the fully saturated state.In principle, standardized gravimetrically wood density measurements at different moisture contents (starting at full saturation) are easily to achieve and the conversion curves obtained are of high value for ecological studies.

Figure 1 :
Figure 1: Relationship between basic wood density and dry wood density in Norway spruce sapwood.Lilac symbols present sapwood specimens without pith tissue (wood beams), pink symbols specimens with pith tissue (stem specimen of young trees).The linear regression is significant at the 0.01 % level (P < 0.0001, n = 94).The equation is as follows: Dry wood density = -42.16+ 1.218 * basic wood density [kg m -³].

Figure 2 :
Figure 2: Relationship between wood density traits and vulnerability to cavitation.P50 is plotted against dry wood density (red symbols) and basic wood density (blue symbols).Quadratic regressions are significant at the 0.01 % level (P < 0.0001, n = 130).

Figure 3 :
Figure 3: Hydraulic vulnerability predicted by different empirical models.P50 predicted for wood density data using either empirical models for basic wood density or dry wood density.Regression lines for the models are given in Fig.2.

Figure 4 :
Figure 4: Volume shrinkage related to dry wood density.The quadratic regression is significant at the 0.01 % level (P < 0.0001, n = 94).Light green dots represent young trees (Austria), green dots young trunkwood (Norway) and dark green dots mature trunkwood (Sweden).

Figure 5 :
Figure 5: Transverse section of Norway spruce trunkwood at different cambial age from a young and an older tree.Anatomy of Norway spruce earlywood and latewood formed in the previous growing season from the trunk of a young small tree from Germany (a) and of a mature tree from southern Sweden (b).P50 of the sample in (a) was -5.27 MPa and in (b) -2.20 MPa.In (a) earlywood is from annual ring three from the third whorl from the top, in (b) from annual ring 19 in breast height.Sections with a thickness of 20 µm were stained with methylene blue; colors were inverted in the digitized images.Reference bars indicate 100 µm.Information on study sites and sampling can be found in Tables1 and 2.

Table 1 :
Plant materialOverview of the sample set and the absolute and relative height (internodes from the top) of the trees where the air seeding pressure (P50), dry wood density and basic wood density were measured.

Table 2 : Data sources Overview
For all trees, dry wood density was assessed (n = 134).P50 was calculated only for six trees harvested in Kalwang; the number of available P50 values was thus 130.For trees originating from Germany no basic wood density dataset existed; basic wood density was calculated with an empirical equation (Equ. 1, in the results part) from the relationship between dry-and basic wood density of 94 specimens.
of the data sources comprising air seeding pressure (P50), dry wood density (DWD) and basic wood density (BWD).