Next Article in Journal
Effects of Stand Density, Age, and Drought on the Size–Growth Relationship in Larix principis-rupprechtii Forests
Next Article in Special Issue
Can Growth Increase of Small Trees after Drought Compensate for Large Trees’ Growth Loss?
Previous Article in Journal
Production of Chestnut Coppice Biomass in a Framework of Low Mechanization and High Expectations to Combat Climate Change and Other Social Expectations
Previous Article in Special Issue
Tracheids vs. Tree Rings as Proxies for Dendroclimatic Reconstruction at High Altitude: The Case of Pinus sibirica Du Tour
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Forest Adaptation to Climate Change: Altitudinal Response and Wood Variation in Natural-Growth Cunninghamia lanceolata in the Context of Climate Change

College of Forestry, Central South University of Forestry and Technology, Changsha 410004, China
*
Author to whom correspondence should be addressed.
Forests 2024, 15(3), 411; https://doi.org/10.3390/f15030411
Submission received: 26 January 2024 / Revised: 18 February 2024 / Accepted: 19 February 2024 / Published: 21 February 2024
(This article belongs to the Special Issue Tree Growth in Relation to Climate Change)

Abstract

:
This research delves into the impact of climate change on the wood traits of Cunninghamia lanceolata across various altitudinal gradients, aiming to understand the influence of altitude and climatic factors like temperature and precipitation on key wood characteristics. Employing a comprehensive approach, samples were collected from different altitudes for detailed phenotypic analysis. Methods included Pearson correlation, principal component analysis, cluster analysis, and random forest analysis. Results revealed significant variations in wood traits such as heartwood ratio, tracheid length, and width across altitudes. Notably, wood traits in lower- and middle-elevation populations exhibited higher variability compared to higher elevations, indicating greater environmental diversity and genetic adaptability at these altitudes. Climatic factors, particularly temperature and precipitation, were found to increasingly influence wood trait variation with altitude. The research concludes that the adaptation of Cunninghamia lanceolata to climate change is significantly influenced by both altitudinal and climatic factors, highlighting their importance in forest genetic breeding and conservation strategies amidst global climate change.

1. Introduction

Cunninghamia lanceolata, endemic to China, is a predominant afforestation species in southern China. It covers one-fourth of the country’s artificial arboreal forests in area and one-third in volume, amounting to 9.87 million hectares and 755 million cubic meters, respectively, leading the nation in both aspects [1]. With a long history of cultivation and a wide range of applications, natural occurrences of Cunninghamia lanceolata have become exceedingly rare due to prolonged human interference and high levels of gene flow, resulting in minimal genetic differentiation among different sources of Cunninghamia lanceolata [2]. This suggests a potential weakening of adaptability against the backdrop of climate change [3]. In 2012, our research team discovered a natural population of Cunninghamia lanceolata in the Xiaoxi National Nature Reserve in Hunan Province, a narrowly distributed ecotype. Locally known as “Iron-Heart Cunninghamia lanceolata”, this variant is distinguished by its dense texture, good rot resistance, high heartwood ratio, and brownish heartwood [4]. Its mechanical properties, such as bending strength, compressive strength, and shrinkage, are significantly superior to other populations of Cunninghamia lanceolata [5]. The quality of wood is a crucial indicator for timber forest cultivation [6]. With the shift in China’s timber market demand, developing precious native timber species is a key step towards modernizing forestry [7]. The superior wood quality of this regional Cunninghamia lanceolata has attracted widespread attention. Extensive research has been conducted on its wood physical properties [4,8], seed traits [9], breeding population construction [10,11], and spatial genetic structure [12]. To further explore the excellent germplasm resources of Cunninghamia lanceolata and enhance the sustainable utilization and management of these resources, our research focuses on wood properties of natural populations of Cunninghamia lanceolata and their correlation with altitudinal climate change.
Genetic diversity, encompassing the genetic variation within biological populations and individuals [13], forms the foundation of species diversity and represents the core of biodiversity. It reflects the genetic richness of species and determines their potential to adapt to environmental changes [14,15]. Phenotypic diversity, as the external manifestation of genetic diversity, results from the interplay of genes and the environment [16]. The phenotypic variation among individuals within a species is influenced not only by specific habitat conditions but also by historical processes and phylogeny [17]. In different environmental settings, populations of the same plant species may undergo varying selective pressures, leading to genetic and phenotypic distinctions among populations [18]. Thus, researching the phenotypic diversity of forest tree populations in diverse ecological environments not only unveils the extent of genetic variation but also aids in understanding the evolutionary level and environmental adaptability of plants [19,20]. Studies have demonstrated that wood traits are significantly influenced by environmental factors such as geographic location (latitude and longitude) [21], altitude [22], light [23], soil [24], and climate [25,26].
According to relevant statistics, by 2020, the global temperature had risen by approximately 0.99 °C compared to the pre-industrial era (1850–1900), leading to an increase in extreme weather events worldwide [27]. Global climate change is accelerating its impact on forest disturbance mechanisms and adaptability [28], with effects on forests being dual-natured; they can be either positive or negative [29]. The adaptive capacity of tree populations to climate change depends on intraspecific genetic variation and phenotypic plasticity [30]. Altitude, a pivotal ecological factor, influences plant growth through gradient changes in temperature, precipitation, light, soil, and other factors [31]. This contributes to adaptive variation in plant phenotypic traits, making altitude an ideal subject for researching genetic variation [32,33]. While previous research has explored wood trait variation across altitudinal gradients from ecological [34,35], physiological [36,37,38], and dendrochronological perspectives [39,40], there remains a gap in understanding the genetic variation of these traits and their impact on the environmental adaptability of trees, particularly from a genetic standpoint. The upper limit of the vertical distribution of Cunninghamia lanceolata often varies with different topographical and climatic conditions. In the high-latitude Dabie Mountains, it grows below an altitude of 700 m, while in the lower latitude region of Dali, Yunnan, it can be found at elevations as high as 2500 m [41]. Investigating the variation of Cunninghamia lanceolata across different altitudinal gradients and the influence of altitudinal climatic factors on wood traits can provide insights into the species’ response to climate change and the role of the environment in shaping traits [42].
Wood traits encompass a range of characteristics, primarily including anatomical traits such as tracheids and microfibril angles, physical traits like wood density, shrinkage and swelling properties, mechanical properties, and chemical traits involving the content of cellulose, hemicellulose, lignin, and various metabolites [43]. Wood density is a pivotal factor affecting plant ecological strategies, influenced by vascular traits and tree growth rates [44,45]. Denser wood is not only harder and more decay-resistant but also exhibits greater wind resistance [46]. In this research, seven representative wood traits were meticulously selected for comprehensive analysis. The annual ring width, indicative of the wood layer formed during a tree’s one-year growth cycle, serves not only as a crucial indicator of growth rate [47] but also mirrors historical climate changes, ecosystem dynamics, and trees’ adaptability to environmental shifts [48]. The heartwood ratio, wood density, and water absorption rate are pivotal in assessing wood quality. The heartwood ratio often signifies wood quality, while wood density and water absorption directly influence wood’s processing performance, strength, and durability [49]. Wood with high water absorption is susceptible to swelling, deformation, and biodegradation [50]. Moreover, variations in wood density, a vital attribute in trees’ carbon storage strategy, emerge as significant indicators for climate change assessment [51,52]. Lastly, the dimensions, shape, and arrangement of tracheids, crucial xylem cells in gymnosperms responsible for water and nutrient transport, not only impact the wood’s physical and mechanical properties [53] but also reflect the plant’s adaptability to its growth environment [54].
In this research, wood core samples were meticulously collected from four distinct altitudinal populations of Cunninghamia lanceolata within its natural habitat. We conducted a thorough analysis of phenotypic variation in seven key wood traits. The objective was to unravel the patterns of phenotypic differentiation and variation in Cunninghamia lanceolata wood traits along altitudinal gradients and to investigate the correlations between these traits and corresponding climatic factors at different altitudes. The findings of this research significantly enhance our understanding of the growth characteristics and adaptive mechanisms of Cunninghamia lanceolata across various altitudes. This is crucial for the discovery and conservation of high-quality germplasm resources of Cunninghamia lanceolata, as well as for bolstering its conservation and sustainable utilization. Furthermore, in the context of global climate change, these insights are pivotal in forecasting future shifts in forest ecosystems and elucidating the mechanisms through which trees adapt to environmental changes.

2. Materials and Methods

2.1. Overview of the Research Area

Xiaoxi National Nature Reserve (hereinafter referred to as Xiaoxi) is located in Xiangxi Tujia and Miao Autonomous Prefecture of Hunan Province (Figure 1), east of the Yunnan–Guizhou Plateau and the central part of the Wuling Mountain Range, which is the intersection of China’s second and third ladders, with the highest altitude of 1327.1 m, and latitude and longitude of 110°6′50″–110°21′35″ E, 28°42′15″–28°53′55″ N. The climate is a subtropical humid monsoon climate, a warm and humid climate with abundant rainfall and an annual average temperature of 13~15 °C, annual precipitation of 1300~1400 mm, and relative humidity of 79%. The soil in Xiaoxi shows distinct zonal characteristics: below 400 m in altitude, there are mountain red soils; between 400 m and 500 m, mountain yellow–red soils; from 500 m to 800 m, mountain yellow soils; and above 800 m, mountain yellow–brown soils. Due to the diversity of habitats in Xiaoxi, the area boasts an exceptionally rich variety of plant species. According to surveys and analyses, there are a total of 2252 species of seed plants in Xiaoxi, including 20 species of gymnosperms and 2232 species of angiosperms. The reserve is home to 43 species of nationally protected plants of Classes I and II, such as Manglietia decidua, Davidia involucrata Baill, Taxus wallichiana var. mairei, and Bretschneidera sinensis [55]. Cunninghamia lanceolata predominantly exists as scattered individuals in primary and secondary forests, mainly distributed at altitudes between 500 m and 1200 m.

2.2. Sample Collection

In November 2022, a comprehensive field survey was conducted on the naturally grown Cunninghamia lanceolata in the Xiaoxi National Nature Reserve. The species was categorized into four altitudinal populations, each separated by an interval of 200 m: population 1 (Pop1) at 550 m, population 2 (Pop2) at 750 m, population 3 (Pop3) at 950 m, and population 4 (Pop4) at 1150 m. For each population, 30 adult individual plants exhibiting normal growth and free from apparent pests and diseases were selected, with a consistent diameter at breast height (DBH, 1.3 m). Due to prior human interference and the species’ altitudinal distribution limits, only 8 sample trees were collected for Pop4. A minimum distance of 100 m was maintained between each plant. Using a 5 mm growth cone, a core sample was extracted through the pith from south to north at a height of 1.3 m from the target tree. The core was then placed into a PVC pipe and numbered for subsequent wood trait analysis. Concurrently, the DBH and height of each sample tree were measured, and GPS coordinates were recorded to locate the sample sites. Detailed information about the sample sites and associated meteorological factors is presented in Table S1.

2.3. Measurement of Traits

Wood properties assessed in this research encompass seven key attributes: average annual ring width (Rb), heartwood ratio (P), basic wood density (WHD), water absorption (Hy), tracheid length (L), tracheid width (D), and the tracheid length-to-width ratio (L/D).
As per the technical guidelines, the wood core was first smoothed with 600# sandpaper. A line was drawn perpendicular to the annual rings in the radial direction, and the total width of the complete annual rings was measured along this line. The total width was precisely measured to 0.01 mm using a ruler, and the number of annual rings within the measurement range was counted. The average width of the annual rings was then calculated by dividing the total width by the number of rings [56]. The heartwood ratio (P) was determined by dividing the disc area, calculated using the radius from the edge of the heartwood to the pith, by the total disc area. This was carried out by measuring the length of the heartwood and the cores, with P calculated as r2/R2 × 100%, where r and R represent the lengths of the heartwood and the core, respectively (Figure 2a). Basic wood density (WHD) was measured using the saturated drainage method. The cores were immersed in water, with water changes every two days, and weighed periodically until a constant saturated water content weight (W1) was achieved. The cores were then dried in an oven at 103 + 2 °C to a constant weight (W2); WHD was calculated as 1(W1/W2 − 0.346), and Hy as (W1 − W2)/W2. The cores were divided radially into three sections (near the pith, middle section, near the bark) and dissociated using a mixture of glacial acetic acid and 30% hydrogen peroxide solution. The dissociated samples were then observed under a bio-digital microscope (model OLYMPUS-BX51) to analyze the tracheids’ morphological characteristics (Figure 2b). The lengths and widths of the tracheids were measured 20 times for each sample using Photoshop 2022 software, with the average of these measurements from the three samples taken as the final value. Finally, the tracheid length-to-width ratio was calculated [57].

2.4. Data Processing and Analysis

Preliminary processing of the measured data was carried out using Excel, and the population mean, standard deviation, intra- and interpopulation coefficients of variation, and Nestle’s analysis of variance (ANOVA) were calculated for each trait using SPSS 27.
(1) Inter- and intrapopulation variation was analyzed based on Nestle’s ANOVA model (significance tests were performed using the Duncan method of multiple comparisons), with a linear model of:
Y i j k = μ + τ i + δ j ( i ) + ε i j k
where: Yijk is the kth observation of the jth monoculture in the ith population, μ is the overall mean, τi is the between-population effect value, δj(i) is the within-population monoculture random effect value, and εijk is the random error.
(2) The coefficient of variation (CV) between populations was calculated using the following formula:
C V = δ / x ¯
where: δ and x represent the standard deviation and mean of the trait, respectively.
The climate data for this research were sourced from the World Climate Database (https://worldclim.org/data/worldclim21.html, accessed on 26 October 2023). This database primarily compiles climate observation records from weather stations globally, spanning from 1950 to 2000. The data were generated through interpolation, utilizing a global climate database with a spatial resolution of 30′ [58]. Utilizing ArcGIS 10.2 software, we obtained estimated meteorological data for each altitudinal sampling point. Subsequently, six representative meteorological factors were extracted from a total of nineteen climatic factors through principal component analysis. These factors included the annual mean temperature (BIO1), the average temperature of the wettest season (BIO8), the average temperature of the driest season (BIO9), the annual precipitation (BIO12), the precipitation of the wettest season (BIO16), and the precipitation of the driest season (BIO17). The specifics of these factors are detailed in Table S1.
Data processing and analysis in this research were conducted within the R language environment. Principal component analysis (PCA) was executed using the FactoMineR package, complemented by graphical presentations created with the factoextra and ggplot2 packages. Cluster analysis employed the Euclidean distance matrix, computed via the ape package, and hierarchical clustering was performed using the Ward.D method, with graphical outputs generated by the RColorBrewer package. The linkET package facilitated correlation analysis, calculating Pearson correlation coefficients and significance p-values, which were further validated through the Mantel test. Data processing and transformation were adeptly handled by the dplyr and magrittr packages, with all graphs being plotted using ggplot2.
To uncover the primary drivers of wood trait variation, four traits—heartwood ratio, tracheid length, tracheid width, and tracheid aspect ratio—showing significant correlations with altitudinal climatic factors, were selected based on correlation analysis results. These traits were then subjected to random forest (RF) analysis [59], utilizing the randomForest package to calculate feature importance. The contribution rate of altitudinal climatic factors to the variation of these four wood traits was determined based on this feature importance. The contributions were then visualized using Origin 2022 software.

3. Results

3.1. Analysis of Variance of Wood Traits of Cunninghamia lanceolata at Different Altitudes

Nested ANOVA conducted on the wood traits of Cunninghamia lanceolata across different altitudes revealed significant variations (refer to Table 1). There were highly significant differences observed among Cunninghamia lanceolata populations in terms of heartwood ratio, tracheid length, and tracheid length-to-width ratio (p < 0.001). Notable differences were also evident in the average annual ring width and tracheid width (p < 0.05). However, no significant differences were found in basic wood density and water absorption rate. The within-population ANOVA indicated a lack of significant differences across all seven wood traits.
Multiple comparisons and comparing the mean values of wood traits at different altitudes (Table 2) showed that the basic wood density, tracheid length, tracheid width, and tracheid length-to-width ratio were the largest in Pop1, reaching 0.42 g/cm3, 2888.06 mm, 45.22 mm, and 67.45, respectively; and they were the smallest in Pop4, with 0.39 g/cm3, 2283.09 mm, 40.99 mm, and 58.20. The mean heartwood ratio and water absorption were the largest in Pop3, reaching 30.51 mm, and the smallest in Pop1, with only 24.35 mm. The average annual ring width was the largest in Pop3, reaching 30.51 mm, and the smallest in Pop1, only 24.35 mm. Heartwood ratio and water absorption were the largest in Pop3, 77.51% and 191.9%, respectively, and the smallest in Pop1, only 53.95% and 174.97%, respectively.
The variation of each wood trait showed (see Table 2) that the coefficient of variation of wood traits was 13.55% on average, and the level of variation of different wood traits varied greatly, with the average variation of each trait ranging from 9.40% to 24.31%, with the smallest variation being in tracheid length (9.40%) and the largest being in the average annual ring width (24.31%). The magnitude of variation of wood traits within populations at different altitudes was also significantly different, with the largest magnitude of variation in Pop3 (9.4%~36.38%) and the smallest magnitude of variation in Pop4 (7.81%~13.38%). In terms of the average coefficient of variation of wood traits within the populations, the largest and smallest coefficients of variation were found in Pop3 (15.14%) and Pop4 (10.27%), respectively.

3.2. Principal Component Analysis and Cluster Analysis of Wood Traits of Cunninghamia lanceolata at Different Altitudes

3.2.1. Principal Component Analysis of Cunninghamia lanceolata Wood Traits at Different Altitudes

The principal component analysis (PCA) conducted on the Cunninghamia lanceolata populations at different altitudes identified three principal components with eigenvalues exceeding 1. The cumulative contribution rate of these first three principal components to the variance in the seven wood traits was 79.79% (as shown in Figure 3, Tables S2 and S3), effectively encapsulating most of the information pertinent to the characteristic variables of the wood traits. The first principal component, accounting for 37.98% of the total variance, included key variables such as water absorption rate, basic wood density, average annual ring width, and the aspect ratio of Cunninghamia lanceolata. The second principal component, encompassing tracheid length and heartwood ratio, explained 23.14% of the total variance. The third principal component, comprising tracheid width and tracheid aspect ratio, accounted for 10.82% of the total variance. The cumulative contributions of these principal components, in descending order of significance, were for tracheid width, basic wood density, water absorption, tracheid aspect ratio, tracheid length, mean annual ring width, and heartwood ratio.

3.2.2. Cluster Analysis of Wood Traits of Cunninghamia lanceolata at Different Altitudes

The clustering analysis of 12 wood traits from four Cunninghamia lanceolata populations at different altitudes revealed distinct clusters, as illustrated in Figure 4. The analysis categorized the 98 samples into four clusters. Cluster I comprised 27 samples, representing 27.55% of the total, predominantly including 16 samples from Pop1 and 7 from Pop2. Cluster II, with the fewest samples, contained 10 samples, accounting for 12.20% of the total, mainly from Pop3. Cluster III encompassed 20 samples, constituting 20.41% of the total, primarily including eight samples from Pop3 and six from Pop4. Cluster IV, the largest population, contained 41 samples, making up 41.84% of the total, with 16 samples from Pop2, 13 from Pop1, and 12 from Pop3.

3.3. Correlation between Cunninghamia lanceolata Wood Traits and Altitude Climate Factors at Different Altitudes

3.3.1. Correlation Analysis between Wood Traits of Cunninghamia lanceolata at Different Altitudes

Pearson’s correlation analysis was conducted among various wood traits of Cunninghamia lanceolata (as depicted in Figure 5 and detailed in the accompanying Table S4). The analysis revealed that the mean annual ring width was significantly negatively correlated with the tracheid aspect ratio (−0.2449 *) and tracheid length (−0.2363 *). It also showed a significant negative correlation with wood density (r = −0.4641 ***) and a significant positive correlation with water absorption rate (r = 0.4668 ***). The heartwood ratio exhibited a significant negative correlation with the tracheid aspect ratio (−0.3011 **) and a highly significant negative correlation with tracheid length (−0.3525 ***). Wood density was found to have a significant positive correlation with the tracheid aspect ratio (0.2182 *) and a highly significant negative correlation with water absorption (−0.9908 ***). The tracheid aspect ratio demonstrated a highly significant negative correlation with water absorption (−0.2280 *) and tracheid width (−0.4466 ***), indicating significant and highly significant negative correlations, respectively. Furthermore, tracheid length showed a highly significant positive correlation with tracheid width (0.3570 ***) and the tracheid length-to-width ratio (0.6573 ***).

3.3.2. Correlation Analysis of Wood Traits of Cunninghamia lanceolata at Different Altitudes with Altitude Climatic Factors

The correlation analysis between wood traits and altitudinal–climatic factors revealed highly significant correlations for three traits: heartwood ratio, tracheid length, and tracheid length–width ratio (Figure 5 and Table S4). The results of detailed analyses are as follows: there was a significant positive correlation between the average annual ring width and altitude (0.2852 *). A highly significant positive correlation was observed between the heartwood ratio and altitude (0.5793 ***), annual precipitation (0.4752 ***), precipitation of the wettest season (0.4153 ***), and precipitation of the driest season (0.4687 ***), along with a highly significant negative correlation with annual average temperature (−0.5019 ***), average temperature of the wettest season (−0.5004 ***), and average temperature of the driest season (−0.5061 ***). Tracheid length demonstrated a highly significant negative correlation with altitude (−0.5525 ***), annual precipitation (−0.4516 ***), precipitation of the wettest season (−0.4497), and precipitation of the driest season (−0.4369 ***) and a highly significant positive correlation with annual average temperature (0.4634 ***), average temperature of the wettest season (0.4626 ***), and average temperature of the driest season (0.4630 ***). Tracheid width showed a significant negative correlation with altitude (−0.2488 *), annual precipitation (−0.2030 *), and precipitation of the driest season (−0.2123 *) and a significant positive correlation with annual average temperature (0.2111 *), average temperature of the wettest season (0.2103 *), and average temperature of the driest season (0.2135 *). The tracheid length–width ratio exhibited a highly significant negative correlation with altitude (−0.3608 ***) and precipitation of the wettest season (−0.3303), a significant negative correlation with annual precipitation (−0.3100 *) and precipitation of the driest season (−0.2960 *), and a significant positive correlation with annual average temperature (0.3095 **), average temperature of the wettest season (0.3095 **), and average temperature of the driest season (0.3080 **). Notably, wood density and water absorption rate showed no significant correlation with altitudinal–climatic factors.

3.3.3. Mantel Test of Cunninghamia lanceolata Populations at Different Altitudes with Altitudinal Climatic Factors

The Mantel test conducted on the wood traits of Cunninghamia lanceolata at different altitudes in relation to altitudinal climatic factors revealed (as shown in Figure 5, Table S5) a highly significant correlation between the wood traits of the Pop3 population and the driest seasonal precipitation (BIO17). Additionally, a highly significant correlation was observed between the wood traits of the Pop4 population and altitude. The wood traits of the Pop3 population exhibited significant or highly significant correlations with altitudinal climate factors. In contrast, the wood traits of Pop1 and Pop4 populations showed partial significant correlations with these factors, whereas the wood traits of the Pop2 population did not display any significant correlation with altitudinal climate factors.

3.4. Random Forest Analysis of Wood Traits of Cunninghamia lanceolata at Different Altitudes

Random forest analysis results (illustrated in Figure 6a) suggested that the altitude factor predominantly influenced the variation of four wood traits: heartwood ratio, tracheid length, tracheid width, and tracheid length-to-width ratio, with contribution rates of 48.81%, 45.71%, 61.94%, and 60.92%, respectively. The variances of these four wood traits were differentially impacted by climatic factors. For instance, the average temperature during the wettest season was the most influential climatic factor for the heartwood ratio (10.18%), tracheid length (9.97%), and tracheid length-to-width ratio (7.13%). In contrast, the precipitation during the driest season had the most significant impact on tracheid width (7.55%). Additionally, the cumulative contributions of temperature factors (BIO1, BIO8, BIO9) to the heartwood ratio, tracheid length, tracheid width, and tracheid length-to-width ratio were 29.72%, 28.22%, 19.53%, and 20.70%, respectively. These were higher than the cumulative contributions of precipitation factors (BIO12, BIO16, and BIO17) to the four wood traits, which were 21.46%, 26.07%, 18.53%, and 18.38%, respectively. This indicates that among the climatic factors, temperature had a greater contribution to the variation of wood traits compared to precipitation.
Further analysis revealed significant differences in the contribution of altitudinal climatic factors to the variability of four wood traits across the altitudinal gradient (as shown in Figure 6b). In the middle- and low-altitude gradients (populations Pop1, Pop2, and Pop3), the altitude factor was the predominant influence on the variation of the four wood traits. However, with increasing altitude, the influence of the altitude factor on the variance of heartwood ratio decreased, while its effects on tracheid length and tracheid aspect ratio initially increased and then decreased. The impact on tracheid width showed an overall decreasing trend. This decrease in the influence of the altitude factor was accompanied by an increased impact of the climatic factor on the four wood traits. Notably, except for the heartwood ratio, the cumulative contribution of climatic factors to tracheid length, tracheid width, and tracheid aspect ratio was greater than that of the altitude factor at higher altitudes (population Pop4), reaching 67.17%, 65%, and 72.16%, respectively. These factors became the primary influences on the wood traits at higher altitudes.

4. Discussion

4.1. Analysis of Variation in Wood Traits of Cunninghamia lanceolata at Different Altitudes

Plant phenotypic traits, resulting from the interplay between genetics and environment, not only mirror genetic diversity but are also pivotal in identifying superior germplasm resources [12]. This research delved into the variation of wood traits in Cunninghamia lanceolata across different altitudes, examining the influence of altitudinal climatic factors. Nested ANOVA results indicated highly significant differences (p < 0.001) in heartwood ratio, tracheid length, and tracheid length-to-width ratio among various altitudes. Mean annual ring width and tracheid width also showed significant differences (p < 0.05). These findings align with Elif Topaloğlu et al.’s research [60] on the anatomical properties of oriental beech wood, underscoring the variability of wood traits across altitudes. However, basic wood density and water absorption rate did not exhibit significant differences across altitudes, suggesting a stronger genetic influence on these traits, with a lesser role for environmental factors [61]. The within-population ANOVA revealed no significant differences in any of the seven wood traits, contrasting with Hongjing Duan et al.’s findings [57], which showed significant within-population variation in Cunninghamia lanceolate wood traits from different seed sources. This discrepancy might stem from factors like similar selection pressures [62], inbreeding [63], gene flow [64], and human interference [65] in Cunninghamia lanceolata populations in the same altitudinal gradient. The variability of Cunninghamia lanceolata wood traits across altitudes ranged from 9.40% to 24.31%, relatively modest compared to the variability (6.43% to 45.55%) in Cunninghamia lanceolate wood traits from different seed sources [57]. The coefficients of variation for middle- and low-altitude populations (Pop1, Pop2, Pop3) were comparable, yet significantly higher than the 10.27% for the high-altitude population Pop4. This could reflect the more complex and diverse habitats at middle and low altitudes, while the high-altitude region might offer more severe and uniform environmental conditions, constraining phenotypic variation [66].
To delve deeper into the contribution of the seven wood traits to overall variance, we conducted a principal component analysis (PCA). This analysis identified the key variables of Cunninghamia lanceolata wood traits at various altitudes and quantified their contributions to total variance. This approach offers valuable insights into the adaptation and growth characteristics of Cunninghamia lanceolata across different altitudinal ranges. The PCA results indicated that the first three principal components accounted for a cumulative contribution of 79.79% to the total variance. The order of contribution of the seven wood traits, from highest to lowest, was tracheid width, basic wood density, water absorption, tracheid length-to-width ratio, tracheid length, average annual ring width, and heartwood ratio. This analysis underscored the principal factors influencing wood traits [67], providing a deeper understanding of Cunninghamia lanceolata’s adaptation and growth at different altitudes [68] and aiding in the selection of superior germplasm resources [69].
Cluster analysis of Cunninghamia lanceolata wood traits across different altitudes revealed that the 98 samples were categorized into four taxa. Interestingly, these clusters did not align strictly with altitude, suggesting that while traits like heartwood ratio, tracheid length, and tracheid length-to-width ratio varied significantly among altitudes, these variations were not solely altitude-dependent. This indicates the potential influence of genetic factors or other environmental elements beyond altitude on wood trait formation. For instance, studies on wood traits in species such as the spruce Picea asperata [70], Japanese larch [71], and poplar [72] have demonstrated that traits like tracheid size and wood density are strongly genetically controlled. Additionally, environmental factors like temperature [73], precipitation [74], soil [75], and their interactions [76] also significantly impact wood traits. Therefore, the variation in Cunninghamia lanceolata wood traits could be attributed to a combination of factors, including but not limited to altitude.

4.2. Variation Patterns of Wood Traits and Ecological Adaptations of Cunninghamia lanceolata at Different Altitudes

Mountains, as unique ecosystems, exhibit high habitat heterogeneity and climatic diversity. Environmental factors such as temperature, light, moisture, and soil exhibit gradient variations with altitude [77]. This gradient provides a distinctive setting for investigating species’ genetic variation and their responses to climate change [78]. As altitude increases, plant growth and development are constrained by temperature and moisture, leading to notable changes in adaptive traits, which in turn influence the wood traits of trees [79,80].
In our Pearson correlation analysis of Cunninghamia lanceolata wood traits at different altitudes, we observed intricate relationships among wood traits and between wood traits and climatic factors at various altitudes. These relationships highlight the significance of quantitative traits in trees, where correlations between two traits might be due to strong linkage disequilibrium or pleiotropic effects of related genes [81]. Our findings revealed that the mean annual ring width was significantly negatively correlated with tracheid length and aspect ratio and highly significantly negatively correlated with wood density. Conversely, it showed highly significant positive correlations with water absorption. These results align with previous studies that underscore the interplay between wood traits and environmental influences [57,82]. Notably, the significant negative correlation between heartwood ratio and tracheid aspect ratio and the highly significant negative correlation with tracheid length suggest a robust connection between heartwood formation and the morphological characteristics of tracheids. Moreover, wood density exhibited a significant positive correlation with tracheid aspect ratio and a highly significant negative correlation with water absorption, underscoring the pivotal role of wood density in determining wood traits.
Further examination of the correlation between wood traits and climatic factors at different altitudes revealed that the heartwood ratio exhibited highly significant positive correlations with altitude, annual precipitation, and precipitation during the wettest and driest seasons and highly significant negative correlations with mean annual air temperature and mean air temperature during the wettest and driest seasons. This pattern suggests that heartwood formation is smaller in environments with higher temperatures and lower precipitation, reflecting the adaptive survival strategies of trees under varying hydrothermal conditions. These observations align with Almeida et al.’s findings [83] regarding the impact of climate on eucalyptus heartwood. Additionally, a significant positive correlation was observed between tracheid length, width, and aspect ratio and altitude and air temperature, indicating an increase in tracheid structural dimensions with rising altitude and air temperature. This trend might be attributed to the influence of climatic conditions on plant growth at higher altitudes [84], as suggested by Fonti et al. [85], where a warming climate is conducive to tracheid growth. Conversely, tracheid structure showed a significant negative correlation with annual precipitation and precipitation during the wettest and driest seasons, implying that excessive moisture may inhibit tracheid growth, highlighting the regulatory role of water conditions on plant cell structure. These insights corroborate previous studies [86,87], underscoring the critical role of environmental factors in shaping plant cell structure. Nonetheless, it is essential to consider the variability across different tree species, geographical locations, and environmental conditions [88,89].

4.3. Analysis of the Contribution Rate of Altitude Climatic Factors to the Variation of Wood Traits of Cunninghamia lanceolata

Random forest analysis results underscored the significant impact of altitude on the variability of wood traits like heartwood ratio and tracheid length, width, and aspect ratio, with respective contribution rates of 48.81%, 45.71%, 61.94%, and 60.92%. This evidence points to the pivotal role of altitude in shaping wood traits, particularly in middle and lower altitudinal gradients where it emerges as the dominant environmental influence. Additionally, the analysis revealed that temperature’s contribution to wood trait variability surpassed that of precipitation, aligning with Castagneri et al. [90] and Dang, Haishan et al. [34], who highlighted temperature’s limiting role in tree growth at higher altitudes.
As altitude increases, the influence of altitude factors on wood trait variation diminishes, while climatic factors become more pronounced. At higher altitudes, climatic factors’ cumulative contribution to tracheid length, width, and aspect ratio surpasses that of altitude, becoming the primary determinants of wood traits. This shift could be attributed to low temperatures and altered precipitation patterns at higher altitudes significantly impacting tree growth and wood traits. For instance, low temperatures may restrict cell division and elongation, influencing tracheid formation and size [91]. Changes in precipitation, on the other hand, can affect trees’ water use efficiency and growth rates [92], thereby altering wood quality and structure. These physiological responses indicate trees’ adaptive adjustments to their environment. At higher altitudes, temperature and precipitation’s limiting effects on tree growth become more evident, making wood traits more sensitive to these climatic factors [93]. Research has shown that at the upper limits of a species’ altitudinal distribution, low temperatures during the growing season can limit metabolic processes like photosynthesis fixation, subsequently affecting cell differentiation and xylem formation [22,90].

5. Conclusions

In this research, the phenotypic variation in wood traits of Cunninghamia lanceolata across diverse altitudinal gradients was meticulously analyzed. Our findings revealed significant differences in traits like heartwood ratio, tracheid length, and tracheid length-to-width ratio among various elevations, with mean annual ring width and tracheid width also exhibiting notable variations. This underscores the substantial variability in Cunninghamia lanceolata wood traits across different elevation zones. Notably, the coefficients of variation for these wood traits were markedly higher in populations at middle and lower elevations compared to those at higher elevations, suggesting a greater influence of diverse environmental conditions and enhanced genetic adaptability at these altitudes.
Principal component analysis identified key variables influencing wood trait variation, with tracheid width, basic wood density, and water absorption emerging as significant contributors to overall variability. These insights are crucial for comprehending the adaptation and growth dynamics of Cunninghamia lanceolata at varying altitudes and are instrumental in identifying superior germplasm resources. Cluster analysis further revealed that wood trait variation is influenced not only by altitude but also by an interplay of genetic and other environmental factors.
Pearson correlation analysis delineated intricate relationships among wood traits and between these traits and climatic factors at different altitudes. These relationships highlight the significance of quantitative traits in trees, where correlations between traits could be attributed to strong linkage disequilibrium or pleiotropic effects of associated genes. Random forest analysis underscored the pivotal role of elevation in influencing wood traits such as heartwood ratio, tracheid length, width, and aspect ratio. As elevation increases, the influence of elevation factors on wood trait variation diminishes, while climatic factors, particularly temperature, play an increasingly limiting role in tree growth at higher elevations.
These findings offer novel perspectives on tree adaptation to climate change and bear significant implications for forest genetic breeding and resource conservation. The pronounced variability observed in populations at lower and middle elevations presents a wealth of genetic resources for breeding programs. Concurrently, preserving high-elevation populations is essential for maintaining their unique adaptive traits. In the context of global climate change, these insights are vital for anticipating future shifts in forest ecosystems.
Future research should delve into the specific mechanisms through which different environmental factors, such as soil and light, impact wood traits in Cunninghamia lanceolata. Understanding how these variations influence the ecological functions and ecosystem services of these trees will enrich our knowledge of forest ecology and genetics, guiding future forest management and conservation strategies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/f15030411/s1, Table S1: Basic information of sampling population of different altitudes of Cunninghamia lanceolata; Table S2: Principal component analysis of wood traits of Cunninghamia lanceolata at different altitudes; Table S3: Contribution of variation in wood traits of Cunninghamia lanceolata at different elevations; Table S4: Correlation coefficients between wood traits of Cunninghamia lanceolata at different elevations and climatic factors at elevation; Table S5: Mantel test for wood traits of Cunninghamia lanceolata at different elevations in relation to climatic factors at elevation.

Author Contributions

Conceptualization, N.Z.; Methodology, G.X.; software, G.X.; investigation, G.X., S.L., T.C.; Data curation, G.X.; Writing—original draft, G.X.; Writing—review and editing, N.Z.; Funding acquisition, N.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a sub-project of the National Key Research and Development Program of China (2021YFD2201303) and a project of Science and Technology Innovation Fund for Postgraduate Students of Central South Forestry University of Science and Technology (2023CX02052).

Data Availability Statement

Data will be made available on request.

Acknowledgments

We thank Yamei Gong for her help and support in the field sampling.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Cui, C.W.; Peng, L.H.; Ma, D.X.; Wang, J.Q.; Jiang, X.Q.; Jiang, X.G.; Ma, X.Q.; Lin, K. Effects of Thinning on Soil Microbial Necromass Carbon in Cunninghamia lanceolata Plantation. Sci. Silvae Sin. 2023, 59, 41–52. [Google Scholar] [CrossRef]
  2. Li, Y.X. Genetic Diversity and Genetic Differentiation of Geographic Seed Sources of Cunninghamia lanceolata. Master’s Thesis, Chinese Academy of Forestry Sciences, Beijing, China, 2015. [Google Scholar]
  3. Li, Y.; Li, M.; Li, C.; Liu, Z. Optimized Maxent Model Predictions of Climate Change Impacts on the Suitable Distribution of Cunninghamia lanceolata in China. Forests 2020, 11, 302. [Google Scholar] [CrossRef]
  4. Zhang, X.; Peng, J.Y.; Shi, J.T.; Xu, Q.Q.; Xu, Z.K.; Huang, F. Wood Anatomical Characteristics and Physical—Mechanical Properties of Dark—Brown Heart Cunningham lanceolate from Hunan. J. Southwest For. Univ. (Nat. Sci.) 2021, 41, 155–160. [Google Scholar]
  5. You, R.; Zhu, N.; Deng, X.; Wang, J.; Liu, F. Variation in Wood Physical Properties and Effects of Climate for Different Geographic Sources of Chinese Fir in Subtropical Area of China. Sci. Rep. 2021, 11, 4664. [Google Scholar] [CrossRef] [PubMed]
  6. Zhou, J.; Yu, C.; Cheng, S.D.; Xu, Q.L.; Ma, X.Q.; Li, M. Comparison of Early Growth and Lateral Branch Development in Characteristics of Cunninghamia lanceolata Clones. J. For. Environ. 2023, 43, 561–568. [Google Scholar] [CrossRef]
  7. Jiang, Y.; Wu, J.Y.; Yao, X.Q.; Lin, Q.; Long, Y.Z. Research and Development Countermeasures on Precious Indigenous Timber Tree Species in Hunan Province. Hunan For. Sci. Technol. 2012, 40, 1–4. [Google Scholar] [CrossRef]
  8. Wang, J. Genetic and Variations of Black Heartwood of Cunninghamia lanceolata. Master’s Thesis, Central South Forestry University of Science and Technology, Changsha, China, 2018. [Google Scholar] [CrossRef]
  9. Hunag, R.C.; Zhu, N.H.; Yang, J.; Yang, X.W. Research on Cone and Seed Qualities among Different Families of Black-heart Wood Cunninghamia lanceolata. Hunan For. Sci. Technol. 2021, 48, 40–44. [Google Scholar] [CrossRef]
  10. Yang, X.W. Construction of Parent Population for Breeding Programme in the Black-Heart Wood Cunninghamia lanceolata. Master’s Thesis, Central South Forestry University of Science and Technology, Changsha, China, 2021. [Google Scholar] [CrossRef]
  11. Yang, X.W.; Zhu, N.H.; Han, Z.Q.; Liu, T.Y.; Li, H.; Yang, J.; Zhu, X.C. Identification of Paternal Parents of Half-sib Progeny Superior Plants of the Black-heart Wood Cunninghamia lanceolata. J. Cent. South Univ. For. Technol. 2021, 41, 26–35. [Google Scholar] [CrossRef]
  12. Zhu, N.; Yang, X.; Han, Z.; Can, X. Research Progress on Cunninghamia lanceolata. In Conifers—Recent Advances; Cristina Gonçalves, A., Fonseca, T., Eds.; IntechOpen: Chongqing, China, 2022. [Google Scholar] [CrossRef]
  13. Hu, S.R.; Xia, M.; Guo, C.Y.; Lu, X.C. Overview of Research Methods on Genetic Diversity of Forest Trees. J. Northeast For. Univ. 2001, 29, 72–75. [Google Scholar]
  14. Wen, Y.F.; Han, W.J.; Wu, S. Plant Genetic Diversity and its Influencing Factors. J. Cent. South Univ. For. Technol. 2010, 30, 80–87. [Google Scholar]
  15. Hughes, A.R.; Inouye, B.D.; Johnson, M.T.J.; Underwood, N.; Vellend, M. Ecological Consequences of Genetic Diversity. Ecol. Lett. 2008, 11, 609–623. [Google Scholar] [CrossRef]
  16. Ma, D.C.; Wang, L.H.; Liang, J. Research Progress on Application of Morphological Markers in Plants. Jiangsu Agric. Sci. 2022, 50, 55–62. [Google Scholar] [CrossRef]
  17. Alcántara-Ayala, O.; Oyama, K.; Ríos-Muñoz, C.A.; Rivas, G.; Ramirez-Barahona, S.; Luna-Vega, I. Morphological Variation of Leaf Traits in the Ternstroemia lineata Species Complex (Ericales: Penthaphylacaceae) in Response to Geographic and Climatic Variation. PeerJ 2020, 8, e8307. [Google Scholar] [CrossRef]
  18. Albarrán-Lara, A.L.; Petit, R.J.; Kremer, A.; Caron, H.; Peñaloza-Ramírez, J.M.; Gugger, P.F.; Dávila-Aranda, P.D.; Oyama, K. Low Genetic Differentiation between Two Morphologically and Ecologically Distinct Giant-Leaved Mexican Oaks. Plant Syst. Evol. 2019, 305, 89–101. [Google Scholar] [CrossRef]
  19. Bai, T.D.; Yu, C.L.; Gan, Z.C.; Lai, H.R.; Yang, Y.C.; Huang, H.C.; Jiang, W.X. Association of Cone and Seed Traits of Pinus yunnanensis var. tenuifolia with Geo-meteorological Factors. Chin. J. Plant Ecol. 2020, 44, 1224–1235. [Google Scholar] [CrossRef]
  20. Zhang, T.J.; Chen, X.H.; Kang, X.K.; Liu, J. Phenotypic Diversity of Leaf Morphologic Traits of Davidia involucrata Natural Populations in Sichuan Province. Chin. J. Ecol. 2019, 38, 35–43. [Google Scholar] [CrossRef]
  21. Yan, Y.; Li, B.J.; He, Z.M.; Liu, Y.H.; Li, M.; Yu, X.T.; Ma, X.Q. Comparative Research of Wood Properties of Cunninghamia lanceolate from Different Seed Sources. For. Res. 2021, 34, 49–57. [Google Scholar] [CrossRef]
  22. Petit, G.; Anfodillo, T.; Carraro, V.; Grani, F.; Carrer, M. Hydraulic Constraints Limit Height Growth in Trees at High Altitude. New Phytol. 2011, 189, 241–252. [Google Scholar] [CrossRef]
  23. Gebauer, R.; Volařík, D.; Urban, J.; Børja, I.; Nagy, N.E.; Eldhuset, T.D.; Krokene, P. Effects of Different Light Conditions on the Xylem Structure of Norway Spruce Needles. Trees 2012, 26, 1079–1089. [Google Scholar] [CrossRef]
  24. D’Orangeville, L.; Côté, B.; Houle, D.; Morin, H.; Duchesne, L. A Three-Year Increase in Soil Temperature and Atmospheric N Deposition Has Minor Effects on the Xylogenesis of Mature Balsam Fir. Trees 2013, 27, 1525–1536. [Google Scholar] [CrossRef]
  25. Dufour, B.; Morin, H. Climatic Control of Tracheid Production of Black Spruce in Dense Mesic Stands of Eastern Canada. Tree Physiol. 2013, 33, 175–186. [Google Scholar] [CrossRef] [PubMed]
  26. Arenas-Navarro, M.; Oyama, K.; García-Oliva, F.; Torres-Miranda, A.; De La Riva, E.G.; Terrazas, T. The Role of Wood Anatomical Traits in the Coexistence of Oak Species along an Environmental Gradient. AoB Plants 2021, 13, plab066. [Google Scholar] [CrossRef] [PubMed]
  27. You, Q.; Jiang, Z.; Yue, X.; Guo, W.; Liu, Y.; Cao, J.; Li, W.; Wu, F.; Cai, Z.; Zhu, H.; et al. Recent Frontiers of Climate Changes in East Asia at Global Warming of 1.5 °C and 2 °C. npj Clim. Atmos. Sci. 2022, 5, 80. [Google Scholar] [CrossRef]
  28. Mary Petritan, A.; Beloiu Schwenke, M. Forest Functioning under Climate Warming and Future Perspectives on Forest Disturbances. Forests 2023, 14, 2302. [Google Scholar] [CrossRef]
  29. Harrington, C.A.; Gould, P.J.; Cronn, R. Site and Provenance Interact to Influence Seasonal Diameter Growth of Pseudotsuga menziesii. Front. For. Glob. Chang. 2023, 6, 1173707. [Google Scholar] [CrossRef]
  30. Harrington, C.A.; Dang, Q.-L.; Man, R.; Inoue, S.; Tedla, B. Editorial: Changing Seasons: How Is Global Warming Affecting Forest Phenology? Front. For. Glob. Chang. 2023, 6, 1257096. [Google Scholar] [CrossRef]
  31. Jiang, A.P.; Jiang, J.M.; Liu, J. Responses of Leaf Traits of Sassafras tsumu (Hemsl.) Hemsl along an Altitudinal Gradient. J. Ecol. 2016, 35, 1467–1474. [Google Scholar] [CrossRef]
  32. Shen, D.; Bo, W.; Xu, F.; Wu, R. Genetic Diversity and Population Structure of the Tibetan Poplar (Populus szechuanica Var. tibetica) along an Altitude Gradient. BMC Genet. 2014, 15, S11. [Google Scholar] [CrossRef]
  33. Reisch, C.; Rosbakh, S. Patterns of Genetic Variation in European Plant Species Depend on Altitude. Divers. Distrib. 2021, 27, 157–163. [Google Scholar] [CrossRef]
  34. Dang, H.; Zhang, Y.; Zhang, K.; Jiang, M.; Zhang, Q. Climate-Growth Relationships of Subalpine Fir (Abies fargesii) across the Altitudinal Range in the Shennongjia Mountains, Central China. Clim. Chang. 2013, 117, 903–917. [Google Scholar] [CrossRef]
  35. Zhang, L.; Jiang, Y.; Zhao, S.; Kang, X.; Zhang, W.; Liu, T. Lingering Response of Radial Growth of Picea Crassifolia to Climate at Different Altitudes in the Qilian Mountains, Northwest China. Trees 2017, 31, 455–465. [Google Scholar] [CrossRef]
  36. Carrer, M.; Castagneri, D.; Prendin, A.L.; Petit, G.; Von Arx, G. Retrospective Analysis of Wood Anatomical Traits Reveals a Recent Extension in Tree Cambial Activity in Two High-Elevation Conifers. Front. Plant Sci. 2017, 8, 737. [Google Scholar] [CrossRef]
  37. Keleş, S.Ö. The Effect of Altitude on the Growth and Development of Trojan fir (Abies nordmanniana subsp. equi-trojani [Asch. & Sint. ex Boiss] Coode & Cullen) Saplings. Cerne 2020, 26, 381–392. [Google Scholar] [CrossRef]
  38. Cocozza, C.; Palombo, C.; Tognetti, R.; La Porta, N.; Anichini, M.; Giovannelli, A.; Emiliani, G. Monitoring Intra-Annual Dynamics of Wood Formation with Microcores and Dendrometers in Picea Abies at Two Different Altitudes. Tree Physiol. 2016, 36, 832–846. [Google Scholar] [CrossRef] [PubMed]
  39. Cai, Q.; Liu, Y.; Qian, H.; Liu, R. Inverse Effects of Recent Warming on Trees Growing at the Low and High Altitudes of the Dabie Mountains, Subtropical China. Dendrochronologia 2020, 59, 125649. [Google Scholar] [CrossRef]
  40. Zhuang, L.; Axmacher, J.C.; Sang, W. Different Radial Growth Responses to Climate Warming by Two Dominant Tree Species at Their Upper Altitudinal Limit on Changbai Mountain. J. For. Res. 2017, 28, 795–804. [Google Scholar] [CrossRef]
  41. Plant Intelligence. Available online: https://www.iplant.cn/info/%E6%9D%89%E6%9C%A8 (accessed on 25 December 2023).
  42. Gong, H.; Yang, M.; Wang, C.; Tian, C. Leaf Phenotypic Variation and Its Response to Environmental Factors in Natural Populations of Eucommia Ulmoides. BMC Plant Biol. 2023, 23, 562. [Google Scholar] [CrossRef]
  43. Xu, M.; Ren, H.Q.; Guo, W. Research Progress of Variation for Wood Properties Cunninghamia lanceolata. J. Northwest Coll. For. 2008, 185–189. [Google Scholar]
  44. Wimmer, R. Wood Anatomical Features in Tree-Rings as Indicators of Environmental Change. Dendrochronologia 2002, 20, 21–36. [Google Scholar] [CrossRef]
  45. Nagavciuc, V.; Mursa, A.; Ionita, M.; Sfeclă, V.; Popa, I.; Roibu, C.-C. An Overview of Extreme Years in Quercus Sp. Tree Ring Records from the Northern Moldavian Plateau. Forests 2023, 14, 894. [Google Scholar] [CrossRef]
  46. Iglesias, C.; Santos, A.; Martínez, J.; Pereira, H.; Anjos, O. Influence of Heartwood on Wood Density and Pulp Properties Explained by Machine Learning Techniques. Forests 2017, 8, 20. [Google Scholar] [CrossRef]
  47. Thybring, E.E.; Fredriksson, M. Wood Modification as a Tool to Understand Moisture in Wood. Forests 2021, 12, 372. [Google Scholar] [CrossRef]
  48. Arsić, J.; Stojanović, M.; Petrovičová, L.; Noyer, E.; Milanović, S.; Světlík, J.; Horáček, P.; Krejza, J. Increased Wood Biomass Growth Is Associated with Lower Wood Density in Quercus petraea (Matt.) Liebl. Saplings Growing under Elevated CO2. PLoS ONE 2021, 16, e0259054. [Google Scholar] [CrossRef]
  49. McDowell, N.G.; Allen, C.D.; Anderson-Teixeira, K.; Aukema, B.H.; Bond-Lamberty, B.; Chini, L.; Clark, J.S.; Dietze, M.; Grossiord, C.; Hanbury-Brown, A.; et al. Pervasive Shifts in Forest Dynamics in a Changing World. Science 2020, 368, eaaz9463. [Google Scholar] [CrossRef] [PubMed]
  50. Ziemińska, K.; Butler, D.W.; Gleason, S.M.; Wright, I.J.; Westoby, M. Fibre Wall and Lumen Fractions Drive Wood Density Variation across 24 Australian Angiosperms. AoB Plants 2013, 5, plt046. [Google Scholar] [CrossRef]
  51. Yeboah, D.; Burton, A.J.; Storer, A.J.; Opuni-Frimpong, E. Variation in Wood Density and Carbon Content of Tropical Plantation Tree Species from Ghana. New For. 2014, 45, 35–52. [Google Scholar] [CrossRef]
  52. Bouriaud, O.; Leban, J.-M.; Bert, D.; Deleuze, C. Intra-Annual Variations in Climate Influence Growth and Wood Density of Norway Spruce. Tree Physiol. 2005, 25, 651–660. [Google Scholar] [CrossRef] [PubMed]
  53. Dahlen, J.; Nabavi, M.; Auty, D.; Schimleck, L.; Eberhardt, T.L. Models for Predicting the Within-Tree and Regional Variation of Tracheid Length and Width for Plantation Loblolly Pine. For. Int. J. For. Res. 2021, 94, 127–140. [Google Scholar] [CrossRef]
  54. Palla, B.; Ladányi, M.; Cseke, K.; Buczkó, K.; Höhn, M. Wood Anatomical Traits Reveal Different Structure of Peat Bog and Lowland Populations of Pinus sylvestris L. in the Carpathian Region. Forests 2021, 12, 494. [Google Scholar] [CrossRef]
  55. Xu, L.; Chen, G.X.; Zhang, D.G.; Li, X.T.; Zhang, Y. Flora of Seed Plants in Xiaoxi National Nature Reserve, Hunan, China. Acta Bot. Boreali-Occident. Sin. 2010, 30, 2307–2316. [Google Scholar]
  56. GB/T 1927314-2022; State Administration of Market Supervision, National Standardization Administration. Test Method for Physical and Mechanical Properties of Wood in Small Specimens without Defects Part 3: Determination of Growth Wheel Width and Latewood Rate. China Standard Press: Beijing, China, 2022.
  57. Duan, H.J.; Cao, S.; Zheng, H.; Hu, D.; Lin, J.; Lin, H.; Hu, R.; Sun, Y.; Li, Y. Variation in the Growth Traits and Wood Properties of Cunninghamia lanceolata from Six Provinces of Southern China. Forests 2016, 7, 192. [Google Scholar] [CrossRef]
  58. Xu, Y.L. Genetic Variation of Natural Population in Pinus yunnanensis Franch. Ph.D. Thesis, Beijing Forestry University, Beijing, China, 2015. [Google Scholar]
  59. Jiao, S.; Chen, W.; Wang, J.; Du, N.; Li, Q.; Wei, G. Soil Microbiomes with Distinct Assemblies through Vertical Soil Profiles Drive the Cycling of Multiple Nutrients in Reforested Ecosystems. Microbiome 2018, 6, 146. [Google Scholar] [CrossRef] [PubMed]
  60. Topaloğlu, E.; Ay, N.; Altun, L.; Serdar, B. Effect of Altitude and Aspect on Various Wood Properties of Oriental Beech (Fagus orientalis Lipsky) Wood. Turk. J. Agric. For. 2016, 40, 397–406. [Google Scholar] [CrossRef]
  61. Rocha, S.M.G.; Vidaurre, G.B.; Pezzopane, J.E.M.; Almeida, M.N.F.; Carneiro, R.L.; Campoe, O.C.; Scolforo, H.F.; Alvares, C.A.; Neves, J.C.L.; Xavier, A.C.; et al. Influence of Climatic Variations on Production, Biomass and Density of Wood in Eucalyptus Clones of Different Species. For. Ecol. Manag. 2020, 473, 118290. [Google Scholar] [CrossRef]
  62. Audigeos, D.; Buonamici, A.; Belkadi, L.; Rymer, P.; Boshier, D.; Scotti-Saintagne, C.; Vendramin, G.G.; Scotti, I. Aesqeaurchaaprtioclerins in the Wild: Natural Genetic Diversity and Selective Pressure in the PIP Gene Family in Five Neotropical Tree Species. BMC Evol. Biol. 2010, 10, 202. [Google Scholar] [CrossRef]
  63. Huang, R.; Zeng, W.; Deng, H.; Hu, D.; Wang, R.; Zheng, H. Inbreeding in Cunninghamia lanceolata: Insight into the Rare Self-Fertilizing Event from a Genetic View. Genes 2022, 13, 2105. [Google Scholar] [CrossRef]
  64. Silva-Montellano, A.; Eguiarte, L.E. Geographic Patterns in the Reproductive Ecology of Agave Lechuguilla (Agavaceae) in the Chihuahuan Desert. II. Genetic Variation, Differentiation, and Inbreeding Estimates. Am. J. Bot. 2003, 90, 700–706. [Google Scholar] [CrossRef]
  65. Jin, Z.X.; Li, J.M. Analysis of the Genetic Diversity of Hepatacodium miconioides at Different Altitude and in Tiantai Mountain in Zhejiang Province and its Relationship with Environmental Factors. J. Zhejiang Univ. (Sci. Ed.) 2005, 452–458. [Google Scholar] [CrossRef]
  66. Schellenberger Costa, D.; Gerschlauer, F.; Kiese, R.; Fischer, M.; Kleyer, M.; Hemp, A. Plant Niche Breadths along Environmental Gradients and Their Relationship to Plant Functional Traits. Divers. Distrib. 2018, 24, 1869–1882. [Google Scholar] [CrossRef]
  67. Zheng, G.B.; Ma, L.; Zhang, Z.M.; Li, M.L.; Liu, S.Y.; Zang, Y.N.; Sun, Y.; Ma, X.Y. Principal Component Analysis and Comprehensive Evaluation of Ningxia Wolfberry Fruits under Water Stress Conditions at Different Fertility Periods. Water Sav. Irrig. 2022, 47–57. [Google Scholar]
  68. Liang, K.L.; Jiang, W.Q.; Zhou, Z.Y.; Guo, X.; Li, X.Z.; Dai, W.A.; Wang, R.; Liu, X.Y. Variation in Main Morphological Characteristics of Amorpha fruticosa Plants in the Qinghai-Tibet Plateau. Acta Ecol. Sin. 2012, 32, 311–318. [Google Scholar] [CrossRef]
  69. Jia, Q.B.; Zhang, H.G.; Zhang, L. Variation in Carbon Conten of Hybrid larch Families and Superior Families Selection. J. Northeast For. Univ. 2016, 44, 1–6. [Google Scholar]
  70. Luo, J.X.; Li, X.Q.; Sun, P.; Wang, L.H.; Hu, G.R.; Wang, F.L.; Zheng, W. Variation Pattern of Tracheid and Wood Density in Natural Population of Picea asperata. J. Beijing For. Univ. 2004, 80–85. [Google Scholar]
  71. Ma, S.X.; Wang, J.H.; Zhang, S.T.; Sun, X.M.; Liang, B.S.; Zhou, D.Y.; Liu, S.M. Genetic Variation of Wood Properties in Japanese larch Clones. For. Res. 2008, 69–73. [Google Scholar]
  72. Yu, H.Y.; Pang, Z.Y.; Yin, C.H.; Ding, C.J.; Wang, F.S.; Zhang, J.W.; Wang, F.L.; Jin, P.L.; Fu, Z.X.; Qu, G.Z.; et al. Variation Analysis on the Growth and Wood Properties Variation of 100 Poplar Clones. J. Northeast For. Univ. 2023. [Google Scholar]
  73. Ramananantoandro, T.; Ramanakoto, M.F.; Rajoelison, G.L.; Randriamboavonjy, J.C.; Rafidimanantsoa, H.P. Influence of Tree Species, Tree Diameter and Soil Types on Wood Density and Its Radial Variation in a Mid-Altitude Rainforest in Madagascar. Ann. For. Sci. 2016, 73, 1113–1124. [Google Scholar] [CrossRef]
  74. Camarero, J.J.; Fernández-Pérez, L.; Kirdyanov, A.V.; Shestakova, T.A.; Knorre, A.A.; Kukarskih, V.V.; Voltas, J. Minimum Wood Density of Conifers Portrays Changes in Early Season Precipitation at Dry and Cold Eurasian Regions. Trees 2017, 31, 1423–1437. [Google Scholar] [CrossRef]
  75. Fang, S.; Sun, D.; Shang, X.; Fu, X.; Yang, W. Variation in Radial Growth and Wood Density of Cyclocarya Paliurus across Its Natural Distribution. New For. 2020, 51, 453–467. [Google Scholar] [CrossRef]
  76. Zhang, H.; Zhang, S.; Chen, S.; Xia, D.; Yang, C.; Zhao, X. Genetic Variation and Superior Provenances Selection for Wood Properties of Larix Olgensis at Four Trials. J. For. Res. 2022, 33, 1867–1879. [Google Scholar] [CrossRef]
  77. Jiang, A.P.; Jiang, J.M.; Liu, J. Genetic Diversity and Genetic Structure of Sassafras tsumu Populations along Altutudinal Gradients in Tianmushan Mountain, China. Chin. J. Appl. Ecol. 2016, 27, 1829–1836. [Google Scholar] [CrossRef]
  78. Zu, K.; Wang, Z. Research progress on the elevational distribution of mountain species in response to climate change. Biodivers. Sci. 2022, 30, 123–137. [Google Scholar] [CrossRef]
  79. Pan, H.L.; Li, M.H.; Cai, S.H.; Wu, J.; Du, Z.; Liu, X.L. Responses of Growth and Ecological of Plant to Altitude. Ecol. Environ. Sci. 2009, 18, 722–730. [Google Scholar] [CrossRef]
  80. Fajardo, A.; Piper, F.I.; García-Cervigón, A.I. The Intraspecific Relationship between Wood Density, Vessel Diameter and Other Traits across Environmental Gradients. Funct. Ecol. 2022, 36, 1585–1598. [Google Scholar] [CrossRef]
  81. Du, Q.; Gong, C.; Wang, Q.; Zhou, D.; Yang, H.; Pan, W.; Li, B.; Zhang, D. Genetic Architecture of Growth Traits in Populus Revealed by Integrated Quantitative Trait Locus (QTL) Analysis and Association Studies. New Phytol. 2016, 209, 1067–1082. [Google Scholar] [CrossRef] [PubMed]
  82. Bouslimi, B.; Koubaa, A.; Bergeron, Y. Intra-Ring Variations and Interrelationships for Selected Wood Anatomical and Physical Properties of Thuja occidentalis L. Forests 2019, 10, 339. [Google Scholar] [CrossRef]
  83. Almeida, M.N.F.D.; Vidaurre, G.B.; Pezzopane, J.E.M.; Lousada, J.L.P.C.; Silva, M.E.C.M.; Câmara, A.P.; Rocha, S.M.G.; Oliveira, J.C.L.D.; Campoe, O.C.; Carneiro, R.L.; et al. Heartwood Variation of Eucalyptus Urophylla Is Influenced by Climatic Conditions. For. Ecol. Manag. 2020, 458, 117743. [Google Scholar] [CrossRef]
  84. Olano, J.M.; Almería, I.; Eugenio, M.; Von Arx, G. Under Pressure: How a M Editerranean High-mountain Forb Coordinates Growth and Hydraulic Xylem Anatomy in Response to Temperature and Water Constraints. Funct. Ecol. 2013, 27, 1295–1303. [Google Scholar] [CrossRef]
  85. Fonti, P.; Bryukhanova, M.V.; Myglan, V.S.; Kirdyanov, A.V.; Naumova, O.V.; Vaganov, E.A. Temperature-induced Responses of Xylem Structure of Larix sibirica (Pinaceae) from the Russian Altay. Am. J. Bot. 2013, 100, 1332–1343. [Google Scholar] [CrossRef] [PubMed]
  86. Zhao, X.P.; Guo, M.H. Response of Wood Tracheid Morphological Features of Pinus koraiensis Plantation to Climate Change. J. Northeast For. Univ. 2009, 37, 1–4. [Google Scholar]
  87. Zhu, L.J.; Li, Z.S.; Wang, X.C. Anatomical Characteristics of Xylem in Tree Rings and its Relationship with Environmental. Chin. J. Plant Ecol. 2017, 41, 238–251. [Google Scholar] [CrossRef]
  88. Xu, J.M.; Zhang, R.; Lv, J.X.; Robert, E. Climatic Response in of Picea crassifolia along Elevation Gradient in Qilian Mountains, Northwestern China. J. Beijing For. Univ. 2015, 37, 102–108. [Google Scholar] [CrossRef]
  89. Xu, J.M.; Lu, J.X.; Bao, F.C.; Huang, R.F.; Zhao, Y.K. Influence of Climatic Factors on Wood Cell Structure. For. Sci. 2011, 47, 151–158. [Google Scholar]
  90. Castagneri, D.; Petit, G.; Carrer, M. Divergent Climate Response on Hydraulic-Related Xylem Anatomical Traits of Picea abies along a 900-m Altitudinal Gradient. Tree Physiol. 2015, 35, 1378–1387. [Google Scholar] [CrossRef]
  91. Nagelmüller, S.; Hiltbrunner, E.; Körner, C. Low Temperature Limits for Root Growth in Alpine Species Are Set by Cell Differentiation. AoB Plants 2017, 9, plx054. [Google Scholar] [CrossRef] [PubMed]
  92. Heilman, K.A.; Trouet, V.M.; Belmecheri, S.; Pederson, N.; Berke, M.A.; McLachlan, J.S. Increased Water Use Efficiency Leads to Decreased Precipitation Sensitivity of Tree Growth, but Is Offset by High Temperatures. Oecologia 2021, 197, 1095–1110. [Google Scholar] [CrossRef] [PubMed]
  93. Maxime, C.; Hendrik, D. Effects of Climate on Diameter Growth of Co-Occurring Fagus sylvatica and Abies alba along an Altitudinal Gradient. Trees 2011, 25, 265–276. [Google Scholar] [CrossRef]
Figure 1. Schematic map of the sampling research area of Cunninghamia lanceolata at different altitudes. (a) Administrative map of China; (b) Administrative map of Hunan; (c) Tujia–Miao autonomous continent of Xiangxi Prefecture; (d) Xiaoxi National Nature Reserve.
Figure 1. Schematic map of the sampling research area of Cunninghamia lanceolata at different altitudes. (a) Administrative map of China; (b) Administrative map of Hunan; (c) Tujia–Miao autonomous continent of Xiangxi Prefecture; (d) Xiaoxi National Nature Reserve.
Forests 15 00411 g001
Figure 2. Schematic diagrams of Cunninghamia lanceolata cores and tracheids at different altitudes; 550 m, 750 m, 950 m, and 1150 m from Pop1, Pop2, Pop3, and Pop4 populations, respectively, representing the average altitude sampled, (a) Cunninghamia lanceolata cores collected from different altitudes; (b) tracheids under the bio-digital microscope.
Figure 2. Schematic diagrams of Cunninghamia lanceolata cores and tracheids at different altitudes; 550 m, 750 m, 950 m, and 1150 m from Pop1, Pop2, Pop3, and Pop4 populations, respectively, representing the average altitude sampled, (a) Cunninghamia lanceolata cores collected from different altitudes; (b) tracheids under the bio-digital microscope.
Forests 15 00411 g002
Figure 3. (a) Fragmentation plot; (b) Principal component analysis 3D plot. Pop1, Pop2, Pop2, Pop4 represent four Cunninghamia lanceolata populations at different elevations.
Figure 3. (a) Fragmentation plot; (b) Principal component analysis 3D plot. Pop1, Pop2, Pop2, Pop4 represent four Cunninghamia lanceolata populations at different elevations.
Forests 15 00411 g003
Figure 4. Cluster analysis diagram. L1–L30 represent 30 samples from Pop1; LH1–LH30 represent 30 samples from Pop2; MH1–MH30 represent 30 samples from Pop1; H1–H8 represent 8 samples from Pop4.
Figure 4. Cluster analysis diagram. L1–L30 represent 30 samples from Pop1; LH1–LH30 represent 30 samples from Pop2; MH1–MH30 represent 30 samples from Pop1; H1–H8 represent 8 samples from Pop4.
Forests 15 00411 g004
Figure 5. Correlation analysis and Mantel test between altitude climatic factors and wood traits. * p < 0.05; ** p < 0.01; *** p < 0.001. BIO1, Mean Annual Temperature; BIO8, Mean Wettest Season Temperature; BIO9, Mean Driest Season Temperature; BIO12, Annual Precipitation; BIO16, Wettest Season Precipitation; BIO17, Driest Season Precipitation. Rb, Average Annual Ring Width; P, Heartwood Ratio; WBD, Basic Wood Density; Hy, Water Absorption Rate; L, Tracheid Length; D, Tracheid Width; L/D, Tracheid Length-to-Width Ratio.
Figure 5. Correlation analysis and Mantel test between altitude climatic factors and wood traits. * p < 0.05; ** p < 0.01; *** p < 0.001. BIO1, Mean Annual Temperature; BIO8, Mean Wettest Season Temperature; BIO9, Mean Driest Season Temperature; BIO12, Annual Precipitation; BIO16, Wettest Season Precipitation; BIO17, Driest Season Precipitation. Rb, Average Annual Ring Width; P, Heartwood Ratio; WBD, Basic Wood Density; Hy, Water Absorption Rate; L, Tracheid Length; D, Tracheid Width; L/D, Tracheid Length-to-Width Ratio.
Forests 15 00411 g005
Figure 6. (a) Stacked diagram of the sources of variation in wood traits. (b) Stacked diagram of the sources of variation in wood traits at different altitudes. Random forest analysis. Values represent the contribution of altitudinal climatic factors to variation in wood traits. BIO1, Mean Annual Temperature; BIO8, Mean Wettest Season Temperature; BIO9, Mean Driest Season Temperature; BIO12, Annual Precipitation; BIO16, Wetter Season Precipitation; BIO17, Driest Season Precipitation; P, Heartwood Ratio; L, Tracheid Length; D, Tracheid Width; L/D, Tracheid Length-to-Width Ratio; 550, 750, 950, 1150 for different altitudes.
Figure 6. (a) Stacked diagram of the sources of variation in wood traits. (b) Stacked diagram of the sources of variation in wood traits at different altitudes. Random forest analysis. Values represent the contribution of altitudinal climatic factors to variation in wood traits. BIO1, Mean Annual Temperature; BIO8, Mean Wettest Season Temperature; BIO9, Mean Driest Season Temperature; BIO12, Annual Precipitation; BIO16, Wetter Season Precipitation; BIO17, Driest Season Precipitation; P, Heartwood Ratio; L, Tracheid Length; D, Tracheid Width; L/D, Tracheid Length-to-Width Ratio; 550, 750, 950, 1150 for different altitudes.
Forests 15 00411 g006
Table 1. Nested ANOVA of wood property traits of natural populations of Cunninghamia lanceolata at different altitudes.
Table 1. Nested ANOVA of wood property traits of natural populations of Cunninghamia lanceolata at different altitudes.
Wood Property TraitsMean SquareF
Among PopulationsWithin PopulationsRandom ErrorAmong PopulationsWithin Populations
Rb, mm207.642453.267267.82183.0616 *0.7854
P1031.1240102.794563.998716.1116 ***1.0662
WBD, gm30.00230.00140.00240.96220.5977
Hy, %763.9594492.2978820.66570.93090.5999
L, um1,076,974.157268,727.261564,557.949016.6823 ***1.0646
D, um59.069821.006420.06512.9439 *1.0469
L/D385.573344.534178.57514.9071 ***0.5668
Note: * p < 0.05; *** p > 0.001. Rb, Average Annual Ring Width; P, Heartwood Ratio; WBD, Basic Wood Density; Hy, Water Absorption Rate; L, Tracheid Length; D, Tracheid Width; L/D, Tracheid Length-to-Width Ratio.
Table 2. Variations of wood property traits of natural populations of Cunninghamia lanceolata at different altitudes.
Table 2. Variations of wood property traits of natural populations of Cunninghamia lanceolata at different altitudes.
Wood Property TraitsPopulation IPopulation Pop2Population Pop3Population Pop4Among Populations
Mean ± SDCV (%)Mean ± SDCV (%)Mean ± SDCV (%)Mean ± SDCV (%)Mean ± SDCV (%)
Rb, mm24.35 ± 6.66b27.3528.01 ± 5.88ab20.9930.51 ± 11.1a36.3829.21 ± 3.61ab12.3628.02 ± 6.8124.31
P53.95 ± 7.58c14.0457.963 ± 8.99c15.5163.58 ± 9.64b15.1577.51 ± 7.95a10.2563.05 ± 8.5413.50
WBD, g/cm30.42 ± 0.05a11.900.42 ± 0.05a11.900.41 ± 0.04a9.760.39 ± 0.04a10.260.41 ± 0.0510.98
Hy, %174.97 ± 27.86a15.92177.64 ± 28.43a16.00180.61 ± 24.25a13.43191.90 ± 25.67a13.38181.28 ± 26.5514.65
L, um2888.06 ± 266.18a9.222589.22 ± 215.29b8.312503.64 ± 296.24b11.832283.09 ± 187.22c8.202629.51 ± 241.239.40
D, um45.22 ± 4.54a10.0444.65 ± 4.92a11.0243.21 ± 4.32ab10.0040.99 ± 3.2b7.8143.52 ± 4.259.75
L/D67.45 ± 10.06a14.9161.08 ± 9b14.7359.14 ± 5.56b9.4058.20 ± 5.60b9.6261.47 ± 7.5612.29
average value-14.77-14.07-15.14-10.27-13.55
Note: SD, standard deviation; CV, coefficient of variation. The underline represents the maximum value and the letters a, b, and c represent the results of multiple comparisons. Rb, Average Annual Ring Width; P, Heartwood Ratio; WBD, Basic Wood Density; Hy, Water Absorption Rate; L, Tracheid Length; D, Tracheid Width; L/D, Tracheid Length-to-Width Ratio.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xie, G.; Liu, S.; Chang, T.; Zhu, N. Forest Adaptation to Climate Change: Altitudinal Response and Wood Variation in Natural-Growth Cunninghamia lanceolata in the Context of Climate Change. Forests 2024, 15, 411. https://doi.org/10.3390/f15030411

AMA Style

Xie G, Liu S, Chang T, Zhu N. Forest Adaptation to Climate Change: Altitudinal Response and Wood Variation in Natural-Growth Cunninghamia lanceolata in the Context of Climate Change. Forests. 2024; 15(3):411. https://doi.org/10.3390/f15030411

Chicago/Turabian Style

Xie, Gongliang, Sen Liu, Ting Chang, and Ninghua Zhu. 2024. "Forest Adaptation to Climate Change: Altitudinal Response and Wood Variation in Natural-Growth Cunninghamia lanceolata in the Context of Climate Change" Forests 15, no. 3: 411. https://doi.org/10.3390/f15030411

APA Style

Xie, G., Liu, S., Chang, T., & Zhu, N. (2024). Forest Adaptation to Climate Change: Altitudinal Response and Wood Variation in Natural-Growth Cunninghamia lanceolata in the Context of Climate Change. Forests, 15(3), 411. https://doi.org/10.3390/f15030411

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