Relative Bark Thickness is Correlated with Tree Species Distributions along a Fire Frequency Gradient

Abstract

The probability of stem survival after fire is strongly influenced by energy allocation to bark because bark thickness affects heat transfer during fire. Greater relative investment in inner bark versus outer bark should also enhance survival because of greater moisture content of inner bark. We measured stem diameter, bark thickness, and habitat preference of five species typical of long-leaf pine savannas, and six species characteristic of adjacent wetlands (pocosins), and calculated relative bark thickness, the inner bark proportion, radial growth, and bark accumulation of each species. We hypothesized that savanna species have thicker bark and greater relative investment in inner bark than pocosin species, because fires occur more frequently in savannas than pocosins. As hypothesized, savanna species have relatively thicker bark than pocosin species. Relative bark thickness and the rate of bark accumulation were correlated with the mean location of a species along the pocosin-to-savanna gradient. However, the inner bark proportion did not differ between savanna and pocosin species. Our results indicate that relative bark thickness is likely the primary bark trait affecting fire-induced topkill and influencing the distribution of species along the pocosin-to-savanna gradient.

Resumen

La probabilidad de sobrevivencia del tronco después de un incendio está fuertemente influenciada por la energía que recibe la corteza, dado que el espesor de la misma afecta la transferencia de calor durante el fuego. Una mayor inversión relativa en las estructuras internas de la corteza por sobre las externas aumenta la sobrevivencia, debido al mayor contenido de humedad que poseen las estructuras internas. Nosotros medimos el diámetro del tronco, la proporción de estructuras internas de la corteza y las preferencias de hábitat de cinco especies típicas de sabanas de pino de hoja larga (o pino palustre) y de seis especies características de humedales adyacentes (pocosins), y calculamos el espesor relativo de la corteza, la proporción de estructuras internas de la corteza, el crecimiento radial, y la acumulación de corteza de cada especie. Hipotetizamos que las especies de sabana tienen una corteza más gruesa e invierten más recursos en la creación de estructuras de la corteza interna, dado que los incendios ocurren más frecuentemente en sabanas que en pocosins. Como fuera hipotetizado, las especies de sabana tienen efectivamente la corteza más gruesa que las especies de pocosins. El ancho relativo de la corteza y su tasa relativa de acumulación fueron correlacionadas con la ubicación media de las especies a lo largo del gradiente desde los pocosin hasta la sabana. Por supuesto, la proporción de la corteza interna no difirió entre las especies de sabana y de pocosins. Nuestros resultados indican que el espesor relativo de la corteza es prácticamente el atributo más importante que afecta la muerte apical inducida por fuego, e influencia la distribución de especies a lo largo del gradiente entre la sabana y los pocosins.

Introduction

In pyrogenic ecosystems, thick bark allows trees to minimize or escape fire damage by protecting the vascular cambium (Gignoux et al. 1997, Hoffmann et al. 2009, Midgley et al. 2010, Lawes et al. 2011a) and xylem (Michaletz et al. 2012). Thick bark reduces heat transfer to the cambium during fire and reduces the likelihood that a lethal cambium temperature is reached (Vines 1968, Uhl and Kauffman 1990, Gashaw et al. 2002, van Mantgem and Schwartz 2003, Lawes et al. 2011a, VanderWeide and Hartnett 2011). Thus, the probability of survival after fire is greater for individuals and species with thicker bark (Harmon 1984, Lawes et al. 2011b, Catry et al. 2012, Hoffmann et al. 2012).

Fire frequency and intensity can select for higher investment in bark thickness (Pausas 2014) due to the importance of thick bark for post-fire survival. In both Brazil and Australia, savannas burn more frequently than adjacent forests, and savanna tree species have relatively thicker bark than forest tree species (Hoffmann et al. 2003, Hoffmann et al. 2009, Lawes et al. 2011a). In North America, oak (Quercus spp.) and pine (Pinus spp.) saplings have thicker outer bark (i.e., rhytidome) in ecosystems that experience frequent, low-intensity fire than in ecosystems that experience less frequent, more intense fire (Jackson et al. 1999). Furthermore, anthropogenic alteration of fire regimes can change selection pressures for energy allocation to bark and species composition. Frequent anthropogenic fires may favor trees with thicker bark (Stephens and Libby 2006), while fire suppression may allow for increased survival of species with thinner bark that are less resistant to fire (VanderWeide and Hartnett 2011).

Although thick bark can protect adult trees from surface fires, saplings generally lack sufficient bark to escape topkill (Hoffmann and Solbrig 2003, Lawes et al. 2011b, Hoffmann et al. 2012, Lawes et al. 2013). Consequently, the ability to accumulate bark quickly should be particularly beneficial in frequently burned environments. Since bark thickness increases as a function of stem diameter as a stem grows, the rate of bark accumulation will depend not only on the relative investment in bark, but also on the rate of stem growth (Midgley et al. 2010, VanderWeide and Hartnett 2011, Hoffmann et al. 2012). Allocation of resources to bark, however, may occur at the expense of growth (Jackson et al. 1999), so it is important to consider both stem growth and relative bark thickness when considering plant strategies for survival in frequently burned environments.

In fire-prone ecosystems, bark thickness is likely the best predictor of species differences in topkill after fire (Lawes et al. 2011b), but bark characteristics such as moisture content also influence heat transfer (Harmon 1984). In fact, bark moisture has been shown to increase fire resistance (Higgins et al. 2012, but see Lawes et al. 2011a). If frequent fire selects for increased bark moisture, it could be manifested as an increase in the ratio of inner bark to total bark thickness. Inner bark, which is living tissue, is moister than outer bark (Vines 1968, van Mantgem and Schwartz 2003, Jones et al. 2004), and several studies have suggested that inner bark provides greater insulation than outer bark (van Mantgem and Schwartz 2003, Scholz et al. 2007).

In this study, we tested the hypotheses that tree species in longleaf pine (Pinus palustris Mill.) savannas have relatively thicker bark and a greater proportion of inner bark than tree species in stream-head pocosins (i.e., wetlands). Specifically, we measured relative bark thickness and the proportion of inner bark of tree species in savanna and pocosin habitats in North Carolina, USA. Upland longleaf pine savannas historically had an average fire return interval of approximately two years (Stambaugh et al. 2011). Lowland stream-head pocosins, which support tree species characteristic of eastern deciduous forests (e.g., Acer rubrum L. and Liriodendron tulipifera L.) and occur within the savanna matrix (Schafale and Weakley 1990, Sorrie et al. 2006), burn every 7 to 50 years (Frost 1995). The pocosin-to-savanna hydrological gradient therefore also represents a fire frequency gradient. Thus, in addition to hypothesizing that savanna species have relatively thicker bark than pocosin species overall, we hypothesized that relative bark thickness and bark accumulation are correlated with the location of species along the pocosin-to-savanna gradient. As fire frequency increases, relatively thicker bark and a greater rate of bark accumulation should be beneficial for protecting the vascular cambium during fire (Gignoux et al. 1997, Hoffmann et al. 2009, Midgley et al. 2010, Lawes et al. 2011a). Because inner bark is moister than outer bark (Vines 1968, van Mantgem and Schwartz 2003, Jones et al. 2004), and higher moisture content may inhibit temperature change during fire, we hypothesized that savanna species have a greater proportion of inner bark than pocosin species.

Methods

This study was conducted at Fort Bragg (in Cumberland and Hoke counties) in the Sandhills region of North Carolina. Mean annual precipitation at Fort Bragg is 1275 mm (Sorrie et al. 2006), and soils range from well drained ultisols to poorly drained entisols (USDA NRCS 2007). The most abundant vegetation type is upland pine-scrub oak sandhill (i.e., savanna), which is dominated by longleaf pine and Quercus spp. (Sorrie et al. 2006). Within the savanna matrix, lowland stream-head pocosin wetlands are abundant. Common pocosin tree species include red maple (Acer rubrum), tulip poplar (Liriodendron tulipifera), and pond pine (Pinus serotina Michx.), and the diversity and density of tree species depends on fire frequency (Sorrie et al. 2006). Similarly, the width of the ecotone between longleaf pine savanna and stream-head pocosin communities is influenced by fire frequency and intensity, which are determined in part by fuel moisture content (Weakley and Schafale 1991). Longleaf and pond pines are often the only tree species in the ecotone, but other pocosin species may be present (Sorrie et al. 2006), particularly in fire-suppressed sites.

Fort Bragg is divided into over 1000 discrete landscape units, which are currently managed with prescribed fire. Each unit is burned approximately every three years. The longleaf pine savannas are highly flammable, but fire spread into the stream-head pocosins is often limited, likely due to high moisture content of fuels (Weakley and Schafale 1991), the fuel packing ratio, and microclimate (Rothermel 1972). Thus, within a burn unit, savanna species are exposed to fire more frequently than pocosin species.

We selected five savanna species (Pinus palustris, Quercus incana W. Bartram, Q. laevis Walter, Q. margarettae [Ashe] Small, and Q. marilandica Münchh.) and six pocosin species (Acer rubrum, Liquidambar styraciflua L., Liriodendron tulipifera, Nyssa biflora Walter, Pinus serotina, and Quercus nigra L.) that are relatively common at Fort Bragg (Table 1). Pinus palustris is best classified as fire resistant; in the grass stage, it is capable of resprouting from the root collar (Boyer 1990), but after bolting, it is capable of resprouting only from apical buds. All other study species are able to resprout basally from the root collar or from belowground after fire. Among oak (i.e., Quercus) species, savanna oaks generally have greater rhizome resprouting potential than forest oaks (Cavender-Bares et al. 2004a). All study species except P. palustris are also likely capable of epicormic resprouting; this occurs most often in A. rubrum, Liquidambar styraciflua, Liriodendron tulipifera, and P serotina (J. Schafer, North Carolina State University, Raleigh, USA, personal observation).

Table 1 Sample sizes, range of stem diameters, and best fit equations for the relationships between stem diameter (in cm) and bark thickness (in mm) and stem diameter and the proportion of inner bark; n.s. indicates a non-significant relationship. Equations were not forced through the origin. Statistics (i.e., R², F, and P) for the best fit equations are given in Figures 1 and 4. Species are in order of their location along the pocosin-to-savanna gradient.

We measured bark thickness and diameter at breast height (DBH; 1.4 m) of 20 to 34 individuals of each species. Individuals were haphazardly selected to cover the range of stem diameters present at Fort Bragg and were located in 34 different sites (i.e., burn units). We used a chisel and mallet to remove three sections of bark (1 cm² to 4 cm²) per stem and measured total bark thickness (at the thickest point) and outer bark (i.e., rhytidome) thickness (with calipers) where the bark was removed. Inner bark thickness was determined as the difference in thickness between total and outer bark. It was not possible to differentiate between outer and inner bark for small individuals (e.g., DBH approximately <3 cm, depending on the species), so only total bark thickness was recorded in these cases. We calculated the mean total bark thickness and the mean proportion of total bark thickness that is inner bark (i.e., the ratio of inner bark thickness to total bark thickness; hereafter referred to as the proportion of inner bark) for each individual.

For each species, we used the regression curve estimation function in SPSS version 19.0 (IBM Corporation, Armonk, New York, USA) to analyze the relationships between stem diameter and bark thickness and between stem diameter and the proportion of inner bark. We calculated relative bark thickness of each individual as the ratio of bark thickness to stem radius (multiplied by 100 %; Hoffmann et al. 2012). Lawes et al. (2013) suggested that relative bark thickness is a reliable measure of fire regime and that relative bark thickness should be determined for trees with smaller diameters because fire is likely to have a greater effect on the survival of these individuals. Thus, we analyzed relative bark thickness and the proportion of inner bark for three sets of individuals: <10 cm diameter, ≥10 cm diameter, and all individuals. Overall differences between savanna and pocosin species in relative bark thickness and the proportion of inner bark were analyzed with independent samples t-tests. Differences among species in relative bark thickness and the proportion of inner bark were analyzed with one-way ANOVAs with species as a fixed factor. Relative bark thickness was Log₁₀-transformed to meet the assumption of normality. Pairwise differences were determined with post-hoc Gabriel tests (Field 2009). Because stem resistance to fire scales as the square of bark thickness (Peterson and Ryan 1986, Lawes et al. 2011a), we also calculated bark thickness squared for a standardized stem size (5 cm diameter) of each species. The overall difference between savanna and pocosin species in bark thickness squared was analyzed with an independent samples t-test.

In 2011, we established 29 10 m wide belt transects that spanned the pocosin-to-savanna gradient. All trees >3 cm DBH in each transect were marked and identified to species. Due to differences in the steepness of the hydrological gradient, transect lengths varied from 24 m to 62 m across sites. To permit a relative measure of position along the gradient, each transect was divided into five equally sized plots and each plot was assigned a value based on location along the pocosin-to-savanna gradient: 0 (pocosin), 0.2, 0.4, 0.6, and 0.8 (savanna). Each individual tree was assigned the value of the plot in which it occurred, allowing us to calculate the mean location of each species along the gradient. Quercus incana did not occur in any transect because it tends to occur farther upland than the uppermost extent of the transects, so it was assigned a value of 1. We used the regression curve estimation function to test for a significant relationship across species between gradient location and mean relative bark thickness (of all individuals; different trees were used in the determination of gradient location and relative bark thickness). The DBH of each marked individual was measured yearly from 2011 to 2013. For the small individuals (<10 cm DBH) that are most susceptible to topkill during fire (Lawes et al. 2011b, VanderWeide and Hartnett 2011, Hoffmann et al. 2012, Lawes et al. 2013), we calculated the bark increment (i.e., increase in bark thickness per year) for each species by multiplying the mean ratio of bark thickness to stem radius (of individuals <10 cm DBH) by the mean increment in stem radius (mm yr⁻¹) of each species. We used the regression curve estimation function to test for a significant relationship across species between gradient location and mean bark increment. Growth data were not available for Q. incana and were available for only one individual of Q. margarettae, so these species were excluded from the analysis of bark increment.

Results

Bark thickness of all species increased with stem diameter (Figure 1). The relationship was curvilinear for most species (Table 1), indicating that relative bark thickness does not remain constant as trees grow. This had little influence on species comparisons, however, because at the species level, the rank order of relative bark thickness was nearly consistent between small (<10 cm DBH) and large (≥10 cm DBH) size classes (Spearman’s rho = 0.927, P < 0.001; Figure 2A and 2B).

Figure 1

Relationship between stem diameter (in cm) and bark thickness (in mm) for all savanna and pocosin species. Best fit regression lines (Table 1) are shown for each species. P < 0.001 for all species.

Figure 2

Mean (+SE) relative bark thickness of individuals <10 cm diameter (A), individuals ≥10 cm diameter (B), and all individuals (C), and the inner bark proportion of individuals <10 cm diameter (D), individuals ≥10 cm diameter (E), and all individuals (F) of each species. Different letters represent significant differences among species at α = 0.05. Species are ordered from left to right according to their location along the pocosin-to-savanna gradient.

Overall, savanna species had higher relative bark thickness than pocosin species (individuals <10 cm: t ₉ = 4.36, P = 0.002; individuals ≥10 cm: t ₉ = 2.29, P = 0.048; all individuals: t ₉ = 4.26, P = 0.002). Within genera, savanna oaks had higher relative bark thickness than the pocosin-dwelling Q. nigra, while P. palustris and P. serotina had similar relative bark thicknesses despite different habitat preferences (Figure 2). For a standardized stem size (5 cm diameter), savanna species had greater insulation (bark thickness squared) than pocosin species (t ₉ = 4.81, P = 0.001); the mean square of bark thickness was approximately 2.75 times higher for savanna than pocosin species (42.85 mm² vs. 15.49 mm² for savanna and pocosin species, respectively; Table 2).

Table 2 Sample sizes, range of stem diameters, mean diameter, and mean relative bark thickness (BT) for individuals <10 cm diameter at breast height (DBH). Because fire resistance scales as the square of bark thickness (BT²), we calculated the bark thickness and the square of bark thickness for a 5 cm DBH individual of each species. Species are in order of their location along the pocosin-to-savanna gradient.

Across all species, there was a positive logarithmic relationship between relative bark thickness (of all individuals) and mean position along the pocosin-to-savanna gradient (Figure 3A). Species that occur in the ecotone had intermediate relative bark thickness. The bark increment for individuals <10 cm DBH increased logarithmically with gradient location (Figure 3B); savanna species accumulated more bark per year than pocosin species. Across species, there was no relationship between growth rate and relative bark thickness (data not shown; R² = 0.015, F_1,7 = 0.108, P = 0.752) or between mean position along the pocosin-to-savanna gradient and growth rate (data not shown; R² = 0.033, F _1,7 = 0.240, P = 0.639).

Figure 3

Relationship between mean location along the pocosin-to-savanna gradient and mean relative bark thickness (A; for all individuals) and bark increment (B; for individuals <10 cm DBH). Each point represents one species. Q. margarettae and Q. incana are not included in panel B.

The proportion of inner bark decreased with stem diameter for four savanna species and three pocosin species (Figure 4; Table 1). Savanna and pocosin species had similar proportions of inner bark (individuals <10 cm: t _7.08 = −0.82, P = 0.440; individuals ≥10 cm: t ₉ = −0.63, P = 0.544; all individuals: t ₉ = −0.58, P = 0.575). At the species level, the rank order of inner bark proportion was similar between small and large size classes (Spearman’s rho = 0.927, P < 0.001; Figure 2D and 2E). The proportion of inner bark in P. palustris and P. serotina was, on average, approximately half that of other species.

Figure 4

Relationship between stem diameter and the proportion of inner bark for all savanna and pocosin species. Best fit regression lines (Table 1) are shown for significant relationships.

Discussion

Overall, savanna species had higher relative bark thickness than pocosin species, as hypothesized. Bark thickness is correlated with the time required for the cambium to reach lethal temperatures (Gashaw et al. 2002, van Mantgem and Schwartz 2003) and, thus, is a strong predictor of stem resistance to fire (Harmon 1984, Lawes et al. 2011a, Lawes et al. 2011b, VanderWeide and Hartnett 2011). Longleaf pine savannas burn more frequently than adjacent stream-head pocosins (Weakley and Schafale 1991, Frost 1995, Stambaugh et al. 2011), so higher relative bark thickness should be beneficial for surviving more frequent fire. In fact, the difference in relative bark thickness between savanna and pocosin species was larger—26 % vs. 15 %—when considering only the smallest individuals (i.e., <10 cm diameter), which are most susceptible to topkill during fire (Lawes et al. 2011b, VanderWeide and Hartnett 2011, Hoffmann et al. 2012, Lawes et al. 2013), rather than all individuals. Given this difference, and the fact that stem resistance to fire scales as the square of bark thickness (Peterson and Ryan 1986, Lawes et al. 2011a), individuals of pocosin species must be approximately 1.6 times larger, on average, than individuals of savanna species to accumulate enough bark to reach the same level of fire resistance. Our results support other studies that have reported higher relative bark thickness in savanna trees compared to forest trees (Hoffmann et al. 2003, Hoffmann et al. 2009, Lawes et al. 2011a).

Using a more continuous characterization of species habitat, we found that relative bark thickness was related to the mean location of species along the pocosin-to-savanna gradient, as hypothesized. The direction of causation, however, is not clear because the distribution of our study species may be strongly influenced by water and nutrient availability (Donovan et al. 2000), as well as by the direct effects of fire. The decline in soil moisture along the pocosin-to-savanna gradient (J. Schafer, unpublished data) selects for drought tolerance, while the increase in fire frequency selects for thick bark. Thus, due to past environmental filtering, there is a correlation between habitat preference and bark thickness. The higher fire frequency in savannas could select for thicker bark in drought-tolerant species, while drought-tolerant species with inherently thicker bark could be more likely to survive frequent fires due to the protective effects of thick bark (Harmon 1984, Lawes et al. 2011b, Catry et al. 2012, Hoffmann et al. 2012). In addition, among small individuals of oak species, those that prefer dry habitats have greater investment in bark thickness than those that prefer wet habitats (Schwilk et al. 2013). Regardless of the mechanisms that control species distributions, greater bark accumulation by ecotone and savanna species indicates that individuals of these species should become fire resistant more quickly than pocosin species.

Contrary to our hypothesis, we found no difference in the proportion of inner bark between savanna and pocosin species overall. We expected savanna species to have a greater relative investment in inner bark because heat is transferred through outer bark faster than through inner bark (van Mantgem and Schwartz 2003) and a higher proportion of inner bark may reduce the susceptibility of topkill. This expectation, however, overlooks tradeoffs and associations between bark traits and functions (Rosell et al. 2013, Poorter et al. 2014) as well as potential differences in the cost of producing and maintaining these two bark components. Specifically, inner bark is a living tissue that is likely to incur substantial maintenance costs, so ultimately there may be little or no net benefit of greater relative investment in inner bark. Furthermore, relative allocation to inner versus outer bark may be limited by phylogenetic constraints because outer bark thickness is a somewhat conserved trait in oaks (Cavender-Bares et al. 2004b). This is supported in our study: the oak species all had high inner bark proportions, while the pines had low inner bark proportions, regardless of habitat.

Despite the lack of a difference in the proportion of inner bark between savanna and pocosin species, the greater insulation provided by inner bark (van Mantgem and Schwartz 2003, Scholz et al. 2007) may be important when considering small individuals, which are more susceptible to topkill (Hoffmann et al. 2012). Within species, there was a widespread trend for small stems to have a higher proportion of inner bark than larger stems. Regardless of whether this trend is the result of size-specific differences in natural selection or simply allometry, the greater relative investment in inner bark should benefit stems when they are most vulnerable to fire.

Stream-head pocosin species tended to have relatively thin bark, which suggests that, for stems of a given diameter, pocosin species will suffer greater topkill during fire than savanna species (Lawes et al. 2011a, Hoffmann et al. 2012). Differences in the proportion of inner bark among species, however, may mediate the influence of bark thickness on stem survival. Nonetheless, our results suggest that species’ distributions along the pocosin-to-savanna gradient may be influenced by relative bark thickness. Pocosin species accumulate bark more slowly than savanna species, but a fire return interval of 7 to 50 years (Frost 1995) may allow individuals of pocosin species to grow large enough, increasing in both height and diameter (Lawes et al. 2011b), to accumulate thick enough bark to survive fire. A longer fire-free interval, however, may lead to higher fuel loads and more intense fires, which can cause greater stem damage (Bova and Dickinson 2005) and mortality (Glitzenstein et al. 1995, Adie et al. 2011) regardless of stem diameter and bark thickness. On the other hand, changes in species composition over a long fire-free interval (Gilliam and Platt 1999) may reduce fire intensity due to changes in vegetation structure and fuel flammability (van Wilgen et al. 1990, Behm et al. 2004). Regardless, fire suppression may contribute to hardwood establishment in savannas (Waldrop et al. 1992, Gilliam and Platt 1999, Van Lear et al. 2005) by allowing pocosin species to reach a diameter (Chatziefstratiou et al. 2013) or accumulate bark thickness sufficient for surviving future fires.