Effects of Harvesting Systems and Bole Moisture Loss on Weight Scaling of Douglas-fir Sawlogs (pseudotsuga Menziesii Var. Glauca Franco)

Characterizing the moisture loss from felled trees is essential for determining weight-to-volume (W-V) relationships in softwood sawlogs. Several factors affect moisture loss, but research to quantify the effects of bole size and harvest method is limited. This study was designed to test whether bole size, harvest method, environmental factors, and the associated changes in stem moisture content of felled Douglas-fir (Pseudotsuga menziesii var. glauca Franco) affected the weight-to-volume relationship of sawlogs. Thirty trees in three size classes (12.7–25.4 cm, 25.5–38.1 cm, 38.2–50.8 cm) were felled and treated with one of two harvesting processing methods. Moisture content was sampled every two days for four weeks. Results showed 6% greater moisture loss in the crowns of stems that retained limbs after felling compared to stems with limbs removed after harvesting. Additionally, moisture loss rate increased as stem size decreased. The smallest size class lost 58% moisture content compared to 34% in the largest size class throughout the study duration. These stem moisture content changes showed a 17% reduction in average sawlog weight within the largest size class, shifting current W-V relationships from 2.33 tons m −3 to 1.94 tons m −3 during the third seasonal quarter for northern Idaho Douglas-fir and potentially altering relationships year-round.


Introduction
Since the 1960s, weight scaling of sawlogs for volume based transactions has developed into a widely applied method used by forest industry across North America [1,2], with interior British Columbia uses weight scaling to process 97% of its harvested timber [3].Developments in timber harvesting techniques and lumber production processes, such as improved supply chain procedures and mechanized machinery, have increased harvesting efficiency and retrievable products from each forest management activity [4].In order to meet the growing demand for increased forest products while reducing scaling costs for mills, a systematic random sample of trucks arriving at the mill is often used to determine truck weight-to-volume (W-V) relationships.The process of weight scaling, however, relies heavily upon regional knowledge of wood characteristics and seasonal trends to accurately estimate harvested volumes [5].
Water content is among the many variables associated with W-V relationships affecting sawlog inventory [6].Stem moisture content in weight scaling systems has been found to be significant in southern hardwoods and not significant in some mid-western softwood sawlogs studies [2,7].Sap flux density and thus log moisture vary spatially within stems, commonly peaking towards the outer edges in the radial direction and declining toward the heartwood, leading to different moisture contents between heartwood and sapwood [8,9].For softwoods, the moisture content of heartwood is typically considerably lower than for sapwood [10][11][12].Additional studies have identified seasonal patterns associated with stem moisture content of pulpwood [13].Previous studies have explored similar methodologies for measuring bole moisture loss of softwoods, with each producing common trends of moisture loss related to season and the presence of tree crowns [14][15][16][17].Subsequent results have shown faster moisture loss immediately following felling, with trees that have intact crowns losing more moisture over the study duration.However, these studies are characterized by smaller sample sizes (8-20 stems) and narrow sample diameter ranges (6-18 cm).
Although decreasingly common in the northern Rocky Mountain region, operational practices that have been the focus of past research, such as sour-felling and floatation of logs for transportation purposes provide insight into changes in wood moisture content for weight scaling [18].For example, sour felling works to encourage moisture loss in harvested stems in an effort to increase the amount of volume that can be transported on each truckload [19,20].Reducing transportations costs to lower biomass moisture content can also provide increased energy yields in forest products burned for energy [20].This process also increases the buoyancy of logs for transportation through floatation [18].However, while increased moisture loss benefits harvest operations based on volume based transactions, it negatively impacts profit margins of operations using weight based transactions by removing a source of profit in the bole moisture content and subsequently total stem weight.Finally, fuel moisture content studies have modeled stem moisture changes and rates based on stem size; limitations of such approaches have included the coarse temporal-resolution of samples (weekly) and that key information such as time since felling (i.e., proportion of dead and green stems) is often not included in the regression models [21,22].Douglas-fir (Pseudotsuga menziesii var.glauca Franco) is one of the predominant species of the forests in the interior Pacific Northwest [23].The water stored in standing timber is progressively reduced by as much as 50% in sapwood during the summer season, relative to the reduced winter draw down or reduction of sapwood moisture content [24,25].The relationship between reduced moisture content in the sapwood of felled timber and the ratio of stem W-V is not well understood [26][27][28].Log buyers and land appraisers are challenged with accurately establishing the value of standing timber; but either stumpage or mill-delivered log values must incorporate changes in order to value the result of felling, processing, and delivery practices.Improved knowledge of the uncertainties in mapping moisture loss within softwood sawlogs with respect to current operational harvesting practices is needed to better understand potential changes to scaled product and ultimately profit margin.
The objective of this study was to determine the effects of bole size, harvesting method, and environmental factors on moisture content of felled Douglas-fir sawlogs and the subsequent effect on sawlog W-V relationships.We hypothesized that stem size and harvest system would affect sawlog moisture content.Specifically, that smaller stems and stem retaining limbs would lose more moisture than larger stems and stems without limbs.In keeping with the concept of seasonal bole water storage and draw-down [13,29], we also hypothesized that stemwood moisture content would be significantly affected by relative humidity, atmospheric temperature, and vapor pressure deficit (VPD).Finally, we believed that including environmental factors such as soil moisture and stem ground contact would not significantly affect changes in moisture content.The change in moisture content observed with respect to the measured variables will help determine the impact of moisture loss on W-V relationships and current weight scaling practices.

Site Description
The research site was selected to represent high potential evapotranspiration using a southern aspect in late summer in order to provide longer radiation exposure to the felled samples.The mature, second growth, closed canopy Douglas-fir (Pseudotsuga menziesii) stand was located in the West Hatter Creek unit (WHC) of the University of Idaho Experimental Forest in north-central Idaho (46°50′05.15″N/116°50′25.19″W).The WHC unit is located 15 km northeast of Moscow, Idaho, on the north slope of the Palouse Range (Figure 1).The Palouse Range is a heterogeneous mosaic of second-growth mixed conifer forest stands distributed across a topographically complex landscape and, consequently, a diverse range of habitat types [30].The stand was predominantly Douglas-fir, but scattered Grand fir (Abies Grandis (Dougl.)Lindl.) was also present in the understory at less than 20 trees per hectare.The sampled stems showed no visible signs of defect or loss of structural integrity.The site reflected a shelterwood stand following the study treatments.Within the sample area, stems were randomly spaced and mixed between size classes, with canopy cover conditions similar for all sample stems.The site featured the expected seasonal and environmental conditions, including low precipitation and average daily temperatures associated with the time of year (Table 1).In addition to gathering gravimetric moisture content, we also collected relative humidity (+/− 0.5%), soil moisture (+/− 0.03 m 3 /m 3 ), soil temperature (+/− 1 °C), air temperature (+/− 1 °C), and precipitation (+/− 0.1mm) throughout the study.The environmental measurements were collected using a Decagon EM50 data logger, Decagon 5-TM soil moisture and temperature sensors, Decagon ECRN-100 high resolution rain gauge, and Hygrochron ™ Ibuttons.Ibuttons are small (1 cm diameter) battery-powered temperature sensors that can be deployed with inexpensive radiation shields to record surface air temperature [31].Data was collected continuously over the duration of the study at 15 minutes intervals.The data was then structured to represent daily means.Through combining the factors of atmospheric temperature and relative humidity both vapor pressure deficit and adjusted degree days were calculated.Vapor Pressure Deficit (VPD) was calculated by taking the difference in saturation vapor pressure and actual vapor pressure (Equation ( 1)).Degree days were determined by assessing the daily high temperature (°C) and calculating the difference from 0 °C (Equation (2) (The inclusion of VPD and degree days offered greater ability to assess variation in daily moisture change during the study.VPD was calculated using the formula expressed in [32].

Tree Measurements
Co-dominant Douglas-fir trees with complete crowns and free of noticeable signs of injury or disease were selected by diameter at breast height (DBH, 1.37 m) to represent the size range of sample selected.A total of 30 trees were chosen and separated by three size classes of 12.7 cm intervals (Table 2).The diameter classes represent commercial sawlog sizes ranging from lower end merchantability standards to the larger limitation of conventional mechanized felling equipment.Samples were felled according to two treatment types representing harvesting systems commonly applied in the northwestern United States.Ground based harvesting systems comprise a large portion of the timber annually harvested in North America.Within the ground based category, independent procedures and operations are adjusted depending on machinery availability and contractor preference.Two of the most prominent methods are cut-to-length harvesting and whole tree harvesting which vary when and where the felled trees are processed for transportation to mills.Samples identified as treatment 1 were felled and immediately processed into sawlogs, which included removing all limbs manually with a chainsaw.This method represented a cut-to-length harvesting system where trees are processed at the stump.Samples identified as treatment 2 were felled and left intact.The process of felling-only and leaving the bole with limbs intact represented a whole-tree harvesting system.Trees were sampled every other day for a period of four weeks.The cores were collected using a 5.15 mm diameter increment borer (Haglöf, Sweden).After felling, all stems were measured for height using a steel tape (Feet, 1/10th) and relative proportion in contact with the ground.Ground contact percentage was measured by dividing the length of the stem in contact with the ground by total stem length.Previous studies have assessed field wood moisture content with samples containing some portion of ground contact [33] and alternative research where samples were suspended above the ground [34].The establishment of ground contact percentage helped to identify the potential effect of soil temperature and moisture on adjusting stem moisture content.
To avoid broken cores or sampling in close proximity to prior cores, grids of each sampling location were identified and each sampling point was separated by 7.62 cm in all directions from other sample points.After the completion of drilling charts on each stem, cores were taken randomly from each sample location throughout the study.In order to better characterize potential vertical variability along the bole, trees were sampled at two locations on each stem.The lower core sample was taken at 1.37 m from the base of the tree established prior to felling, and the upper sample was taken at the mid-point of the crown.The upper stem sample locations averaged diameters of 15.7, 18.5, and 24.1 cm, respectively, for the 12.7-25.4,25.5-38.1,and 38.2-50.8cm size classes.For each sample removed, heartwood and sapwood were separated using a razor blade.During analysis, total stem moisture content was derived by accounting for the heartwood to sapwood ratio of from each sample location.Samples were immediately placed into air-tight, pre-weighed plastic bags and placed in an environmental container.The heartwood-sapwood barrier was determined visually.In the laboratory, samples were weighed within three hours of collection and an initial weight was recorded with the samples still in bags.The predetermined bag weight was then subtracted and samples were dried according the ASTM International Standards for drying wood to obtain moisture content [35].After drying, samples were reweighed and the percentage of moisture content of each sample was calculated as the difference between the original sample mass and oven-dry mass divided by the oven-dry mass, often considered to determine moisture content on a green basis (Equation ( 3)) [36].Similar processes using Equation (3), or variations of it, have been implemented in previous studies exploring mositure content of surface fuels in fire models [33,34].

Statistical Analysis
All statistical analysis was completed using R Statistical Package [37].A series of one way ANOVA tests were used to evaluate the significance of the variables days since harvest, treatment type, size class, and sample stem location.To further understand the influence of size class, treatment type, and sample stem location, Tukey's HSD test was used to evaluate interactions between these variables (α = 0.05).
From the ANOVA analysis a series of linear mixed-effects model were developed using the significant variables and also days since harvest, VPD, RH, soil moisture content, soil temperature, atmospheric temperature, and percentage of ground contact for each felled log segment to test the Ho that size class, and environmental factors have no effect on stem moisture content of felled sawlogs.Within the model analysis, sample trees were treated as random variables to reduce the magnitude of temporal effects and limit the influence of standard error across the study days.All models were tested using a similar form to: Y( mc) = β0 + β1 Days Since Harvest + β2 Size Class + β3 VPD + β4 (1|tree) + u.Akaike's Information Criterion (AIC) test was used to assess which model performed the best in predicting stem moisture content over the study duration.AIC evaluates the statistical quality of a model in comparison to others by balancing goodness of fit against model complexity [38], where individual AIC values are unit-less and non-meaning without other model values to compare it against.When comparing between model AICs, lower values indicate better fitting models [39].In the case of negative reported values, it is the highest negative value that identifies the best model.The delta AIC (Δi) statistic is commonly used to assess the statistical difference between competing models and is calculated as [40,41] (Equation ( 4)).
where AICi is the value for an individual competing model and min(AIC) is the minimum AIC value among tested models.The AIC and R 2 values from the model were derived from conditional R 2 and AIC procedures established in Nakagawa and Schielzeth, 2013 [42].

Moisture Loss
The ANOVA analysis indicated that only size class by itself exhibited significant influence on moisture loss.In testing the interactions between the variables, Tukey's HSD identified significant differences in moisture loss when considering both the interaction between size class and treatment type, and the interaction between sample location and treatment type (Table 3).Within the smallest size class (12.7-25.4cm), 58% of the stem moisture content was lost over the four weeks compared to 34% in the largest size class (38.2-50.8cm) (Figure 2).The interaction between whole-tree harvesting with size class accelerated moisture loss over the rates observed in the cut-to-length harvest system (Table 3).The effect of treatment type on stem moisture loss was also significant when looking at sample location within the stem (Figure 3 and Table 3), with upper stem locations treated with cut-to-length harvesting losing moisture significantly slower than stems treated with whole-tree harvesting.The change in stem moisture content slowed as average atmospheric relative humidity (RH) and stem moisture content drew closer later in the study (Figure 4).Although stem heartwood moisture content was lower to begin with, mean heartwood moisture content experienced less loss than the mean sapwood moisture content over the study period, losing an average of 42% compared to 62% (Table 4).

Optimal Model
Different combination of size class, treatment type, VPD, RH, days since harvest, soil moisture content, soil temperature, atmospheric temperature, and percentage of ground contact for each felled log segment, and the significant interaction terms were used within the linear mixed-effects model analysis.This rigorous analysis approach has not been previously applied to predicting moisture loss of harvested sawlogs, where each sample tree is utilized as a random variable to account for influences in standard error between sample days on the same tree.The application of a linear mixed-effects model recognized the interaction between sample days through using each sample tree as a random variable.Model evaluation using the AIC showed the strongest model (i.e., smallest integer or highest negative) to be Model 3: Y(mc) = β0 + β1 Days Since Harvest + β2 Size Class + β3 VPD + β4 (1|tree) + u (Table 5).The inclusion of either treatment type or RH by themselves did not improve the models' fit enough to overcome the added model complexity (Table 5).Model 3 used the covariates days since harvest, harvest treatment, and VPD combined to explain 82% (R 2 = 0.819) of the total variability in the bole moisture content of the sampled stems (Table 6).The identification of these variables supported the hypothesis that size class and environmental factors would significantly affect the prediction of stem moisture content, which declined with increasing time after harvest.In contrast though, when including soil temperature, soil moisture, and stem percentage in contact with the ground, there was no significant influence on stem moisture content of the harvested sawlogs.

ΔAIC Analysis: Secondary Model
The ΔAIC analysis identified that Model 5 was moderately comparable (i.e., ΔAIC < 10) and accounted for a higher amount of the model error [39][40][41]43].Model 5 explained a higher proportion of variation in the response variable based on R 2 (0.849 Model 5 vs. 0.819 Model 3).While Model 3 provided the best balance between fit and complexity of the models compared, Model 5 could be argued to have greater implications for management prescriptions as it included the interaction between size class and treatment types.This, as well as the operational significance of the current management shift towards younger harvesting rotations and increased commercial thinning where stems contain greater crown ratio, should be kept in mind.The individual effects of days since harvest, size class * treatment, and VPD are shown in Table 7 using size class 1, treatment 1, and their interaction as the base treatment in the model.While the moisture content loss was significant when comparing stem size regardless of treatment, the impact of the loss was much greater when viewed within the context of active logging operations.The results showed whole tree harvesting stems averaged lower moisture content over the duration of the study than cut-to-length.

Table 7.
The regression values of a generalized linear mixed effects model assessing moisture content of felled stems using Size Class 1, TRT 1, and their interaction as the base treatment.
*Denotes interaction term between two variables.

Discussion
In support of our hypothesis that stem size would affect the moisture content in felled Douglas-fir, we saw significant correlation between stem moisture content and stem size class.The greatest amount of moisture loss was observed in the days immediately following harvesting.In addition, the majority of the moisture lost was contained in the sapwood compared to the heartwood, which makes sense when considering that the majority of each tree's pre-treatment moisture was located primarily in the sapwood.The best model identified by the AIC test included the variables days since harvest, size class, and VPD.Although the AIC model with the smallest number should be considered the best overall, limited differences between the best two models (<10) supports using the second best model when operational data is available.In accordance with ΔAIC comparisons [39], the second best AIC model was also evaluated, which included an interaction term between size class and treatment type.Studies looking at seasonal changes in live tree stem moisture have reported a range of results, including increased sapwood moisture content during the growing season [27,44].In contrast, several studies show increased stem moisture content in winter months [7,13].While the standing, seasonal moisture content of conifers fluctuates by region and species, the moisture content of felled trees exhibits more variability.The large quantity of water lost suggests that even under alternative seasonal conditions, the effect of harvest system and time between harvest and delivery could prove significant when working in a weight scaling system.Tree age and size can affect wood characteristics and in turn water storage capacity [45].
The effect of stem size on moisture content over time from this study is consistent with conventional fuel moisture models used to characterize 1000 and 10,000 hours fuel moisture content.These rates classify the time needed for moisture to change 1/3 of the way towards equilibrium moisture content.While these rates prove efficient for fire predictions and fuel modeling, they contain general practices that reduce their effectiveness when applied to weight scaling.For example, they do not account for the time a given species has been dead or downed.In addition, the percentage of bark cover is not identified and base moisture contents for 1000 hours fuels are averaged over 7-day periods.
We found that harvesting and processing of sawlogs affects log moisture content, and that the results depend on sawlog diameter.Because moisture content directly affects sawlog weight, the relationship of W-V, and thus log value, would be affected when using weight scaling systems that do not account for these variables.As identified in the results, stems harvested using whole-tree harvesting system with limbs remaining on the stem after harvesting and prior to processing showed greater moisture loss than stems harvested using cut-to-length harvesting, where processing occurs at the stump and logs are limbed and bucked to length immediately.While the harvest method was important in our study, significant moisture loss occurs in each scenario.In practice, operators may modify practices to minimize loss, depending on which system is being deployed and what current environmental conditions are at the time of harvest.For example, in the current study, 38.2-50.8cm.sawlogs lost an average of 34% moisture content from each log over a one month period after harvest, regardless of harvest system.Assuming a truckload contains 24,500 kilograms, roughly matching state averages and legal highway limits, a 34% moisture content loss would result in 4.16 metric tons of lost water weight or 17% loss in total truckload value if paid using a weight basis system.Current industry weight-to-volume relationships for Intermountain West Douglas-fir during the third seasonal quarter are 5.51 tons/mbf or 2.33 tons m −3 .A 17% loss in weight would shift the ratio to 4.57 tons/mbf or 1.94 tons m −3 .This decrease in weight-to-volume ratio could impact the financial productivity of landowners and contractors.Converting between Scribner board foot to cubic volume was derived from Spelter, 2004 [46].While the very dry, warm conditions during the study may have increased moisture loss, moisture loss was present in all felled trees and should be considered operationally.
Altering harvesting practices to change bole moisture content and improve profit margins greatly depends on the scaling unit employed.Losses in moisture content are viewed negatively when associated with transactions based on weight.However, loss in bole moisture content has long been sought out when product transactions are based on mbf or cubic meter volume.Methods such as "sour felling", which consists of whole tree harvesting followed by transpirational drying to allow moisture loss have been applied in efforts to increase loaded truck net volume per load [16].While effective, these practices encourage increased defect in harvested product in the form of checking, insect infestation, and staining [17].The increase in volume per load profit can be counteracted in the subsequent loss of usable product; retaining the idea that sour felling is best applied when working with pulpwood or woody biomass where profit margins are limited and defect is less important [14,16,17,47].Sour felling practices are also applicable in harvesting systems where floating logs is used for transportation or storage purposes.Decreases in the bole moisture content result in increased buoyancy for periods of storage or transportation.Log storage in water, when properly executed, reduces checking, insect activity, and staining of logs [48].Ultimately, the perception of harvesting method and bole size on moisture loss in sawlogs should be considered in direct relation to the unit of volume measure.With western United States scaling practices shifting towards increased transactions through weight scaling, long standing notions of moisture loss in felled sawlogs must be adjusted accordingly in order to maximize profit margin.
Beyond environmental factors present during drying, bark cover and thickness is often assumed to influence drying rates [49,50].However, when investigating moisture content of branches, Boddy, 1983 [51] did not find drying to be significantly faster in bark-less branches.In contrast, Nicholls and Brackley, 2008 [52], identified bark as a barrier to moisture loss.This is supported by Defo and Brunette, 2006 [53], who found that the stem drying rate of aspen is proportional to the percentage of bark absent.Bark is reported to make up 10-25 percent of the over-bark volume and weight of a tree [54,55].Therefore, depending on location and season, bark can prove to be beneficial or costly to the forester or contractors [56].However, felled trees lose moisture through either diffusion or transpiration.While the presence of bark would affect the diffusion process, transpiration would be unaffected.In contrast, remaining tree foliage would be affected by transpiration.Therefore, if stem moisture loss is dominated by transpiration then the effects of bark coverage may go unnoticed.Further investigation concerning the pathway of stem moisture loss could help quantify the importance of bark on stem moisture loss.The simulated processing techniques for the current study were completed using chainsaws as opposed to machinery, leaving the majority of the stem bark intact.This alteration from expected seasonal and operational bark loss from machine processing limits the understanding of bark influence on moisture loss in this study.
While the species characteristics are relevant in mapping stem moisture loss rates, local environmental conditions drive changes to solid wood moisture content.The combination of low relative humidity, high mean daily temperatures, and limited precipitation on the southern aspect provided an extreme setting for moisture loss.The long-term effect of seasonal relative humidity in a region will cause non-living wood to fluctuate with a delayed lag time until equilibrium with the environment is reached.The rate of drying increases with the difference between relative humidity and stem moisture content at the time of harvesting.In order for foresters and contractors to understand and adapt to the variability in sawlog moisture content when working with weight scaling systems in mountainous terrain, sawlog size, harvesting method, and present atmospheric conditions must be considered.

Conclusions
The maximum moisture loss occurred in both the smallest size class and stems in the whole-tree harvesting system.The greatest moisture content variation and subsequent moisture loss also took place in the sapwood portion of the stem.With these results, we reject the null hypothesis that harvested stem moisture loss was unrelated to bole size and harvest system of felled Douglas-fir trees.Trees felled and left intact with limbs proceeded to show increased moisture loss vs. stems that were felled and immediately processed into sawlogs.Additionally, we predicted that the moisture content during late summer would not be significantly affected by stem ground contact percentage, bark coverage, or soil characteristics, and the results supported that prediction.Consistent with the concept of evapotranspiration, stem moisture content is best predicted by a combination of VPD, stem size, and harvest system.The probable reason for changes in the moisture content variation among treatment types is the increased surface area present of trees containing limbs and foliage.The effect of stem size on moisture content was also likely driven by the reduced amount of external energy required to extract moisture from a smaller object.These findings suggest that understanding the implications of seasonal harvests exposing harvested timber to strong climatic and environmental conditions can significantly alter expected results from areas utilizing weight scaling practices.

Figure 1 .
Figure 1.The study unit located Northeast of Moscow, Idaho in the West Hatter Creek Unit of the University of Idaho Experimental Forest (46°50′05.15″N/116°50′25.19″W).

Figure 2 .
Figure 2. Moisture content (MC) of stem wood, dry weight basis, over 28 days following harvesting; mean MC values for trees on a given sampling date are shown as solid shapes averaged across both harvest systems, and error bars show one standard error (SE).

Figure 3 .
Figure3.The mean moisture content of stem wood separated by treatment type and stem location in the time following harvesting; Treatment 1 represented a cut-to-length harvesting system and Treatment 2 a whole-tree system; these results show the MC in the upper portions of the stems.

Figure 4 .
Figure 4.The changes in stem wood moisture content and site relative humidity are shown over the duration of the study; mean moisture contents and relative humidity are shown in solid shapes with error bars showing one standard error.

Table 2 .
Characteristics of the Douglas-fir trees sampled for stem moisture content loss rate.

Table 3 .
Results of Tukey's HSD test to determine ad hoc differences in bole moisture content between size classes, stem sample location, and treatment type.
* Denotes interaction term between two variables.

Table 4 .
The beginning and ending stem moisture contents of all sample stems regardless of size separated by heartwood and sapwood.

Table 5 .
The results of Akaike's Information Criterion (AIC) testing to determine the significant model variables.

Table 6 .
The regression values of a generalized linear mixed effects model assessing moisture content of felled stems using Treatment 1 as the base treatment.
Model: Y