Multi-Factorial Load Analysis of Pine Sawlogs in Transport to Sawmill

: This study analyzed the variability of truckloads of large-sized pine logs transported to a furniture manufacturing mill, depending on the origin of the timber, delivery period, length of the transported logs, volume of the load and its mass. It was assumed that the volume and mass of the transported timber load depends on the season of the year and the origin of the timber, which has a signiﬁcant impact on the mass of one cubic meter of the load. The analysis of the wood origin parameters (location of its growth) took into account the type of forest habitat, soil type and stand abundance. The characteristics of the tree parameters from which the roundwood was obtained took into account age, diameter at breast height (DBH) and height of the stand. The origin of the timber was determined for 1063 timber loads from 40 forest districts transported over a 12 month period. We obtained an average timber load volume of 29.34 m 3 which includes single shipments di ﬀ ering signiﬁcantly in their volume of between 24.21 and 36.51 m 3 . The origin of the timber, habitat conditions and stand parameters of the harvesting site inﬂuence load mass and the mass of 1 m 3 of load. The delivery date has an impact on the studied elements of the sawlogs load, with two similar periods—spring and summer. The volume of the transported timber is already at a stabilized level in accordance with legal regulations, which is no longer dependent on other factors and only on the accepted conversion factors for roundwood density and the mass of the empty vehicle set. This has not inﬂuenced a reduction in the mass of the transported timber.


Introduction
Based on the Statistical Yearbook of Forestry 2019, the total forest area in Poland is 9242 thousand ha. The area of all coniferous trees is 6309 thousand ha (68.4% of total forest area). The main tree species is Scots pine (Pinus Sylvestris) with an area of 5366 thousand ha (58.2% of total forest area) [1]. The natural potential and proper management of forests allows the increase of the harvest of raw material, to secure increasing demand from wood industry. In 2018, State Forestssold 47.7 million m 3 of total timber volume, of which sawmill coniferous timber (softwood saw timber) constituted 14.5 million m 3 [1]. In terms of raw material, large-size pine wood is one of the basic assortments sold by SF and used by the wood industry in Poland.
There are only 29 companies, which buy over 100 thousand m 3 of wood. Twenty processed up to 500 thousand m 3 . Two companies have processing capacity at the level of 500-1000 thousand m 3 . Four companies process over 1 million m 3 [2]. In 2018, 30 of the largest sawmill companies in Poland processed 6440 thousand m 3 , including 7 sawmills processing over 250 thousand m 3 [3]. However, 70% of the State Forests' contractors, who buy 9% of the total sold volume, are companies processing less than 2000 m 3 of wood annually [2]. In 2011, about 309,000 persons were employed in about 67,000 wood sector companies in Poland [4]. Among the entities involved in processing large-sized pine logs assortment originates. The butt end is characterized by the highest density at 816 kg·m −3 , and the lowest values are recorded from the top parts of the stem at 707 kg·m −3 [32].
In order to unequivocally determine the mass of transported raw wood material, a resolution providing a normative wood density for each species was introduced in Poland in 2012, and the mass of a load is defined as result of multiplying the volume and the determined wood density for inspections by the relevant public road services [33]. In fact, this resulted in increasing the average load of large-sized pine wood for a single transport from 25.57 m 3 in 2012 to 29.14 m 3 in 2016 [34].
The current availability of information allows a more extensive analysis of this issue by also adding elements of the taxonomic description of the forest from which the timber originates. Habitat elements, soil conditions, age, height or site index, which may have a potential impact on the delivery of large-sized pine logs to the large processing plant belonging to an international furniture manufacturer.
We wanted to present in the article an analysis of various factors affecting the volume and mass of the transported pine sawlogs. The scope of research included seasons of the year and parameters of the forest stand, such as habitat, soil type, height, dbh, and so forth, which has a significant impact on the mass of one cubic meter of the load (kg·m −3 ).
The scope of work to achieve the research aim and assumptions formulated in this way included: • analyzing the timber delivery to the recipient at different times of the year, • analyzing the volume of timber in a unit of load (truck load), • determining the mass of the load, • analyzing the number of logs in a unit of load, • determining the parameters of the forest habitat from which the timber was harvested, • determining the parameters of the trees from which the timber was harvested.

Material and Methods
Appropriate analyses were performed on the deliveries of the raw material to one of the largest recipients of timber in Poland, a mill purchasing sawlogs pine timber, both tree-length and shortwood logs. The analyses were conducted from 1 April 2016 to 29 March 2017 on the premises of the plant processing sawlogs pine timber in the Warmińsko-Mazurskie Voivodeship. During the analyzed period, 310,207 m 3 of sawlogs pine timber was delivered to the plant in 10,800 deliveries of wood. Randomly selected deliveries of wood were analyzed; the unit load on vehicle of the deliveries was obtained from a single location, separated. The total number of observations was 1514 deliveries. Nevertheless, a description of the stand could be obtained for 1063 deliveries, where the supplied wood came from one harvest operation in the same place. The remaining 451 deliveries were combined loads from more than one place. In the conducted research, we analyzed 1514 truckloads of sawlogs pine timber from 40 forest districts, which were transported by 59 different vehicle units (e.g. truck with trailer, truck with semi-trailer) in specific months of the year. The timber from the forest was delivered by 11 private transport companies contracted by the consignee (wood processing mill). In the analyzed transports from January (100 observations) and in subsequent months (observations in February-39, March-431, April-176, May-101, June-118, July-10, August-239, September-25, October-270, November-83) until in December there were 12 observations, and 44,429 m 3 of wood was delivered to the mill.
Large-size wood (L), is wood with a thin end diameter of 14 cm (excluding bark), calculated in single pieces. In terms of quality and size, large-size wood is divided into four classes-WA, WB, WC, WD and into 2 sub-classes-general purpose wood (0) and special purpose wood (1). The large-size general purpose wood is comparable to the assortment defined as sawmill wood. The markings (WA, WB, WC, WD) refer to the quality of the raw material, where WA is the best class and WD the worst. [35].
The timber loads were transported by external companies contracted by the mill. The delivery date was determined from the recipient's (mill's) documents as well as the shipping documents issued by the forest inspectorate. Information given in the Wood Delivery Note, verified in a sawmill during wood receptionnumber of logs in the load (pieces), length of logs in the load (m), load volume (m 3 ). Calculated-mean volume of one log (m 3 ), mass of 1 m 3 of load (kg). The number of pieces of sawlogs and tree-length logs as well as their length per unit of load was determined on the basis of the shipping documents describing each piece of the timber load, providing its length, diameter and volume. The volume (load volume) of the transported timber in cubic meters was determined on the basis of the shipping documents, which form the basis for accounting for the timber delivered from supplier to consignee after verification. At the same time, it should be noted that when timber is sold in Poland and subsequently transported in Poland, the reported volumes do not include bark. The share of bark is significant and may be as high as several dozen percent [36,37] of the pine stem, affecting the mass of the transported cargo. Therefore, in the analyses, we discuss the mass of the whole load (wood, bark, moisture) or the mass of one cubic meter of the load (kg·m −3 ).
The actual mass of each load (kg) was determined from the difference between the gross vehicle mass plus load (GVW) and the empty vehicle mass after unloading (tare). The mass of the vehicle unit was obtained by weighing it on the stationary scale of the mill at the time of delivery of the raw material, and then after unloading, the empty vehicle unit was weighed again (tare).
The analysis of the parameters of the timber's origin (place where it grew) took into account such factors as forest habitat type, soil type and stand abundance. The characteristics of the parameters of the trees from which the timber was harvested took into account the age, diameter at breast height (DBH) and stand height. The forest inspectorate specifies the harvesting site on the shipping documents by providing the forest address from the State Forests Information System (SILP). The forest address is unambiguously specified by the code of the regional directorate, forest district, precinct, forester area, division up to the smallest area unit-plot. Having precise data on the origin of the harvested timber, the characteristics of the forest habitat and stand were determined on the basis of the Forest Data Bank (BDL: https://www.bdl.lasy.gov.pl/portal/mapy-en), which contains forest survey descriptions from the management operations of the forest inspectorates. The factors characterizing the origin of the timber in the loads are stand parameters such as age, DBH and height. To perform the analyses, the results were grouped for stand age (years), DBH (cm), tree height (m) and abundance (m 3 ·ha −1 ) into categories: • Some of the distinguished stand groups were not included in the statistical analyses due to the small number of their transports, less than 20.
Based on the literature [11,12,26], we know that the volume and mass of loads in the transport of freshly felled wood differs, therefore the mass of one cubic meter of the load, calculated as the quotient of the load mass (kg) and its volume (m 3 ), was used for the comparative analyses of the effect of timber origin.
The obtained results were statistically analyzed using the STATISTICA 12 package. In all analyzed periods, the distributions of the variables for all parameters deviated from the normal distribution. Therefore, the significance of the differences was determined mainly by using the Kruskal-Wallis and Dunn tests. The Mann-Whitney test was conducted to examine the significance of selected characteristics depending on the procedure in which the wood was harvested (pre-final cut, final cut). Due to the heterogeneity of variance, analyses were performed with the Kruskal-Wallis test on the significance of differences in load volume (m 3 ), load mass (kg), one cubic meter of load (kg·m −3 ) and the number of logs (pieces) in a timber load for individual months of the year and depending on forest habitat and soil. Additionally, the multiple comparison of mean ranks test (Dunn tests) was conducted, and the analyses were performed at a significance level of 0.05. The dependencies analyzed were for the date of delivery in the month, the tested habitat and stand characteristics, as well as for specific seasons (in accordance with the calendar). The need to purchase wood from many forest districts results in a significant diversity of timber origin, therefore we also checked the significance of habitat characteristics and timber origin depending on delivery date. For the statistical tests Kruskal-Wallis, Dunn and Manna-Whitney were not taken for some of the observations groups (e.g. of the month (July i December), stand age) due to the low number of results in a specific group (less than 20).
In order to determine the strength of the relationship (correlation) between the analyzed parameters (load volume, load mass, mass of 1 m 3 load) and parameters describing forest stand (stand age, DBH, tree height as well as habitat abundance and habitat), a Spearman rank correlation test was performed.

Characteristics of the Origin of the Timber
The delivered sawlogs pine timber came from 40 forest districts and the following forest habitat types: The most common soils in the analyzed habitats from which the timber was harvested were RDb-rusty podzolic soils (436 plots), RDw-(proper) podzolic soils (332), RDbr-rusty brown soils (160), Bw-(proper) podzolic soils (66) and Bgw-(proper) gley-podzolic soils (17). Only a few observations (from 3 to 7) were recorded of various other soil types for the remaining 17 transports, thus further analyses were performed only on the transported timber from the soil types enumerated above.
In the case of 873 transports, the timber was cut based on final harvesting, the remaining cases were pre-final cuttings.
The characteristics of the harvested stands are presented in Table 1. The average age of the trees was 87.79 years with a span of 30 to 191 years (SD = 24.33). The average height was 26.20 m (SD = 3.17), with the lowest stand being 14 m high and the highest − 36 m. The average DBH of the trees was 33.42 cm (determined on the basis of the forest survey descriptions). The forest areas from which the timber was harvested were characterized by an average abundance of 327.10 m 3 ·ha −1 with SD = 94.19. The statistical analysis with the Kruskal-Wallis test indicated a significant difference (p < 0.05) between the characteristics of the wood's origin (Abundance and DBH p = 0.0000, Age p = 0.0008, Tree height p = 0.0004) from 40 forest districts depending on the month of delivery.
The multiple comparison of mean ranks test showed differences in stand abundance for deliveries between the months of January with June March with April, June and April with June and May with November and November with April, May, June. The timber deliveries in February, from July to October, and in December do not differ in terms of the types of forest habitats from which the timber originated. Tree heights are very similar and the difference is only between the April and October deliveries. The DBH of harvested trees, despite the difference shown, is only statistically significant for deliveries between March and October as well as April and October. In general, significant differences for tree age were also found (p = 0.0008), however, the multiple comparison of mean ranks test showed statistically significant differences between deliveries in December and March, April, May, June, August, October and November. On this basis, it can be concluded that the wood of the analyzed deliveries originates from sites having very similar habitat and stand conditions.

Characteristics of the Sawlogs Pine Timber Loads
The mill acquires shortwood or tree-length pine logs. The observed transports included shortwood logs of: 3.7 m (12 transports), 4.0 m (6), 4.4 m (588), 5.0 m (816) and tree-length timber of 10 m (3), 13 m (7), 14 m (9). Tree-length logs were delivered in 47 transports, while 23 transports were recorded during the study with mixed log lengths and tree-length logs (described hereinafter as a mix). Currently, the mill is receiving mainly 5.0 m logs. The average values of the parameters characterizing timber loads are presented in Table 2. The analyzed timber loads were characterized by an average volume of 29.34 m 3 with very small differences (SD = 1.53 and a median of 29.14), but individual shipments occurred with significantly different volumes, from a minimum of 24.21 m 3 to a maximum of 36.51 m 3 . The transported cubic meters of timber determined the resulting load mass, where the average was 30,330 kg, broadly ranging from 21,750 kg to 38,910 kg, which with the same tree species as well as assortment (sawlogs logs) indicates that other factors are influencing the load mass than just volume. A reliable indicator for analyses and comparisons of variable loads (volume and mass) can be the mass of 1 m 3 of wood in the load, which amounted to 1,033.87 kg·m −3 , with similar results in the population at SD = 57.87, and a median of 1,031.61 kg·m −3 . There was an average of 146 logs (shortwood or tree-length) in one load, and the range of the results was from a minimum of 28 logs (the transport of only tree-length pieces) to 360 logs (3.7 m log lengths from a pre-final cut).

Analysis of the Factors Influencing the Load Parameters and Their Significance
A Mann-Whitney test was used to analyze the significance of selected characteristics of the loads depending on the type of cut from which the wood was obtained (pre-final and final cuts). The statistical analysis at the assumed significance level of 0.05 shows no significant differences between the loads only for the volume of transported timber (p = 0.65108). The lack of differences in load volume may be the result of the introduced wood conversion factors [30], where drivers, knowing the mass of an empty vehicle unit, take such a load volume that after applying the conversion factor of 740 kg·m −3 , the total mass of the vehicle unit does not exceed the legally permissible 40,000 kg.
We analyzed the significance of the differences in the tested characteristics of the load depending on the length of the transported timber. Due to the small number of observations (from 3 to 9) for  We analyzed the significance of the differences in the tested characteristics of the load depending on the length of the transported timber. Due to the small number of observations (from 3 to 9) for log lengths of 4.0, 10.0, 13.0 and 14.0 m in the loads, they were excluded from the Kruskal-Wallis test. Statistically significant differences were found for load volume (p = 0.0000), mass of 1 m 3 of load (p = 0.0001), number of logs (pieces) in the load (obtained p = 0.000) and load mass (p = 0.0059) ( Figures  1a, 1b, 1c, 1d). The diagrams for given load volumes (Figure 1a), mass of 1 m 3 of wood in a load ( Figure 1c ) and especially for the number of logs (pieces) (Figure 1d In accordance with the methodology, analyses were performed of the significance of the examined load characteristics depending on the delivery date. The comparison of the values obtained for particular characteristics is presented in Figures 2a to 2d. The Kruskal-Wallis test showed statistically significant differences in the tested characteristics of the load depending on the delivery date, with a statistical value at 0.0000 ≤ p ≤ 0.0281 (0.0281 for the volume of timber in the load). Comparing the load volume in individual months with the multiple comparison of mean ranks test, no differences were found in relation to the month of transporting the load (Figure 2a), which is explained by the use of the wood conversion factors when the timber is issued by the foresters, and transports are conducted by the same vehicle units. The results for July and December were not included in the statistics due to the small (10-12) observation samples. Statistically significant differences for load mass and 1 m 3 of load mass were mainly found between particular groups of In accordance with the methodology, analyses were performed of the significance of the examined load characteristics depending on the delivery date. The comparison of the values obtained for particular characteristics is presented in Figure 2a-d. The Kruskal-Wallis test showed statistically significant differences in the tested characteristics of the load depending on the delivery date, with a statistical value at 0.0000 ≤ p ≤ 0.0281 (0.0281 for the volume of timber in the load). Comparing the load volume in individual months with the multiple comparison of mean ranks test, no differences were found in relation to the month of transporting the load (Figure 2a), which is explained by the use of the wood conversion factors when the timber is issued by the foresters, and transports are conducted by the same vehicle units. The results for July and December were not included in the statistics due to the small (10-12) observation samples. Statistically significant differences for load mass and 1 m 3 of load mass were mainly found between particular groups of months: January, February, March and November (the winter period) and April, May, June, August, September (spring, summer) (Figure 2b,c). For the number of sawlogs or tree-length logs in one delivery, the transport date is significant (statistically significant differences at p = 0.0000). Differences also occur between single pairs of months: March and June as well as June with January and March (Figure 2d).  (Figures 2b and 2c). For the number of sawlogs or tree-length logs in one delivery, the transport date is significant (statistically significant differences at p = 0.0000). Differences also occur between single pairs of months: March and June as well as June with January and March (Figure 2d).  The literature contains analyses of timber transports by season [12,29,30]. We analyzed changes in the values of the characteristics of timber loads depending on the season. Statistically significant differences (p = 0.0000) occur between the value of load mass as well as the mass of 1 m 3 of load and the number of logs (pieces) in the load. The Kruskal-Wallis test showed no differences for load volume (p = 0.1807) (Figure 3a). The load mass (Figure 3b) and 1 m 3 mass (Figure 3c) do not differ statistically only between spring and summer (multiple comparison of mean ranks test). The number of logs in a transport differs only between the data for spring and winter (Figure 3d). The literature contains analyses of timber transports by season [12,29,30]. We analyzed changes in the values of the characteristics of timber loads depending on the season. Statistically significant differences (p = 0.0000) occur between the value of load mass as well as the mass of 1 m 3 of load and the number of logs (pieces) in the load. The Kruskal-Wallis test showed no differences for load volume (p = 0.1807) (Figure 3a). The load mass (Figure 3b) and 1 m 3 mass (Figure 3c) do not differ statistically only between spring and summer (multiple comparison of mean ranks test). The number of logs in a transport differs only between the data for spring and winter (Figure 3d).
The origin of timber loads from different types of forest habitats influences the results obtained characterizing the loads except for load volume -p = 0.3516 and load mass -p = 0.0201 (Figure 4a). Statistically significant differences occur for the load mass of timber from Bśw compared to BMw and LMśw (Figure 4b). Similar differences occur for the mass of 1 m 3 of the load from Bśw compared to BMśw, LMśw and Lśw and between BMśw and LMśw (Figure 4c). The number of logs (pieces) in the load also statistically differs between timber originating mainly from Bśw (different from LMśw and BMw) and BMśw (different from BMw and LMśw) (Figure 4d). The origin of timber loads from different types of forest habitats influences the results obtained characterizing the loads except for load volume -p = 0.3516 and load mass -p = 0.0201 (Figure 4a). Statistically significant differences occur for the load mass of timber from Bśw compared to BMw and LMśw (Figure 4b). Similar differences occur for the mass of 1 m 3 of the load from Bśw compared to BMśw, LMśw and Lśw and between BMśw and LMśw (Figure 4c). The number of logs (pieces) in the load also statistically differs between timber originating mainly from Bśw (different from LMśw and BMw) and BMśw (different from BMw and LMśw) (Figure 4d).  Analyses were performed of the influence of forest soil types occurring in the areas from which the timber loads originated on the examined characteristics. Loads from the five most frequent forest soil types (Bw, Bgw, RDb, RDbr, RDw) were analyzed ( Figure 5), and no statistically significant differences were found for the results of load volume (p = 0.4726) and load mass (p = 0.634) depending on soil type. Figures 5a, 5b, 5c, 5d show small changes in median values for the examined characteristics and for the median results ± standard error. Statistically significant differences were found for the results of the mass of 1 m 3 of load and the number of logs (pieces) between RDbr soils and RDb, RDw as well as Bw soils (Figures 5c, 5d). Analyses were performed of the influence of forest soil types occurring in the areas from which the timber loads originated on the examined characteristics. Loads from the five most frequent forest soil types (Bw, Bgw, RDb, RDbr, RDw) were analyzed ( Figure 5), and no statistically significant differences were found for the results of load volume (p = 0.4726) and load mass (p = 0.634) depending on soil type. Figure 5a-d show small changes in median values for the examined characteristics and for the median results ± standard error. Statistically significant differences were found for the results of the mass of 1 m 3 of load and the number of logs (pieces) between RDbr soils and RDb, RDw as well as Bw soils (Figure 5c,d). The last factors characterizing the origin of the timber in the loads are stand parameters such as age, DBH and height. Some of the distinguished stand abundance groups (10-50, 51-100 or 551-600) were not included in the statistical analyses due to the small number of their transports, less than 20. The statistical analysis (Kruskal-Wallis test at p < 0.05) showed that there were no statistically significant differences between the results of load volume and the abundance value of stands from which the timber was harvested for the remaining characteristics. The results of the assessment of the significance of the differences are presented in Table 3. Table 3. The results of the assessment of the significance of the differences of the selected load characteristics depending on the parameters describing the stand.

Measure
Timber load volume (m 3 )  The last factors characterizing the origin of the timber in the loads are stand parameters such as age, DBH and height. Some of the distinguished stand abundance groups (10-50, 51-100 or 551-600) were not included in the statistical analyses due to the small number of their transports, less than 20. The statistical analysis (Kruskal-Wallis test at p < 0.05) showed that there were no statistically significant differences between the results of load volume and the abundance value of stands from which the timber was harvested for the remaining characteristics. The results of the assessment of the significance of the differences are presented in Table 3. Table 3. The results of the assessment of the significance of the differences of the selected load characteristics depending on the parameters describing the stand. The multiple comparison of mean ranks test for the value of timber depending on stand height as well as for the number of logs in the load depending on abundance showed no statistically significant differences between the results. Statistically significant differences in the mass of 1 m 3 of load occurred depending on stand age in the category of 91-100 years and those from 60 to 90 years, as well as between 111-120 years and 71-80 and 121-130 years (Figure 6a). The mass of 1 m 3 of load also differs depending on the DBH of the tree between groups of 21-30 cm and 41-50 cm (Figure 6b) and in all height groups (Figure 6c). Due to stand abundance, the mass of 1 m 3 of load varies between the 401-450 value group and those with values of 151 to 350 m 3 ·ha −1 (3 groups) as well as the 451-500 value group with those of 101-400 m 3 ·ha −1 (6 groups) (Figure 6d). The stand age from which the timber originated also influences the differences in load volume between the groups aged 81-90 years and 101-110 years (this is the only statistically significant difference), as well as the load mass between the group aged 71-80 years and those in the 91-100 and 111-120 age groups. Figure 6 shows the distribution of the obtained results of the mass of 1 m 3 (kg) of load depending on the parameters of the stands from which the timber was obtained.  (Figure 6d). The stand age from which the timber originated also influences the differences in load volume between the groups aged 81-90 years and 101-110 years (this is the only statistically significant difference), as well as the load mass between the group aged 71-80 years and those in the 91-100 and 111-120 age groups. Figure 6 shows the distribution of the obtained results of the mass of 1 m 3 (kg) of load depending on the parameters of the stands from which the timber was obtained. Statistically significant relationships (correlations) were obtained between the analyzed parameters, load mass, mass of 1m 3 , and parameters describing the origin of the wood, except the age of the trees (p < 0.05). For the volume of load, only significant dependence on the number of pieces of logs in the load was obtained ( Table 4). The mass of the load depends mainly on volume (0.617103). The number of logs in the load is correlated with the parameters of the tree from which the timber Statistically significant relationships (correlations) were obtained between the analyzed parameters, load mass, mass of 1m 3 , and parameters describing the origin of the wood, except the age of the trees (p < 0.05). For the volume of load, only significant dependence on the number of pieces of logs in the load was obtained ( Table 4). The mass of the load depends mainly on volume (0.617103). The number of logs in the load is correlated with the parameters of the tree from which the timber assortment was obtained, but these are not significant dependencies (0.278936-0.429014).

Load Mass as a Function of Load Volume
The main problem in organizing the transport of timber is that it is not possible to unequivocally determine its mass. The research conducted concerned the same tree species and a largely limited assortment (in terms of length and diameter) as well as timber origin (see Section 3.1. Statistical analysis of the characteristics of the origin of the timber). A regression model was developed for load mass M (kg) as a function of load volume V (m 3 ). The assessment of the model parameters is presented in Table 4. For the developed model, the standard error of the estimate is 1689.2 and the coefficient of determination r 2 is equal to 0.4126 ( Table 5). The developed model took the following form: where: M − load mass (kg); V − load volume (m 3 ) Table 5. Assessment of model parameters on load mass. We compared the calculated mass of the load according to the presented relationship. Calculated mass as multiplying the volume of the load and the factor of 740 kg·m −3 (legally determined [33]) and real mass obtained from weighing in the plant. We received a statistically significant linear relationship (p < 0.05) for the mass of the load depending on volume (Figure 7). As is presented in Figure 7, the calculated load mass as the product of multiplying the volume of the load and the factor of 740 kg⋅m −3 (legally determined [33]) is significantly understated in relation to the real load masses. The obtained values of load mass on the basis of formula (1) are similar to the real load mass, where the average are 30,329 kg (median 30,145) and 30,330 kg (median 30,200), respectively. The Mann-Whitney test analysis showed no statistically significant differences between the real mass of the load obtained from weighing and calculated on the basis of the obtained formula (p = 0.245116). Of course, the mass calculated on the basis of the conversion factor 740 is very different from the others.
The weight of the load depends mainly on its volume, but also on the characteristics of the origin of the wood. The greater strength of the relationship we get for one cubic meter of load, which is still dependent on the type of forest habitat and the age of the trees. This is confirmed by other authors dealing with the impact of forest stand conditions on tree parameters [36,37,[39][40][41]. The number of logs in the load mainly depends on the size of the transported wood, but it is a derivative of the age, diameter and height of the tree from which the load was obtained. So from older and larger trees we get a greater assortment (which is understandable) and the result is a negative, significant correlation (−0.429014).
Differences in the mass of the load in individual months, with indicated lack of differences in the volume of the load of wood occuring mainly between the winter (heaviest loads) and early As is presented in Figure 7, the calculated load mass as the product of multiplying the volume of the load and the factor of 740 kg·m −3 (legally determined [33]) is significantly understated in relation to the real load masses. The obtained values of load mass on the basis of formula (1) are similar to the real load mass, where the average are 30,329 kg (median 30,145) and 30,330 kg (median 30,200), respectively. The Mann-Whitney test analysis showed no statistically significant differences between the real mass of the load obtained from weighing and calculated on the basis of the obtained formula (p = 0.245116). Of course, the mass calculated on the basis of the conversion factor 740 is very different from the others.
The weight of the load depends mainly on its volume, but also on the characteristics of the origin of the wood. The greater strength of the relationship we get for one cubic meter of load, which is still dependent on the type of forest habitat and the age of the trees. This is confirmed by other authors dealing with the impact of forest stand conditions on tree parameters [36,37,[39][40][41]. The number of logs in the load mainly depends on the size of the transported wood, but it is a derivative of the age, diameter and height of the tree from which the load was obtained. So from older and larger trees we get a greater assortment (which is understandable) and the result is a negative, significant correlation (−0.429014).
Differences in the mass of the load in individual months, with indicated lack of differences in the volume of the load of wood occuring mainly between the winter (heaviest loads) and early autumn and summer, which is caused by vegetation of trees. Comparing the load masses depending on the month of delivery (Dunn's test) and the mass of one cubic meter we get the same differences, this indicates that conversion factors can be develop for wood from a specific area.
When transporting timber, the high variability of species, as well as assortments and moisture content of the wood, impact the determination of the mass of transported timber and GVW-gross vehicle weight [11,26,28], which largely depends on the load volume and its density [11,15,29]. Our study shows that, currently, the volume of the wood is similar in the entire analyzed period (no statistically significant differences) and does not depend on season of the year or origin of the wood. The lack of statistically significant differences in load volume when analyzing transports using semi-trailer trucks was demonstrated by Trzciński and others [12], at the same time indicating significant differences (however of little practical significance) between the results from autumn and the remaining seasons of the year for vehicle units with a trailer. The load volume slightly differed for transports in the autumn and summer [30]. The obtained average volume of 29.34 m 3 (SD = 1.53) is similar to the research results of other authors of 29.47 m 3 [30] or 29.07-30.29 m 3 [12]. The value of the average load volume obtained from the study of 1514 transports is slightly higher (0.72 m 3 ) than the average of 28.72 m 3 for all 10,800 deliveries transported during the year to the mill, which provided 310,207 m 3 of timber. The lack of variability in load volume shown in our study results from the current legal regulations in Poland [33], where a mass of one cubic meter of wood converter was introduced to regulate timber transports, which is 740 kg·m −3 for pine. A forester releases an amount of wood to the load to ensure that the total mass of a vehicle unit (GVW) does not exceed 40,000 kg, which, with a limited number of vehicle types delivering wood to a specific recipient (known set mass − tare) determines the variability of load volume [12,[16][17][18]30]. The empty mass (tare), ranging from 13,800 to 23,700 kg, was observed for the truck and semi-trailer, trailer 20,224 kg, platform 14,941 kg and dolly 20,071 kg [12]. Introducing mass of one cubic meter of wood coefficients [33] resulted in an increase in the average volume of timber loads [34]. It is very practical in wood transport, where an increase in the demand for wood in a sawmill by 40% caused an increase in the number of courses by only 23% [34]. At the same time, overloading the car may contribute to the risk of traffic safety, damage to the car and administrative penalties resulting from legal provisions [42].
One of the main problems in this field is to determine the mass of 1 m 3 of the load, which would allow transport companies to decide on the volume of the timber load using as much DCM as possible and the capacity of the vehicle units, when timber is delivered from a wide range of different areas to consignees. Wood is a natural material which differs greatly in origin [39,43], density [44], moisture content [43][44][45], but it also varies depending on the part of the tree it comes from or seasonal variations [45][46][47]. Another factor influencing load mass is the density of 1 m 3 of the load and not the density of the wood, because in the study, the load also includes bark-which is not quantified by the forester in the load mass-as well as the moisture content of the wood or even at times, other factors, such as soil, snow, or ice on the logs. The share of bark in the biomass of pine stems depends on the age and origin of the wood and is defined at the significant level of 9.3% [36], or even 13.5% or 17% [48][49][50]. Load density (mass of 1 m 3 ) is strongly influenced by the parameters [38,39,42,[44][45][46] of the timber. The mass of 1 m 3 of load was on average 1033.87 kg at SD = 57.87 kg and is consistent with previously published studies [30] and with the 979 kg reported by Tomczak and others [37]. In many cases, significant differences in the mass of 1 m 3 of load were found, depending on the date of transport or the origin of the wood, due to the forest habitat or stand age and its parameters, which corresponds to the dependencies presented in the literature [15,29,43]. The ratio of the actual mass of one cubic meter (1033.87 kg·m 3 ) to the mass provided in the tables [33] (740 kg·m 3 ) amounted to 1.4:1, is similar to previous published studies [34]. The real mass of the timber was about 40% higher than the mass, which can be estimated on the basis of the wood density tables developed for the road transport needs. The results of the mass of 1 m 3 of load in the range from 763.15-1,343.07 kg·m 3 (  (Table 2), gives the level of overloading of 12,796-19,298 kg. These results are similar to other authors [12,16,22,26] and therefore DCM has been raised in some countries [17,18,20]. In this situation, the presented formula can be used in practice, where with the known mass of the empty set (truck units weigh in many sawmills), during loading, a safe amount of loaded wood can be calculated to not overload the truck units. This pattern obviously works well for loads of wood from the analyzed State Forest Districts, and further research with wood from another region of Poland can verify it.
The origin of wood from various regions of the country and forest habitats has a huge impact on its characteristics, which is confirmed by studies of the physical properties of wood [39][40][41][43][44][45][46][47]50].
However, while no variability was found for load volume of the timber transports, load mass is varied, having an average of 30,330 kg with a range of 21,750 to 38,910 kg (SD = 2203) and this depends on the season of the year or the origin of the timber [29,43]. Because of the actual load mass and its variability, transport vehicles are often loaded with more than the allowable total mass [11,15,16,26,30,51], and introducing the legally prescribed conversion factors for individual tree species [33] resulting in an increased load mass and gross vehicle weight [30], which our research confirms. This is the result of not taking into account the content of water and bark in wood (as described above) and the origin of wood from various habitats in the wood mass conversion [33].
Drivers who know the mass of a vehicle unit to be 15,000-20,000 kg [12,30] take such a volume of timber that after using the conversion factor of 740 kg·m −3 , the total mass of the transport does not exceed the legally permissible 40,000 kg. They are not guided by the actual conditions and parameters of the wood. The difference in the mass of 1 m 3 of load based on our research is an average of 293 kg·m −3 (39.5%), which results in an average GVW of 55,800 kg. The large variation in the volume, as well as the mass of 1m3 of the load in relation to the normative conversion factor (740 kg·m −3 ) contributes to real GVW increased by 17% [12,30].

Conclusions
The research conducted enabled the factors influencing load volume and load weight mass of sawlogs pine wood to be characterized. These factors include the date of delivery, origin of the timber, mainly the type of forest habitat, as well as the stand age and height of trees from which the wood was harvested.
In the analyzed period, the average mass of delivered timber loads, determined on the basis of weighing, was 30,330 kg, oscillating in the range of 21,750-38,910 kg. The relatively lighter loads were transported in the early spring months, while the heaviest were in the autumn and winter season. The decisive factor was the date when the wood was harvested and transported. The limit mass is highly correlated (0.617) with its volume, and the developed linear relationship allows to calculate the limit mass that does not differ statistically from reality.
Parameters determining the origin of the wood-habitat abundance, soil type and stand parameters (BDH, tree height) from which the wood was obtained has an impact on the mass of the load. The mass of one cubic meter of load is additionally dependent on the type of forest habitat. These are statistically significant relationships, but the strength of their relationship is not very high (from 0.0627 to 0.2987).
The volume of transported wood is at a stable level with an average of 29.34 m 3 (SD = 1.53 m 3 ), which does not depend on the factors characterizing the origin of the wood and the delivery date. This is the effect of implemented conversion factors for wood density (740 kg·m −3 ). We can state with a high degree of probability that the real mass of one cubic meter of a load is between 996 and 1070 kg·m −3 and is greater than the conversion factor of 740 kg·m −3 implemented in the legislation. This may result in an average underestimation of the load weight of 86,114 kg. This indicates the need to verify the adopted values of the kg·m −3 conversion factor, and an incorrectly selected indicator contributes to overloading the truck units. verify the adopted values of the kg⋅m conversion factor, and an incorr contributes to overloading the truck units.