Abstract: This article investigates the relationships between soil characteristics (physical and chemical) and wood properties of Persian oak in three different elevation sites. For this purpose, 27 trees were randomly chosen and cut in Zagros forests in western Iran. The test samples were prepared at the stem (breast height) to examine physical and biometrical properties. For each elevation site, four soil samples were obtained at a 0-20 cm soil depth under the canopy of each tree to measure soil properties, including clay, silt soil, sand soil, electrical conductivity, pH, nitrogen, phosphorus, potassium, and organic matter content. Then, the relationship of soil and wood properties was determined by principal component analysis. Results specified that there are a positive correlation between wood density and volumetric swelling with clay and available potassium. Moreover, the results revealed a positive correlation between fiber length, cell wall thickness, and fiber diameter with electrical conductivity, sand percentage, and total nitrogen content, respectively.
Keywords: Quercus brantii, soil properties, wood density, wood fiber, wood quality.
ARTÍCULO
Influence of soil physicochemical properties on biometrical and physical features of persian oak wood
Huijun Dong
Mohammad Dahmardeh Ghalehno
Mohsen Bahmani mohsen.bahmani@sku.ac.ir
Elham Ghehsareh Ardestani
Leila Fathi
Received: 22 September 2021
Accepted: 12 September 2022
Zagros forest in western Iran covers more than five million ha (Nazari et al. 2021). This forests consists sub-Mediterranean, semiarid temperate climate, mainly consisting of deciduous, broad-leaved trees (Olfat and Pourtahmasi 2010). Crategus spp. and Pyrus spp., Quercus spp. (oaks), and Pistacia mutica (wild pistachio) (Ghazanfari et al. 2004) are the three main wood tree species in this district. The ring-porous tree species Quercus brantii (covering more than 50 % of the Zagros forest area) has dominated tree woody species in the Zagros forests growing at 450 m to 2700 m altitudes (Sagheb-Talebi 2005). Native to western Asia and with more than 50 % coverage of the Zagros forests in Iran, it is one of the ecologically and economically most important tree species (Safari et al. 2022). They play an important role in water protection, water regulation, and soil protection, as well as perform sanitary and recreational functions. Persian oak is a big tree with height of 20 m and a big spherical crown and it has generally ovoid leaves with serrated margins. The fruit of Q. brantii is elongated and oval in velvety white bowl and conical (Alikhani et al. 2014). It has been reported that climate and environmental factors affect the quality of wood (Zobel and van Buijtenen 1989, Nazari et al. 2020). Climatic factors (temperature and rainfall), geographic features (latitude, elevation, and slope), biotic factors (animals, plants, humans, fungi, bacteria), and physical and chemical properties of soil can impact wood properties (Topaloglu et al. 2016, Nazari et al. 2021). Although numerous studies investigated the relationships between ecological indexes and wood properties (Rossi 2015, Topaloglu et al. 2016, Sousa et al. 2018), our knowledge is still limited regarding the relation between soil properties and variations in wood quality.
Research on the physical and chemical properties of the soil that determine the properties of wood is complicated because it is affected by many other factors, such as the age of the tree, heredity, etc., and the relationship between parameters. Although the physical and chemical properties of the soil have a key role in the growth of trees and wood quality, their effects on wood properties are still rarely investigated. (Marini et al. 2021). Rigatto et al. (2004) indicated that there is a significant correlation between soil physicochemical properties and wood quality. Cutter et al. (2004) found that sites with favorable soil properties significantly influenced the wood density. Tufekcioglu et al. (2005) studied the influence of soil characteristics on hybrid poplar growth in Turkey and concluded that the content of clay and Mg+2 was negatively correlated with the pH of the soil, while the range of phosphorus and sand of the soil was positively correlated with hybrid poplar growth. Bektas et al. (2003) investigated the site effect on heartwood and sapwood portions of Turkish Calabrian pine. They stated that the differences in heartwood and sapwood might be related to ecological indicators such as elevation, lime content, soil organic matter, and soil type. Yáñez-Espinosa et al. (2001) indicate that soil texture and water salinity are closely associated with the anatomical characteristics of Mexican mangrove populations. Kiaei (2014) pointed out that there is a significant correlation between soil clay and sand with wood properties of Eldar pine (Pinus eldarica), in contrast a negative relationship between wood properties and soil chemical properties was found. In Iran, there is no study related to the influence of soil physical and chemical properties on the wood properties of Persian oak.
Principal Component Analysis (PCA) is an algorithm in biometrics. It is a statistical technique that makes use of orthogonal transformation to transform a set of data of probably linked variables into a set of values of unrelated linear variables. It is also a tool that reduces multidimensional data to lower dimensions while retaining most of the information. (Karamizadeh et al. 2013). PCA is a sorting method that develops theoretical variables and reduces the total residual sum of squares after fitting a straight line to the species data. The best wood properties values for the sites are selected by PCA (Kooch et al. 2008).
Considering the importance of the chemical and physical properties of the soil to the wood quality and the value of Persian oak to Iranian forests, this research aimed was to investigate the influence of soil properties on physical properties (dry density and volumetric swelling) and wood fiber properties (fiber length, cell wall thickness, and fiber diameter) in three different elevations and to inspect the relationship between wood properties with soil. The relationship of mentioned variables was determined by PCA.
This study was carried out in trees of Persian oak (Quercus brantii Lindl) collected from tree sites (three forest sample stands along an altitudinal gradient, below 1900, 1900-2000 and above 2200 m a.s.l.) in the oak forests of Bazoft. Due to their remote location and the absence of evidence of human impact, it is assumed that all stands have been developed under the influence of natural impacts and disturbances. In each site, nine plots were implemented, accordingly the 27 plots were selected in study area. In each plot, all live trees of at least 7,5 cm diameter at breast height (DBH 1,3 m above the root collar) were identified, and their diameter at breast height, height, crown length and perpendicular diameters were recorded within 0,1-hectare plot. The caliper, Vertex and diameter tape were used to measurements of trees diameters, height and crown diameter, respectively. Then stand-level variables such as stand basal area (square meter per hectare), stand density (the number of trees per hectare), stand quadratic mean diameter) QMD) (cm) and tree-level variables such as stem basal area (m2), stem volume (m3), crown diameter (m) and crown basal area (m2) in each site were calculated. Then in each plot, those dominant trees that healthy with the largest diameters at breast height (DBH) without any defects and reaction wood were sampled and one disc was taken from the tree trunk at DBH for the determination of wood properties. Then, the tree age at breast height (ABH) was obtained by counting the annual rings of the sampled disks (Figure 1). Finally, the mean annual diameter increment (MADI) was obtained by dividing DBH by the number of annual rings. In total, 27 healthy trees (9 trees× 3 sites) were selected. The main dimensions of the tree and site properties are presented in Table 1.
Physical properties, namely oven-dry density and volumetric swelling were measured from 5 cm thick disks cut down from each tree. The disks were cut from the stem diameter at breast height (DBH = 1,30 m aboveground). Sample tests were prepared following ISO 13061-14 (2016). The wood sample size was 3×2×2 cm3. Two hundred seventy samples were prepared from various portions of the disks. The size of samples in both green and oven-dry situations was measured with a slide caliper. An electronic balance was used to measure the oven-dry density. Dimensional changes from the green to dry conditions were considered to calculate the volumetric swelling. The mentioned physical properties were determined in accordance with the formula (Equation 1 and Equation 2 ):
Where: D0 - oven-dry density (g.cm-3), M0 - oven-dry mass (g), V0 - oven-dry volume (cm3), αv - volumetric swelling (%), Vs. - volume in state of saturate (cm3).
Franklin (1945) was applied to separate the individual wood fiber. In detail, Saturation of wood samples (15×10×2 mm3) in a mixture (1:1) of oxygenized water and acetic acid were performed, and then the samples were maintained inside an oven for 48 hours at 65 ºC ± 5 °C. After maceration, the samples were washed in distilled water. Finally, they were submerged to distilled water, shacked, and fiber parameters such as fiber length, cell wall thickness, and fiber diameter were evaluated under a light microscope. From each slice, at least 50 fibres were used for the measurements.
To calculate soil properties, four soil samples were obtained at a 0-20 cm soil depth under the canopy of each tree and mixed. Physicochemical properties of soil were measured, including the percentage of clay, percentage of silt, percentage of sand, electrical conductivity (Ec), soil reaction (pH), total nitrogen (N), available phosphorus (P), available potassium (K), and percentage of organic matter.
Soil treatments were air-dried, then passed through a 2 mm sieve. Analysis of soil samples was carried out at the laboratory of Agricultural and Natural Resources Research Center of Isfahan Province, Isfahan, Iran. Percentage of clay, percentage of silt, percentage of sand were calculated by the hydrometer method (Bouyoucos 1962). The pH and EC were determined using a pH/EC meter. The total nitrogen was calculated using the Kjeldahl method (Zarinkafsh 1993). Phosphorus content was calculated by the Olsen method (Nelson and Sommers 1996). Available potassium was measured by Flame photometry (Zarinkafsh 1993). Organic matter content was estimated as in Walkley and Black (1934).
In this study, the influence of different sites was studied on the wood properties. Analysis of variance (ANOVA) and Duncanʼs test were applied to differentiate the significance of means and their grouping, respectively. PCA was also considered for determining the relationships of the site conditions and soil properties with the wood properties.
Table 1 gives descriptive statistics analysis of Q. brantii population in three studied sites.
According to the ANOVA results, there are no significant differences in dry wood density between different sites. However, the sites affected the volumetric swelling (VS) (P < 0,01) (Table 3), with the highest and lowest values measured at sites 2 and 1, respectively. Nevertheless, no significant differences were found between sites 2 and 3. The average values of VS were 17,34 % ± 4,13 % (Table 4).
The analysis of variance shows that, in addition to the regular changes of fiber length at different sites, the influence of sites on fiber length was highly significant (p < 0,01 Table 2). In addition, the value of fiber length declined with raising altitudes (Table 3). The average values of fiber length were 0,87 mm ± 0,08 mm. ANOVA results indicated significant differences in fiber diameter among the sites (P < 0,01; Table 2), with the maximum and minimum values which are measured at sites 1 and 2, respectively. It was found that there are no significant differences between sites 2 and 3. The mean fiber diameter values were 20,50 μm ± 1,25 μm (Table 3). Analysis showed no statistically significant differences between the sites for cell wall thickness.
PCA discoveries hypothetical components that account for most of the variance in multivariate data (Davis 1986, Harper 1999). For plotting goals, PCA can reduce the data set to only two components (the first two variables). It can also be assumed that the most critical components are related to other main variables. In other words, it is necessary to discover the best axis between each axis before PCA can obtain a real axis. PCA scree plot of the studied variables is illustrated in Figure 2. The 'Scree plot' can also indicate the number of essential components. The plot showed a line plot of principal components in multivariate statistics. It can indicate the number of variables to retain in main components to keep in PCA.
In Table 4, the eigenvalues and proportion of variance showed using the component axes. The first six principal components (PCs) have eigenvalues > 1 and contribute ~89,09 % of the total variation in the set of data. Therefore, these PCs provide a decent approximation of the variation in the 21 axes (Figure 2). Jolliffe and Cadima (2016) showed that the threshold of cumulative ~70 % variation is appropriate to keep the PCs for PCA. Although the first four PCs have eigenvalues > 1 and contribute ~75,42 %, these PCs will be complicated to visualize at once, and these PCs need to implement pairwise visualization (Table 4).
In Table 4, the first and second PCs are the major to describe variance across the factors (soil and wood factors). The percentage of eigenvalue for the two PCs is 34,31 % and 18,95 % respectively. The place of soil variables, wood features, and sites in different parts of Figure 2 and Figure 3 is based on the correlation coefficient between factors.
As shown in Table 5, eigenvectors include scores that show the weight of each factor (soil and wood factors) on the two PCs (-1< eigenvectors<+1); Kent and Coker (1992) demonstrated that if the score of the eigenvector is close to ±1 for each factor, the factor is more important to weigh on the PCs. The eigenvectors for fiber length, volumetric swelling, percentage of clay, percentage of organic matter, and Ec for the first PC, the value of DBH, height, stem basal area, and stem volume for the second PC is more significant relative to other eigenvectors factors. These factors perform the main role in explaining the first and second PCs of alteration.
As shown in Figure 3, the first PC describes 34,31 % of the total alteration. This PC is most greatly correlated with the soil properties available potassium, pH, percentage of sand, percentage of clay, percentage of organic matter, Ec; and wood properties crown basal area, crown diameter, fiber length, and volumetric swelling. The second PC explains 18,95 % of the alteration highly influenced by the soil properties total nitrogen, available phosphorus, percentage of silt, and wood properties DBH, height, stem basal area, stem volume, annual diameter increment, wood dry density, and fiber diameter.
In Figure 4, we see that sites were placed separately from each other. The PCA analysis of the soil physical and chemical factors and wood physical and biometric parameters drive an independent distribution of the samples that seem grouped when they relate to the same site (Figure 4). Therefore, PCA displays apparent relationships of wood features fiber length and fiber diameter and soil properties Ec, pH, and percentage of sand in site 1 (bottom right). The PCA data also shows that site 2 (upper center) has the highest relationship with the crown basal area, crown diameter, cell wall thickness and percentage of silt, total nitrogen, and available phosphorus. Site 3 (left) has the best relationship with volumetric swelling, stem basal area, DBH, and stem volume, and percentage of organic matter, percentage of clay, and available potassium (Figure 4).
The width of the annual rings will be closely related to the wood density, which is important in ring-porous species. The density of wood affect the main features of wood quality. Therefore, wood quality can change with tree genetic, site conditions e.g. the soil components, light and water, and other factors (Vintoniv et al. 2007, Grześkiewicz 2007, Hemery et al. 2008, Sopushynskyy and Teischinger 2013). The density of wood is generally lower on more fertile soils, which are conductive to faster tree growth. This lead to an increase in the percentage of wood vessels with thinner fiber walls on fertile soils (Lima et al. 2010, Mevanarivo et al. 2020).
According to PCA results, wood density has a negatively correlate with silt and phosphorus and a weak positive correlation with clay, sand, and potassium content. Similar results were also observed by Kiaei (2014). Parsa (1993) stated that sandy soils reduce the density of the wood. Maharani and Fernandes (2015) demonstrated that mineral elements (N, P, K) had a major correlation with fiber length and wood density on Shorea leprosula, and these elements had a notable correlation with wood density on S. parvifolia. In the meantime, these elements had a non- considerable correlation with fiber length.
A number of studies have reported how components in the soil, especially the nitrogen content, affect wood properties. Variation of wood density with components in the soil might vary according to the amount of the nutrients present or applied, and on the species (Downes et al. 2014, Eufrade-Junior et al. 2017, Mevanarivo et al. 2020).
Volumetric swelling is positively correlated with available potassium, organic matter content, clay, and total nitrogen. Meanwhile, it is negatively correlated with Ec. Moya and Perez (2008) reported that normal radial shrinkage and normal tangential shrinkage were the best-correlated factors with soil properties, even though that the less correlated factors were normal volumetric shrinkage and specific gravity.
In the current study, fiber length positively correlated with the Ec and negatively correlated with the available potassium, while fiber diameter positively correlated with soil sand percentage. It was found that the fiber diameter had no significant correlation with elements of N, P, and K. Cell wall thickness is positively correlated with the total nitrogen and available potassium. In contrast, available phosphorus and percentage of organic matter presented a weak relationship with cell wall thickness. Larson (1994) demonstrated that the mineral elements have direct effects on vascular cambium. For the particular case of phosphorous, this element increases cell division in the vascular cambium, and this element allows improved performance in plants. The increase in cambial activity is followed by changes in the structure, mainly because of pores of larger size, fibers with thinner cell walls, and more significant, presence of parenchyma cells. Moya and Perez (2008) found such anatomical structures make a diminution in the normal tangential shrinkage.
In this study, the relationship between soil physicochemical properties and some wood properties of Persian oak was investigated. Soil characteristics (physical and chemical) had a significant influence on oak wood properties. It can be stated that soil physical properties correlated with fiber diameter, while soil chemical properties associated with cell wall thickness. Meanwhile, soil physical and chemical properties correlated with wood dry density, volumetric swelling, and fiber length. Zagros forests are particularly important for protecting soil and water supply and they are known as a second renewable source of cellulose in Iran that are threatened by environmental factors such as climate change, disease and pests, and human activities. To forests rehabilitate with oak species, soil with a high content of potassium, nitrogen, and organic matter is suggested because these factors play a crucial role in the wood properties and wood quality.
We are grateful for financial support from the University of Zabol (Grant No. UOZ-GR-9618-111) and Shahrekord Universty.
♠Corresponding author: mohsen.bahmani@sku.ac.ir
