Next Article in Journal
Disentangling the Roles of Topography, Patch, and Land Use on Conservation Trait Status of Specialist Birds in Marginal Forest Land Use Types
Next Article in Special Issue
An ERF Transcription Factor Gene from Malus baccata (L.) Borkh, MbERF11, Affects Cold and Salt Stress Tolerance in Arabidopsis
Previous Article in Journal
A Framework for Characterizing and Regulating Ecosystem Services in a Management Planning Context
Previous Article in Special Issue
Tree Growth and Water-Use Efficiency Do Not React in the Short Term to Artificially Increased Nitrogen Deposition
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 11
Issue 1
10.3390/f11010105
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Role of Climate Niche, Geofloristic History, Habitat Preference, and Allometry on Wood Density within a California Plant Community

by
Rebecca A. Nelson
1,*,
Emily J. Francis
2,
Joseph A. Berry
3,
William K. Cornwell
4 and
Leander D. L. Anderegg
3,5
1
Department of Biology, School of Humanities and Sciences, Stanford University, Stanford, CA 94305, USA
2
Department of Integrative Biology, College of Natural Sciences, University of Texas at Austin, Austin, TX 78712, USA
3
Department of Global Ecology, Carnegie Institution for Science, Stanford, CA 94305, USA
4
School of Biological, Earth and Environmental Sciences, University of New South Wales, Sydney, NSW 2052, Australia
5
Department of Integrative Biology, College of Letters & Science, University of California, Berkeley, CA 94720, USA
*
Author to whom correspondence should be addressed.
Forests 2020, 11(1), 105; https://0-doi-org.brum.beds.ac.uk/10.3390/f11010105
Submission received: 1 December 2019 / Revised: 21 December 2019 / Accepted: 7 January 2020 / Published: 14 January 2020
(This article belongs to the Special Issue Relationship between Forest Ecophysiology and Environment)
Browse Figures
Review Reports Versions Notes

Abstract

:
Research Highlights: To better understand within-community variation in wood density, our study demonstrated that a more nuanced approach is required beyond the climate–wood density correlations used in global analyses. Background and Objectives: Global meta-analyses have shown higher wood density is associated with higher temperatures and lower rainfall, while site-specific studies have explained variation in wood density with structural constraints and allometry. On a regional scale, uncertainty exists as to what extent climate and structural demands explain patterns in wood density. We explored the role of species climate niche, geofloristic history, habitat specialization, and allometry on wood density variation within a California forest/chaparral community. Materials and Methods: We collected data on species wood density, climate niche, geofloristic history, and riparian habitat specialization for 20 species of trees and shrubs in a California forest. Results: We found a negative relationship between wood density and basal diameter to height ratio for riparian species and no relationship for non-riparian species. In contrast to previous studies, we found that climate signals had weak relationships with wood density, except for a positive relationship between wood density and the dryness of a species’ wet range edge (species with drier wet range margins have higher wood density). Wood density, however, did not correlate with the aridity of species’ dry range margins. Geofloristic history had no direct effect on wood density or climate niche for modern California plant communities. Conclusions: Within a California plant community, allometry influences wood density for riparian specialists, but non-riparian plants are ‘overbuilt’ such that wood density is not related to canopy structure. Meanwhile, the relationship of wood density to species’ aridity niches challenges our classic assumptions about the adaptive significance of high wood density as a drought tolerance trait.

1. Introduction

Wood density (the mass per unit volume of wood) is considered an essential functional trait associated with mechanical support, carbon and nutrient storage, drought tolerance, water transport, and pathogen defense [1,2,3]. Variation in wood functional traits often reflect ecological tradeoffs that are influenced by allometric, biogeographic, and phylogenetic factors [1]. Understanding wood density patterns can help elucidate spatial trends in aboveground biomass [4]. Moreover, wood density is a useful predictor of growth and mortality rates in the tropics [5] and is related to growth and performance in Mediterranean climates [6]. Finally, models of carbon stocks often rely on wood density, making estimates of wood density operationally as well as ecologically important [7,8].
Wood density often relates to mechanical stability [9,10,11] with crown architecture, diameter, and height correlating with wood density at small scales [8,12]. Tree height is a key trait for understanding patterns in wood density and allometry in Mediterranean communities [6]. Allometric equations that relate wood density to height and diameter have explained site-specific patterns across tropical ecosystems [8,13]. The relationship of wood density to diameter–height allometry at an individual site level has varied across studies. Whereas some studies have found that tree height at a given diameter is unrelated to wood density [14,15], Iida et al. (2012) found that, at a constant tree height, higher wood density is associated with smaller tree diameters [16], and wood density is positively related to diameter-corrected tree height in rainforests [17]. Simultaneously, in a California mixed evergreen forest, wood density correlated strongly with vessel traits and microsite soil moisture [18,19], indicating that moisture availability can be an environmental filter for wood density at very small scales. The variation in the wood density–allometry relationship across different climate regimes and site conditions may be partially explained by the availability of water at particular microsites, with some individuals and species at or close to drainages being able to access consistent water throughout the year, while other species experience strong seasonal fluctuations in water availability. In other words, riparian habitat specialization could influence the degree to which structural constraints and climate niche relate to wood density.
Global and regional meta-analyses suggest that higher wood density correlates with lower precipitation and higher mean annual temperature [6,11,20]. Across a latitudinal gradient in Central America, patterns in wood density more strongly correlated to precipitation and aridity than temperature [21]. Wood density is positively associated with aridity in New Caledonia, and a global meta-analysis [22,23] suggests that high wood density is adaptive in dry environments. However, precipitation and temperature explain variation in tree height–diameter relationships across North America as well, possibly influencing wood density indirectly through structural constraints [24]. The findings that rainfall amount can modify wood density–allometry relationships at a regional level in water-stressed climates suggests at least some direct effect of climate on wood density [25], and denser wood is correlated with lower drought mortality [26]. Variations in shade tolerance and altitude add further nuance to relationships between wood density and rainfall patterns in tropical and Mediterranean climates [12,27].
While many studies have examined the relationship between climate and wood density across a broad geographic scale, and a number of tropical studies have highlighted the importance of structural or allometric constraints on wood density at a site level, few studies have assessed both structural and climatic constraints on wood density in temperate communities. Moreover, phylogenetic signal, biogeography, and geofloristic history (plant assemblages recorded in the fossil record) can strongly structure the composition and functional traits of temperate plant communities [11,22,27,28]. As a result, community composition that is the contingent outcome of past climates could drive functional trait patterns that are not mechanistically linked to the contemporary climate [28]. In water-stressed ecosystems, uncertainty exists as to whether mechanical demands, habitat specialization, geofloristic history, or macroclimatic niche best explain variation in wood density.
The California floristic province, with its high plant diversity and endemism, is vulnerable to increasing drought frequency and intensity exacerbated by anthropogenic climate change. Microsite variation in soil moisture influences the wood density and tree height of the local community across microsites in northern California [18,19]. This suggests environmental filtering for higher wood density species in drier sites. The current species composition of California mixed evergreen forests resulted from the confluence of several geofloras (plant community assemblages identified in the fossil record), especially the temperate Arcto-Tertiary geoflora and the subtropical Madro-Tertiary geoflora [29,30]. The Madro-Tertiary and Arcto-Tertiary geofloras, both originating in the Tertiary Period, are distinguished by their biogeographic origin and differing community compositions. The Madro-Tertiary geoflora, consisting mostly of subtropical, sclerophyll species, is associated with the development of a Mediterranean climate in California and is the ancestor of many extant chaparral specialists [29,30]. The Arcto-Tertiary geoflora, compromised of mixed deciduous forest species, advanced into California from more northerly climates [29,30].
We investigated whether structure, climatic niche, geofloristic history, and habitat specialization related to wood density variation among members of a regional woody plant community, sampling 20 species from a northern California plant community across multiple sites. We tested four hypotheses: (1) We hypothesize that the structural demand of supporting a plant canopy will constrain wood density values such that species with tall, thin stems must have higher wood density to support their canopies. We predict that ‘taper’ (the ratio of plant basal diameter to plant height) will be negatively related to wood density. (2) We predict that, within a community, more arid adapted species will have higher wood density. Thus, we predict that wood density and metrics of water availability will be negatively related for the dry and wet edges of species’ aridity niche, with stronger relationships at the dry edge. (3) We hypothesize that species with a Madro-Tertiary geofloristic history would have stronger wood density–climate relationships than species with an Arcto-Tertiary geofloristic history because the Madro-Tertiary geoflora arose in conjunction with California’s dry Mediterranean climate. (4) Finally, we hypothesize that riparian habitat specialization may greatly influence the relationship between wood density and structure as well as wood density and climate niche.

2. Materials and Methods

2.1. Study Site

We sampled naturally occurring vegetation at three sites near Stanford, CA: Jasper Ridge Biological Preserve (37.40311, −122.24428), a 481-ha preserve in the eastern foothills of the Santa Cruz Mountains; the Stanford Quarry (37.39999, −122.16154); and the Stanford main campus (37.42875, −122.17919) (Figure 1). The climate is Mediterranean with a mean annual temperature of 14 °C, mean annual precipitation of 570 mm, and most rainfall occurring from November to April. The vegetation communities present at all three sites included chaparral and broadleaf, mixed evergreen forest. All three sites have loam to gravely clay loam soils derived from similar mixed sedimentary and metamorphic residuum and alluvium [31].

2.2. Species Selection

We selected woody species native to northern California. In the spring and summer of 2017, we sampled five individuals from 13 native species growing on the Stanford campus or the Stanford Quarry. We sampled reasonably abundant species (five or more individuals in a ~200 m search area) and focused sampling on the largest individuals of each species present based on a visual survey. We included comparable data of wood density and plant allometry from Cornwell and Ackerly (2009) collected at the Jasper Ridge Biological Preserve, making for a total of 20 species of California trees and shrubs [19]. A range of habitats including riparian forests and chaparral/oak woodlands were sampled at all three sites to capture the dominant woody species across multiple microsites. We only sampled species that grew naturally at our sites. All three sites have similar climate and native plant communities. See Table S1 for a full list of species sampled and which datasets these species occurred in. For data on the heights, diameters, wood density values, and other traits measured of the study species, see Tables S2 and S3.

2.3. Trait Measurement

We obtained three branches per individual for at least five individuals sampled per species to use for measuring branch wood density via a volumetric method. We randomly selected mid- to upper-canopy, sun-exposed branches from each individual using hand pruners or a 4.5 m pole pruner. Branches of ~1 cm diameter and ~3 cm length were peeled of bark, their wet volume measured via mass displacement, and their dry weight measured following 72+ h drying at 65 °C. We measured individual height with an inclinometer. We obtained basal diameter for larger trees as diameter at breast height. We then calculated ‘taper’ as the basal diameter to height ratio, back calculating the basal diameter from diameter at breast height measurements assuming a conic stem where necessary. We averaged branch wood density per individual, and then averaged wood density and taper of individuals per species. We used branch wood density in order make our results comparable to wood density values obtained by Cornwell and Ackerly (2009) [19], who used data from three-year old branches from five individuals per species (see Preston et al. (2006) [18] for details of wood density measurements).

2.4. Habitat Preferences

We used information in the National Wildlife Federation and Audubon field guides as well as Calscape.org to determine whether species in our study were riparian specialists [32,33,34]. If the literature described a species as preferring rivers, streambanks, or riparian areas, we recorded these species as riparian specialists. Our dataset contained 11 riparian species and 9 non-riparian species with a mix of riparian and non-riparian species sampled at each of the three sites (Figure S1). Table S2 provides information on whether species were riparian and the sample size per species.

2.5. Paleohistory

We conducted a literature search in order to determine whether each species had ancestors in the Madro-Tertiary or Arcto-Tertiary geoflora. We assigned a geoflora category based off the most taxonomically specific level possible; for most species, we obtained geofloristic information at the genus level. For chaparral species, we placed them in the Madro-Tertiary group based on literature stating that the current chaparral community evolved from the Madro-Tertiary geoflora.

2.6. Climate Data

We determined climatic parameters describing each species climate niche using the WorldClim database and geo-referenced specimens from North America in the Global Biodiversity Information Facility (GBIF) database [35]. We used the full climate range for each species. We obtained 1960–1990 long-term averages of the following annual climatic parameters at 2.5 arc-minutes from the WorldClim 1.4 database [36]: Minimum temperature (Tmin), maximum temperature (Tmax), mean annual temperature (MAT), mean annual precipitation (MAP), and precipitation seasonality (coefficient of variation of precipitation across months, Pseason).We also obtained aridity data from the Global Aridity and PET Database [37], which is derived from WorldClim input data, including potential evapotranspiration (PET) and climate moisture deficit (CMD). CMD was calculated as the amount of water by which PET exceeds precipitation. Climate data were extracted for all GBIF specimen locations for each species. The extent of each species’ climate niche was calculated based on the 5th and 95th percentile for each of the above climate variables.

2.7. Data Analysis

We averaged all wood density and taper values per species (taper did not show variation with plant height within-species) and constructed a series of linear models for the data, using R version 3.3.3. First, we constructed linear models relating wood density to taper, habitat preference, the interaction between the two, and the interaction without the habitat preference main effect. We selected the most parsimonious model based on the second-order Aikike information criterion (AICc, or AIC for small sample sizes) using the ‘stats’ package and the ‘MuMIn’ [38]. We also tested for univariate relationships between wood density and climate niche parameters related to water availability (CMD, MAP, PET) for the driest edge and wettest edge of each species (based on the 5th and 95th percentile of their species range), and then selected the best single climate variable by built linear models for niche parameters that were initially found to be significant using AICc. Finally, we combined the best structural model and the best climate model to build a single full model and all potential nested models, and then used AICc to select the most parsimonious of these models. The equation for our final model was as follows: WD ~5th% CMD + Riparian Specialization × Taper, where WD is wood density and CMD is climate moisture deficit. Finally, we verified that this inference was not purely driven by phylogenetic non-independence by refitting the final model using a phylogenetic generalized least squares (PGLS) model. We used the Smith and Brown 2018 [39] phylogeny (which contained all species except Baccharis pilularis, which was replaced with a congener to estimate branch lengths) and a Brownian motion covariance structure [40] using the ape [41], geiger [42,43], and nlme [44] R packages. See Figure S2 for the results of the phylogenetic analysis.

3. Results

A higher taper (basal diameter to height ratio) refers to woody plants that are wider in basal diameter at a given height (Figure 2). For riparian habitat specialist species, branch wood density decreased with increasing taper (Figure 2). For non-riparian species, no correlation existed between branch wood density and taper (Figure 2). The taper model with the lowest AICc had a main effect of taper and a taper x riparian interaction (but no riparian main effect), though this model was only a marginal improvement over a model with only a riparian specialization main effect (ΔAICc = 1.1, likelihood ratio test p = 0.06017) (Table 1).
Few correlations existed between wood density and species driest range edge, on the contrary to the existing relationships between wood density and site water availability in global meta-analyses (Figure 3). However, positive relationships between wood density and the dryness of a species’ wettest range boundary were found for all metrics of moisture availability (5th percentile CMD, PET, and MAT; 95th percentile MAP), with 5th percentile CMD being the best climate predictor of wood density (Table 2). Riparian specialists had generally lower wood density than non-specialists for the same wet edge aridity, but showed a similar positive relationship between wood density and wet edge aridity (Figure 3). Adding a main effect of habitat specialization improved the best climate model (ΔAICc = 6.86, Table 2), but adding an interaction between climate and habitat specialization did not improve the model. These results were qualitatively similar if climate niches were constructed from only California occurrence records in the California Consortium of Herbaria [45].
Although geoflora suggested small differences between the climate niches of Madro-Tertiary and Arcto-Tertiary species, especially for the minimum temperature range (Figure 4 and Figure 5), geoflora did not directly affect wood density ( ∆ AICc = 0.62918, likelihood ratio test p = 0.2794) or strongly differentiate the climate niches of modern California plant communities (Figure 4 and Figure 5).
Combining all evidence, aridity of the wet range edge and, for riparian specialists, taper explained almost three quarters of the variation in wood density across this California plant community. The most parsimonious full model overall included a main effect of 5th percentile CMD and a riparian specialization x taper interaction (Table 3, Figure 6) and had an R2 of 0.74. On the basis of a decomposition of variance analysis (assessing the R2 of nested climate and structural models), structure and habitat specialization explained more of the variance than climate, with species 5th percentile CMD explaining 32% of the total variance and the taper × riparian specialization interaction explaining 52% of the variance, with 10% of the variance shared between them. All model effects remained significant and qualitatively similar when controlling for phylogenetic non-independence using a PGLS model [39], though visual inspection of the congeneric species in the dataset suggested that the lack of relationship between taper and wood density for non-riparian species may be driven by the high wood density and high taper values of the four species of oaks in the dataset (Figure S2).

4. Discussion

In a Mediterranean woody plant community, we found that canopy allometry was only related to wood density in riparian specialists, that geofloristic provenance had little influence on wood density, and that species climate niche had unexpected relationships with wood density. These results contrast with both global analyses of wood density and site-specific investigations of wood density variation within individual communities. We discuss the implications of our results for the context dependence of the climatic and structural controls on wood density, and the potential for disconnects between small-scale environmental filtering at microsites and large-scale environmental filtering of a regional species pool.
We found that allometry only relates to wood density in riparian contexts. In the California mixed evergreen forest, the relationship between wood density and allometry, however, does not appear as consistent as the wood density–allometry structural equations proposed by Chave et al. (2005) for tropical forests [13]. This lack of consistency may be explained by the differing functional strategies of riparian and non-riparian species. We found that riparian specialists had a negative relationship between wood density and taper consistent with Chave et al. (2005) [13], indicating that species with high wood density are taller at a given diameter than species with low wood density. This result is also consistent with Iida et al. (2012), who found that species with higher wood density have smaller diameters at a given height in Malaysia [16]. However, these findings contrast with those of Aiba & Nakashizuka (2009) as well as King, Davies, Tan, & Noor (2006), who found no relationship between wood density and height for a given diameter in tropical forests [14,15]. Predictions from engineering theory suggest that high wood density does not necessarily increase the mechanical stability of taller trees at a given diameter, which further contrasts with our results [46]. This negative relationship between wood density and taper is consistent with previous work at one of our field sites that showed strong relationships between wood density, allometry, and microsite soil moisture [18,19]. This allometric relationship with wood density for riparian plants may be similar to allometry–wood density relationships in the wet tropics, because both ecosystems are not water-limited.
Differences in allometry–wood density relationships between riparian and non-riparian species may reflect evolutionary processes. Chave et al. (2009) suggested that wood density is related to water transport efficiency [1]. Stronger selection on wood density–allometry relationships may, therefore, occur in wet climates where efficient transportation of water would be advantageous. Non-riparian species, however, lack this wood density–allometry relationship, suggesting that, in water-limited ecosystems, these species are ‘overbuilt’ for other reasons [47,48]. Riparian forests tend to be more light-limited than immediately adjacent non-riparian forests owing to their closed, tall canopies, so selection may occur in riparian habitats to favor taller, competitive trees, giving wood density more of a structural role in riparian than non-riparian ecosystems, but not in systems with less light competition. The high taper, non-riparian species were all oaks, including species with deep roots that can penetrate groundwater [49,50]. A large dataset covering more sites would be needed to determine if the oaks have a phylogenetic signal.
The California floristic province formed from the confluence of a diverse range of geofloras [29], but this history did not explain the variation in wood density. The Madro-Tertiary geoflora has been associated with the modern chaparral species in California and associated with the development of a Mediterranean climate in California [30]. Species with Madro-Tertiary ancestry did not significantly differ in climate niche than species with temperate Arcto-Tertiary history, suggesting that paleoclimatic influences do not clearly shape modern climate niches in this community. There were neither strong influences of paleoclimatic history on wood density directly nor on the effects of climate or taper on wood density. A more comprehensive study of the traits and climate niches of Madro-Tertiary versus Arcto-Tertiary geofloras could reveal more subtle legacies of their biogeographic origins, but we found little preliminary evidence of geofloristic legacies in this community.
Surprisingly, the wet edge of climate moisture deficit best predicted wood density when compared with other climate variables. This contrasts with meta-analyses conducted at regional and global scales where mean annual temperature and mean annual precipitation of the sampling site explained patterns in wood density across many species [11,20,21] and higher wood density was associated with decreased drought mortality [26]. The positive relationship between wood density and climate moisture deficit is consistent with previous literature, but it is surprising that this trend does not appear for the dry edge. For a Mediterranean plant community, Camarero (2019) found that the relationship between April precipitation and tree growth was the strongest in the species with the highest wood density and that the wet season particularly influenced radial growth rates [6]. In contrast, for a Mediterranean climate in South Africa, xylem density positively correlated with seasonal water stress, suggesting that the dry edge of climate niche was more of an influence [10]. For tropical and subtropical forests, wood density positively related to aridity across sites [22,25,27]. Meanwhile, in a tropical forest, deciduous versus evergreen leaf habit mediated the degree of variation in wood density in response to rainfall [51]. A water-stressed Mediterranean climate may select for lower wood density in wet environments rather than higher wood density in dry environments. Structural and competitive constraints may have more of an influence on wood density in wetter ecosystems, where it is more advantageous to grow tall and fast. Selection for extremely efficient water transport may occur in wetter climates where increased conductivity has an advantageous impact on plant fitness.
In order to better understand landscape-level variation in wood density, a more nuanced approach is required beyond the species-mean climate–wood density correlations used in global analyses. Moreover, examining relationships between climate and wood density alone may overlook the influence of structural constraints and habitat specialization on wood density. Structure influenced wood density only for riparian species, suggesting that the relationship between wood density and allometry is highly context-specific, and that differences in habitat specialization across a landscape should be accounted for. A disconnect often exists between the climate and wood density correlations observed at regional levels and the microsite influences on wood density at localized scales [6,11,18,19,20]. Our findings caution against the overgeneralization of how structural constraints influence wood density across ecosystems. Further research is needed to explore the complexity of factors, such as climate, allometry, and habitat specialization, that shape wood density variation at a community level. Understanding how these factors drive functional tradeoffs in wood traits can provide more nuanced insights into the biological mechanisms by which these tradeoffs occur. Elucidating the extent to which climate and local habitat explain wood functional traits is relevant in the context of increasing water stress resulting from anthropogenic climate change [26,52,53]. Furthermore, fire greatly shapes the structure and composition of California ecosystems and relates to wood density in tropical forests, yet the extent to which wood density relates to species’ responses to fire remains unknown [54]. Finally, an exploration of within-species patterns of wood density and allometry between individuals might elucidate some of the fundamental structural constraints driving the patterns we found in riparian specialists in this comparative study.

5. Conclusions

In summary, we found that a negative relationship occurred between the basal diameter to height ratio and wood density only for riparian habitat specialists. Geofloristic history had little influence on relationships between wood density and climate in our woody plant community. The wet edge of climate moisture deficit further explained variation in wood density. These results suggest a more nuanced approach is required in order disentangle the drivers of wood density variation at a community level.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1999-4907/11/1/105/s1, Tables S1–S3 and Figures S1 and S2.

Author Contributions

Conceptualization, R.A.N., E.J.F., J.A.B., W.K.C., and L.D.L.A.; methodology and sampling design and data collection, R.A.N. and L.D.L.A.; analysis, R.A.N. and L.D.L.A.; data curation, R.A.N., W.K.C., and L.D.L.A.; writing—Original draft preparation, R.A.N., E.J.F., and L.D.L.A.; writing—Review and editing, R.A.N., E.J.F., J.A.B., W.K.C., and L.D.L.A.; visualization, R.A.N. and L.D.L.A.; supervision, L.D.L.A.; project administration, R.A.N. and L.D.L.A.; funding acquisition, R.A.N. and L.D.L.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a Stanford Vice Provost for Undergraduate Education small grant. L.D.L.A. received funding from the National Science Foundation (DBI-1711243) and a Climate and Global Change postdoctoral fellowship from the National Oceanic and Atmospheric Administration. Any opinions, findings, and conclusions or recommendations expressed in this material are those of the author(s) and do not necessarily reflect the views of the National Science Foundation.

Acknowledgments

Thank you to Rodolfo Dirzo for providing initial advice on this project.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Chave, J.; Coomes, D.; Jansen, S.; Lewis, S.L.; Swenson, N.G.; Zanne, A.E. Towards a worldwide wood economics spectrum. Ecol. Lett. 2009, 12, 351–366. [Google Scholar] [CrossRef]
  2. Ogle, K.; Pathikonda, S.; Sartor, K.; Lichstein, J.W.; Osnas, J.L.D.; Pacala, S.W. A model-based meta-analysis for estimating species-specific wood density and identifying potential sources of variation. J. Ecol. 2014, 102, 194–208. [Google Scholar] [CrossRef] [Green Version]
  3. Steffenrem, A.; Solheim, H.; Skrøppa, T. Genetic parameters for wood quality traits and resistance to the pathogens Heterobasidion parviporum and Endoconidiophora polonica in a Norway spruce breeding population. Eur. J. For. Res. 2016, 135, 815–825. [Google Scholar] [CrossRef] [Green Version]
  4. Baker, T.R.; Phillips, O.L.; Malhi, Y.; Almeida, S.; Arroyo, L.; Di Fiore, A.; Erwin, T.; Killeen, T.J.; Laurance, S.G.; Laurance, W.F.; et al. Variation in wood density determines spatial patterns in Amazonian forest biomass. Glob. Chang. Biol. 2004, 10, 545–562. [Google Scholar] [CrossRef]
  5. Wright, S.J.; Kitajima, K.; Kraft, N.J.B.; Reich, P.B.; Wright, I.J.; Bunker, D.E.; Condit, R.; Dalling, J.W.; Davies, S.J.; Díaz, S.; et al. Functional traits and the growth–mortality trade-off in tropical trees. Ecology 2010, 91, 3664–3674. [Google Scholar] [CrossRef]
  6. Camarero, J.J. Linking functional traits and climate-growth relationships in Mediterranean species through wood density. IAWA J. 2019, 40, 215–240. [Google Scholar] [CrossRef] [Green Version]
  7. Sarmiento, C.; Patiño, S.; Timothy Paine, C.E.; Beauchêne, J.; Thibaut, A.; Baraloto, C. Within-individual variation of trunk and branch xylem density in tropical trees. Am. J. Bot. 2011, 98, 140–149. [Google Scholar] [CrossRef]
  8. Fayolle, A.; Doucet, J.L.; Gillet, J.F.; Bourland, N.; Lejeune, P. Tree allometry in Central Africa: Testing the validity of pantropical multi-species allometric equations for estimating biomass and carbon stocks. For. Ecol. Manag. 2013, 305, 29–37. [Google Scholar] [CrossRef]
  9. Pratt, R.B.; Jacobsen, A.L.; Ewers, F.W.; Davis, S.D. Relationships among xylem transport, biomechanics and storage in stems and roots of nine Rhamnaceae species of the California chaparral. New Phytol. 2007, 174, 787–798. [Google Scholar] [CrossRef]
  10. Jacobsen, A.L.; Agenbag, L.; Esler, K.J.; Pratt, R.B.; Ewers, F.W.; Davis, S.D. Xylem density, biomechanics and anatomical traits correlate with water stress in 17 evergreen shrub species of the Mediterranean-type climate region of South Africa. J. Ecol. 2007, 95, 171–183. [Google Scholar] [CrossRef]
  11. Zhang, S.B.; Slik, J.W.F.; Zhang, J.L.; Cao, K.F. Spatial patterns of wood traits in China are controlled by phylogeny and the environment. Glob. Ecol. Biogeogr. 2011, 20, 241–250. [Google Scholar] [CrossRef]
  12. Poorter, L.; Lianes, E.; Moreno-de las Heras, M.; Zavala, M.A. Architecture of Iberian canopy tree species in relation to wood density, shade tolerance and climate. Plant Ecol. 2012, 213, 707–722. [Google Scholar] [CrossRef] [Green Version]
  13. Chave, J.; Andalo, C.; Brown, S.; Cairns, M.A.; Chambers, J.Q.; Eamus, D.; Fölster, H.; Fromard, F.; Higuchi, N.; Kira, T.; et al. Tree allometry and improved estimation of carbon stocks and balance in tropical forests. Oecologia 2005, 145, 87–99. [Google Scholar] [CrossRef]
  14. King, D.A.; Davies, S.J.; Tan, S.; Noor, N.S.M. The role of wood density and stem support costs in the growth and mortality of tropical trees. J. Ecol. 2006, 94, 670–680. [Google Scholar] [CrossRef]
  15. Aiba, M.; Nakashizuka, T. Architectural differences associated with adult stature and wood density in 30 temperate tree species. Funct. Ecol. 2009, 23, 265–273. [Google Scholar] [CrossRef]
  16. Iida, Y.; Poorter, L.; Sterck, F.J.; Kassim, A.R.; Kubo, T.; Potts, M.D.; Kohyama, T.S. Wood density explains architectural differentiation across 145 co-occurring tropical tree species. Funct. Ecol. 2012, 26, 274–282. [Google Scholar] [CrossRef]
  17. Read, J.; Evans, R.; Sanson, G.D.; Kerr, S.; Jaffré, T. Wood properties and trunk allometry of co-occurring rainforest canopy trees in a cyclone-prone environment. Am. J. Bot. 2011, 98, 1762–1772. [Google Scholar] [CrossRef]
  18. Preston, K.A.; Cornwell, W.K.; DeNoyer, J.L. Wood density and vessel traits as distinct correlates of ecological strategy in 51 California coast range angiosperms. New Phytol. 2006, 170, 807–818. [Google Scholar] [CrossRef]
  19. Cornwell, W.K.; Ackerly, D.D. Community assembly and shifts in plant trait distributions across an environmental gradient in coastal California. Ecol. Monogr. 2009, 79, 109–126. [Google Scholar] [CrossRef] [Green Version]
  20. Swenson, N.G.; Enquist, B.J. Ecological and evolutionary determinants of a key plant functional trait: Wood density and its community-wide variation across latitude and elevation. Am. J. Bot. 2007, 94, 451–459. [Google Scholar] [CrossRef]
  21. Martínez-Cabrera, H.I.; Jones, C.S.; Espino, S.; Jochen Schenk, H. Wood anatomy and wood density in shrubs: Responses to varying aridity along transcontinental transects. Am. J. Bot. 2009, 96, 1388–1398. [Google Scholar] [CrossRef] [PubMed]
  22. Ibanez, T.; Chave, J.; Barrabé, L.; Elodie, B.; Boutreux, T.; Trueba, S.; Vandrot, H.; Birnbaum, P. Community variation in wood density along a bioclimatic gradient on a hyper-diverse tropical island. J. Veg. Sci. 2017, 28, 19–33. [Google Scholar] [CrossRef]
  23. Nabais, C.; Hansen, J.K.; David-Schwartz, R.; Klisz, M.; López, R.; Rozenberg, P. The effect of climate on wood density: What provenance trials tell us? For. Ecol. Manag. 2018, 408, 148–156. [Google Scholar] [CrossRef]
  24. Hulshof, C.M.; Swenson, N.G.; Weiser, M.D. Tree height-diameter allometry across the United States. Ecol. Evol. 2015, 5, 1193–1204. [Google Scholar] [CrossRef] [PubMed]
  25. Onoda, Y.; Richards, A.E.; Westoby, M. The relationship between stem biomechanics and wood density is modified by rainfall in 32 Australian woody plant species. New Phytol. 2010, 185, 493–501. [Google Scholar] [CrossRef]
  26. Greenwood, S.; Ruiz-Benito, P.; Martínez-Vilalta, J.; Lloret, F.; Kitzberger, T.; Allen, C.D.; Fensham, R.; Laughlin, D.C.; Kattge, J.; Bönisch, G.; et al. Tree mortality across biomes is promoted by drought intensity, lower wood density and higher specific leaf area. Ecol. Lett. 2017, 20, 539–553. [Google Scholar] [CrossRef]
  27. Chave, J.; Muller-Landau, H.C.; Baker, T.R.; Easdale, T.A.; Hans Steege, T.E.R.; Webb, C.O. Regional and phylogenetic variation of wood density across 2456 neotropical tree species. Ecol. Appl. 2006, 16, 2356–2367. [Google Scholar] [CrossRef] [Green Version]
  28. Ackerly, D.D. Adaptation, niche conservatism, and convergence: Comparative studies of leaf evolution in the California chaparral. Am. Nat. 2004, 163, 654–671. [Google Scholar] [CrossRef]
  29. Axelrod, D.I. Evolution of the madro-tertiary geoflora. Bot. Rev. 1958, 24, 433–509. [Google Scholar] [CrossRef]
  30. Axelrod, D.I. Axelrod. Bull. Geol. Soc. Am. 1957, 68, 19–46. [Google Scholar] [CrossRef]
  31. Soil Survey Staff, Natural Resources Conservation Service, United States Department of Agriculture. Web Soil Survey. Available online: https://websoilsurvey.sc.egov.usda.gov/ (accessed on 16 December 2019).
  32. Kershner, B.; Tufts, C.; Nelson, G. National Wildlife Federation Field Guide to Trees of North America; Sterling Publishing Company: New York, NY, USA, 2008. [Google Scholar]
  33. Little, E. The Audubon Society Field Guide to North American Trees: Western Region; Knopf: New York, NY, USA, 1988. [Google Scholar]
  34. Calscape. Available online: https://calscape.org/ (accessed on 21 December 2019).
  35. GBIF Home Page. Available online: GBIF.org.
  36. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  37. Zomer, R.J.; Trabucco, A.; Bossio, D.A.; Verchot, L.V. Climate change mitigation: A spatial analysis of global land suitability for clean development mechanism afforestation and reforestation. Agric. Ecosyst. Environ. 2008, 126, 67–80. [Google Scholar] [CrossRef]
  38. Barton, K. Package “MuMIn”. R Package Version 2015, 1, 18. [Google Scholar]
  39. Smith, S.A.; Brown, J.W. Constructing a broadly inclusive seed plant phylogeny. Am. J. Bot. 2018, 105, 302–314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Slater, G.J.; Harmon, L.J.; Wegmann, D.; Joyce, P.; Revell, L.J.; Alfaro, M.E. Fitting models of continuous trait evolution to incompletely sampled comparative data using approximate Bayesian computation. Evolution 2012, 66, 752–762. [Google Scholar] [CrossRef] [PubMed]
  41. Paradis, E.; Schliep, K. Ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 2018, 35, 526–528. [Google Scholar] [CrossRef]
  42. Pennell, M.W.; Eastman, J.M.; Slater, G.J.; Brown, J.W.; Uyeda, J.C.; FitzJohn, R.G.; Alfaro, M.E.; Harmon, L.J. Geiger v2.0: An expanded suite of methods for fitting macroevolutionary models to phylogenetic trees. Bioinformatics 2014, 30, 2216–2218. [Google Scholar] [CrossRef]
  43. Harmon, L.J.; Weir, J.T.; Brock, C.D.; Glor, R.E.; Challenger, W. GEIGER: Investigating evolutionary radiations. Bioinformatics 2007, 24, 129–131. [Google Scholar] [CrossRef] [Green Version]
  44. Pinheiro, J.; Bates, D.; DebRoy, S.; Sarkar, D.; Team, R.C. Nlme: Linear and Nonlinear Mixed Effects Models. R Package Version 2019, 3, 1–140. [Google Scholar]
  45. Baldwin, B.G.; Thornhill, A.H.; Freyman, W.A.; Ackerly, D.D.; Kling, M.M.; Morueta-Holme, N.; Mishler, B.D. Species richness and endemism in the native flora of California. Am. J. Bot. 2017, 104, 487–501. [Google Scholar] [CrossRef] [Green Version]
  46. Anten, N.P.R.; Schieving, F. The role of wood mass density and mechanical constraints in the economy of tree architecture. Am. Nat. 2010, 175, 250–260. [Google Scholar] [CrossRef] [PubMed]
  47. Hacke, U.G.; Sperry, J.S.; Pockman, W.T.; Davis, S.D.; McCulloh, K.A. Trends in wood density and structure are linked to prevention of xylem implosion by negative pressure. Oecologia 2001, 126, 457–461. [Google Scholar] [CrossRef] [PubMed]
  48. Lachenbruch, B.; McCulloh, K.A. Traits, properties, and performance: How woody plants combine hydraulic and mechanical functions in a cell, tissue, or whole plant. New Phytol. 2014, 204, 747–764. [Google Scholar] [CrossRef] [PubMed]
  49. Mahall, B.E.; Tyler, C.M.; Cole, E.S.; Mata, C. A comparative study of oak (Quercus, Fagaceae) seedling physiology during summer drought in southern California. Am. J. Bot. 2009, 96, 751–761. [Google Scholar] [CrossRef]
  50. Mclaughlin, B.C.; Zavaleta, E.S. Predicting species responses to climate change: Demography and climate microrefugia in California valley oak (Quercus lobata). Glob. Chang. Biol. 2012, 18, 2301–2312. [Google Scholar] [CrossRef]
  51. Tarelkin, Y.; Hufkens, K.; Hahn, S.; Van den Bulcke, J.; Bastin, J.F.; Ilondea, B.A.; Debeir, O.; Van Acker, J.; Beeckman, H.; De Cannière, C. Wood anatomy variability under contrasted environmental conditions of common deciduous and evergreen species from central African forests. Trees 2019, 33, 893–909. [Google Scholar] [CrossRef]
  52. Pretzsch, H.; Biber, P.; Schütze, G.; Kemmerer, J.; Uhl, E. Wood density reduced while wood volume growth accelerated in Central European forests since 1870. For. Ecol. Manag. 2018, 429, 589–616. [Google Scholar] [CrossRef]
  53. Sapes, G.; Serra-Diaz, J.M.; Lloret, F. Species climatic niche explains drought-induced die-off in a Mediterranean woody community. Ecosphere 2017, 8, e01833. [Google Scholar] [CrossRef]
  54. Brando, P.M.; Nepstad, D.C.; Balch, J.K.; Bolker, B.; Christman, M.C.; Coe, M.; Putz, F.E. Fire-induced tree mortality in a neotropical forest: The roles of bark traits, tree size, wood density and fire behavior. Glob. Chang. Biol. 2012, 18, 630–641. [Google Scholar] [CrossRef]
Forests 11 00105 g001 550
Figure 1. Map of study sites. Green markers represent the geographic location of each site.
Figure 1. Map of study sites. Green markers represent the geographic location of each site.
Forests 11 00105 g001
Forests 11 00105 g002 550
Figure 2. Branch wood density (g/cm3) vs basal diameter to height ratio for California trees and shrubs. Trees with a low basal diameter to height ratio are tall and skinny, while trees with a high basal diameter to height ratio are short and wide, as represented by tree symbols. The arrow represents increasing basal diameter to height ratio. Red points/lines indicate riparian specialists and black points/lines non-riparian specialists. Dotted lines show the model fit for the habitat only linear model (Table 1), while solid lines show the model fit for the taper × habitat model.
Figure 2. Branch wood density (g/cm3) vs basal diameter to height ratio for California trees and shrubs. Trees with a low basal diameter to height ratio are tall and skinny, while trees with a high basal diameter to height ratio are short and wide, as represented by tree symbols. The arrow represents increasing basal diameter to height ratio. Red points/lines indicate riparian specialists and black points/lines non-riparian specialists. Dotted lines show the model fit for the habitat only linear model (Table 1), while solid lines show the model fit for the taper × habitat model.
Forests 11 00105 g002
Forests 11 00105 g003 550
Figure 3. Branch wood density vs. (a) species 5% and 95% climate moisture deficit (CMD) and (b) species 5% CMD for riparian vs. non-riparian California trees and shrubs. Black dots represent non-riparian specialist species, while red dots represent species that are riparian habitat specialists. The black line indicates a significant relationship between wood density and CMD at the 5th percentile, but no significant effect at the 95th percentile of CMD. PET, potential evapotranspiration; MAP, mean annual precipitation.
Figure 3. Branch wood density vs. (a) species 5% and 95% climate moisture deficit (CMD) and (b) species 5% CMD for riparian vs. non-riparian California trees and shrubs. Black dots represent non-riparian specialist species, while red dots represent species that are riparian habitat specialists. The black line indicates a significant relationship between wood density and CMD at the 5th percentile, but no significant effect at the 95th percentile of CMD. PET, potential evapotranspiration; MAP, mean annual precipitation.
Forests 11 00105 g003
Forests 11 00105 g004 550
Figure 4. Branch wood density g/cm3 vs. (a) 5th percentile climate moisture deficit (mm) and (b) 5th percentile Tmin (°C) for California trees and shrubs. Taper vs. (c) mean climate moisture deficit and (d) 5th percentile Tmin for California trees and shrubs.
Figure 4. Branch wood density g/cm3 vs. (a) 5th percentile climate moisture deficit (mm) and (b) 5th percentile Tmin (°C) for California trees and shrubs. Taper vs. (c) mean climate moisture deficit and (d) 5th percentile Tmin for California trees and shrubs.
Forests 11 00105 g004
Forests 11 00105 g005 550
Figure 5. (a) Boxplot of wood density values and (b) taper (basal diameter to height ratio) values per geoflora type. NS denotes non-significance.
Figure 5. (a) Boxplot of wood density values and (b) taper (basal diameter to height ratio) values per geoflora type. NS denotes non-significance.
Forests 11 00105 g005
Forests 11 00105 g006 550
Figure 6. Branch wood density vs basal diameter to height ratio (taper) for (a) riparian and (b) non-riparian species, colored by the 5th percentile climate moisture deficit shown for each species (drier in lighter blue). Trend lines are only for the wood density (WD)~taper relationships.
Figure 6. Branch wood density vs basal diameter to height ratio (taper) for (a) riparian and (b) non-riparian species, colored by the 5th percentile climate moisture deficit shown for each species (drier in lighter blue). Trend lines are only for the wood density (WD)~taper relationships.
Forests 11 00105 g006
Table
Table 1. Model selection for structure and habitat specialization. WD refers to branch wood density. Taper is basal diameter to height ratio. Habitat refers to whether or not species are riparian specialists. AIC, Aikike information criterion.
Table 1. Model selection for structure and habitat specialization. WD refers to branch wood density. Taper is basal diameter to height ratio. Habitat refers to whether or not species are riparian specialists. AIC, Aikike information criterion.
ModelDegrees of FreedomAICc
Null (WD~1)2−24.22
WD~taper3−22.29
WD~habitat3−31.84
WD~taper + habitat4−29.22
WD~taper + habitat + taper/habitat5−29.49
WD~taper + taper/habitat4−32.95
Table
Table 2. Model selection for wet edge climate effects. WD refers to branch wood density. Habitat refers to whether or not species are riparian specialists. CMD 0.05 is the 5th percentile climate moisture deficit. MAP 0.05 is the 5th percentile for mean annual precipitation. MAT 0.05 is the 5th percentile for mean annual temperature. PET 0.05 is the 5th percentile for potential evapotranspiration.
Table 2. Model selection for wet edge climate effects. WD refers to branch wood density. Habitat refers to whether or not species are riparian specialists. CMD 0.05 is the 5th percentile climate moisture deficit. MAP 0.05 is the 5th percentile for mean annual precipitation. MAT 0.05 is the 5th percentile for mean annual temperature. PET 0.05 is the 5th percentile for potential evapotranspiration.
ModelDegrees of FreedomAICc
Null2−21.20
WD~CMD 0.053−29.21
WD~MAP 0.053−22.04
WD~PET 0.053−28.84
WD~MAT 0.053−25.75
WD~CMD 0.05 + habitat4−36.07
WD~CMD 0.05 + habitat/CMD 0.054−27.80
WD~CMD 0.05 + habitat + habitat/CMD 0.055−32.58
Table
Table 3. Model selection for most parsimonious model overall. WD refers to branch wood density. Taper is basal diameter to height ratio. Habitat refers to whether or not species are riparian specialists. CMD 0.05 is the 5th percentile climate moisture deficit.
Table 3. Model selection for most parsimonious model overall. WD refers to branch wood density. Taper is basal diameter to height ratio. Habitat refers to whether or not species are riparian specialists. CMD 0.05 is the 5th percentile climate moisture deficit.
ModelDegrees of FreedomAICc
Null (WD~1)2−21.20
WD~CMD 0.053−29.21
WD~CMD 0.05 + taper + habitat + taper/habitat6−37.76
WD~CMD 0.05 + taper + taper/habitat5−41.74
WD~taper + habitat + taper/habitat5−29.49
WD~CMD 0.05 + habitat4−36.07
WD~habitat3−31.84
WD~taper3−22.29

Share and Cite

MDPI and ACS Style

Nelson, R.A.; Francis, E.J.; Berry, J.A.; Cornwell, W.K.; Anderegg, L.D.L. The Role of Climate Niche, Geofloristic History, Habitat Preference, and Allometry on Wood Density within a California Plant Community. Forests 2020, 11, 105. https://0-doi-org.brum.beds.ac.uk/10.3390/f11010105

AMA Style

Nelson RA, Francis EJ, Berry JA, Cornwell WK, Anderegg LDL. The Role of Climate Niche, Geofloristic History, Habitat Preference, and Allometry on Wood Density within a California Plant Community. Forests. 2020; 11(1):105. https://0-doi-org.brum.beds.ac.uk/10.3390/f11010105

Chicago/Turabian Style

Nelson, Rebecca A., Emily J. Francis, Joseph A. Berry, William K. Cornwell, and Leander D. L. Anderegg. 2020. "The Role of Climate Niche, Geofloristic History, Habitat Preference, and Allometry on Wood Density within a California Plant Community" Forests 11, no. 1: 105. https://0-doi-org.brum.beds.ac.uk/10.3390/f11010105

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop