Forests 2022, 13, x. https://doi.org/10.3390/xxxxx www.mdpi.com/journal/forests Article Influence of Timber Moisture Content on the Ultrasonic Wave Velocity Measurement of Tectona grandis L.F. and Cupressus lusitanica M from Costa Rica Viviana Paniagua 1,*, Javier Corrales 2, Cindy Torres 2,3 and Beatriz González-Rodrigo 4,* 1 Architecture Department, University of Costa Rica; viviana.paniaguahernandez@ucr.ac.cr 2 Chemical Engineering Department, University of Costa Rica; javier.corralesvargas@ucr.ac.cr (J.C.); cindy.torres@ucr.ac.cr (C.T.) 3 Material Science and Engineering Research Center 4 Civil Engineering School, Universidad Politécnica de Madrid, Spain * Correspondence: viviana.paniaguahernandez@ucr.ac.cr (V.P.); beatriz.gonzalez.rodrigo@upm.es (B.G.-R.) Abstract: Ultrasound is a reliable non-destructive method commonly used to evaluate the state of a piece of wood. The effect of the moisture content (MC) on the timber wave velocity, which is differ- ent depending on the wood species, must be considered. MC adjustment factors are important if accurate and comparable results are to be achieved. Thus, the goal of this study is to define a model that allows obtaining the adjustment factor to predict the standard velocity (12% of MC (V12) from velocity obtained at different MC for two Costa Rican commercial species: Tectona grandis L.F (teak) and Cupressus lusitanica M. (cypress). This effect was studied on small clean specimens during the desorption stages, from the fiber saturation point to the oven-dry state, controlling the specimen’s mass and MC on 62 specimens. With this data, the rate of change in ultrasound velocity per MC was modeled. Thus, the applicability of already published moisture adjustment models for conifers and hardwood tropical species was proved. The results showed that the proposed model coefficients adjust better than the ones obtained from the wood science literature, which makes them suitable to describe ultrasound velocity in different moisture conditions (VH). Keywords: cypress; teak; dynamic modulus of elasticity; moisture content; ultrasonic velocity 1. Introduction Ultrasonic testing is one of the preferred non-destructive testing (NDT) methods for estimating the anisotropic elastic properties of wood materials, like young’s modulus, which is one of the most important variables for wood. This NDT is a suitable tool to analyze the wood’s mechanical behavior due to lower costs of equipment and staff, com- bining speed, ease, and accuracy in the results acquired (Mascarenhas et al., 2020). This method estimates the dynamic modulus of elasticity (MOEdyn) using the measurement of wave velocity and bulk density (ρ) of wood pieces (Bachtiar et al., 2015). Therefore, con- sidering that MOEdyn exhibits a nearly perfect linear relationship with the static modulus (MOEapp) (Chauhan and Sethy, 2016), ultrasound testing represents a valuable alternative for wood testing. It is used effectively to predict lumber grade and strength properties, assess heritage structural elements, detect delamination, and monitor the in-use perfor- mance of timber during service life (Casado et al., 2012; Marhenke et al., 2018; Neuenschwander et al., 2013; Sobue, 1993; Wang et al., 2003; Wang, 2000; Yamasaki et al., 2017). Additionally, it could be used to determine outdoor degradation or to evaluate wood structures subject to the constant action of weathering and high humidity, such as outdoor decks and other structures (Mascarenhas et al., 2020). Citation: Paniagua, V.; Corrales, J.; Torres, C.; González-Rodrigo, B. Influence of Timber Moisture Content on the Ultrasonic Wave Velocity Measurement of Tectona grandis L.F. and Cupressus lusitanica M from Costa Rica. Forests 2022, 13, x. https://doi.org/10.3390/xxxxx Academic Editor(s):Joris Van Acker Received: 7 July 2022 Accepted: 5 August 2022 Published: date Publisher’s Note: MDPI stays neu- tral with regard to jurisdictional claims in published maps and institu- tional affiliations. Copyright: © 2022 by the authors. Submitted for possible open access publication under the terms and con- ditions of the Creative Commons At- tribution (CC BY) license (https://creativecommons.org/licen- ses/by/4.0/). Forests 2022, 13, x FOR PEER REVIEW 2 of 11 Empirical correction factors have been developed to predict standard velocity (12% of MC (V12)) from the obtained velocity at different MC (VH) (Sanabria et al., 2011; Sandoz, 1989a). For these models, it must be considered that the VH is influenced when timber contains only bound water, below the FSP, or when it also contains free water, above FSP (Sandoz, 1993). This author estimated that MC influence is 8 times stronger below the FSP. Hence, the relation VH and MC can be represented using a bilinear model with a change of the slope at the FSP (Sandoz, 1993; Yamasaki et al., 2017). The pattern of VH vs. average MC during desorption varies within species (Mishiro, 1996a, 1996b). Consequently, several researchers have performed corresponding adjust- ment factors for different species and moisture levels (Gerhards, 1982; Wang, 2008; Unterwieser and Schickhofer, 2011a; Nocetti et al., 2015; Gonçalves et al., 2017). Gonçalves, et al. (2017) tested three different tropical hardwoods, obtaining three different slopes (45, 55, and 66) ms−1. In general, most of the studies on timber ultrasound wave velocity that are available in the literature are focused on softwood species, although hardwood species have been studied and incorporated into Brazilian standards (NBR 15521, 2007). There is no evidence of the existence of a higher influence of MC and the ultrasound velocity comparing hard- woods and softwood species. While Llana et al. (2020) observed that hardwoods pre- sented a lower slope when moisture decreased, other authors have compared the slope of the ultrasound velocity with changes in the MC, obtaining almost the same slope for soft- wood and hardwoods. Goncalves, R. and Leme da Costa (2008) found almost the same variation of ultrasound velocity with MC, 25 m s−1 for softwood and 24 ms−1 for hardwood. Similar values were obtained by Saadat-Nia et al. (2011), respectively, (33 and 31) ms−1. However, Peng et al. (2016) reported higher differences among three hardwoods from (24 to 60) ms−1 than among three softwoods from (36 to 42) ms−1. In Costa Rica, many efforts have been made to increase structural wood consump- tion, which is promoted by the country’s commitments to be carbon-neutral. Different studies have focused on characterizing native species of timber proceeding from local plantations (Moya-Roque et al., 2010). Wood degradation in tropical environments and acid atmospheres from volcanic emanations have also been studied. The widespread use of wood in this environment is critical for the expansion of commercial utilization of these species in volcanic areas. The detailed methodology proposed for the evaluation of this degradation in tropical areas has been analyzed (Paniagua-Hernández et al., 2021). Con- sidering that the ultrasound method allows determining the timber degradation in the same piece over time, the MC variations must be considered to adjust VH to compare field data to continue that investigation. The species of interest Tectona grandis L.F (teak) and Cupressus lusitanica M. (cypress), two representative species commonly used in construction in Costa Rica, proceed from local forest plantations. As a hardwood, young teak is the most commercialized in Costa Rica and the one that has been exported for the longest time, whereas cypress is the most widespread co- niferous plant in Costa Rica (Barrantes, Ugalde, 2019). The use of these wood species as structural material raises the need for standard parameters to determine the degradation of the wood over time and in different moisture scenarios. Regarding teak, studies of its MOEdyn are available (Chauhan and Sethy, 2016; Ilic, 2003; Karlinasari et al., 2018; Mascarenhas et al., 2020; Solorzano et al., 2012) but no evi- dence has been found on the correlation factors between VH and MC, nor have they been found for cypress. The goal of this study is to determine appropriate moisture adjustment functions that allow comparing mechanical properties considering those obtained at different MCs for these two Costa Rican commercial species. This research contributes to the expansion of their commercial use and the increase of the knowledge of best practices databases in the country. Forests 2022, 13, x FOR PEER REVIEW 3 of 11 2. Materials and Methods 2.1. Wood Material This investigation was carried out with 62 small clean specimens (30 teak, 32 cypress). The samples were supplied in green condition (freshly cut), as solid wood from planta- tions obtained from local sawmills located in La Rita, Guápiles, Limón, (teak) and Paso Llano, San José de la Montaña, Heredia (cypress). The plantation’s ages were 14 (teak) and 40 years (cypress), respectively, due to the maturity representation at which they are traded locally. The specimens dimensions were 1 × 1 × 16 (25(b) × 25(h) × 410(l) mm), ac- cording to (ASTM-D143, 2014). 2.2. General Procedure The specimens were taken to the laboratory with a MC greater than FSP. Then, they were gently dried, and the initial ultrasonic wave velocity (considered the first valid meas- urement) was taken when the species reached 28% ± 1% MC  the FSP. In literature, values between 28% and 32% as the FSP have been used by Sandoz (1993) and Steiger (1995b). In this study, the FSP was considered as recommended by Unterwieser and Schickhofer (2011) at 28% for ultrasonic runtime measurement, and not 32%, which is recommended for the natural frequency measurement method. Specimens were subjected to systematic monitoring during the desorption stages from FSP to 5% of MC. They were kiln-dried at a dry-bulb temperature of 103 C and measurements of ToF, weight, and dimensions were taken every hour until they reached 5% of moisture, with at least 6 measurements for each specimen. Temperature control was considered, which means cooling the specimens after coming out of the oven in order not to affect the mass measurement. The bulk density ρ was calculated from weight and vol- ume. The longitudinal wave propagation velocity of each specimen was measured repeat- edly, three times, and, then averaged, as recommended by Yamasaki et al. (2017). To avoid errors in the readings due to the pressure in the transducers, the recommendations by Carrasco et al. (2021) were considered. The MC was obtained according to the secondary oven drying (ASTM-D4442-16, 2016). From each specimen, a 5-cm sample was removed. The sample piece’s mass meas- ured wet and dried allows obtaining the moisture content of the test specimen at various test stages (%) by applying Equation (1). 𝑀𝐶 = [(𝑚𝑢 − 𝑚0)/𝑚0] × 100 (1) where MC = moisture content of the test specimen at various test stages (%); mu = test specimen mass at various test stages (kg), and m0 = the oven-dried mass of the specimen (kg). Bulk density (kgm−3) was determined. 2.3. Wave Velocity Measurement Procedure Wave propagation time measurements were performed each hour following the de- scribed method. Ultrasonic velocity VH was measured with a portable ultrasound device (Sylvatest duo-Lausanne, Switzerland), and based on the determination of the time of flight (ToF) of the ultrasonic wave between transmitter and receiver, which is generated by piezoelectric transducers with a frequency of about 54 kHz. Measurements were ob- tained in a direct contact position with the transducers placed facing each other on oppo- site edges of specimens at the center of the section in the longitudinal direction (Figure 1) since this is the most consistent parameter for the global evaluation of structural members. Three consecutive measurements were taken to verify each value, using the average. The ultrasound velocity was estimated from Equation (2). 𝑉𝐻 = 𝐿/𝑇𝑜𝐹 (2) where VH = ultrasound velocity (ms−1), L = specimen’s length (m), ToF = time of flight (s). Forests 2022, 13, x FOR PEER REVIEW 4 of 11 The final consideration is the constant temperature (room temperature 20 C) in the functional range of the model since according to Kang and Booker (2002), the frequency spectrum of the transmitted signal is lightly affected by temperature. This is also sup- ported by Bucur et al. (2006), who found a dominant effect of MC over temperature on wave propagation by studying the differences in the trends of VH across the moisture range, along with the higher decreasing rates per unit moisture. Figure 1. Representation of ultrasound velocity TOF measurement in the longitudinal direction and explaining wood structure in different scales to visualize the phenomenon of wave propagation and the affectation by water molecules. Central figure adapted from (Nishimura et al., 2018). The percentage of the velocity variation for every 1% MC variation, VH from the V12 relation is presented in Equation (3) (Arriaga et al., 2017; Gonçalves et al., 2017; Sandoz, 1993; Steiger, 1995; Unterwieser and Schickhofer, 2011) 𝑉𝐻/𝑉12 = 1 + 𝑘𝑢 (𝑀𝐶 − 12) (3) where V12 = velocity at 12% MC, VH = velocity obtained in wood with MC below FSP, ku = adjustment factors for each timber species. 2.4. Statistical Analysis The VH/V12 ratio was obtained for each of the measurements made. A Strubb analysis was performed at p = 0.05, 95% confidence intervals. All the data were studied using a non-linear regression analysis for each specie and the mean squared error and coefficient of determina- tion were obtained. A nonlinear regression (nlinfit) was performed using Matlab® software. Table 1. Symbols in the present paper. MC Moisture content (%) FSP Fiber saturation point (at 28% MC) VH Ultrasound velocity at certain moisture conditions (ms−1) VH/V12 The ratio of ultrasound velocity at certain MC (VH) to ultrasound velocity at 12% MC (V12) change rate V12 Ultrasound velocity at 12% MC (ms−1) Vfsp Ultrasound velocity at FSP 28% MC (ms−1) ΔU Moisture content variation ρ Density (kgm−3) ku Adjustment factors for each timber specie ToF Time of flight MOEdyn Dynamic modulus of elasticity (GPa) Forests 2022, 13, x FOR PEER REVIEW 5 of 11 ρV2 MOEdyn by calculation of VH and density (GPa) NDT Non-destructive testing 3. Results and Discussion The set of all determined properties is shown in Table 2. It must be considered that teak specimens come from the thinning of forest plantations, as is common in the Costa Rican internal market, since the wood from the final felling is exported. On the other hand, the cypress consumed in this country usually comes from older plantations. Therefore, it should be considered that in this study, teak is mainly sapwood, and for cypress, even though the percentage of heartwood is considerable (because of its age), it is difficult to visually determine the percentage because it does not show a noticeable color change be- tween sapwood and heartwood. This explains the fact that both species tested in this stu- dio have very similar densities, although teak has higher V12 than cypress, related with a higher modulus of elasticity (Table 2). Costa Rican teak timber increases in density with age, as reflected in literature, from values of 580 kgm−3 (17 years) to 610 kgm−3 (28 years) (Castro and Raigosa, 2000). The average density for mature teak wood is 800 kgm-3 and 600 kgm−3 for cypress (Moya-Roque et al., 2010). Considering the aforementioned charac- teristics, the obtained values are in the expected range with a relative variation of less than 10% among each species replication. Table 2. Set of all determining properties at 12% MC (ρ, Vfsp, V12) for cypress and teak. Species/Age Density ρ (kgm−3) Vfsp (ms−1) V12 (ms−1) Cypress/40 years 533.37 ± 52.73 3860.76 ± 397.38 4655.16 ± 400.27 Teak/14 years 531.30 ± 49.22 4463.80 ± 360.69 5152.39 ± 371.37 Figures 2 and 3 show the relation between MC and ultrasound velocity from FSP to oven dry conditions for teak and cypress, respectively. These scatter plots found that ul- trasound velocity increases as MC decreases, showing an inverse relationship with a linear tendency. Then, a linear adjustment function describes the correlation between the varia- bles below FSP, which was calculated considering the overall values presented by speci- mens (Table 2). Figure 2. Correlation between MC and relative ultrasound velocity measurement (5%–28%) trend line for Teak. The different slope ku were determined for both species. Using Equation (3), the cor- rection factors for some NDT results to a standard 12% MC were calculated. The proposed model for teak is presented in Equations (4) and (5) describes the proposed cypress model. Forests 2022, 13, x FOR PEER REVIEW 6 of 11 𝑉𝐻/𝑉12 = 1 − 0.0083 ∙ (𝑀𝐶 − 12) (4) 𝑉𝐻/𝑉12 = 1 − 0.0106 ∙ (𝑀𝐶 − 12) (5) Figure 3. Correlation between MC and relative ultrasound (5%–28%) trend line for Cypress. The experimental adjustment parameters are shown in Equations (4) and (5). These were found using the nlinfit tool, where the factors were obtained using least squares with a 95% confidence interval, considering that the ku parameter was k1 = 1. The performance of the model found was evaluated with two statistical parameters: the mean squared error (MSE) and the coefficient of determination (R2). In the case of the model for the cypress, it presented an R2 of 0.8480 and an MSE of 6.33 × 10−4, while for teak it showed an R2 of 0.8486 and an MSE of 4.15 × 10−4. Therefore, the results of the relation- ship between ultrasound velocity and MC confirm the reliability of the coefficients for prediction purposes following the literature on the subject (Sandoz, 1990 and 1993; Gon- çalves, et al., 2017; Íñiguez et al. 2015; Unterwieser and Schickhofer 2011). A Matlab routine was programmed for testing the performance of available models in the literature and the developed models with the experimental datasets obtained in this study. Models for hardwoods were tested with the teak dataset, while the cypress dataset was used to test softwood models. Table 3 shows the statistical performance of the different models for the experimental data, distinguishing models for hardwoods and softwoods. In Figures 4 and 5, the adjust- ment according to the different parameters can be observed graphically applied to the experimental data for teak and cypress. Table 3. Adjustment (R2 and MSE) of the proposed models versus literature models tested with the experimental teak and cypress data for the estimation of references velocity at 12% MC (V12) from velocity (VH) with different moisture content. Equation Specie Author ku R2 MSE Hardwood models/Dataset teak (14 years) T Teak (Proposed model) 0.0083 0.8486 4.14 × 10−4 H1 Eucalyptus grandis Gonçalves, et al. (2017) 0.0066 0.7893 5.78 × 10−4 H2 Corymbia citridora Gonçalves, et al. (2017) 0.0083 0.8485 4.16 × 10−4 H3 Eucalyptus pellita Gonçalves, et al. (2017) 0.0067 0.7893 5.78 × 10−4 Softwood models/Dataset cypress (40 years) C Cypress (Proposed model) 0.0106 0.8479 6.34 × 10−4 Forests 2022, 13, x FOR PEER REVIEW 7 of 11 C1 - Íñiguez, et al. (2015) 0.0100 0.8449 6.47 × 10−4 C2 Fir and Spruce Sandoz (1990, 1993) 0.0053 0.5052 0.0021 C3 Spruce and Scots pine Unterwieser Schickho- fer (2011) 0.0060 0.5914 0.0017 ku = adjustment factors for each timber species. Note: all Equations are for 5% < MC < 28% except C1 for 5% < MC < 30% and C2 ≤ 32%). In the equations proposed by different authors to infer V12 from a VH estimated in MC below FSP, the ku varied from 0.53% to 1%. In this research, values of ku were 0.83% for teak (Equations (4) and (T) in Table 3), 1.06% for cypress (Equations (5) and (C) also in Table 3). These values are within the interval obtained by the authors cited in Table 3, using a VH derived for MC below FSP. The obtained results show a suitable fit in the models using the developed parame- ters, with high determination coefficients and low mean square errors, which means a good statistical fit in the order of other studies in the literature, taking into account that they were fitted with small clear specimens. It should be considered that MC influence is not totally independent of dimensions. According to Llana et al. (2020) dimensions can also influence the measurement of MC during drying studies. The MC inside small clear specimens is more homogenous than it is in structural sizes, where a MC gradient can exist from the inside to the outside, and slenderness (length/width) may also affect low- frequency ultrasound devices (from 23 to 54 kHz) (Llana et al., 2020). Assuming the previous considerations, the studied teak shows coefficients in the or- der of those obtained for larger pieces (metric scale) of Corymbia citriodora, being different from those of Eucalyptus grandis and Eucalyptus pellita. As is summarized in Table 3, teak experimental data were evaluated with the models for these three hardwoods species (H1), (H2), and (H3) developed by (Gonçalves et al., 2017). It must be considered that these species presented greater average densities (E. grandis (508 kg m−3), E. pellita (737 kg m−3), C. citriodora (813 kg m−3)) than teak (531 kg m−3). Thus, it should be considered that E. grandis and E. pellita, although different in den- sity, present similar behavior related to velocity variation as a function of MC, whereas C. citriodora presents a greater variation, with a slope similar to the tested teak (Figure 4). Therefore, it should be considered that variation in ultrasound velocity by moisture influ- ence is intrinsic to the species. Figure 4. Comparison of adjustment equations to obtain longitudinal references velocity at 12% MC (V12) for hardwood with the experimental data teak. Forests 2022, 13, x FOR PEER REVIEW 8 of 11 Similar behavior is observed when the experimental data for cypress were analyzed with softwood models (C1), (C2), and (C3), with densities around 450 kgm−3, which are lower than the 533 kgm−3 obtained for cypress. It was observed that the cypress proposed model (C) is similar to (C1) proposed by (Íñiguez et al., 2015), which was used for four Spanish (radiata, scots, salzmann, and maritime) pines (Llana et al., 2020). It is interesting to have found similarities only with one of the three models, as is shown in Table 3. To understand the cause of this, additional research is required. Considering the spe- cie’s anatomical and chemical characteristics. Such as porosity, the geometry of its trache- ids or vessels (arrangement, size, frequency, contents), other factors associated with the composition of its fibers (cellulose, hemicellulose, and lignin) and hence its hydration po- tential, and climatic factors influencing variations in the species characteristics should also be considered (Ferreira, 2013). Figure 5. Comparison of adjustment equations to obtain longitudinal references velocity at 12% MC (V12) for softwood with the experimental data cypress. For the two species studied, the ultrasound velocity increased as the MC decreased. Using the relation 1% of MC decrease, the VH increased 43.9 ms−1 for teak and 50.5 ms−1 for cypress. These data contrast with what was obtained in the study carried out by Llana et al. (2020), where it was determined that the ultrasound velocity was influenced by the MC more in hardwoods than in softwoods. It must be considered that cypress has ranges of velocity variation (ToF) of 1466.45 ms−1 and moisture content variation (ΔU) of 23.69%, whereas teak has ranges of 1464.98 ms−1 and 24.46%, respectively. The species tested, teak (531.30 kgm−3), and cypress (533.37 kgm−3), although pos- sessing similar densities, presented a different behavior related to the ultrasound velocity variation as a function of MC, with cypress showing a greater variation (slope), which is attributable to the higher percentage of the heartwood due to its age. Thus, further re- search can deepen in the relationship between volumetric contractions (shrinkage and swelling) and the anatomy of wood cells and their relationships with bound water move- ment. The contribution of this work lies in highlighting the need to evaluate with experi- mental data the models available in the literature. As shown in this research, the variation between models is substantial, and it is the reason why some models fit better than others. Therefore, empirical models and comparisons are necessary. Forests 2022, 13, x FOR PEER REVIEW 9 of 11 4. Conclusions The present research provides adjustment functions to obtain the velocity at 12% MC of Costa Rican young teak (14 years) and mature cypress (40 years) timber at different MCs employing ultrasound measurements. The coefficient between the velocities ob- tained at different MC and the reference velocity (V12) shows a linear relationship with MC, with a determination coefficient above 85% for both species studied. During FSP desorption at the oven-dry condition, the ultrasound velocity in the lon- gitudinal direction can be corrected, for every 1% in MC by 43.9 ms−1 for teak and 50.5 ms−1 for cypress. The influence of MC on the ultrasound velocity can be considered slightly higher for cypress (conifer) than for teak (hardwood). Experimental data on teak and cypress have been used to test proposed models, along with those previously published by researchers for coniferous and hardwood spe- cies. Teak model presents similar adjustment factors to the hardwood model for Brazilian Corymbia citrodora and the cypress model showed agreement with the model of Spanish pines wood. These models can be applied to moisture correction in NDT for teak and cypress in Costa Rica. Obtained results open the way for monitoring the in-use performance of these timber over its lifespan and determining outdoor degradation, as it is being done in ex- treme acid environments. It will contribute to the expansion of their commercial use and knowledge of best practice databases in the country. Author Contributions: Funding: Acknowledgments: Roger Moya-Roque, Forestry Innovation Research Center, Instituto Tecno- lógico de Costa Rica (TEC), Cartago, Costa Rica, for his advice in NDT methodology. Guillermo Íñiguez-González from the Timber Construction Research Group, School of Forest Engineering and Natural Resources, Universidad Politécnica de Madrid, Spain for the advice in literature and meth- ods. National Laboratory of Materials and Structural Models (LANAMME) of Universidad de Costa Rica. San José, Costa Rica, for their support in the destructive tests. Conflicts of Interest: The authors declare no conflict of interest. References 1. Arriaga, F.; Llana, D.F.; Esteban, M.; Íñiguez-González, G. Influence of length and sensor positioning on acoustic time-of-flight (ToF) measurement in structural timber. Holzforschung 2017, 71, 713–723. https://doi.org/https://doi.org/10.1515/hf-2016-0214. 2. ASTM-D143. ASTM-D143 Standard Test Methods for Small Clear Specimens of Timber. ASTM International: West Conshohocken, PA, USA, 2014. 3. ASTM-D4442-16. Standard Test Methods For Direct Moisture Content Measurement of Wood And Wood-Based Materials. ASTM International: West Conshohocken, PA, USA, 2016. 4. Bachtiar, E.V.; Sanabria, S.J.; Niemz, P. Elastic Characteristics of Wood by Means of Ultrasonic Waves and Mechanical Test. Dans FPInnovations (Éd.). In Proceedings of the 5th ISCHP International Scientific Conference on Hardwood Processing, QC, Canada, 2015; pp. 12–19. 5. Barrantes, A.; Ugalde, S. Usos y aportes de la madera en Costa Rica estadísticas 2019 precios 2020. Oficina Nacional Forestal ONF. (2019, août). 6. Bucur, V.; Eberhardsteiner, J.; Mang, H.; Waubke, H. Acoustics of Wood, 2nd ed.; Publisher: New York, NY, USA, 2006. (s.n.). 7. Carrasco, E.V.M.; Alves, R.C.; Smits, M.A.; Pizzol, V.D.; Oliveira, A.L.C.; Mantilla, J.N.R. Influence of the applied pressure of the transducer on the propagation speed of the ultrasonic wave in wood. Holzforschung 2021, 75, 1097–1103. https://doi.org/10.1515/HF-2020-0272. 8. Casado, M.; Acuña, L.; Basterra, L.A.; Ramón-Cueto, G.; Vecilla, D. Grading of structural timber of Populus × euramericana clone I-214. Holzforschung 2012, 66, 633-638. https://doi.org/10.1515/HF-2011-0153. 9. Castro, F.; Raigosa, J. Crecimiento y propiedades físico-mecánicas de la madea de teca (Tectona grandis) de 17 años de edad en San Joaquín, de Abangares, Costa Rica. Agron. Costarric. 2000, 24, 7–23. 10. Chauhan, S.; Sethy, A. Differences in dynamic modulus of elasticity determined by three vibration methods and their relationship with static modulus of elasticity. Maderas. Cienc. Y Tecnol. 2016, 18, 373–382. https://doi.org/10.4067/S0718- 221X2016005000034. Forests 2022, 13, x FOR PEER REVIEW 10 of 11 11. Ferreira, A.T.B. Avaliação da Estrutura Anatômica e da Densidade Aparente do Lenho e do Carvão de Árvores de Eucalyptus sp. e de Corymbia sp. Ph.D. Thesis, Universidade de São Paulo, São Paulo, Brazil, 2013. (Doctoral dissertation). 12. Gerhards, C.C. Effect of moisture content and temperature on the mechanical properties of wood: An analysis of immediate effects. Wood Fiber Sci. 1982, 14, 4–36. 13. Goncalves, R.; da Costa, O.A.L. Influence of moisture content on longitudinal, radial, and tangential ultrasonic velocity for two Brazilian wood species. Wood Fiber Sci. 2008, 40, 580–586. 14. Gonçalves, R.; Lorensani, R.; Pedroso, C.B. Moisture-related adjustment factor to obtain a reference ultrasonic velocity in structural lumber of plantation hardwood. Wood Mater. Sci. Eng. 2017, 13, 254–261. https://doi.org/10.1080/17480272.2017.1313312. 15. Ilic, J. Dynamic MOE of 55 species using small wood beams. Holz Als Roh-Und Werkst. 2003, 61, 167–172. https://doi.org/10.1007/s00107-003-0367-8. 16. Íñiguez-González, G.; Arriaga, F.; Esteban, M.; Llana, D.F. Reference conditions and modification factors for the standardization of nondestructive variables used in the evaluation of existing timber structures. Constr. Build. Mater. 2015, 101, 1166–1171. https://doi.org/10.1016/j.conbuildmat.2015.05.128. 17. Kang, H.; Booker, R.E. Variation of stress wave velocity with MC and temperature. Wood Sci. Technol. 2002, 36, 41–54. https://doi.org/10.1007/s00226-001-0129-x. 18. Karlinasari, L.; Azmi, M.I.; Priadi, T. The changes in color and dynamic modulus of elasticity of five important Indonesian tropical wood species after 10 months of outdoor exposure. J. Indian Acad. Wood Sci. 2018, 15, 149–157. https://doi.org/https://doi.org/10.1007/s13196-018-0220-7. 19. Llana, D.F.; Íñiguez-González, G.; Esteban, M.; Hermoso, E.; Arriaga, F. Timber moisture content adjustment factors for nondestructive testing (NDT): Acoustic, vibration and probing techniques. Holzforschung 2020, 74, 817–827. https://doi.org/10.1515/hf-2019-0187. 20. Marhenke, T.; Neuenschwander, J.; Furrer, R.; Twiefel, J.; Hasener, J.; Niemz, P.; Sanabria, S.J. Modeling of delamination detection utilizing air-coupled ultrasound in wood-based composites. NDT E Int. 2018, 99, 1–12. https://doi.org/10.1016/j.ndteint.2018.05.012. 21. Mascarenhas, A.R.P.; de Melo, R.R.; Pimenta, A.S.; Stangerlin, D.M.; de Oliveira Corrêa, F.L.; Sccoti, M.S.V.; de Oliveira Paula, E.A. Ultrasound to estimate the physical-mechanical properties of tropical wood species grown in an agroforestry system. Holzforschung 2020, 75, 879–891. 22. Mishiro, A. Ultrasonic Velocity and Moisture Content in Wood II. Ultrasonic velocity and average moisture content in wood during desorption (1) ; Moisture content below the fiber saturation point. Mokuzai Gakkaishi 1996, 42, 612–617. 23. Mishiro, A. Ultrasonic Velocity and Moisture Content in Wood III. Ultrasonic velocity and average moisture content in wood during desorption (2) : During desorption from a water-saturated condition. Mokuzai Gakkaishi 1996, 42, 930–936. 24. Moya-Roque, R.; Muñoz-Acosta, F.; Salas-Garita, C.; Berrocal-Jiménez, A.; Leandro-Zúñiga, L.; Esquivel-Segura, E. Tecnología de madera de plantaciones forestales: Fichas técnicas. Rev. For. Mesoam. Kurù 2010, 7, 18–19. 25. NBR 15521. NBR 15521 Non-Destructive Testing-Ultrasonic Testing—Mechanical Classification of Dicotyledonous Sawn Wood. Publisher: Sao Pablo, Brazil, 2007. 26. Neuenschwander, J.; Sanabria, S.J.; Schuetz, P.; Widmann, R.; Vogel, M. Delamination detection in a 90-year-old glulam block with scanning dry point-contact ultrasound. Holzforschung 2013, 67, 949–957. https://doi.org/10.1515/HF-2012-0202. 27. Nishimura, H.; Kamiya, A.; Nagata, T.; Katahira, M.; Watanabe, T. Direct evidence for α ether linkage between lignin and carbohydrates in wood cell walls. Sci. Rep. 2018, 8, 1–11. 28. Nocetti, M.; Brunetti, M.; Bacher, M. Effect of moisture content on the flexural properties and dynamic modulus of elasticity of dimension chestnut timber. Eur. J. Wood Wood Prod. 2015, 73, 51–60. https://doi.org/10.1007/s00107-014-0861-1. 29. Paniagua-Hernández, V.; Torres, C.; González-Rodrigo, B.; Sasa, Z. Methodology to Analyze the Degradation of Structural Timber in Acidic Atmospheres. In Proceedings of the Dans World Conference of Timber Engineeringerence on Timber Engineering, Santiago de Chile, Chile, 2021. 30. Peng, H.; Jiang, J.; Zhan, T.; Lü, J. (2016). Influence of density and moisture content on ultrasound velocities along the longitudinal direction in wood. Sci. Silvae Sin. 2016, 52, 117–124. 31. Saadat-Nia, M.A.; Brancheriau, L.; Gallet, P.; Enayati, A.A.; Pourtahmasi, K.; Honarvar, F. Ultrasonic wave parameter changes during propagation through poplar and spruce reaction wood. BioResources 2011, 6, 1172–1185. 32. Sanabria, S.J.; Furrer, R.; Neuenschwander, J.; Niemz, P.; Sennhauser, U. Air-coupled ultrasound inspection of glued laminated timber. Holzforschung 2011, 65, 377–387. https://doi.org/10.3929/ethz-b-000036154. 33. Sandoz, J.L. Grading of construction timber by ultrasound. Wood Sci. Technol. 1989, 23, 95–108. 34. Sandoz, J.L. Moisture content and temperature effect on ultrasound timber grading. Wood Sci. Technol. 1993, 27, 373–380. 35. Sobue, N. Simulation study on stress wave velocity in wood above fiber saturation point. J. Jpn. Wood Res. Soc. (Jpn.) 1993, 39, 271–276. 36. Solorzano, S.; Moya, R.; Murillo, O. Early prediction of basic density, shrinking, presence of growth stress, and dynamic elastic modulus based on the morphological tree parameters of Tectona grandis. J. Wood Sci. 2012, 58, 290-299. https://doi.org/10.1007/s10086-012-1261-y. 37. Steiger, R. Versuche an Fichten-Kanthölzern: Biegemoment-Normalkraft-Interaktion. Institut Für Baustatik und Konstruktion (IBK), Eidgenössische Technische Hochschule Zürich (ETH). ETH Zurich: Zïrich, Swiss, 1995. https://doi.org/10.3929/ethz-a-001473586. Forests 2022, 13, x FOR PEER REVIEW 11 of 11 38. Unterwieser, H.; Schickhofer, G. Influence of moisture content of wood on sound velocity and dynamic MOE of natural frequency- and ultrasonic runtime measurement. Eur. J. Wood Wood Prod. 2011, 69, 171-181. https://doi.org/10.1007/s00107-010- 0417-y. 39. Wang, S.Y.C. Experimental data correlation of the dynamic elastic moduli, velocity, and density of solid wood as a function of moisture content above the fiber saturation point. Holzforschung 2000, 54, 309–314. 40. Wang, S.Y.; Lin, C.J.; Chiu, C.M. The adjusted dynamic modulus of elasticity above the fiber saturation point in Taiwania plantation wood by ultrasonic-wave measurement. Holzforschung 2003, 57, 547–552. https://doi.org/10.1515/HF.2003.081. 41. Wang, X. Effects of Size and Moisture on Stress Wave E-rating of Structural Lumber. In Proceedings of the Dans 10th World Conference on Timber Engineering, Miyazaki, Japan, 2008; pp. 1–8. 42. Yamasaki, M.; Tsuzuki, C.; Sasaki, Y.; Onishi, Y. Influence of moisture content on estimating Young’s modulus of full-scale timber using stress wave velocity. J. Wood Sci. 2017, 63, 225–235. https://doi.org/10.1007/S10086-017-1624-5.