Next Article in Journal
Short-Term Effects of Bamboo Biochar and Oyster Shell Powder on Soil Organic Carbon Fraction, Microbial Respiration, and Enzymatic Stoichiometry in a Lei Bamboo Plantation
Next Article in Special Issue
Forest Gaps Slow the Humification Process of Fir (Abies faxoniana Rehder & E.H.Wilson) Twig Litter during Eight Years of Decomposition in an Alpine Forest
Previous Article in Journal
Identifying Risk Factors and Evaluating Occupational Safety in South Korean Forestry Sector
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Classification of Log Decay Classes and an Analysis of Their Physical and Chemical Characteristics Based on Artificial Neural Networks and K-Means Clustering

1
College of Forestry, Northwest A&F University, Yangling 712100, China
2
College of Mechanical and Electronic Engineering, Northwest A&F University, Yangling 712100, China
3
Shaanxi Forestry Survey and Planning Institute, Xi’an 710082, China
4
Qinling National Forest Ecosystem Research Station, Ningshan 711600, China
*
Author to whom correspondence should be addressed.
Forests 2023, 14(4), 852; https://doi.org/10.3390/f14040852
Submission received: 25 February 2023 / Revised: 14 April 2023 / Accepted: 19 April 2023 / Published: 21 April 2023
(This article belongs to the Special Issue Forest Succession and Leaf Litter Decomposition)

Abstract

:
Most existing methods for determining log decay levels normally use variations in log surface characteristics, and the results are subject to human subjectivity, which is uncertain and inaccurate. In order to investigate a novel method for the quantitative determination of log decay levels, we randomly selected log samples from four species (Pinus tabulaeformis, Larix principis-ruprechtii, Betula albosinensis and Quercus aliena var. acuteserrata) with different levels of decay and determined their basic physicochemical characteristics in the laboratory. An artificial neural network (ANN) model was used to predict the hardness values of the log samples with different levels of decay at different moisture contents. The hardness was then used as a clustering factor to quantify the decay levels of the log via K-means clustering analysis. The variations in and correlations between the basic physicochemical factors of the log specimens were investigated between the different decay classes and between the different tree species, and then ANOVA and correlation analysis were used to verify the reliability of the clustering results. The results showed that the prediction of the hardness of the decayed log by the ANN was very effective and that the highly significant variability in the dry matter content, basic density and some basic chemical element contents between the log samples that were classified into different decay grades confirmed the reliability of the clustering results. This study explores an innovative method for the quantitative determination of log decay classes.

1. Introduction

Logs are important components of forest ecosystems and have received increasing amounts of attention because of their importance in ground maintenance, biodiversity conservation and carbon and nutrient cycling [1,2]. A log is defined as a log fragment with a diameter of 10 cm or more at the thickest end, a length of more than 1 m, and an inclination of 45° or more from the vertical [3]. Logs are coarse woody residue, as are standing dead trees, root stumps and dead branches. These forms of woody residue can be created by competitive exclusion during tree growth, the natural mortality of old trees, other natural factors (such as rain, snow and fire) and human disturbance (such as logging) [1]. Typically, the decomposition of a downed log occurs through a combination of respiration, leaching and natural shattering, which is a very complex process that does not always occur continuously [1,4,5]. The decay level of the log is generally classified from I to V. At different decay levels, the physicochemical properties of the log (i.e., density, hardness, moisture content, cellulose content, nitrogen (N) and phosphorus (P) content, etc.) change accordingly.
Log decay includes the process of carbon that is fixed in the log and released into nature [6], which is an important link within the material cycle of forest ecosystems that can maintain the stability and biodiversity of forest ecosystems, provide necessary living environments, and nutrients for decomposers (such as insects and invertebrates) and promote the rapid renewal of forests [7]. At the same time, different kinds of logs have different decay dynamics. There are many factors that affect the decomposition of logs, and the interactions between them are complex. Additionally, because log decomposition takes a long time (i.e., decades or even centuries), the factors of log decomposition have not been studied enough and need to be further explored in depth [8]. The density of the log is the main physical indicator of the log decay process; meanwhile, the hardness of the log is another indicator that is used to evaluate decay. The determination and analysis of the density and hardness of logs with different decay processes can help judge the decay class of logs more accurately.
Commonly used log decay class systems mainly rely on morphological characteristics that can be easily observed and recognized in the field [9], such as the presence of bark and branches, log form and color, the shape and integrity of the main trunk, etc. These characteristics have been used as a basis for constructing grading index systems, and other characteristics of the log have been combined to create the final systems. The earliest decay class system was the five-level classification system proposed by Logel in 1972. Then, thanks to continuous improvement from many other scholars, the decay class system for logs that is currently in use was formed [10,11,12]. This classification system is the most widely used and is the preferred system for many scholars for the classification of log decay classes.
This traditional grading system can quickly and roughly determine the decay class of a log, but it also has certain limitations. In general, the determination of the decay class of a log is usually based on the human judgment of experienced workers, meaning some other subjective factors can be introduced that may affect the results of the rest of the experiments based on the determined decay class [13]. In addition, there is a lack of parameterized indicators to describe the decomposition dynamics of the log [14].
The relative density of the log is species-specific and is one of the most important indicators for judging the level of decay [15,16]. The density of the log can very closely be related to species, the growing environment, microbial colonization, and the time of death of the log, with a less dense log being more susceptible to decomposition and exhibiting a higher respiration rate [17]. There is a strong relationship between log density and decay class [18], with log density appearing the greatest at the time of death and decreasing throughout the respiration [1,19], leaching and natural comminution process until all nutrients have entered the soil [20].
As the decay process proceeds, the hardness of the log changes, and in the early stages of decay, this hardness may even decrease faster than the density [14]. Additionally, some logs may become hard and non-decayable on their surface due to environmental factors [1], thus affecting the determination of the decay class. In addition, the hardness of the log gradually decreases with the loss of mass during decay, and the rate of decline is related to the species [21], meaning that the change in hardness is not uniform and the decaying log tends to show local variations in hardness under the action of decomposers [22].
Logs are mainly composed of lignin, cellulose and hemicellulose [23], the content of which also affects the decay process of the log to some extent [24,25]. All three have differences in degradability and structure: cellulose and hemicellulose are single chains of glucose molecules that have simple structures [26], and so can be more easily degraded by various organisms in the early and middle stages of decomposition and are usually considered as the main sources of C in forests [27]; lignin is an aromatic compound with a more complex structure that can only be degraded by specific soil microorganisms and is more resistant to decay than cellulose and hemicellulose [14,28,29,30,31]. Lignin is commonly found in higher plants, and studies have shown that lignin content has a greater effect on the decomposition of logs than other chemical contents [32], meaning that lignin may decompose the fastest during early log decomposition [33]. Cellulose and hemicellulose are the main components of plant cell walls and are important energy sources for microorganisms that are involved in the decomposition of logs [34], so the degree of decomposition for cellulose and hemicellulose can be correlated with the decay of the log and there is a relationship between the decay grade of the log and cellulose and hemicellulose contents.
There is no systematic and quantitative method to determine the decay class of logs, so the judgment of log decay classes in the field mainly relies on some visible characteristics (i.e., the structural integrity of leaves, the decay of heartwood and sapwood, the decay of the crown and branches, etc.) and human qualitative judgment, which can be influenced by human subjective factors and the environment. Hardness is the easiest physical factor to measure, and it can quantitatively reflect changes in log decay in the field, but there is no standard for determining the decay class according to the hardness value. Therefore, it is important to develop a set of simple, efficient, scientifically sound and applicable criteria for determining log decay classes in the field.
In this study, log samples from four typical tree species in the Huoditang forest area of the Qinling Mountains were used as research objects for the following objectives: (1) to determine the hardness, density, elemental contents and environmental factors of decayed log and systematically analyze the relationships between each index; (2) to conduct quantitative grading based on the results of cluster analysis and construct an artificial neural network model to determine the decay class of the log from the measured hardness values. In this study, we innovatively propose a quantitative determination method for the decay grade of logs, which could provide a scientific basis for the determination of the decay class of logs.

2. Materials and Methods

2.1. Site

The present study was carried out in the 33°25′–33°29′ N and 108°25′–108°30′ E area of the southern slope of the middle part of the Qinling Mountains. The altitude of the study area was 800–2500 m above sea level, with a mean annual precipitation of 900–1200 mm and evaporation of 800–950 mm. The average annual temperature was 8–10 °C, with 1100–1300 h of sunshine. The growing season last for 6 months. The terrain in the area is diverse, with broken slopes and steep mountains. The vegetation in the study area was rich and diverse. The main tree species were P. tabulaeformis, B. albosinensis, and Pinus. armandii and Q. aliena var. acuteserrata.

2.2. Materials

In this study, 46 log samples from four major forest trees (P. tabulaeformis, L. principis rupprechtii, B. albosinensis, and Q. aliena var. acuteserrata) with different levels of decay from the Huoditang area of the Qinling Mountains were collected using the quadrat method.

2.3. Experimental Procedures

2.3.1. Determination of Basic Factors

The samples were collected from the field, and for larger, mildly decayed logs, we utilized a hacksaw to cut 5 cm thick discs. For heavily decayed logs, we extracted partial samples using a knife and packed them into an aluminum box of known volume. Then, these were placed in an oven at 75–85 °C and dried to a constant weight. The dry weight m1 of each sample was recorded. The dried samples were sprayed with water for 2 h, 4 h, 6 h, 8 h and 10 h. The wet weights m2 of each sample was then also measured. Following treatment, spraying was continued until the specimens reached water saturation, at which point the water saturation weight m3 of each specimen was recorded. The volume V0 of each specimen under water saturation was determined using the drainage method, and the basic density and dry matter content (DMC) of each specimen were also calculated. The National Renewable Energy Laboratory (NREL) method [35] was used to determine the acid-soluble lignin, acid-insoluble lignin, cellulose, hemicellulose, glucose and xylose contents of the samples.
B a s i c   d e n s i t y = m 1 V 0  
D r y   m a t t e r   c o n t e n t       D M C = m 1 m 3
m o i s t u r e   c o n t e n t = m 2 m 1 m 2 × 100 %
where m1 is the dry weight of a specimen, m2 is the wet weight of a specimen, V0 is the wet volume of the sample under saturated water absorption, and m3 is the weight of the specimen under saturated water absorption.

2.3.2. Determination of Hardness

To determine the hardness of the specimens, their moisture contents were calculated, and the hardness of the cross-sections was determined in accordance with the established method (i.e., the “Method for Determination of Hardness of Log”) using a Universal Mechanical Experiment Instruments.

2.3.3. Data Analysis

  • Artificial Neural Network
Since each indicator had different dimensions and dimensional units, this situation affected the results of data analysis, and in order to eliminate the dimensional impact between indicators, data normalization was required to solve the comparability between data indicators. At the same time, in order to speed up the convergence of the model, we chose to normalize the data first [36]. The normalization of the indicator data could be expressed as follows:
T i = i = 1 N i
i ^ = i T i
where i represents the measured value of the experimental sample for each indicator, Ti represents the measured total value of all experimental samples for each indicator, and i ^   is the normalized data for each experimental sample measurement for each indicator.
We developed a three-hidden layer artificial neural network using a Python open-source library (PyTorch), in which the first and second layers used a rule as the activation function and the third layer used sigmoid as the activation function (Figure 1). Considering that tree species are categorical variables and cannot directly participate in ANN calculations, they need to be encoded. We used One-hot Encoding to convert all the categories of each multiple-category feature, which had m categories, into m binary features. To be specific, we encoded P. tabulaeformis, L. principis rupprechtii, B. albosinensis and Q. aliena var. acuteserrata as 1000, 0100, 0010, and 0001, respectively. This approach ensured that each feature corresponded to a specific category and effectively captured their characteristics. Our neural network structure consisted of one input layer, three hidden layers, and one output layer. The input layer comprised 13 neurons corresponding to 13 feature inputs, including varieties of trees 1, varieties of trees 2, varieties of trees 3, varieties of trees 4, moisture content, DMC, basic density, acid-insoluble lignin content, acid-soluble lignin content, cellulose content, hemicellulose content, glucose content, and xylose content, with 184 sets of data (175 sets as training data (95%) and 9 sets as the test data (5%)). We set the learning rate to 0.0001, the number of epochs to 1000, and the loss function to MSE loss.
The mean absolute percentage error (MAPE), root mean square error (RMSE), and determination coefficient (R2) was utilized to evaluate the performance of the ANN. The MAPE and RMSE values were considered the most important performance criteria, while R2 was used to measure how much of the original data could be explained by the established model. When the RMSE and MAPE values approached 0 and R2 approached 1, ANN predictions were optimum [37].
M A P E = 1 N { i = 1 N [ | x i y x i | ] } × 100
R M S E = i = 1 N ( y i x i ) 2 N
x ¯ = 1 N i = 1 N x i
R 2 = 1 i = 1 N ( x i y ) 2 i = 1 N ( x i x ¯ ) 2
where yi represents the experimental output, xi represents the predicted output, N represents the total number of samples and x ¯   represents the mean of the predicted outputs.
  • Hardness clustering
In this study, the K-means clustering method was used to cluster the hardness values for different moisture contents. Due to the limited indicators considered in the clustering process, we did not take the dimensional effect on experimental results into account and used raw data to visually present the clustering results. We also analyzed the differences in the hardness values among different tree species and found no significant variations among the selected tree species. Thus, the tree species factor was not considered in this study.
  • Normal distribution test and ANOVA
In different clustering results, there were significant differences in the number of different tree species among different decay classes. To ensure consistency in the number of samples in each set of experiments, we randomly selected 60 sets of data for normal analysis testing and ANOVA. The analysis of variance (ANOVA) was performed by classifying the differences between the categories into random differences, between-category differences and between-case differences. Only between-class differences and between-case differences needed to be considered and the homogeneity of variance ensured that the between-case differences between the categories were consistent. Additionally, between-category differences were expressed as the mean of each category. Firstly, a normal distribution test was performed on the other determination factors based on the results of the cluster analysis. A kurtosis absolute value of less than 10 and a skewness absolute value of less than 3 were accepted as the normal distribution. ANOVA was performed on the other factors based on the clustering results, and the results of the normal distribution test. A homogeneity of the variance test was performed using Levene’s test, and a significance level of greater than 0.05 was considered to meet the condition of the homogeneity of variance. ANOVA and Pearson correlation analysis was used to compare the differences between the tree species, decay classes, and the content of each physicochemical characteristic to test whether the differences were significant and whether there were correlations between the factors. The graphics were drawn using Origin 2022 drawing software.
Based on the clustering results, we attempted to judge the decay classes of the decayed log samples by calculating the Euclidean metric of the hardness of the decayed logs from the center of each cluster.
We calculated the Euclidean metric from the row data to the cluster centers as follows:
E u c l i d e a n   m e t r i c = ( h ^ h i ˜ ) 2
where h ^   is the hardness data and h i ˜   is the hardness data of the clustering centers, respectively.
This specific method is shown in the following flowchart(Figure 2):

3. Results

3.1. Prediction of the Hardness Value of Log under Different Levels of Decay and Different Moisture Content Conditions Using an ANN

The MAPE was determined to be 0.02069, 0.02134 and 0.01716 for the prediction of the hardness values for the training, test and overall datasets (Table 1). The training set RMSE was 1.0264, the test set RMSE was 0.6187, and the overall dataset RMSE was 0.9756, indicating that the ANN could predict changes in the hardness of the log under different levels of decay and moisture content conditions using the MAPE and RMSE criteria.
The R2 values for the training, test and overall datasets were 0.9925, 0.9964 and 0.9931, respectively. They were all above 0.99, indicating that the model fitted the regression well. The testing set R2 values was higher than the training set R2 values, indicating that the model was not overfitted. In both the training and test sets (Figure 3), the predicted hardness values converged into a straight line, indicating that they were close to the actual distribution of the hardness values and that the artificial neural network model was very accurate in predicting the hardness values of the logs under different levels of decay and different moisture content conditions.

3.2. K-Means Clustering Analysis and the Classification of Log Decay Classes Based on the Hardness of Log

With reference to the hardness values under different moisture content conditions, a total of five clustering centers were classified (Figure 4). The hardness value of the level 1 cluster center was 41.57 N/mm2, the hardness value of the level 2 cluster center was 32.17 N/mm2, the hardness value of the level 3 cluster center was 22.50 N/mm2, the hardness value of the level 4 cluster center was 14.36 N/mm2, and the hardness value of the level 5 cluster center was 7.64 N/mm2. Most specimens had hardness values concentrated in the level 2 cluster (54 samples), while the smallest number of specimens was concentrated in the level 1 cluster (21 samples).
The hardness values of the wood gradually decreased from level 1 to level 5 cluster centers. The range of variation in the hardness values was insignificant between the level 1 and level 3 clusters, all of which were above 31 N/mm2. From the level 3 to level 4 clusters, a significant change in hardness was observed, but between the level 4 and level 5 cluster centers, the change in hardness was insignificant. The results of the ANOVA analysis show that there were extremely significant differences in DMC, BD, Ce, Glu, and H among different decay classes, which proved the robustness of the clustering results (Table 2).
Therefore, we hypothesized that the basic physicochemical characteristics of wood could be used as inputs to predict the hardness values of logs using our ANN. We input the predicted hardness values and calculated the Euclidean distance between the hardness values and each hardness cluster center. Then, the cluster center corresponding to the minimum value of the distance was found to be the decay grade of the wood.

3.3. Physicochemical Characteristics of Logs with Different Decay Grades

DMC, basic density and hardness showed a highly significant (p < 0.01) negative correlation among the different decay classes of the same species, i.e., the DMC, basic density and hardness values of each species decreased as the decay level increased (Figure 5). The cellulose content of P. tabulaeformis, L. principis rupprechtii and B. albosinensis all showed a highly significant (p < 0.01) negative correlation with the decay classes of the log. L. principis rupprechtii showed a significantly (p < 0.05) negative correlation with the acid-insoluble lignin content and highly significant (p < 0.01) positive and negative correlations with acid-soluble lignin and cellulose contents, respectively, between the different levels of decay. B. albosinensis showed a highly significant (p < 0.01) positive correlation with acid-insoluble lignin content and highly significant (p < 0.01) negative correlations with acid-soluble lignin and hemicellulose content. Q. aliena var. acuteserrata showed a significant (p < 0.05) positive correlation with the acid-insoluble lignin and acid-soluble lignin contents and a significant (p < 0.05) negative correlation with the hemicellulose content. The rest of the factors were not significantly correlated. The variations in glucose and xylose contents among the different decay classes of the same tree species were consistent with the variations in the cellulose and hemicellulose contents, respectively.
There were highly significant (p < 0.01) differences in DMC, basic density, acid-insoluble lignin content and hardness among the decay classes of the same tree species (Figure 6). The acid-insoluble lignin contents of B. albosinensis and Q. aliena var. acuteserrata, the cellulose contents of B. albosinensis, L. principis rupprechtii and Q. aliena var. acuteserrata, and the hemicellulose contents of L. principis rupprechtii and Q. aliena var. acuteserrata all showed highly significant (p < 0.01) differences between the decay classes. All the factors for all the species showed decreasing trends with increasing decay levels, except for acid-soluble lignin and acid-insoluble lignin contents, which showed increasing trends.
There were highly significant (p < 0.01) differences in acid-insoluble lignin, acid-soluble lignin and hemicellulose contents between different species with the same decay class and highly significant (p < 0.01) differences in the basic density and cellulose content between species with a decaying class of 1, 2 or 5. However, there were no significant differences in the hardness values between the species with each decay class (p > 0.05).

4. Discussion

In this paper, we present a novel method to quantitively determine log decay grades according to changes in log hardness. In the process of collecting the samples, we judged and recorded the qualitative decay classes of the collected samples with reference to the qualitative classification method proposed by Yan et al. [9]. The existing qualitative decay classification methods are susceptible to human subjectivity, and the results can be inaccurate [13]. At the same time, the number of indicators used in existing qualitative decay classification methods is large; the processes are complex, the types of species investigated, and the criteria used can vary, thereby preventing comparisons. By contrast, quantitative judgments can determine log decay grades more easily, more accurately, and faster, especially in the field.
The physicochemical factors of log decay that were used in this study were found to be closely related to other factors of the log, such as density, moisture content, cellulose content and lignin content, in many existing studies [38,39]. Wang et al. employed the lignin and cellulose content as indicators to assess the decay level of Pinus. koraiensis and Juglans mandshurica Maxim., and elucidated the mechanisms underlying the changes in lignin, cellulose, and hemicellulose content during log decomposition [40]. Oberle et al. used an osmometer to measure the mechanical parameters related to the hardness of the decaying log and found that the relationship between log density and hardness varied with the coarse woody debris (CWD) values of three key factors: moisture content, species and the degree of decay [14]. This confirmed the usefulness of estimating the density of decaying logs using hardness values. However, they only observed samples from 1–3-year-old logs, which is not representative of the complete log decay process. Due to the long-term nature of the log decay process, only randomly selected log specimens with different levels of decay from the understory were used for our experimental research analysis in this study.
Existing studies have shown that logs have a good ratio of mechanical strength to density [41]. The present study innovatively introduced hardness as a physical characteristic for the classification of log decay classes and demonstrated the reliability of the clustering results through differential variations in representative physicochemical characteristics during the decay process. However, due to the long-term nature of the understory decay process and the random nature of the sampling process, the effects of the causes of decay, the diameter of the decayed log, heartwood and sapwood, the elemental contents and the environment were not taken into account in this study. Therefore, more factors should be considered and investigated in future studies to determine the effects of changes in hardness during the log decay process. Additionally, due to environmental factors during the decay process and human factors during sample collection, the effects of bark, sapwood and core wood on the lignin, cellulose and hemicellulose contents were also not considered in this study.
A complete set of hardness measurements for decaying logs has not yet emerged. Compared to standing timber, the hardness of fallen timber tends to decrease as the decay grade increases and the porosity within the fallen timber increases. At a certain moisture content, the decaying log has a softer texture. Some existing log hardness values have traditionally been known as Janka and Brinell hardness, but they are only suitable for laboratory conditions or for the determination of standing logs [42]. Oberle et al. measured the hardness of decaying logs using a similar instrument to the penetrometer and examined the mechanical properties of logs to illustrate how their mechanical properties change over time as the log decays [14]. It has also been shown that as dead logs are exposed to sunlight and extreme temperatures for long periods of time, they gradually develop a hard decay-resistant crust. To address this particular situation and the problem that hardness cannot be measured in the late stages of decay when specimens have been dissected, we developed an artificial neural network to predict the changes in the hardness of logs during the decay process using some basic physicochemical characteristics. Van Nguyen et al. proposed a novel experimental approach for the prediction of changes in the hardness of the log using artificial neural networks [43]. In comparison, training and validation data values were similar when the artificial neural network was used to predict the missing hardness values in this study, subjectively indicating that the predictions were reliable.
As the level of decay increased, the dry matter content, basic density, cellulose content, hemicellulose content and hardness of the log samples from the four species all showed downward trends, which was in agreement with the results of other studies [44,45,46,47]. The main reason for this was that, with the passage of time, the degree of decay increased, and the log was subjected to leaching, external forces and microbial decomposition, resulting in a gradual decrease in the material content of the log [20]. The acid-soluble and acid-insoluble lignin contents both increased with the increase in decay grade. This was possible because lignin degradation was slower due to the preferential degradation of cellulose by fungi; thus, the lignin contents increased in the decaying log [48,49,50].
It is well known that the hardness of the log depends on other characteristics of log chemistry and anatomy in addition to some basic physical characteristics, including the spatial distribution of internal cellulose micro-angles and voids, which vary considerably between species [51]. However, in this study, we did not consider the internal material pore structure of the decaying log from each species, which included both broadleaf and coniferous species from different elevations in the Qinling Mountains, although there were significant differences in the material contents between species. In the future, a more suitable method to determine the hardness of the decaying log should be developed to take into account the differences between tree species in different climatic zones in order to improve existing quantitative hardness classification systems.

5. Conclusions

The main objective of this study was to develop a classification system that was suitable for the quantitative classification of the decay grade of logs. The three basic indices of MAPE, RMSE and R2 were used to verify the prediction accuracy of log hardness values under different decay degrees and moisture contents using an artificial neural network model. By analyzing the correlations between and variations in the basic physicochemical characteristics of log samples with different decay grades, the feasibility of quantitatively judging the log decay grade using the log hardness value as the main clustering factor was supported. With the increase in the decay grade of the log, the basic density, dry matter content, cellulose content, hemicellulose content and hardness value of the log showed downward trends, while the acid-soluble and acid-insoluble lignin contents showed upward trends. The quantitative log decay grade classification system proposed in this study was more reliable than the current subjective classification system. In the future, more decaying samples from different tree species should be collected and examined to identify more factors that could be used to explore changes in the internal pore structure of logs with different degrees of decay so as to improve the applicability and reliability of this proposed system.

Author Contributions

Conceptualization, J.Y.; Data curation, W.W., W.Z., S.H., H.H. and X.W.; Investigation, W.Z., H.H., H.Q., W.W. and N.R.; Methodology, W.W., S.H. and W.Z.; Project administration, J.Y.; Software, W.Z. and S.H.; Writing—original draft, W.W. and W.Z.; Writing—review and editing, J.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the projects “Effect of downed log on forest soil respiration and its regulatory mechanism in the Qinling Mountains”(32271647) and “Respiration flux of downed log and its influence mechanism in the Qinling Mountains” (31800372) from the National Natural Science Foundation of China, and “Study on the key technology of forest carbon sink measurement and monitoring in the Qinling Mountains”(2023-YBSF-307) from the Key Research and Development Program in Shaanxi Province of China.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Harmon, M.E.; Franklin, J.F.; Swanson, F.J.; Sollins, P.; Gregory, S.V.; Lattin, J.D.; Anderson, N.H.; Cline, S.P.; Aumen, N.G.; Sedell, J.R.; et al. Ecology of coarse woody debris in temperate ecosystems. Adv. Ecol. Res. 1986, 15, 122–302. [Google Scholar]
  2. Brassard, B.W.; Chen, H.Y.H. Stand structural dynamics of North American boreal forests. Crit. Rev. Plant Sci. 2006, 25, 115–137. [Google Scholar] [CrossRef]
  3. Harmon, M.E.; Chen, H. Coarse woody debris dynamics in two old-growth ecosystems. Bioscience 1991, 41, 604–610. [Google Scholar] [CrossRef]
  4. Garrett, L.; Davis, M.; Oliver, G. Decomposition of coarse woody debris, and methods for determining decay rates. N. Z. J. For. Sci. 2007, 37, 227–240. [Google Scholar]
  5. Mackensen, J.; Bauhus, J.; Webber, E.J. Decomposition rates of coarse woody debris: A review with particular emphasis on Australian tree species. Aust. J. Bot. 2003, 51, 27–37. [Google Scholar] [CrossRef]
  6. Xu, H.D.; Wei, J.D.; Di, Y.N.; Qin, R.X.; Zhen, Z.L. Changes in Carbon Sequestration in Wood During Decay by Brown- and White-Rot Fungi. J. Biobased Mater. Bioenergy 2020, 14, 414–419. [Google Scholar] [CrossRef]
  7. Brown, S. Measuring carbon in forests: Current status and future challenges. Environ. Pollut. 2002, 116, 363–372. [Google Scholar] [CrossRef]
  8. Yuan, J. Study on Coarse Woody Residues of Main Forest Types in Huoditang Forest Area of Qinling Mountains. Ph.D. Thesis, Northwest A&F University, Xianyang, China, 2016. [Google Scholar]
  9. Yan, E.; Wang, X.; Huang, J. Concept and classification of coarse woody debris in forest ecosystems. Acta Ecol. Sin. 2005, 25, 158–167. [Google Scholar] [CrossRef]
  10. Maser, C.; Anderson, R.G.; Cromack, K.; Williams, J.T.; Martin, R.E. Dead and Down Woody Material. In Wildlife Habitats in Managed Forests of the Blue Mountains of Oregon and Washington; USDA Forest Service: Washington, DC, USA, 1979. [Google Scholar]
  11. Sollins, P. Input and decay of coarse woody debris in coniferous stands in western Oregon and Washington. Can. J. For. Res. 1982, 12, 18–28. [Google Scholar] [CrossRef]
  12. Spies, T.A.; Franklin, J.F.; Thomas, T.B. Coarse woody debris in Douglas-fir forest of western Oregon and Washington. Ecology 1988, 69, 1689–1702. [Google Scholar] [CrossRef]
  13. Larjavaara, M.; Muller-Landau, H. Comparison of decay classification, knife test, and two penetrometers for estimating wood density of coarse woody debris. Can. J. For. Res. 2010, 40, 2313–2321. [Google Scholar] [CrossRef]
  14. Oberle, B.; Dunham, K.; Milo, A.; Walton, M.; Young, D.; Zanne, A. Progressive, idiosyncratic changes in wood hardness during decay: Implications for dead wood inventory and cycling. For. Ecol. Manag. 2014, 323, 1–9. [Google Scholar] [CrossRef]
  15. Spänhoff, B. Simple Method for Rating the Decay Stages of Submerged Woody Debris. J. N. Am. Benthol. Soc. 2001, 20, 385–394. [Google Scholar] [CrossRef]
  16. Vlad, R.; Zhiyanski, M.; Dinca, L.; Sidor, C.; Constandache, C.; Gheorghe, P.; Ispravnic, A.; Blaga, T. Assessment of the density of wood with stem decay of Norway spruce trees using drill resistance. Comptes Rendus L’académie Bulg. Sci. Sci. Mathématiques Nat. 2018, 71, 1502–1510. [Google Scholar]
  17. Chambers, J.Q.; Schimel, J.P.; Nobre, A.D. Respiration from coarse wood litter in central Amazon forests. Biogeochemistry 2001, 52, 115–131. [Google Scholar] [CrossRef]
  18. Lutes, D.C.; Hardy, C.C. Lodgepole Pine Bole Wood Density and Decay Rate 1, 11, and 22 Years After Felling in Central Montana, United States. Front. For. Glob. Chang. 2021, 4, 687567. [Google Scholar] [CrossRef]
  19. Maser, C.; Trappe, J.M. The Seen and Unseen World of the Fallen Tree; Pacific Northwest Forest and Range Experiment Station, US Department of Agriculture, Forest Service: Portland, OR, USA, 1984.
  20. Yuan, J.; Zhang, S.X. Study on the variation characteristics of fallen wood density and water content of natural secondary Pinus tabulaeformis forest in Huoditang, Qinling Mountains. J. Cent. South Univ. For. Technol. 2012, 32, 105–109. [Google Scholar]
  21. Fujihira, M.; Nakamura, Y.; Isoda, N.; Hikita, Y. Relationship between decay resistance and changes of bending strength of structural lumber of wood-framing construction by fungal attack. Mokuzai Gakkaishi 1997, 43, 589–594. [Google Scholar]
  22. Maeda, K.; Ohta, M.; Momohara, I. Relationship between the mass profile and the strength property profile of decayed wood. Wood Sci. Technol. 2015, 49, 331–344. [Google Scholar] [CrossRef]
  23. Berg, B.; McClaugherty, C. Plant Litter: Decomposition, Humus Formation, Carbon Sequestration; Springer: Berlin, Germany, 2014. [Google Scholar]
  24. Talbot, J.M.; Treseder, K.K. Interactions among lignin, cellulose, and nitrogen drive litter chemistry–decay relationships. Ecology 2012, 93, 345–354. [Google Scholar] [CrossRef]
  25. Rutigliano, F.A.; De Santo, A.V.; Berg, B.; Alfani, A.; Fioretto, A. Lignin decomposition in decaying leaves of Fagus sylvatica L. and needles of Abies alba Mill. Soil Biol. Biochem. 1996, 28, 101–106. [Google Scholar] [CrossRef]
  26. Wu, A.M.; Yin, R.; Xu, Z.F.; Zhang, L.; You, C.M.; Liu, Y.; Li, H.; Wang, L.X.; Liu, S.N.; Zhang, Y.Q.; et al. Forest gaps slow lignin and cellulose degradation of fir (Abies faxoniana) twig litter in an alpine forest. Geoderma 2022, 424, 116010. [Google Scholar] [CrossRef]
  27. He, W.; Fu, Z.W.; Zhang, D.J.; Yang, W.Q.; Tan, B.; Zhao, Y.Y.; Wu, Q.Q. The effects of forest gaps on cellulose degradation in the foliar litter of two shrub species in an alpine fir forest. Plant Soil 2015, 303, 109–122. [Google Scholar] [CrossRef]
  28. He, W.; Wu, F.Z.; Yang, W.Q.; Tan, B.; Zhao, Y.Y.; He, M. Lignin Degradation in Foliar Litter of Two Shrub Species from the Gap Center to the Closed Canopy in an Alpine Fir Forest. Ecosystems 2016, 19, 115–128. [Google Scholar] [CrossRef]
  29. Tan, B.; Yin, R.; Yang, W.Q.; Zhang, J.; Xu, Z.F.; Liu, Y.; He, S.Q.; Zhou, W.; Zhang, L.; Li, H.; et al. Soil fauna show different degradation patterns of lignin and cellulose along an elevational gradient. Appl. Soil Ecol. 2020, 155, 103673. [Google Scholar] [CrossRef]
  30. Means, J.E.; Cromack, K.J.; MacMillan, P.C. Comparison of decomposition models using wood density of Douglas-fir logs. Can. J. For. Res.-Rev. Can. Rech. For.—Can. J. Forest Res. 1985, 15, 1092–1098. [Google Scholar] [CrossRef]
  31. Ganjegunte, G.K.; Condron, L.M.; Clinton, P.W.; Davis, M.R.; Mahieu, N. Decomposition and nutrient release from radiata pine (Pinus radiata) coarse woody debris. For. Ecol. Manag. 2004, 187, 197–211. [Google Scholar] [CrossRef]
  32. Rahman, M.M.; Tsukamoto, J.; Rahman, M.M.; Yoneyama, A.; Mostafa, K.M. Lignin and its effects on litter decomposition in forest ecosystems. Chem. Ecol. 2013, 29, 540–553. [Google Scholar] [CrossRef]
  33. Klotzbücher, T.; Kaiser, K.; Guggenberger, G.; Gatzek, C.; Kalbitz, K. A new conceptual model for the fate of lignin in decomposing plant litter. Ecology 2011, 92, 1052–1062. [Google Scholar] [CrossRef]
  34. Chen, Y.M.; Liu, Y.; Zhang, J.; Yang, W.Q.; Deng, C.C.; He, R.L. Cumulative cellulolytic enzyme activities and initial litter quality in prediction of cellulose degradation in an alpine meadow of the eastern Tibetan Plateau. J. Plant Ecol. 2020, 13, 51–58. [Google Scholar] [CrossRef]
  35. Sluiter, A.; Hames, B.; Ruiz, R.O.; Scarlata, C.; Sluiter, J.; Templeton, D. Determination of structural carbohydrates and lignin in biomass. Biomass Anal. Technol. Team Lab Anal. Proced. 2004, 2011, 1–14. [Google Scholar]
  36. Yi, Y.G.P. Property analysis of linear dimensionless methods. Stat. Res. 2008, 196, 93–100. [Google Scholar] [CrossRef]
  37. Haghbakhsh, R.; Adib, H.; Keshavarz, P.; Koolivand, M.; Keshtkari, S. Development of an artificial neural network model for the prediction of hydrocarbon density at high-pressure, high-temperature conditions. Thermochim. Acta 2013, 551, 124–130. [Google Scholar] [CrossRef]
  38. Han, Y.; Zhang, Y.; Jin, G. Effects of decay class and diameter class on moisture content and wood density in a typical mixed broadleaf-Korean pine forest. J. Nanjing For. Univ. Nat. Sci. Ed. 2020, 44, 133–140. [Google Scholar]
  39. Mu, Z.B.; Pan, H.; Wen, X.H.; Chao, Y.S.; You, W.B.; Liu, J.S.; Cai, C.T.; He, D.J. Effects of decay grades and diameter classes on physicochemical properties of woody debris in Tsuga longibracteata forest in Tianbaoyan National Nature Reserve. Shengtaixue Zazhi 2022, 41, 246–255. [Google Scholar] [CrossRef]
  40. Wang, L.H.; Zhang, X.L.; Xu, G.Q.; Xu, H.D.; Wu, J.Z. Using Lignin Content, Cellulose Content, and Cellulose Crystallinity as Indicators of Wood Decay in Juglans mandshurica Maxim. and Pinus koraiensis. Bioresources 2014, 9, 6205–6213. [Google Scholar] [CrossRef]
  41. Larson, P.R. Wood Formation and the Concept of Wood Quality; Yale University, School of Forestry: New Haven, CT, USA, 1969; pp. 1–54. [Google Scholar]
  42. Doyle, J.; Walker, J.C.F. Indentation of wood by wedges. Wood Sci. Technol. 1985, 19, 47–55. [Google Scholar] [CrossRef]
  43. Van Nguyen, T.H.; Nguyen, T.T.; Ji, X.; Lanh Do, K.T.; Guo, M. Using Artificial Neural Networks (ANN) for Modeling Predicting Hardness Change of Wood during Heat Treatment. IOP Conf. Ser. Mater. Sci. Eng. 2018, 394, 032044. [Google Scholar] [CrossRef]
  44. Shorohova, E.; Kapitsa, E.; Kuznetsov, A.; Kuznetsova, S.; de Gerenyu, V.L.; Kaganov, V.; Kurganova, I. Coarse woody debris density and carbon concentration by decay classes in mixed montane wet tropical forests. Biotropica 2022, 54, 635–644. [Google Scholar] [CrossRef]
  45. Di Lella, S.; Tognetti, R.; La Porta, N.; Lombardi, F.; Nardin, T.; Larcher, R. Characterization of Silver fir Wood Decay Classes Using Sugar Metabolites Detected with Ion Chromatography. J. Wood Chem. Technol. 2019, 39, 90–110. [Google Scholar] [CrossRef]
  46. Russell, M.B.; Fraver, S.; Aakala, T.; Gove, J.H.; Woodall, C.W.; D’Amato, A.W.; Ducey, M.J. Quantifying carbon stores and decomposition in dead wood: A review. For. Ecol. Manag. 2015, 350, 107–128. [Google Scholar] [CrossRef]
  47. Yuan, J.; Cheng, F.; Zhu, X.; Li, J.X.; Zhang, S.X. Respiration of downed logs in pine and oak forests in the Qinling Mountains, China. Soil Biol. Biochem. 2018, 127, 1–9. [Google Scholar] [CrossRef]
  48. Harmon, M.E.; Fasth, B.; Woodall, C.W.; Sexton, J. Carbon concentration of standing and downed woody detritus: Effects of tree taxa, decay class, position, and tissue type. For. Ecol. Manag. 2013, 291, 259–267. [Google Scholar] [CrossRef]
  49. Lombardi, F.; Cherubini, P.; Tognetti, R.; Cocozza, C.; Lasserre, B.; Marchetti, M. Investigating biochemical processes to assess deadwood decay of beech and silver fir in Mediterranean mountain forests. Ann. For. Sci. 2013, 70, 101–111. [Google Scholar] [CrossRef]
  50. Petrillo, M.; Cherubini, P.; Sartori, G.; Abiven, S.; Ascher, J.; Bertoldi, D.; Camin, F.; Barbero, A.; Larcher, R.; Egli, M. Decomposition of Norway spruce and European larch coarse woody debris (CWD) in relation to different elevation and exposure in an Alpine setting. iForest-Biogeosci. For. 2016, 9, 154–164. [Google Scholar] [CrossRef]
  51. Tze, W.T.Y.; Wang, S.; Rials, T.G.; Pharr, G.M.; Kelley, S.S. Nanoindentation of wood cell walls: Continuous stiffness and hardness measurements. Compos. Part A Appl. Sci. Manuf. 2007, 38, 945–953. [Google Scholar] [CrossRef]
Figure 1. The artificial neural network structure.
Figure 1. The artificial neural network structure.
Forests 14 00852 g001
Figure 2. A method of quantitatively dividing the level of log decay by hardness value.
Figure 2. A method of quantitatively dividing the level of log decay by hardness value.
Forests 14 00852 g002
Figure 3. The Q–Q plots of the relationships between the actual and predicted hardness values of the log samples.
Figure 3. The Q–Q plots of the relationships between the actual and predicted hardness values of the log samples.
Forests 14 00852 g003aForests 14 00852 g003b
Figure 4. The projections of the different ranges of hardness values among the different decay classes. The total number of samples involved in clustering was 184.
Figure 4. The projections of the different ranges of hardness values among the different decay classes. The total number of samples involved in clustering was 184.
Forests 14 00852 g004
Figure 5. The correlation analysis between the various indicators between the decay classes and tree species. * p < 0.05; ** p < 0.01; TS, tree species; DC, decay class; BD, basic density; ASL, acid-soluble lignin; AIL, acid-insoluble lignin; Ce, cellulose; He, hemicellulose; Glu, glucose; Xyl, xylose; H, hardness. Red labels indicate positive correlations and blue labels indicate negative correlations.
Figure 5. The correlation analysis between the various indicators between the decay classes and tree species. * p < 0.05; ** p < 0.01; TS, tree species; DC, decay class; BD, basic density; ASL, acid-soluble lignin; AIL, acid-insoluble lignin; Ce, cellulose; He, hemicellulose; Glu, glucose; Xyl, xylose; H, hardness. Red labels indicate positive correlations and blue labels indicate negative correlations.
Forests 14 00852 g005
Figure 6. Variations in the physicochemical factors of different tree species and different decay classes. BD, basic density; ASL, acid-soluble lignin; AIL, acid-insoluble lignin; Ce, cellulose; He, hemicellulose; Glu, glucose; Xyl, xylose; H, hardness. P, P. tabulaeformis; L, L. principis rupprechtii; B, B. albosinensis; Q, Q. aliena var. acuteserrata. Different small letters within the same decay class or capital letters within the same tree species are significantly different (p < 0.05).
Figure 6. Variations in the physicochemical factors of different tree species and different decay classes. BD, basic density; ASL, acid-soluble lignin; AIL, acid-insoluble lignin; Ce, cellulose; He, hemicellulose; Glu, glucose; Xyl, xylose; H, hardness. P, P. tabulaeformis; L, L. principis rupprechtii; B, B. albosinensis; Q, Q. aliena var. acuteserrata. Different small letters within the same decay class or capital letters within the same tree species are significantly different (p < 0.05).
Forests 14 00852 g006
Table 1. The evaluation results for the predicted changes in the hardness of logs under different levels of decay and different moisture content conditions.
Table 1. The evaluation results for the predicted changes in the hardness of logs under different levels of decay and different moisture content conditions.
Performance CriteriaData Sets
Training DataTesting DataAll
MAPE0.020690.021340.01716
RMSE1.02640.61870.9576
R20.99250.99640.9931
Table 2. ANOVA test for the differences in physicochemical properties among different tree species and decay classes (hardness clustering results). (p-value) TS, tree species; DC, decay class; BD, basic density; ASL, acid-soluble lignin; AIL, acid-insoluble lignin; Ce, cellulose; He, hemicellulose; Glu, glucose; Xyl, xylose; H, hardness.
Table 2. ANOVA test for the differences in physicochemical properties among different tree species and decay classes (hardness clustering results). (p-value) TS, tree species; DC, decay class; BD, basic density; ASL, acid-soluble lignin; AIL, acid-insoluble lignin; Ce, cellulose; He, hemicellulose; Glu, glucose; Xyl, xylose; H, hardness.
Indictors DMCBDASLAILCeHeGluXylH
p-Value
TS0.4960.4960.0410.0000.0000.0000.0000.0000.971
DC0.0000.0000.6680.7890.0000.2950.0000.0320.000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wen, W.; Zhang, W.; He, S.; Hu, H.; Qiao, H.; Wang, X.; Rao, N.; Yuan, J. The Classification of Log Decay Classes and an Analysis of Their Physical and Chemical Characteristics Based on Artificial Neural Networks and K-Means Clustering. Forests 2023, 14, 852. https://doi.org/10.3390/f14040852

AMA Style

Wen W, Zhang W, He S, Hu H, Qiao H, Wang X, Rao N, Yuan J. The Classification of Log Decay Classes and an Analysis of Their Physical and Chemical Characteristics Based on Artificial Neural Networks and K-Means Clustering. Forests. 2023; 14(4):852. https://doi.org/10.3390/f14040852

Chicago/Turabian Style

Wen, Wen, Wenjun Zhang, Shirong He, Haitao Hu, Hailiang Qiao, Xiao Wang, Nan Rao, and Jie Yuan. 2023. "The Classification of Log Decay Classes and an Analysis of Their Physical and Chemical Characteristics Based on Artificial Neural Networks and K-Means Clustering" Forests 14, no. 4: 852. https://doi.org/10.3390/f14040852

APA Style

Wen, W., Zhang, W., He, S., Hu, H., Qiao, H., Wang, X., Rao, N., & Yuan, J. (2023). The Classification of Log Decay Classes and an Analysis of Their Physical and Chemical Characteristics Based on Artificial Neural Networks and K-Means Clustering. Forests, 14(4), 852. https://doi.org/10.3390/f14040852

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