Kinetic analysis of color changes in keyaki (Zelkova serrata) and sugi (Cryptomeria japonica) wood during heat treatment

The kinetics of color changes in keyaki (Zelkova serrata Makino) and sugi (Cryptomeria japonica D. Don) wood during heat treatment were examined. The color of wood specimens treated at 90, 120, 150, and 180 °C was measured by an imaging spectrophotometer and expressed using CIELAB color parameters. At any treatment temperature, values for L* and $ \Updelta E_{ab}^{*} $ decreased and increased in both wood species, respectively, with increased treatment time. Changes in a* and b* varied depending on wood species and treatment temperature. The color changes were successfully analyzed using the kinetic approach applying time–temperature superposition method. This approach elucidated and accurately predicted color changes during heat treatment.     

Causes for the small scale variability of nitrate concentration in seepage water of an N saturated mature spruce stand

Context

In N-saturated forests nitrate concentrations in seepage water ( $ {\text{N}}{{\text{O}}_3}{^{ - }_{\text{seepage}}} $ ) regularly show high spatial variability even within homogeneous stands. Up to now the reasons of this variability are not fully understood.

Aims

The main objective was to identify the crucial parameters that control spatial variability of $ {\text{N}}{{\text{O}}_3}{^{ - }_{\text{seepage}}} $ at the H?glwald site.

Methods

We investigated a multitude of parameters (e.g. N turnover, root biomass, soil chemistry, soil physics, stand parameters) and related them to $ {\text{N}}{{\text{O}}_3}{^{ - }_{\text{seepage}}} $ , measured in 40?cm depth with suction cups.

Results

A small number of biological parameters (net N mineralization, root distribution, and stand density) explained up to 93?% of the variability of $ {\text{N}}{{\text{O}}_3}{^{ - }_{\text{seepage}}} $ in linear regression models. Net N-mineralization rates in the humus layer and fine root biomass in the upper mineral soil influenced $ {\text{N}}{{\text{O}}_3}{^{ - }_{\text{seepage}}} $ positively. Fine root biomass in deeper soil layers (30?C40?cm depth) and stand density had a negative influence.

Conclusion

The rate of net N mineralization in the organic layer is decisive for the nitrate production in the soil. Roots in the upper mineral soil increase $ {\text{N}}{{\text{O}}_3}{^{ - }_{\text{seepage}}} $ by intensive water uptake but excluding nitrate at the same time. The variation of these two parameters is responsible for most of the small-scale variability of $ {\text{N}}{{\text{O}}_3}{^{ - }_{\text{seepage}}} $ .     

Evaluation of the wave attenuation function of a coastal black pine Pinus thunbergii forest using the individual-based dynamic vegetation model SEIB-DGVM

The wave attenuation function of a Japanese black pine forest was evaluated based on its growth at different initial planting densities (P _ini) using the spatially explicit, individual-based, dynamic global vegetation model. The forest dynamics were simulated for 150 years utilizing datasets for tree density and stem diameter at different stand ages obtained in the field. To elucidate the ability of the forest to reduce the wave height $ \eta_{t} $ (m), a long linear wave that propagates on dry ground was assumed. The attenuation of $ \eta_{t} $ (m) was expressed as follows: $ \eta_{t} = \eta_{t0} \exp ( - k_{\text{i}} x) $ , where $ \eta_{t0} $ , x, and k _i are the initial wave height (m), the distance (m), and the wave attenuation coefficient (m^?1), respectively. The tree destruction caused by the waves was considered in order to estimate k _i. The model suggested that there was a peak age that maximized k _i and was dependent on $ \eta_{t} $ , and that the maximum k _i attained decreased with increasing $ \eta_{t} $ . When P _ini was varied widely from 0.5 to 4 m^?2, the maximum k _i for a relatively low wave height (≤3 m) changed dramatically. For example, when $ \eta_{t} = 2{\text{ m}},$ the maximum k _i ranged from 0.008 to 0.031 m^?1, depending on P _ini. Thus, utilizing a relatively low P _ini would be an efficient way of quickly creating a forest capable of sufficient wave attenuation in areas where a relatively high wave height (≥4 m) is expected. It was concluded that regular harvesting and planting would be required to realize the full potential of the coastal forests to attenuate waves, and that tailoring P _ini is one of the management options that could be used to establish a wave prevention forest.     

A new composite k-tree estimator of stem density

This study presents a generally applicable and robust k-tree composite estimator of density. We propose to estimate stem density by a weighted average $ \left( {\hat{\lambda }_{\text{aic}} } \right) $ of 16 individual density estimators. The weights given to individual estimators are inversely proportional to the relative fit (Akaike’s corrected information criterion) of each estimator to the assumed distribution of observed k-tree distances. The performance of the proposed estimator is evaluated in simulated simple random sampling with k?=?3 and 6 in 58 forest stands (54 actual and 4 simulated) and 600 replications. Sample sizes were 15 and 30 locations per stand. Eleven estimators were novel, including three designed for regular spatial patterns. Absolute stand-level bias with k?=?6 varied from 0.1 to 8.1% (mean 1.8%), and a bias larger than 6% was limited to 3 stands with either pronounced density gradients or a strong clustering of stem locations. Root mean squared errors were approximately 16% (k?=?6 and n?=?15) versus 12% for sampling with comparable fixed-area plots. Coverage of computed 95% confidence intervals ranged from 0.72 to 0.99 (median?=?0.98 with n?=?15 and 0.95 with n?=?30), with 98% of all intervals achieving a coverage of 0.85 or better. In seven stands used in an assessment of a novel spatial point pattern reconstruction k-tree density estimator (RDE) by Nothdurft et al. (Can J For Res 40:953–967, 2010), the average absolute bias of $ \hat{\lambda }_{\text{aic}} $ with k?=?6 was 1.5 versus 0.7% for $ \hat{\lambda }_{\text{RDE}} $ .     

Genetics of wood quality attributes in Western Larch

Context

Wood quality traits are important to balance the negative decline of wood quality associated with selection for growth attributes in gymnosperm breeding programs. Obtaining wood quality estimates quickly is crucial for successful incorporation in breeding programs.

Aims

The aims of this paper are to: (1) Estimate genetic and phenotypic correlations between growth and wood quality attributes, (2) Estimate heritability of the studied traits, and (3) Assess the accuracy of in situ non-destructive tools as a representative of actual wood density.

Methods

Wood density (X-ray densitometry), tree height, diameter, volume, resistance drilling, acoustic velocity, and dynamic modulus of elasticity were estimated, along with their genetic parameters, for 1,200, 20-year-old trees from 25 open-pollinated families.

Results

Individual tree level heritabilities for non-destructive evaluation attributes were moderate ( $ {\widehat{h}}_i^2=0.37-0.42 $ ), wood density and growth traits were lower ( $ {\widehat{h}}_i^2=0.23-0.35 $ ). Favorable genetic and phenotypic correlations between growth traits, wood density, and non-destructive evaluation traits were observed. A perfect genetic correlation was found between resistance drilling and wood density (r _G ?=?1.00?±?0.07), while acoustic velocity and dynamic modulus of elasticity showed weaker genetic correlations with wood density (r _G ?=?0.25?±?0.24;?0.46?±?0.21, respectively).

Conclusion

This study confirmed that resistance drilling is a reliable predictor of wood density in western larch, while the weak genetic correlations displayed by acoustic velocity and dynamic modulus of elasticity suggest limited dependability for their use as fast in situ wood density assessment methods in this species.     

Near-infrared spectroscopy prediction of southern pine No. 2 lumber physical and mechanical properties

This study investigated near-infrared spectroscopy (NIRS) to rapidly estimate physical and mechanical properties of No. 2 2 × 4 southern pine lumber. A total of 718 lumber samples were acquired from six mills across the Southeast and destructively tested in bending. From each piece of lumber, a 25-mm-length block was cut and diffuse reflectance NIR spectra were collected from the transverse face using a FOSS 5000 scanning spectrometer. Calibrations were created using partial least squares (PLS) regression and their performance checked with a prediction set. Overall moderate predictive ability was found between NIRS and the properties for the calibration and prediction sets: block specific gravity (SG) (R ² = 0.66 and R _p ²  = 0.63), lumber SG (0.54 and 0.53), modulus of elasticity (MOE) (0.54 and 0.58), and modulus of rupture (MOR) (0.5 and 0.4). Model performance for MOE (R _p ²  = 0.70) and MOR (R _p ²  = 0.50) improved when performing PLS regression on a matrix containing lumber SG and NIR spectra. Overall NIRS predicted MOE better than linear models using lumber SG (R ² = 0.46), whereas lumber SG (R ² = 0.51) predicted MOR better than NIRS. Overall NIRS has reasonably good predictive ability considering the small volume of wood that is scanned with the instrument.     

The toughness contribution of bamboo node to the Mode I interlaminar fracture toughness of bamboo

Bamboo is a kind of biological composite reinforced by unidirectional long fibers. The cleavage strength along grain of bamboo internode is low; however, the existence of bamboo node can hinder the propagation of interlaminar crack to make up for the defect of weak opening mode fracture toughness along interlamination. In this article, the double cantilever beam method was applied to test the Mode I interlaminar fracture toughness of Moso bamboo internode specimens and specimens with node to study the difference of the Mode I interlaminar fracture toughness between Moso bamboo internode specimens and specimens with node. The results are shown as follows: the Mode I interlaminar fracture toughness of Moso bamboo internode specimens was \( G_{{{\text{I}}C}}^{\text{Internode}} \)  = 498 J/m² (SD = 65 J/m²); the Mode I interlaminar fracture toughness of Moso bamboo specimens with node was \( G_{{{\text{I}}C}}^{\text{Node}}\)  = 1,431 J/m² (SD = 198 J/m²). It can be seen that the Mode I interlaminar fracture toughness of bamboo specimens with node was higher than that of bamboo internode specimens, and the toughness contribution of node to bamboo Mode I interlaminar fracture toughness was 1.87 times. The conclusion was drawn that bamboo node can contribute a lot to hinder the interlaminar fracture of bamboo.     

Within-tree variations in the surface free energy of wood assessed by contact angle analysis

Although the heterogeneity of wood could preclude assessment of the performances of adhesion, the within-tree variations in its thermodynamic properties still remain unexplored. This study analyzed the surface free energy of wood (γ _W ) calculated by the geometric mean approach along the radius and at different sampling heights in black spruce and balsam fir with a statistical model, employing apparent contact angles of three testing liquids. Sampling height and species affected the apparent contact angles, and consequently the polar ( $ \gamma_{W}^{P} $ ) and dispersed ( $ \gamma_{W}^{D} $ ) components of γ _W , but no significant difference was observed along the radius. On average, γ _W was 43.77?mJ?m^?2 at the stem base and gradually decreased, attaining 27.19?mJ?m^?2 at 14?m of height. $ \gamma_{W}^{P} $ was markedly higher than $ \gamma_{W}^{D} $ with an opposite pattern along the tree height. These findings demonstrated the huge variability in γ _W along the stem, which could be related to the structural or chemical features of wood.     

Genetic structure of Tunisian natural carob tree (Ceratonia siliqua L.) populations inferred from RAPD markers

? Seven RAPD markers were used to assess the genetic diversity and structure of ten Tunisian natural Ceratonia siliqua L. populations from different geographic and bioclimatic zones.

? The species maintain a high diversity within population as estimated by the percentage of polymorphic loci and Shannon’s index (P% = 76.31, $\bar H'_{pop} = 0.569$ ). The range of variation between populations was large. Populations from the upper semi-arid bioclimates, with more continuous distribution area showed the highest level of variation.

? A high genetic differentiation among populations (Φ_ST = 0.250 and $\bar G_{ST} = 0.347$ ), as a result of population isolation was revealed. Nevertheless, the genetic structure is in accordance with bioclimate indicating that ecological factors also should influence differentiation. Populations from the sub-humid, upper semi-arid and mean semi-arid zones clustered together and were distinct from those of the lower semi-arid ones.

? Conservation strategy should be made according to the level of polymorphism within population and bioclimate.

     

Surface characterization of a series of cured polyurethane adhesives by inverse gas chromatography and adhesion strength testing of plywood

Inverse gas chromatography (IGC) was used to determine the surface properties of a series of polyurethane adhesives, A ₁, A ₂, A ₃ and A ₄ cured with water. The weight percentages of isocyanate group (NCO) in polyurethane adhesives were 5.3, 7.0, 13.0 and 19.6%, respectively. Four n-alkanes, C₆, C₇, C₈ and C₉ were chosen as non-polar probes to characterize the dispersive component of surface free energy, $ \gamma_{s}^{d} . $ Acetone, tetrahydrofuran and ethyl acetate were chosen as polar probes to detect the Lewis acid–base parameters, K _a and K _b. The trend of $ \gamma_{s}^{d} $ and K _a and K _b of the series of cured adhesives was also shown in this paper. The cured adhesives were all amphoteric, but predominantly Lewis basic. The adhesion strength of the poplar plywood bonded with polyurethane adhesives was tested. The results showed that the adhesion strength of plywood increased with increasing the NCO content of the adhesives and K _b/K _a of the cured adhesives.     

Allometric equations for estimating aboveground biomass of Coffea arabica L. grown in the Rift Valley escarpment of Ethiopia

Coffee, Coffea arabica L., which is native to Ethiopia, is the world’s most widely traded tropical agricultural commodity. While much is known about the productivity and management of coffee for coffee beans little attention has been given to the plants overall biomass production and carbon sequestration. The objective of this study was to develop and evaluate allometric equations for estimating the aboveground biomass of C. arabica plants growing in indigenous agroforestry system in the Rift Valley escarpment of south-eastern Ethiopia. Coffee plays an important role in providing income and in sustaining these productive systems. Biomass harvesting of 31 plants with 54 stems was carried out in a 40 km² area varying in elevation from 1,500 to 1,900 m. The stem accounted for most (56 %) of plant biomass, followed by branches (39 %) and twigs plus foliage (5 %). Plant mean biomass was 22.9 ± 15.8 kg. Power equations using stem diameter measured at either 40 cm (d ₄₀) or at breast height (d, 1.3 m) with and without stem height (h) were evaluated. The square power equation, $ Y \; = \; b_{ 1} d_{ 40}^{ 2} $ , was found to be the best (highest ranked using goodness-of-fit statistics) for predicting total and component biomass. The reliability of the prediction decreased in the order: stem > branches > twigs plus foliage. A cross-validation procedure showed that equation parameterization was stable and coefficients reliable. Our parameterized square power equation for total aboveground biomass was also found to be better than the equations parameterized by Hairiah et al. (Carbon stocks of tropical land use systems as part of the global C balance: effects of forest conversion and options for clean development activities, International Centre for Research in Agroforestry, Bogor, 2001) and Segura et al. (Agroforest Syst 68:143–150, 2006) for C. arabica grown in agroforestry systems, confirming the importance of parameterization of allometric equations with site specific data when possible.     

A simple theory to link bole surface area,stem density and average tree dimensions in a forest stand

A geometrical model of a forest stand has been analyzed. A forest stand has been modeled as a population of cones which was described by the change of total bole surface area with density \(\hat{S}(N)\) , relation between density and a horizontal dimension (radius r) r(N), and the relation between vertical dimension (generatrix l) and radius l(r). It has been shown that there are close relationships between \(\hat{S}(N)\) , l(r) and r(N). In case of \(\hat{S}(N) = const\) , power exponent of l(r) can be predicted from the power exponent of r(N) and vice versa. A comparison of the model analysis with the data available on Scots pine (Pinus sylvestris L.) stands has been performed. In spite of the model simplicity, its inferences proved to be workable in many cases where the data can be interpreted as a dynamics of an even-aged forest stand. In particular, if the estimation of total bole surface area is constant, the power exponent in the relation of diameter and stand density DBH(SD) can be calculated on the basis of the power exponent in the relation of height and diameter H(DBH) and vice versa. Possible limitations and the meaning of the analysis are discussed.     

Ranking larch genotypes with the Rigidimeter: relationships between modulus of elasticity of standing trees and of sawn timber

? Direct assessment of modulus of elasticity (MOE) on standing trees is attractive for breeders to evaluate genotypes prior to selection: this can be done using the Rigidimeter, a bending-based measurement device.

? In this study, we tested its reliability to properly rank genotypes by relating trunk MOE with MOEs estimated with a vibrating analysis system (Bing) on different types of conditioned wood specimens from the same trees (boards and standardised 2×2×30 cm-clear-wood specimens). One hundred and ten trees from different genotypes of hybrid larch (Larix × eurolepis) were tested.

? Mean trunk MOE was 7 300 MPa with a similar value obtained for sawn boards. Clear-wood specimens MOE increased from pith to bark from less than 6 000 MPa to nearly 9 000 MPa. Moderate correlations (r = 0.48–0.61) were found at the individual tree level between trunk MOE and MOE of wood samples.

? Single specimen MOE was shown to be strongly related to a linear combination of trunk MOE and sample position.

? At the genotype mean level, trunk MOE was highly correlated with wood samples MOE (r = 0.80–0.91). Ranking of genotypes based on trunk MOE was mostly consistent with that based on standardised specimens.

? It was concluded that besides other operational advantages which are discussed, the Rigidimeter is a valuable tool for breeders to routinely evaluate and rank genotypes for stiffness prior to further selection.

     

Performance of Pinus elliottii, Pinus caribaea, their F1, F2 and backcross hybrids and Pinus taeda to 10 years in the Mesopotamia region, Argentina

Performances of Pinus taxa were studied to 10 years of age in two trials in each of Misiones and Entre Ríos provinces across the Mesopotamia region of Argentina. Taxa comprised 22 populations from sources in Argentina, Australia, Brazil and Zimbabwe including Pinus elliottii var. elliottii (Pee), Pinus caribaea var. hondurensis (Pch), their four, inter-specific hybrids (F₁, F₂ and backcrosses from F₁ to Pch and to Pee—all as broadly based bulks); other Pee and Pinus taeda (Pt) comprised narrower or unspecified bulks. Variable numbers of taxa were missing at each site. Mean survival across sites at age 10 years ranged 53.2–91.3% averaging 74.2%. Analysis of variance of plot means indicated population effect was statistically significant (p < 0.05) for all or most growth and quality traits at all sites. However, significant differences from the nominated check seedlot at the Entre Ríos sites (Pee, Australia) were extremely rare, while quite common at the northern, Misiones sites (check seedlot a Pt population). In the Misiones trials, F₁, F₂ and both backcross hybrids showed better stem straightness than Pee and Pt from Argentina, generally with statistically significant differences (p < 0.05). Pt showed lowest forking scores (desirable). Taxon × environment interaction was statistically significant (p < 0.01) for growth traits only (p > 0.05). However, this interaction contributed an average of only 34.1% of the taxon variance suggesting a lack of practical importance. Taxa most suitable for deployment in the Mesopotamia region, Argentina are suggested.     

Evaluation of the visual stress grading standard on French Spruce (Picea excelsa) and Douglas-fir (Pseudotsuga menziesii) sawn timber

? In this paper an evaluation of the visual grading standard for softwood sawn timber was made.

? In order to do so, visual grading according to EN 518 and theoretical grading according to EN 338 and EN 384 (measurements of MOE, MOR and density) were applied to lumber. Two batches of 111 and 102 French boards were graded, respectively, of Spruce (Picea excelsa) and Douglas-fir (Pseudotsuga menziesii). For the visual grading the most discriminant criterion was noted: knots, cracks, wane, etc.

? Finally, the results of the two grading methods were compared, and it was shown that the visual stress grading gave quite low results for our two French species.

     

Within-stem variations in mechanical properties of <Emphasis Type="Italic">Melia azedarach</Emphasis> planted in northern Vietnam

Within-stem variations in the mechanical properties of 17–19-year-old Melia azedarach planted in two sites in northern Vietnam were examined by destructive and nondestructive methods. Wood samples were collected from 10, 50, and 90% of the radial length from pith on both sides (North and South) at 0.3, 1.3, 3.3, 5.3, and 7.3 m heights above the ground. The mean values in whole trees of wood density (WD), modulus of rupture (MOR), modulus of elasticity (MOE), and dynamic modulus of elasticity (E_d) at 12% moisture content were 0.51 g/cm³, 78.58 MPa, 9.26 GPa, and 10.93 GPa, respectively. Within the stem, the radial position was a highly (p?<?0.001) significant source of variation in mechanical properties. MOR, MOE, and E_d increased from pith to bark. WD had a strong positive linear relationship with both MOR (r?=?0.85, p?<?0.001) and MOE (r?=?0.73, p?<?0.001). This suggests that it is potentially possible to improve mechanical properties through controlling WD. MOR had also a strong linear relationship with E_d (r?=?0.84, p?<?0.001). This indicates that E_d is a good indicator to predicting the strength of wood if the density of measured element is known. Besides, the stress wave method used in this study provides relatively accurate information for determining the stiffness of Melia azedarach planted in northern Vietnam.     

Bayesian inference for multi-environment spatial individual-tree models with additive and full-sib family genetic effects for large forest genetic trials

Context

The gain in accuracy of breeding values with the use of single trial spatial analysis is well known in forestry. However, spatial analyses methodology for single forest genetic trials must be adapted for use with combined analyses of forest genetic trials across sites.

Aims

This paper extends a methodology for spatial analysis of single forest genetic trial to a multi-environment trial (MET) setting.

Methods

A two-stage spatial MET approach using an individual-tree model with additive and full-sib family genetic effects was developed. Dispersion parameters were estimated using Bayesian techniques via Gibbs sampling. The procedure is illustrated using height growth data at age 10 from eight large Tsuga heterophylla (Raf.) Sarg. second-generation full-sib progeny trials from two series established across seven sites in British Columbia (Canada) and on one in Washington (USA).

Results

The proposed multi-environment spatial mixed model displayed a consistent reduction of the posterior mean and an increase in the precision of error variances $ \left( {\sigma _{e}^{2}} \right) $ than the model with ??sets in replicates?? or incomplete block alpha designs. Also, the multi-environment spatial model provided an average increase in the posterior means of the narrow- and broad-sense individual-tree heritabilities (h _N ² and h _B ² , respectively). No consistent changes were observed in the posterior means of additive genetic correlations (r _Ajj??).

Conclusion

Although computationally demanding, all dispersion parameters were successfully estimated from the proposed multi-environment spatial individual-tree model using Bayesian techniques via Gibbs sampling. The proposed two-stage spatial MET approach produced better results than the commonly used nonspatial MET analysis.     

Spectral image database for observing the quality of Nordic sawn timbers

In this study, a comprehensive spectral image database of Nordic sawn timbers for public use was measured. Economically significant Finnish wood species birch (Betula sp.), Norway spruce (Picea abies) and Scots pine (Pinus sylvestris) were chosen for inclusion in the database. The total of samples was 107 containing heartwood, sapwood, decayed wood, blue stain, mold, resin, early wood, late wood, knots, cracks, pith, reaction wood and bark. Board and crosscut samples were measured in frozen, melted and room-dried conditions. The reflectance of samples was measured over a 300- to 2,500-nm wavelength range. Additionally, the photoluminescence of samples excited by an ultraviolet B light source was measured. The spot size used was 250 μm with an 80 mm \(\times \) 200 mm imaging area, and produced all in all 44 million spectra. In this paper, examples of the possibilities of this spectral image database as a means of detection of the spatial distribution of aromatic lignin and the moisture content (MC) of nonfrozen timber were introduced and provided. From the results, it was found that it was possible to detect the lignin distribution from spectral images, and simple and robust methods for wood MC estimations were also introduced.     

Application of VIS/NIR spectroscopy for estimating chemical,physical and mechanical properties of cork stoppers

The aim of this study was to evaluate the potential of visible and near infrared spectroscopy (Vis/NIRS) in predicting the chemical, physical and mechanical behavior of single-piece natural corks stoppers used for sealing wine bottles. Two training sets of 90 and 150 cork stoppers were used to obtain four spectra per sample in different positions: two of the stopper bases (transversal section) and two of the stopper sides (tangential section and radial section). The samples were scanned in the range of 400–2,500 nm using a Foss-NIRSystems 6500 SY II spectrophotometer equipped with a remote reflectance fiber-optic probe. On each training set, two-thirds of the samples were used to develop modified partial least square (MPLS) calibration equations, and the remaining one-third of the sample for the external validation of these MPLS equations. The best equations were obtained for the transversal section, which is the recommended one when applying Vis/NIRS technology to cork. The best results for the chemical composition were obtained for waxes and total polyphenols, showing coefficient of determination of the cross validation (r _cv ² ) values of 0.64 and 0.56 and coefficient of determination of the external validation (r _EV ² ) values of 0.53 and 0.55, respectively. The best equation for the physical and mechanical parameters was obtained for moisture content (r _cv ²  = 0.86 and r _EV ²  = 0.85), with somewhat lower results for density, compression force and extraction force (r _cv ²  = 0.66, 0.72, 0.52 and r _EV ²  = 0.52, 0.49, 0.51, respectively). The SECV (standard error of cross validation) and SEP (standard error of external validation) were similar for all the physical and mechanical parameters, thus confirming the robustness of the equations. MPLS model for moisture content fulfills the requirements for screening (RPD >2.5), but MPLS models obtained for waxes, total polyphenols, density, compression force and extraction force are not good enough for routine analysis or quality control. The results obtained from the MPLS models based on Vis/NIRS technology would permit the continuous quality control of humidity in the production line as well as obtaining information about certain chemical components (extractives contents) and some physical and mechanical parameters (density, extraction force and compression force).