Next Article in Journal
Responses to Water Deficit and Salt Stress in Silver Fir (Abies alba Mill.) Seedlings
Previous Article in Journal
Refined Systems of National Accounts and Experimental Ecosystem Accounting Versus the Simplified Agroforestry Accounting System: Testing in Andalusian Holm Oak Open Woodlands
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mean Annual Wood Density Variations of Larix gmelinii (Rupr.), Quercus mongolica Fisch. ex Ledeb., and Pinus tabulaeformis Carr. at Two Different Stem Heights

1
Chair of Forest Growth, Albert-Ludwigs-University Freiburg, Tennenbacher Str. 4, 79106 Freiburg, Germany
2
Research Institute of Forest Policy and Information, Chinese Academy of Forestry, No. 1 Dongxiaofu, Haidian District, Beijing 100091, China
*
Author to whom correspondence should be addressed.
Forests 2020, 11(4), 394; https://doi.org/10.3390/f11040394
Submission received: 12 March 2020 / Revised: 24 March 2020 / Accepted: 26 March 2020 / Published: 2 April 2020
(This article belongs to the Section Wood Science and Forest Products)

Abstract

:
Forests are a large carbon sink with an additional substitution effect in the merchantable timber compartment of harvested trees, where carbon stored within the same volume of wood varies depending on wood density. Here, we investigated mean annual air-dry wood density variations depending on cambial age, annual radial increment, and two different stem heights of Larix gmelinii (Rupr.), Quercus mongolica Fisch. ex. Ledeb., and Pinus tabulaeformis Carr. from a first climatic region (Mulan Forest) and exclusively of P. tabulaeformis from a second climatic region (Zhongtiaoshan Forest) in the temperate zone of China. We applied linear mixed-effects models with partly transformed variables and estimated marginal means for pairwise comparisons. Results showed that mean wood density was not significantly different between L. gmelinii (0.626 g cm−3) and Q. mongolica (0.596 g cm−3), but significantly different between P. tabulaeformis from the two different climatic regions (0.445 g cm−3 in Mulan Forest and 0.521 g cm−3 in Zhongtiaoshan Forest). Mean annual wood density within trees except for P. tabulaeformis from Mulan Forest was initially increasing until an intermediate cambial age, after which it decreased again to lower values. These findings showed that tree age had to be considered in assessing carbon sequestration in wood. It also could play an important role in decision making for forest management in Mulan Forest and show the benefit of the wood properties and carbon storage potential of the faster growing L. gmelinii compared to Q. mongolica. Furthermore, these findings gave an indication that intermediate old forest stands for some tree species accumulated more carbon per year within their woody biomass than young stands or old growth forests. Our results may have an impact on the planning of rotation lengths and of tree species composition for forest stands in Mulan Forest and Zhongtiaoshan Forest.

1. Introduction

Forest ecosystems are some of the largest sinks of terrestrial carbon on the globe [1]; they absorb billions of tons of atmospheric CO2 yearly and contribute to climate change mitigation [2]. Within a forest, carbon is stored mainly in the living aboveground biomass and in the organic matter of the mineral soil [3]. The living aboveground biomass can be divided into several compartments, where the merchantable timber compartment has an additional carbon storage function when accounting for the duration of its use as wood for human use and for its substitution effect [4]. The substitution effect can be increased if the used wood fiber is reused in different products in a wood cascade chain [5].
The amount of carbon stored within the same volume of wood varies between tree species depending on wood density. For many tree species, average wood density values are known and accessible (e.g., Global Wood Density Database [6]). Next to its importance for biomass production and the carbon storage potential, wood density also has an effect on the quality of the manufactured end-product as it is linked to several physical and material characteristics such as stiffness and strength [7]. Denser and therefore stronger wood supports and protects the plant against gravity, wind, snow load, and other external factors [8]. Density varies between trees of the same species as it depends on multiple factors such as local climate [9], genetic provenance [10], and age [11]. Furthermore, wood density also varies within an individual tree both horizontally between tree rings and vertically between different stem heights [12].
In artificially regenerated forests, often multiple tree species—both coniferous and broadleaved—can potentially grow on the same site if they are adapted to the present site conditions. The amount of produced biomass and stored carbon in a forest therefore can be influenced by choosing different tree species for regeneration. The analysis of wood density variations of different tree species in the same forest could provide results with direct implications for forest management and silvicultural treatments in order to optimize carbon storage.
Wood density variations have been described for various tree species (e.g., Fagus sylvatica L. [11], Picea abies (L.) H. Karst. [12], Pinus radiata D. Don. [13], Tsuga heterophylla (Raf.) Sarg. [14]), where wood density varies based on tree ring width, cambial age, climate, latewood percentage, radial distance from the pith, and stem height. For economically important tree species in China, wood density variations have also been investigated (e.g., Cunninghamia lanceolata (Lamb.) Hook [15], Pinus massoniana Lamb. [16], and Pinus tabulaeformis Carr. [17,18,19].
Larix gmelinii (Rupr.) is a coniferous tree species with a natural distribution in east Russia, Mongolia, north-east China, and North Korea [20]. Its average air-dry wood density is between 0.528 and 0.599 g cm−3 according to Zanne et al. [6] and 0.684 g cm−3 according to Huang et al. [21]. Jiang et al. [22] found that wood density decreased gradually with increasing stem height.
The tree species P. tabulaeformis occurs naturally in a belt from central towards north-eastern China [23]. Average wood density values are 0.360 g cm−3 [6], 0.432 g cm−3 after 15 days of air drying [17], and 0.551 g cm−3 [18]. Xu [19] found that specific gravity (i.e., the density of wood divided by the density of water) of the wood is highest close to the pith and decreased until the cambial age of 10, after which it increased again towards the bark.
Quercus mongolica Fisch. Ex. Ledeb has its natural range of distribution in east Russia, east Mongolia, north Japan, North Korea, and north-east China [24]. Its average wood density is between 0.603 g cm−3 and 0.660 g cm−3 [6] and 0.748 g cm−3 after 15 days of air drying [24].
For inter- and intra-annual wood density measurements, many different approaches are known and have been used, such as X-ray densitometry [25], video image analysis [26], light reflectance spectroscopy (i.e., near-infrared reflectance [27], reflected-light image analysis [28], THz-imaging [29], and blue reflectance [30]), measurement of drill resistance [31], high-frequency (HF) densitometry [32], microwave scanning [33], and neutron imaging [34].
In the majority of scientific publications related to inter- and intra-annual wood density measurements, X-ray densitometry has been the established method (e.g., Franceschini et al. [35]; Jyske et al. [12]; Bouriaud et al. [11]). However, X-ray densitometry has the disadvantage of having high equipment costs, relying on a time-consuming preparation and measurement of wood samples [32]. A faster and inexpensive alternative is the HF densitometry, which measures relative density values [36,37]. Wassenberg et al. [37] calibrated these relative density values to true volumetric mass density and compared the results to volumetrically determined true wood density and density profiles obtained from image analysis of wood cell structure. Results showed that with linear calibration functions, measured values of HF densitometry can be calibrated to true volumetric mass density values. Boden et al. [36] studied the resolution abilities of HF densitometry and showed that, depending on the geometrical dimensions of the micro-electrode measuring system, it is possible to distinguish variations in wood density up to a resolution of 1 µm.
Here, we investigated mean annual air-dry wood density ( d ¯ ) variations depending on cambial age (CA), annual radial increment (ARI), and stem heights (SH) at 1.3 m (breast height (BH)) and between 10 m and 11 m (upper height (UH)) of L. gmelinii and Q. mongolica and of P. tabulaeformis from two different climatic regions. Air-dry wood density was determined by dividing the mass of air-dry wood by the volume of air-dry wood in the same state (other than for example basic wood density, where the mass of air-dry wood is divided by the volume of water saturated wood). We hypothesized that:
  • Averaged over both stem heights and all observed CA, d ¯ was highest for Q. mongolica and lowest for P. tabulaeformis.
  • For P. tabulaeformis, d ¯ had no significant difference between the two different climatic regions.
  • With increasing CA, d ¯ increased for the two coniferous tree species L. gmelinii and P. tabulaeformis.
  • For Q. mongolica, we expected d ¯ to decrease with increasing CA.
  • For the effect of SH, we expected to observe significantly higher d ¯ at BH compared to at UH at the same CA for Q. mongolica and significantly lower d ¯ for the coniferous tree species.

2. Materials and Methods

2.1. Research Area

The research area of L. gmelinii, Q. mongolica, and P. tabulaeformis was located in Mulan Forest of Hebei Province (referred to as Mulan) and for the second location of P. tabulaeformis in Zhongtiaoshan Forest of Shanxi Province (referred to as Zhongtiaoshan) in Northeast China (Figure 1). Both locations are dominated by temperate coniferous forests, which consist mainly of P. tabulaeformis at an elevation between 950 m and 1650 m above sea level.
The climate at both locations can be classified as a temperate climate, which is influenced by the East Asian Summer Monsoon (EASM) with distinct hot and semi-humid summers and cold and dry winters [39]. Annual precipitation in Mulan ranges from 240 mm to 680 mm. Due to the EASM, about 70% of the total precipitation occurs in the wet season between May and July. A mean annual air temperature of 5.1 °C, ranging from −12.8 °C to 21.0 °C, characterizes this region (meteorological station of Weichang County, 844 m a.s.l., Hebei Province, 1951–2013). Annual precipitation in Zhongtiaoshan ranges from 331 mm to 885 mm, also with a strong influence of EASM. Mean annual temperature is 11.5 °C, ranging from −4.8 °C to 24.4 °C (aggregated climate data from Qinshui meteorological station, 887 m a.s.l., Yangcheng meteorological station, 661 m a.s.l., Gaoping meteorological station, 837 m a.s.l., and Lingchuan meteorological station, 1313 m a.s.l.). According to the Köppen climate classification system, the climate in Mulan is classified as Monsoon-influenced warm-summer humid continental climate (Dwb) and the climate in Zhongtiaoshan as Monsoon-influenced hot-summer humid continental climate (Dwa) [40].

2.2. Collection of Cross-Sections

Trees selected for stem disc collection were taken exclusively from the predominant and dominant canopy layer [41], as these trees experienced less influence from competitors on their growth, and their survival rate was higher compared to dominated individuals in a stand. In addition, sample trees were selected from closed single-species stands, which were not affected by gaps, skid trails, forest roads, and admixed tree species. Furthermore, selected sample trees had healthy and evenly shaped crowns and no visible injuries, damages, or infections (see Table 1 for additional information of sample trees). Moreover, selected forest stands were as old as possible to maximize the observation range of cambial age.
Cross-sections were taken from felled trees at 1.3 m (BH) and between 10 m and 11 m stem height (UH), respectively. The collected cross-sections were then cut along one radius to facilitate a single controlled crack towards the pith while air drying for several weeks. Sampling took place during the growing season of 2014 in Mulan and during the growing season of 2018 in Zhongtiaoshan.

2.3. Density Measurement

Of each cross-section, four radii towards the four cardinal points were cut out into bars and glued on a glass plate. However, if a wood deformation due to an ingrown branch or injury was located directly on a radius towards a cardinal point and caused therefore a variation in tree ring structure and wood density that was not due to annual growth patterns, then the closest alternative radius without any deformation was measured instead.
For surface preparation of the wood bars, a single-point diamond fly cutter equipped with air bearings was used [42]. Here, a diamond blade that moved at a constant rotation speed along a circular track produced smooth surfaces by cutting small chips of the cell wall off the wood bars. Wassenberg et al. [43] showed that the use of the single-point diamond fly cutter with air bearings for wood surface preparation was superior compared to the preparation with a core microtome and sanding with P600 grit.
In this study, the HF densitometry was used to measure relative intra-annual wood density variations [32]. This method utilizes the dielectric properties of wood, which are closely related to wood density. Measurements took place via a micro-electrode system, where a slit-shaped probe had surface contact with a wood bar. The probe consisted of a high-frequency transmitting electrode and a receiving electrode, which were both hermetically separated by a thin metallic foil and had an integration width of 175 µm. The induced electromagnetic stray field propagated through the wood bar, and the detected voltage was higher in denser wood due to its higher local dielectric constant and was linked to the different ratios of cell wall to air-filled cell lumen. If temperature, moisture content, orientation of wood-fibers, and excitation frequency of the electric field are known, then the dielectric constant is linearly related to wood density [37,44].
The wood bars were mounted on a motorized linear travel stage that moved along a cross-sectional radius. The slit of the probe was oriented perpendicular to the measured radius. For HF densitometry, it is crucial that the surface preparation of the wood bars fulfills a high standard, ensuring that no air gaps at any point during measurement are present and that the probe has always full contact with the wood bar. Otherwise, the gap of air between the probe and the wood bar will lead to a decrease in the output voltage and therefore implies lower density [32,43].
Both the single-point diamond fly cutter equipped with air bearings and a prototype of the HF densitometer were developed and set up at the Chair of Forest Growth of the Albert-Ludwigs-University Freiburg, Germany. The utilized software was developed at the Chair of Forest Growth of the Albert-Ludwigs-University Freiburg, Germany, which enabled superimposition of wood density profiles to tree ring series.

2.4. Data Processing

HF measurement of all cross-sectional radii took place over several days. To normalize all measurements to a date where a calibration function was calculated, a reference radius was measured at each day of measurement. The mean density of the total daily reference radius ( x ¯ R i ) was then calculated and divided by the mean density of the total reference radius from the day when the calibration function was calculated ( x ¯ R c a l ). Each measured relative density value was multiplied by the corresponding daily quotient:
d n o r m i = x ¯ R i x ¯ R c a l d i ,
where d n o r m i is the normalized relative volumetric density value at the day of HF measurement i , x ¯ R i is the mean relative volumetric density value of the total reference radius at the day of HF measurement i , x ¯ R c a l is the mean relative volumetric density value of the total reference radius at the day of HF measurement when a calibration function was calculated, and d i is the relative volumetric density value at the day of HF measurement i . The calculated daily quotients ranged from 0.950 to 1.027.
Calibration functions were repeatedly determined to convert measured relative volumetric density from HF densitometry to true volumetric mass density. These calibration functions were based on the linear relationship between the measured relative volumetric density of 6 tropical wood specimens and their true mass density, which was known by determining the weight and volume of the air-dried wood specimen. By use of the calibration functions, the values of voltage can be calculated into values of density:
y = a x + b
where y is the density of air-dry wood in g cm−3, a is the parameter for the slope of the linear relationship, b is the parameter for the intercept at HF output of 0, and x is the measured value of relative density from HF densitometry in volts.
For each density profile, the inter-year transient was removed, as it was an artifact from the integration width of the slit-shaped probe. Specifically, within the data points of a single year, all data points up to the first local minimum and after the last local maximum were removed (see Figure 2 for an example).
Of the cross-sectional radii, mean annual density was calculated as the arithmetic mean of all measured values. d ¯ was only analyzed from the ring of CA 2, as the innermost tree ring of the cut wood bars was in many cases only fragmented due to the width of the saw, and therefore, no true mean values of annual gravimetric density could be obtained. Furthermore, sampling of stem discs took place during the growing season, and the last tree ring of the current calendar year was also not analyzed, as it could never be assured that wood formation processes of the current year had already finished.

2.5. Data Analysis

All statistical methods and data analysis were computed in the R programming environment using the Graphical User Interface R Studio [45,46].
d ¯ was analyzed by using linear mixed-effects models (LMEM), which are especially useful for analysis of non-independent data from repeated measurements with sporadically missing values. The LMEMs were calculated for species subsets using the lmer function of the package lme4 [47]. As scatter plots of mean annual air-dry wood density for each tree species plotted over CA and ARI revealed partly parabolic patterns and partly uneven variability (Figure 3), LMEMs were applied for subsets of species and location, and variables were partly transformed. The models were generally formulated as:
D i , j , k , l = μ + C A i + A R I j + S H k + C A × A R I i , j + C A × S H i , k + A R I × S H i , k + C A × A R I × S H i , j , k + I l + ε i , j , k , l ,
where D i , j , k , l is the response variable mean annual air-dry wood density ( d ¯ ) of the corresponding species and location (L. gmelinii, Q. mongolica, and P. tabulaeformis from Mulan and P. tabulaeformis from Zhongtiaoshan), μ is the general experimental mean, C A i is the continuous fixed-effect C A of cambial age i , A R I j is the continuous fixed-effect A R I of annual radial increment j , S H k is the categorical fixed-effect S H of stem height k (BH and UH), C A × A R I i , j , C A × S H i , k , A R I × S H j , k , and C A × A R I × S H i , j , k are the interaction terms, I l is the random-effect I of the individual tree l , and ε i , j , k , l is the residual error term. To obtain more normal residuals, for L. gmelinii, D i , j , k , l was applied with a power transformation using the exponent of 0.5, and C A i was formulated as a second order polynomial term of CA; for Q. mongolica, C A i and A R I j were formulated as second order polynomial functions of CA and ARI; for P. tabulaeformis from Mulan, D i , j , k , l was applied with a power transformation using the exponent of 0.5; and for P. tabulaeformis from Zhongtiaoshan, D i , j , k , l was applied with a power transformation using the exponent of 0.75, and C A i was formulated as a second order polynomial term of CA. Diagnostic plots of the residuals of each LMEM for species subsets are shown in Appendix A Figure A1, Figure A2, Figure A3 and Figure A4 and were satisfying.
To calculate ANOVA (analysis of variance) tables with p -values and F -statistics, the lmerTest package was used, where the anova and summary functions for objects of the lmerMod class were overloaded to provide p -values for fixed-effects [48].
Post-hoc tests between subsets were performed as pairwise comparisons using the function emmeans of the package emmeans [49]. To avoid higher ranking interactions and therefore potentially misleading results, the model was formulated for the whole and untransformed dataset with the subset ( S ) of species and location as fixed-effect and cambial age ( C A ), stem height ( S H ) and individual tree ( I ) as random-effects. Diagnostic plots of the residuals are shown in Appendix B and were satisfying. For visualization of interaction effects with continuous predictor variables of the LMEMs, plots were produced based on effect displays obtained by the function Effects in the package effects [50].

3. Results

3.1. Effect of Tree Species

A significant effect of subsets ( p < 0.001 ) was caused by significantly lower d ¯ for P. tabulaeformis from Mulan (0.445 g cm−3) compared to all other tree species (L. gmelinii, 0.626 g cm−3, p < 0.001 ; Q. mongolica, 0.596 g cm−3, p < 0.001 ; P. tabulaeformis from Zhongtiaoshan, 0.521 g cm−3, p = 0.050 ). d ¯ of P. tabulaeformis from Zhongtiaoshan was also significantly lower than d ¯ of L. gmelinii ( p = 0.005 ), but no significant difference of d ¯ between Q. mongolica and L. gmelinii and of d ¯ between Q. mongolica and P. tabulaeformis from Zhongtiaoshan was observed (Figure 4).

3.2. Effect of Stem Height

According to the LMEMs for each subset, a significant effect for SH was only found for L. gmelinii and P. tabulaeformis from Mulan (Table 2). For L. gmelinii, d ¯ at BH was on average 0.656 g cm−3 and at UH on average 0.625 g cm−3. For P. tabulaeformis from Mulan, d ¯ at BH was on average 0.463 g cm−3 and at UH on average 0.436 g cm−3.

3.3. Effect of Cambial Age

The fixed-effect of CA was significant for all subsets (Table 2). As can be seen from the effect plot of d ¯ over CA for all subsets, no clear linear increasing or decreasing age trends were visible, but partially parabolic patterns were present in the data (Figure 5). The interaction of CA with SH for subsets was significant for all groups, but Q. mongolica (Table 2). The effect plots for the predictor variables CA and SH of subset models showed similar patterns for both SH of L. gmelinii and P. tabulaeformis from Mulan (Figure 6a,b). The significant interaction effect of CA and SH of P. tabulaeformis from Zhongtiaoshan showed similar density values for the first 40 to 50 years. However, the fitted model indicated a convex and overall increasing pattern for d ¯ at UH, while at BH d ¯ , was overall decreasing in a concave pattern (Figure 6c).

3.4. Effect of Annual Radial Increment

ARI had a significant effect on d ¯ for all subsets, but P. tabulaeformis from Zhongtiaoshan (Table 2). The effect plot of d ¯ over ARI revealed an increasing linear pattern for d ¯ with increasing ARI for L. gmelinii and a decreasing linear pattern for P. tabulaeformis from Mulan. For Q. mongolica, increasing d ¯ with increasing ARI was observed for ARI below 3 mm, after which, it leveled off and decreased again (Figure 7).

3.5. Effect of Cambial Age, Annual Radial Increment, and Stem Height

The three-way interaction effect of CA, ARI, and SH was significant for all subsets (Table 2). Effect plots of the predictor variables CA and SH for multiple levels of ARI showed for L. gmelinii that the parabolic pattern at BH was diminishing with increasing ARI and that a close to linear increase was present. At UH, the curve of d ¯ over CA at an ARI of 1 mm had an overall decreasing shape between 15 years and 72 years. At ARI of 2 mm and 3 mm, a parabolic shape was present and shifting towards higher d ¯ with increasing ARI. At ARI of 4 mm and 5 mm, an increase in d ¯ for observed CA was present (Figure 8a–f).
For Q. mongolica, the corresponding effect plots revealed a different development. At BH, the parabolic pattern of d ¯ with an increase in the first 20 to 25 years remained over all levels of ARI, but was increasing on average. At UH, the curve at ARI of 1 mm was convex, whereas for ARI of 2 mm, the shape for the observed range was concave with a steep slope (Figure 9a–c).
Effect plots for P. tabulaeformis from Mulan showed a linear decrease of d ¯ with increasing CA at BH for tree rings with ARI of 1 mm and almost constant values of d ¯ at UH. With an ARI of 3 mm, the relationship changed to linearly increasing for both SH, which became stronger for higher values of ARI and which was also always more pronounced at BH than at UH (Figure 10a–e).
For P. tabulaeformis from Zhongtiaoshan, the corresponding effect plots revealed a concave parabola at BH that was especially visible at ARI of 2 mm and 3 mm. At UH, modelled d ¯ shifted from a convex pattern for ARI below 3 mm towards a more linear increase at ARI of 3 mm and 4 mm (Figure 11a–e).

4. Discussion

In this study, we used linear mixed-effects models and pairwise comparisons of estimated marginal means to analyze variations of mean annual air-dry wood density ( d ¯ ) of L. gmelinii, Q. mongolica, and P. tabulaeformis from Mulan and P. tabulaeformis from Zhongtiaoshan over CA and ARI and between breast height (BH) (1.3 m) and upper stem height (UH) (between 10 m and 11 m).
One of our main findings was that d ¯ was significantly higher for L. gmelinii compared to P. tabulaeformis from both locations and that Q. mongolica was only significantly higher than P. tabulaeformis from Mulan, but not significantly different from L. gmelinii and P. tabulaeformis from Zhongtiaoshan (Figure 4), which was not fully in line with our first hypothesis. We formulated our first hypothesis according to Zanne et al. [6], where average wood density for Q. mongolica was between 0.603 g cm−3 and 0.660 g cm−3, average wood density for L. gmelinii between 0.528 g cm−3 and 0.599 g cm−3, and average wood density of P. tabulaeformis 0.360 g cm−3. However, the values of d ¯ for L. gmelinii found in this study were higher than most values from the literature [6,51], while those values of d ¯ for Q. mongolica were lower than most values from the literature [6,24]. Wassenberg et al. [52] showed for multiple tree species that, as wood density values of a single tree can vary within the stem, estimated mean values also vary based on the sampling strategy. It is therefore possible that under a different sampling strategy, mean wood density values would show higher similarity to cited values.
Other than hypothesized, the means of d ¯ were significantly higher for P. tabulaeformis from Zhongtiaoshan than those of P. tabulaeformis from Mulan. Figure 3g,h shows that different ARI for the different locations could not be the cause as d ¯ remained systematically higher at Zhongtiaoshan over the whole observed range of ARI. Furthermore, the longer observation period of CA in Mulan of 100 years was not an explanation for the observed differences of d ¯ as also for this predictor, d ¯ remained relatively stable (Figure 5). A possible explanation could be differences in climatic conditions, as meteorological data indicated a warmer and moister climate at Zhongtiaoshan, which in turn may had led to adaptions to different site conditions. In fact, when looking at average tree height at 50 years of tree age, trees in Mulan were 15.9 m high and trees in Zhongtiaoshan 12.2 m high, indicating different site productivities. Gindl et al. [53] found that lignin content in secondary cell wall layers of terminal latewood tracheids of P. abies positively correlated with the air temperature of August and September. Rossi et al. [54] suggested that cell production of conifers in the early growing season allowed tracheid differentiation to be completed and cell wall lignification to be adequate. As the growing season in Zhongtiaoshan was longer and started earlier than in Mulan with higher monthly air temperature, a higher amount of lignification and therefore denser wood could be the result. Larson [55] found on Pinus elliottii Engelm. in the south-east United States that the percentage of latewood within a tree ring increased with higher precipitation in June and July. Additionally, the evolution of different provenances could also explain differences in d ¯ between the two locations, leading to small differences in wood anatomy and growth responses. Follow-up research with the consideration of climatic parameters and an analysis of growth response to climatic conditions could provide explanations of different wood density values for P. tabulaeformis. Furthermore, additional research on wood density depending on different site productivity classes could lead to explanations of our observed results.
Contradicting to our third hypothesis, we did not always observe obvious increasing patterns of d ¯ with increasing CA for the coniferous tree species (Figure 5); only at UH of P. tabulaeformis from Zhongtiaoshan, an overall increase was present (Figure 6c). The expected reason was that with increasing CA, ARI decreased, while at the same time, the ratio of latewood to earlywood within a tree ring increased [56]. The majority of cells in the xylem of coniferous tree species are tracheids, which have a higher ratio of cell wall to lumen in latewood, causing a higher wood density compared to earlywood. However, for P. tabulaeformis from Mulan, decreasing d ¯ with increasing ARI could be observed, which was expected to be linked with increasing d ¯ over increasing CA (Figure 7). This finding did not correspond with the results of Xu [19], who observed an increase of wood density after the 10th tree ring for P. tabulaeformis. Possibly an increase in sample size and a follow up study could lead to clearer results.
An expected decrease of d ¯ for increasing CA of Q. mongolica could be observed only after an initial increase (Figure 5). A Similar decrease of wood density with increasing CA was shown for other oak species such as Quercus suber L. [57] and Quercus petraea Liebl. [58]. As a ring-porous broadleaved tree species, Q. mongolica produces large diameter tracheas for water transportation in the beginning of the growing season. With narrower tree rings, the ratio of tracheas within the xylem increases, resulting in a decrease of wood density [56]. As expected, d ¯ increases with increasing ARI, which was expected to be linked with decreasing d ¯ over CA (Figure 7). The same relationship of increasing wood density with increasing tree ring width was found by Zhang et al. [59] for Q. petraea and Quercus robur L. This led to the conclusion that for the investigated trees, no clear linear relation between CA and ARI was present. Especially at UH with an average observation period of 28 years, long term growth may not be detectable. However, a parabolic pattern in specific gravity over CA was observed for the closely related tree species Quercus liaotungensis Koidz., which supported our results [60].
For the effect of SH, we expected to observe significantly higher d ¯ at BH compared to at UH at the same CA for Q. mongolica and significantly lower d ¯ for the coniferous tree species. For annual growth during the same calendar year, different studies for multiple tree species showed that no differences in ARI were present, explained solely by different stem heights [61,62], and that ARI was highly similar at different stem heights [63]. We assumed therefore that ARI at younger CA at UH was less than ARI at younger CA at BH. For tree rings with less ARI, we expected less d ¯ for Q. mongolica and higher d ¯ for both coniferous tree species, as explained already in relation to tree ring anatomy.
A significant effect of SH was only present for L. gmelinii and P. tabulaeformis from Mulan. Figure 3 indicates that d ¯ was lower at UH of both subsets. Repola [64] found for P. sylvestris in southern Finland also a decrease of wood density with increasing SH. Furthermore, the pattern of decreasing wood density with increasing SH was observed for Pinus massoniana Lamb. and Pinus koraiensis Siebold & Zucc. in China [16,65]. However, Jyske et al. [12] found for P. abies an increasing pattern of wood density with increasing SH in a given CA, which showed that the results in the literature were not consistent for the family Pinaceae.
A contrasting development of d ¯ with increasing CA was observed for the two different SH of P. tabulaeformis from Zhongtiaoshan; at BH, d ¯ overall decreased in a concave pattern, while at UH, d ¯ increased in a convex pattern (Figure 6c). Additionally, Figure 11 reveals that the initial convex pattern at UH changed towards a more linear increasing pattern for higher values of ARI. Molteberg and Høibø [66] found a similar convex pattern for basic wood density over cambial age in 28-year-old P. abies stems. Furthermore, Petty et al. [67] found for P. abies and Picea sitchensis (Bong.) Carr. an initial decrease of wood density followed by a gradual increase. However, for all curves at UH, data were present for relatively short observation periods, leading to some degree of uncertainty in the modelled values.

5. Conclusions

Our study showed that mean wood density was not significantly different between L. gmelinii and Q. mongolica. This finding could play an important role in decision making for forest management in Mulan and showed the benefit of the wood properties and carbon storage potential of the faster growing L. gmelinii compared to Q. mongolica, as indicated by overall narrower ARI for Q. mongolica and overall higher sample trees of L. gmelinii (Figure 3e,f, Table 1). Furthermore, we found that wood density within a tree was not monotonically increasing with cambial age, but for some subsets, levelling off eventually and even decreasing to lower annual values. This finding gave an indication that intermediate old forest stands possibly accumulate more carbon per year within their woody biomass than young stands or old growth forests, while at the same time, the time of maximum mean annual carbon sequestration would be delayed. This conclusion may have an impact on the planning of rotation lengths for forest stands in Mulan and Zhongtiaoshan. However, it also showed that tree age needs to be considered in assessing carbon sequestration in wood.
Furthermore, we stress that due to the small sample size of five trees per subset, the external validity of our results was not given, and the results were only valid for pre-dominant and dominant trees of the same tree species in the same geographical region, Mulan and Zhongtiaoshan. Nevertheless, our work served well as an exploratory study, indicating interesting and novel results and interactions of effects on wood density variations. We therefore see high benefit in follow-up research on wood density variations within stems of tree species in China with increased sample size and possibly additional levels of the effect of stem height. Furthermore, additional attention needs to be given to a climate analysis of growth reactions and the effect of site productivity and tree age on wood density.

Author Contributions

Conceptualization, L.S., S.W., and H.S.; data curation, L.S. and S.-M.H.; formal analysis, L.S., Z.C., and S.-M.H.; funding acquisition, L.S., S.W., and H.S.; investigation, L.S.; methodology, L.S.; project administration, S.W. and H.S.; resources, L.S., S.-M.H., S.W., and H.S.; software, L.S.; supervision, S.W. and H.S.; validation, L.S., Z.C., and H.S.; visualization, L.S.; writing, original draft, L.S.; writing, review and editing, Z.C., S.-M.H., S.W., and H.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the German Federal Ministry of Education and Research (BMBF) within the project Lin4Future, Grant Number 01LZ1701A, the German Federal Ministry of Food and Agriculture (BMEL) within the project Chinawald, Grant Number 28I-010-01, and the National Natural Science Foundation of China (NSFC), Support Code 31270681. The article processing charge was funded by the German Research Foundation (DFG) and the University of Freiburg in the funding program Open Access Publishing.

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; nor in the decision to publish the results.

Appendix A

Figure A1. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of L. gmelinii.
Figure A1. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of L. gmelinii.
Forests 11 00394 g0a1
Figure A2. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of Q. mongolica.
Figure A2. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of Q. mongolica.
Forests 11 00394 g0a2
Figure A3. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of P. tabulaeformis from Mulan.
Figure A3. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of P. tabulaeformis from Mulan.
Forests 11 00394 g0a3
Figure A4. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of P. tabulaeformis from Zhongtiaoshan.
Figure A4. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model of P. tabulaeformis from Zhongtiaoshan.
Forests 11 00394 g0a4

Appendix B

Figure A5. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model for pairwise comparisons of estimated marginal means between subsets of tree species and location.
Figure A5. Plot of model residuals against fitted values (a), histogram of model residuals (b), and QQ-plot of model residuals (c) of the linear mixed-effects model for pairwise comparisons of estimated marginal means between subsets of tree species and location.
Forests 11 00394 g0a5

References

  1. Pan, Y.; Birdsey, R.A.; Fang, J.; Houghton, R.; Kauppi, P.E.; Kurz, W.A.; Phillips, O.L.; Shvidenko, A.; Lewis, S.L.; Canadell, J.G.; et al. A Large and Persistent Carbon Sink in the World’s Forests. Science 2011, 333, 988–993. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Canadell, J.G.; Raupach, M.R. Managing Forests for Climate Change Mitigation. Science 2008, 320, 1456–1457. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Fahey, T.J.; Woodbury, P.B.; Battles, J.J.; Goodale, C.L.; Hamburg, S.P.; Ollinger, S.V.; Woodall, C.W. Forest carbon storage: Ecology, management, and policy. Front. Ecol. Environ. 2010, 8, 245–252. [Google Scholar] [CrossRef] [Green Version]
  4. Sathre, R.; O’Connor, J. Meta-analysis of greenhouse gas displacement factors of wood product substitution. Environ. Sci. Policy 2010, 13, 104–114. [Google Scholar] [CrossRef]
  5. Sathre, R.; Gustavsson, L. Energy and carbon balances of wood cascade chains. Resour. Conserv. Recycl. 2006, 47, 332–355. [Google Scholar] [CrossRef]
  6. Zanne, A.E.; Lopez-Gonzalez, G.; Coomes, D.A.; Ilic, J.; Jansen, S.; Lewis, S.L.; Miller, R.B.; Swenson, N.G.; Wiemann, M.C.; Chave, J. Data from: Towards a worldwide wood economics spectrum. Dry. Digit. Repos. 2009. [Google Scholar] [CrossRef]
  7. Niklas, K.J.; Spatz, H.-C. Worldwide correlations of mechanical properties and green wood density. Am. J. Bot. 2010, 97, 1587–1594. [Google Scholar] [CrossRef]
  8. 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]
  9. Briffa, K.R.; Osborn, T.J.; Schweingruber, F.H.; Jones, P.D.; Shiyatov, S.G.; Vaganov, E.A. Tree-ring width and density data around the Northern Hemisphere: Part 1, local and regional climate signals. Holocene 2002, 12, 737–757. [Google Scholar] [CrossRef]
  10. Hannrup, B.; Cahalan, C.; Chantre, G.; Grabner, M.; Karlsson, B.; Bayon, I.L.; Jones, G.L.; Müller, U.; Pereira, H.; Rodrigues, J.C.; et al. Genetic parameters of growth and wood quality traits in Picea abies. Scand. J. For. Res. 2004, 19, 14–29. [Google Scholar] [CrossRef]
  11. Bouriaud, O.; Bréda, N.; Le Moguédec, G.; Nepveu, G. Modelling variability of wood density in beech as affected by ring age, radial growth and climate. Trees 2004, 18, 264–276. [Google Scholar] [CrossRef]
  12. Jyske, T.; Mäkinen, H.; Saranpää, P. Wood density within Norway spruce stems. Silva Fenn. 2008, 42, 439–455. [Google Scholar] [CrossRef] [Green Version]
  13. Ivković, M.; Gapare, W.; Wu, H.; Espinoza, S.; Rozenberg, P. Influence of cambial age and climate on ring width and wood density in Pinus radiata families. Ann. For. Sci. 2013, 70, 525–534. [Google Scholar] [CrossRef] [Green Version]
  14. DeBell, D.S.; Singleton, R.; Gartner, B.L.; Marshall, D.D. Wood density of young-growth western hemlock: Relation to ring age, radial growth, stand density, and site quality. Can. J. For. Res. 2004, 34, 2433–2442. [Google Scholar] [CrossRef]
  15. Hu, D.; Ruan, Z.; Qian, Z.; Huang, Y.; Huang, X.; Chen, Z. The Genetics and Variation of Wood Density in Chinese Fir Clones and Its Interaction to the Growth Characters. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFD2004&filename=ZNLB200405006&v=MTY4MzkxRnJDVVI3cWZZK1p1Rnlqa1ZyL0xQeVBIYkxHNEh0WE1xbzlGWW9SOGVYMUx1eFlTN0RoMVQzcVRyV00 (accessed on 11 March 2020).
  16. Liu, Q.; Zhang, R.; Jin, G.; Chu, D.; Zhou, Z. Variation of Ring Width and Wood Basic Density and Early Selection of Pinus massoniana Provenances. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFD2010&filename=LYKE201005010&v=MTYzMTdyM1BLVFRBYTdHNEg5SE1xbzlFWklSOGVYMUx1eFlTN0RoMVQzcVRyV00xRnJDVVI3cWZZK1p1RnlqbFU (accessed on 11 March 2020).
  17. Hong, P. Pinus tabulaeformis. In Enzyklopädie der Holzgewächse: Handbuch und Atlas der Dendrologie; Wiley-VCH: Weinheim, Germany, 2004; ISBN 978-3-527-32141-4. [Google Scholar]
  18. Yin, M.; Yang, L.; Yin, W.; Lei, Q.; Tan, X.; Zhang, Y.; Li, Z. Dynamic changes of carbon storage in trunk volume of Pinus tabulaeformis, Robinia pseudoacacia and Populus euramericana. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFD2011&filename=BJLY201105013&v=MjczNzhSOGVYMUx1eFlTN0RoMVQzcVRyV00xRnJDVVI3cWZZK1p1RnlqbFU3L0FKeWZIZDdHNEg5RE1xbzlFWjQ (accessed on 11 March 2020).
  19. Xu, Y. Specific Gravity Variation in Pinus Tabulaeformis. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFD9093&filename=HZNY199103011&v=MzAyMzc3dkJMVGZGZDdLeEY5RE1ySTlFWllSOGVYMUx1eFlTN0RoMVQzcVRyV00xRnJDVVI3cWZZK1p1RnlqbFU (accessed on 11 March 2020).
  20. Debreczy, Z.; Rácz, I.; Musial, K. Conifers around the World: Conifers of the Temperate Zones and Adjacent Regions. Vol. 1; DendroPress: Budapest, Hungary, 2011; ISBN 978-963-219-063-1. [Google Scholar]
  21. Huang, P.; Zhuo, L.; Li, G.; Wei, C.; Zhou, W.; Li, Y.; Chen, T.; Wu, J. Studies on Physical and Mechanical Properties of Dahurian Larch-Wood of Two Wood Colours. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFD9495&filename=DBLY403.003&v=MjAxMDU3ZTRIYy9NcjR3cUY1NE9mZ2c1emhBVTRqaDRPWDZUckgwM2ViQ1VSYjZmWnVkdUZ5dmdWUT09SVMvSGQ (accessed on 11 March 2020).
  22. Jiang, L.; Liu, M.; Liu, Y. Variation of Wood Basic Density and Early Selection of Dahurian Larch and Mongolian Pine. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFD2013&filename=BJLY201301003&v=MjY0ODV4WVM3RGgxVDNxVHJXTTFGckNVUjdxZlkrWnVGeWpsVWIzT0p5ZkhkN0c0SDlMTXJvOUZaNFI4ZVgxTHU (accessed on 11 March 2020).
  23. Farjon, A. Pinaceae: Drawings and Descriptions of the Genera Abies, Cedrus, Pseudolarix, Keteleeria, Nothotsuga, Tsuga, Cathaya, Pseudotsuga, Larix and Picea; Regnum vegetabile; Koeltz: Königstein, Germany, 1990; ISBN 978-3-87429-298-6. [Google Scholar]
  24. Puhua, H. Quercus mongolica. In Enzyklopädie der Holzgewächse: Handbuch und Atlas der Dendrologie; Wiley-VCH: Weinheim, Germany, 2004; ISBN 978-3-527-32141-4. [Google Scholar]
  25. Schweingruber, F.H.; Fritts, H.C.; Bräker, O.U.; Drew, L.G.; Schär, E. The X-Ray Technique as Applied to Dendroclimatology. Tree-Ring Bul. 1978, 38, 61–91. [Google Scholar]
  26. Park, W.-K.; Telewski, F.W. Measuring maximum latewood density by image analysis at the cellular level. Wood Fibre Sci. 2007, 25, 326–332. [Google Scholar]
  27. Hoffmeyer, P.; Pedersen, J.G. Evaluation of density and strength of Norway spruce wood by near infrared reflectance spectroscopy. Holz Roh- Werkst. 1995, 53, 165–170. [Google Scholar] [CrossRef]
  28. Sheppard, P.R.; Graumlich, L.J.; Conkey, L.E. Reflected-light image analysis of conifer tree rings for reconstructing climate. Holocene 1996, 6, 62–68. [Google Scholar] [CrossRef]
  29. Koch, M.; Hunsche, S.; Schumacher, P.; Nuss, M.C.; Feldmann, J.; Fromm, J. THz-imaging: A new method for density mapping of wood. Wood Sci. Technol. 1998, 32, 421–427. [Google Scholar] [CrossRef]
  30. McCarroll, D.; Pettigrew, E.; Luckman, A.; Guibal, F.; Edouard, J.-L. Blue Reflectance Provides a Surrogate for Latewood Density of High-latitude Pine Tree Rings. Arct. Antarct. Alp. Res. 2002, 34, 450–453. [Google Scholar] [CrossRef]
  31. Rinn, F.; Schweingruber, F.-H.; Schär, E. RESISTOGRAPH and X-Ray Density Charts of Wood. Comparative Evaluation of Drill Resistance Profiles and X-ray Density Charts of Different Wood Species. Holzforschung 1996, 50, 303–311. [Google Scholar] [CrossRef]
  32. Schinker, M.G.; Hansen, N.; Spiecker, H. High-Frequency Densitometry-A New Method for the Rapid Evaluation of Wood Density Variations. IAWA J. 2003, 24, 231–239. [Google Scholar] [CrossRef]
  33. Johansson, J.; Hagman, O.; Fjellner, B.-A. Predicting moisture content and density distribution of Scots pine by microwave scanning of sawn timber. J. Wood Sci. 2003, 49, 312–316. [Google Scholar] [CrossRef]
  34. Mannes, D.; Lehmann, E.; Cherubini, P.; Niemz, P. Neutron imaging versus standard X-ray densitometry as method to measure tree-ring wood density. Trees 2007, 21, 605–612. [Google Scholar] [CrossRef] [Green Version]
  35. Franceschini, T.; Bontemps, J.-D.; Gelhaye, P.; Rittie, D.; Herve, J.-C.; Gegout, J.-C.; Leban, J.-M. Decreasing trend and fluctuations in the mean ring density of Norway spruce through the twentieth century. Ann. For. Sci. 2010, 67, 816. [Google Scholar] [CrossRef] [Green Version]
  36. Boden, S.; Schinker, M.G.; Duncker, P.; Spiecker, H. Resolution abilities and measuring depth of High-Frequency densitometry on wood samples. Measurement 2012, 45, 1913–1921. [Google Scholar] [CrossRef]
  37. Wassenberg, M.; Montwé, D.; Kahle, H.-P.; Spiecker, H. Exploring high frequency densitometry calibration functions for different tree species. Dendrochronologia 2014, 32, 273–281. [Google Scholar] [CrossRef]
  38. East Asia Free Map, Free Blank Map, Free Outline Map, Free Base Map Hydrography, States, Main Cities, Names. Available online: https://d-maps.com/carte.php?num_car=28780&lang=en (accessed on 23 January 2020).
  39. Li, Q.; Liu, Y.; Song, H.; Cai, Q.; Yang, Y. Long-term variation of temperature over North China and its links with large-scale atmospheric circulation. Quat. Int. 2013, 283, 11–20. [Google Scholar] [CrossRef]
  40. Köppen, W. Klassifikation der Klimate nach Temperatur, Niederschlag und Jahreslauf. Pet. Mitt. 1918, 64, 193–203. [Google Scholar]
  41. Kraft, G. Beiträge zur Lehre von den Durchforstungen, Schlagstellungen und Lichtungshieben; Klindworth: Hannover, Germany, 1884. [Google Scholar]
  42. Spiecker, H.; Schinker, M.G.; Hansen, J.; Park, Y.-I.; Ebding, T.; Döll, W. Cell structure in tree rings: Novel methods for preparation and image analysis of large cross sections. IAWA J. 2000, 21, 361–373. [Google Scholar] [CrossRef]
  43. Wassenberg, M.; Schinker, M.; Spiecker, H. Technical aspects of applying high frequency densitometry: Probe-sample contact, sample surface preparation and integration width of different dielectric probes. Dendrochronologia 2015, 34, 10–18. [Google Scholar] [CrossRef]
  44. Torgovnikov, G.I. Dielectric Properties of Wood and Wood-Based Materials; Springer: Berlin, Germany, 2012; ISBN 978-3-642-77455-3. [Google Scholar]
  45. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2019. [Google Scholar]
  46. RStudio Team. RStudio: Integrated Development Environment for R; RStudio, Inc.: Boston, MA, USA, 2019. [Google Scholar]
  47. Bates, D.; Mächler, M.; Bolker, B.; Walker, S. Fitting Linear Mixed-Effects Models Using lme4. J. Stat. Softw. 2015, 67, 1–48. [Google Scholar] [CrossRef]
  48. Kuznetsova, A.; Brockhoff, P.B.; Christensen, R.H.B. lmerTest Package: Tests in Linear Mixed Effects Models. J. Stat. Softw. 2017, 82, 1–26. [Google Scholar] [CrossRef] [Green Version]
  49. Lenth, R. emmeans: Estimated Marginal Means, aka Least-Squares Means. Available online: https://CRAN.R-project.org/package=emmeans (accessed on 29 May 2019).
  50. Fox, J.; Weisberg, S. An R Companion to Applied Regression, 3rd ed.; SAGE: Los Angeles, CA, USA, 2019; ISBN 978-1-5443-3647-3. [Google Scholar]
  51. Liu, M.Y.; Liu, Y.B.; Jiang, L.C. Variation of Wood Density in Larch in Northeast of China. Adv. Mater. Res. 2012, 487, 38–42. [Google Scholar] [CrossRef]
  52. Wassenberg, M.; Chiu, H.-S.; Guo, W.; Spiecker, H. Analysis of wood density profiles of tree stems: Incorporating vertical variations to optimize wood sampling strategies for density and biomass estimations. Trees 2015, 29, 551–561. [Google Scholar] [CrossRef]
  53. Gindl, W.; Grabner, M.; Wimmer, R. The influence of temperature on latewood lignin content in treeline Norway spruce compared with maximum density and ring width. Trees 2000, 14, 409–414. [Google Scholar] [CrossRef]
  54. Rossi, S.; Deslauriers, A.; Anfodillo, T.; Morin, H.; Saracino, A.; Motta, R.; Borghetti, M. Conifers in cold environments synchronize maximum growth rate of tree-ring formation with day length. New Phytol. 2006, 170, 301–310. [Google Scholar] [CrossRef] [Green Version]
  55. Larson, P.R. Effect of Environment on the Percentage of Summerwood and Specific Gravity of Slash Pine; Yale School of Forestry & Environmental Studies Bulletin Series; Yale University: New Haven, CT, USA, 1957. [Google Scholar]
  56. Fritts, H.C. Tree Rings and Climate; Academic Press: London, UK, 1976; ISBN 978-0-12-268450-0. [Google Scholar]
  57. Knapic, S.; Louzada, J.L.; Leal, S.; Pereira, H. Radial variation of wood density components and ring width in cork oak trees. Ann. For. Sci. 2007, 64, 211–218. [Google Scholar] [CrossRef] [Green Version]
  58. Bergès, L.; Dupouey, J.-L.; Franc, A. Long-term changes in wood density and radial growth of Quercus petraea Liebl. in northern France since the middle of the nineteenth century. Trees 2000, 14, 398–408. [Google Scholar] [CrossRef]
  59. Zhang, S.-Y.; Owoundi, R.E.; Nepveu, G.; Mothe, F.; Dhôte, J.-F. Modelling wood density in European oak (Quercus petraea and Quercus robur) and simulating the silvicultural influence. Can. J. For. Res. 1993, 23, 2587–2593. [Google Scholar] [CrossRef]
  60. Zhang, S.-Y.; Zhong, Y. Effect of growth rate on specific gravity of East-Liaoning oak (Quercus liaotungensis) wood. Can. J. For. Res. 1991, 21, 255–260. [Google Scholar] [CrossRef]
  61. Bouriaud, O.; Bréda, N.; Dupouey, J.-L.; Granier, A. Is ring width a reliable proxy for stem-biomass increment? A case study in European beech. Can. J. For. Res. 2005, 35, 2920–2933. [Google Scholar] [CrossRef]
  62. Sprengel, L.; Stangler, D.F.; Sheppard, J.; Morhart, C.; Spiecker, H. Comparative Analysis of the Effects of Stem Height and Artificial Pruning on Seasonal Radial Growth Dynamics of Wild Cherry (Prunus avium L.) and Sycamore (Acer pseudoplatanus L.) in a Widely Spaced System. Forests 2018, 9, 174. [Google Scholar] [CrossRef] [Green Version]
  63. van der Maaten-Theunissen, M.; Bouriaud, O. Climate–growth relationships at different stem heights in silver fir and Norway spruce. Can. J. For. Res. 2012, 42, 958–969. [Google Scholar] [CrossRef]
  64. Repola, J. Models for vertical wood density of Scots pine, Norway spruce and birch stems, and their application to determine average wood density. Silva Fenn. 2006, 40. [Google Scholar] [CrossRef] [Green Version]
  65. Liang, D.; Jin, Y.; Zhao, G.; Dong, Y.; Leng, W.; Chen, C.; Wang, H.; Zhao, X. Variance Analyses of Growth and Wood Characteristics of 50 Pinus Koraiensis Clones. Available online: https://kns.cnki.net/KCMS/detail/detail.aspx?dbcode=CJFQ&dbname=CJFDLAST2016&filename=BJLY201606007&v=MTc4ODRkN0c0SDlmTXFZOUZZNFI4ZVgxTHV4WVM3RGgxVDNxVHJXTTFGckNVUjdxZlkrWnVGeWpsVnIzQkp5Zkg (accessed on 11 March 2020).
  66. Molteberg, D.; Høibø, O. Development and variation of wood density, kraft pulp yield and fibre dimensions in young Norway spruce (Picea abies). Wood Sci. Technol. 2006, 40, 173–189. [Google Scholar] [CrossRef]
  67. Petty, J.A.; Macmillan, D.C.; Steward, C.M. Variation of Density and Growth Ring Width in Stems of Sitka and Norway Spruce. Forestry 1990, 63, 39–49. [Google Scholar] [CrossRef]
Figure 1. Map of north-east China indicating the research area of Mulan with the upper red circle and of Zhongtiaoshan with the lower red circle [38] (modified).
Figure 1. Map of north-east China indicating the research area of Mulan with the upper red circle and of Zhongtiaoshan with the lower red circle [38] (modified).
Forests 11 00394 g001
Figure 2. Exemplary relative density profile of LM4 at breast height (BH), radius towards north. (a) shows the total relative density profile from 1923 to 2013, the red box indicates the last five years; (b) shows the last five years of the relative density profile with the inter-year transients being highlighted with red hatching; (c) shows the processed last five years of the relative density profile, where inter-year transients are removed. All data points are in relative density values as the output of high-frequency (HF) densitometry in volts.
Figure 2. Exemplary relative density profile of LM4 at breast height (BH), radius towards north. (a) shows the total relative density profile from 1923 to 2013, the red box indicates the last five years; (b) shows the last five years of the relative density profile with the inter-year transients being highlighted with red hatching; (c) shows the processed last five years of the relative density profile, where inter-year transients are removed. All data points are in relative density values as the output of high-frequency (HF) densitometry in volts.
Forests 11 00394 g002
Figure 3. Scatter plots of measured mean annual air-dry wood density over cambial age (CA) (top row) and over annual radial increment (ARI) (bottom row) for all subsets L. gmelinii, Q. mongolica, and P. tabulaeformis from Mulan and P. tabulaeformis from Zhongtiaoshan. Darker scatters show mean annual air-dry density values at BH and lighter scatters at upper height (UH).
Figure 3. Scatter plots of measured mean annual air-dry wood density over cambial age (CA) (top row) and over annual radial increment (ARI) (bottom row) for all subsets L. gmelinii, Q. mongolica, and P. tabulaeformis from Mulan and P. tabulaeformis from Zhongtiaoshan. Darker scatters show mean annual air-dry density values at BH and lighter scatters at upper height (UH).
Forests 11 00394 g003
Figure 4. Estimated marginal means of d ¯ for all subsets L. gmelinii, Q. mongolica, and P. tabulaeformis from Mulan and P. tabulaeformis from Zhongtiaoshan. Error bars display the standard errors of the estimated marginal means. Different lowercase letters indicate significant differences with p-values < 0.05, α-level for pairwise comparisons.
Figure 4. Estimated marginal means of d ¯ for all subsets L. gmelinii, Q. mongolica, and P. tabulaeformis from Mulan and P. tabulaeformis from Zhongtiaoshan. Error bars display the standard errors of the estimated marginal means. Different lowercase letters indicate significant differences with p-values < 0.05, α-level for pairwise comparisons.
Forests 11 00394 g004
Figure 5. Effect plot of the significant effect of the predictor variable CA on d ¯ . Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Figure 5. Effect plot of the significant effect of the predictor variable CA on d ¯ . Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Forests 11 00394 g005
Figure 6. Effect plots of significant interaction effects of the predictor variables SH and CA on d ¯ . Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Figure 6. Effect plots of significant interaction effects of the predictor variables SH and CA on d ¯ . Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Forests 11 00394 g006
Figure 7. Effect plot of the significant effect of the predictor variable ARI on d ¯ . Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Figure 7. Effect plot of the significant effect of the predictor variable ARI on d ¯ . Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Forests 11 00394 g007
Figure 8. Effect plots of the predictor variables stem height (SH) and CA on d ¯ for ARI of 1, 2, …, 6 mm of L. gmelinii. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Figure 8. Effect plots of the predictor variables stem height (SH) and CA on d ¯ for ARI of 1, 2, …, 6 mm of L. gmelinii. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Forests 11 00394 g008
Figure 9. Effect plots of the predictor variables SH and CA on d ¯ for ARI of 1, 2, and 3 mm of Q. mongolica. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Figure 9. Effect plots of the predictor variables SH and CA on d ¯ for ARI of 1, 2, and 3 mm of Q. mongolica. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Forests 11 00394 g009
Figure 10. Effect plots of the predictor variables SH and CA on d ¯ for ARI of 1, 2, …, 5 mm of P. tabulaeformis from Mulan. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Figure 10. Effect plots of the predictor variables SH and CA on d ¯ for ARI of 1, 2, …, 5 mm of P. tabulaeformis from Mulan. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Forests 11 00394 g010
Figure 11. Effect plots of the predictor variables SH and CA on d ¯ for ARI of 1, 2, …, 5 mm of P. tabulaeformis from Zhongtiaoshan. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Figure 11. Effect plots of the predictor variables SH and CA on d ¯ for ARI of 1, 2, …, 5 mm of P. tabulaeformis from Zhongtiaoshan. Solid lines are modeled d ¯ values, and dashed lines are the standard errors of modeled values.
Forests 11 00394 g011
Table 1. Sample trees with tree species, location, age at 1.3 m stem height in years, tree height in m, diameter at breast height (DBH) in cm, elevation in m, and slope exposition.
Table 1. Sample trees with tree species, location, age at 1.3 m stem height in years, tree height in m, diameter at breast height (DBH) in cm, elevation in m, and slope exposition.
Tree No.Species LocationAge at 1.3 m Stem Height (y)Tree Height (m)DBH (cm)Elevation (m)Exposition
LM1L. gmelinii
Mulan
8027.0053.01450NW
LM27927.1658.01430NW
LM38027.6448.81430NW
LM49425.6552.41293NE
LM59325.3052.71338E
QM1Q. mongolica
Mulan
5511.7527.51352NA
QM26615.3028.61130W
QM36615.4228.61352S
QM46614.0629.21347S
QM57013.6528.21185N
PM1P. tabulaeformis
Mulan
6719.8542.81162NW
PM27024.6041.5954NE
PM310125.8050.6961NE
PM49824.7042.21020NW
PM510023.7446.81013NW
PZ1P. tabulaeformis
Zhongtiaoshan
7111.9821.21223NW
PZ28715.4034.51190NE
PZ37516.1441.41614NW
PZ47216.5063.01211NW
PZ55517.9034.51578N
Table 2. ANOVA table of the linear mixed-effects models for subsets of species and location. p-values < 0.05 α-level are highlighted in bold type.
Table 2. ANOVA table of the linear mixed-effects models for subsets of species and location. p-values < 0.05 α-level are highlighted in bold type.
Fixed-Effect(s)SubsetF-Valuep-ValueSubsetF-Valuep-Value
C A L. gmelinii29.69<0.001P. tabulaeformis Mulan19.39<0.001
A R I 11.14<0.00156.26<0.001
S H 87.36<0.00126.17<0.001
C A × A R I 6.350.00217.80<0.001
C A × S H 16.28<0.0017.560.006
A R I × S H 34.29<0.00120.52<0.001
C A × A R I × S H 14.25<0.0015.950.015
C A Q. mongolica12.38<0.001P. tabulaeformis Zhongtiaoshan4.500.012
A R I 12.08<0.0010.610.436
S H 0.220.6420.400.529
C A × A R I 7.71<0.00110.62<0.001
C A × S H 1.960.1424.590.011
A R I × S H 0.020.9780.490.482
C A × A R I × S H 3.260.0123.050.048

Share and Cite

MDPI and ACS Style

Sprengel, L.; Cheng, Z.; Hipler, S.-M.; Wu, S.; Spiecker, H. Mean Annual Wood Density Variations of Larix gmelinii (Rupr.), Quercus mongolica Fisch. ex Ledeb., and Pinus tabulaeformis Carr. at Two Different Stem Heights. Forests 2020, 11, 394. https://doi.org/10.3390/f11040394

AMA Style

Sprengel L, Cheng Z, Hipler S-M, Wu S, Spiecker H. Mean Annual Wood Density Variations of Larix gmelinii (Rupr.), Quercus mongolica Fisch. ex Ledeb., and Pinus tabulaeformis Carr. at Two Different Stem Heights. Forests. 2020; 11(4):394. https://doi.org/10.3390/f11040394

Chicago/Turabian Style

Sprengel, Lars, Zhongqian Cheng, Sandra-Maria Hipler, Shuirong Wu, and Heinrich Spiecker. 2020. "Mean Annual Wood Density Variations of Larix gmelinii (Rupr.), Quercus mongolica Fisch. ex Ledeb., and Pinus tabulaeformis Carr. at Two Different Stem Heights" Forests 11, no. 4: 394. https://doi.org/10.3390/f11040394

APA Style

Sprengel, L., Cheng, Z., Hipler, S. -M., Wu, S., & Spiecker, H. (2020). Mean Annual Wood Density Variations of Larix gmelinii (Rupr.), Quercus mongolica Fisch. ex Ledeb., and Pinus tabulaeformis Carr. at Two Different Stem Heights. Forests, 11(4), 394. https://doi.org/10.3390/f11040394

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