Bottomland hardwoods are among the most diverse and productive forest ecosystems in the southeastern United States and are critically important for the provision of timber and non-timber ecosystem services. Red oaks, the dominant species in this group of forests, are of high ecological and economic value. Stem profile models are essential for accurately estimating the merchantable volume of oak trees, which is also closely indicative of total tree biomass and other ecosystem services given their allometric relationships. This study aims to develop and compare stem profiles among three red oak species in an 18-year old plantation forest using destructive sampling. Sixty trees randomly selected from an oak restoration plantation in the Arkansas Delta were felled for measuring the diameter-outside-bark (DOB) and diameter-inside-bark (DIB) at different stem heights. These sample composed of twenty trees from each of three species: cherry bark oak (CBO - Quercus pagoda Raf), Nuttall oak (NUT - Quercus texana Buckley), and Shumard oak (SHU - Quercus shumardii Buckl). Multiple models, including the segmented-profile model, form-class profile model, and second-and third-order polynomial models were fitted and compared. Results demonstrate that the form-class profile model was the best fitted for CBO and NUT, whereas the third-order polynomial model was the best for SHU. CBO tends to grow taller and has a higher wood density than NUT and SHU. These findings will inform restoration and management decisions of bottomland hardwood forests, especially red oaks in the region.
Cherry Bark OakNuttall OakShumard OakTaper ModelsWood DensitySoutheastern United StatesForest InventoryIntroduction
Bottomland hardwoods are among the most diverse and productive forest ecosystems in the southeastern US. They are critically important for providing many benefits including timber, water regulation, wildlife habitats and biodiversity, natural sceneries, and carbon sequestration and storage, among others (Saikku 1996, Kellison & Young 1997, EPA 2016). This group of forests once covered nearly 12 million ha in the region; however, sixty percent of bottomland hardwood forests have been lost in the past two hundred years mainly due to agricultural expansion (EPA 2016). Extensive efforts have been undertaken to restore bottomland hardwood forests in recent decades, mainly through plantations on cropland that was originally converted from forests (Allen 1997). Among the most important tree species in this ecosystem are oaks, particularly red oaks (Allen 1997, Kellison & Young 1997, Lockhart et al. 2010), which provide food for a variety of wildlife and lumber and veneer for the production of high quality furniture, flooring, and other products. Yet, the availability of accurate growth and yield models for the forests in general and oaks in particular is limited compared to other commercial species like pines (Rauscher et al. 2000). Information on young oak stand development is particularly limited. Such limitations have hindered the development and implementation of effective forest restoration and management decisions to enhance the provisions of timber and non-timber ecosystem services from this group of forests (Matney & Parker 1992, Rupsys & Petrauskas 2010).
Stem profile studies can help improve the accuracy of growth and yield estimations and projections while providing valuable information about the characteristics of trees and forests. Tree profile equations estimated statistically using empirical data predict the diameter at various heights along the bole of a tree (Methol 2001, Ounekham 2009), improving volume and carbon estimates. Profile equations allow for predicting volumes to merchantable top diameters, a fundamental basis for forest managers to compute volumes to any top diameter as well as diameter/height at a specific height/diameter along the stem. However, there is a lack of published profile equations for oaks in the hardwood forests. Given the ecological importance of oaks in the forests and the high value of oak wood products, the development of profile equations will aid in their volume and value estimation.
Profile modeling is an active area of research in forestry given the flexibility and utility afforded by the various existing models (Westfall & Scott 2010, Tian et al. 2017). Taper equations represent the mathematical expression of stem diameter change with respect to height on the basis of species, age, and stand condition (Husch et al. 1972, Brooks et al. 2008). A number of taper functions have been developed for different tree species with various forms including simple taper functions, variable-form taper functions, and segmented polynomial taper functions (Methol 2001, Jiang et al. 2005). Simple taper functions mainly define tree profiles with a single continuous equation for the whole bole (Bruce et al. 1968, Kozak et al. 1969, Ormerod 1973, Hilt 1980, Gordon et al. 1996, Zakrzewski & MacFarlane 2006). Though simple, they were unable to precisely describe the whole bole profile especially the base of the stem (Jiang et al. 2005). By contrast, the variable-form taper functions assume that the stem form/shape continuously changes with the height (Lee et al. 2003); therefore, those equations describe the stem profile with an exponent variable to characterize the neiloid, paraboloid, and conic (Newnham 1988, Kozak 2004). The advantages of the variable-form taper functions are that they describe the stem form through a continuous function and can predict/estimate the upper-stem diameters more accurately. However, the disadvantage of these functions are their inability to be logically associated to estimating total volumes.
Segmented polynomial taper functions describe stem shape through fitting a distinct equation for each segment and then binding multiple equations together to produce and overall segmented structure estimation (Max & Burkhart 1976, Cao et al. 1980, Clark et al. 1991, Jiang et al. 2005). Max & Burkhart (1976) developed a well-known profile equation for loblolly pine (Pinus taeda) involving a segmented polynomial regression approach, which delineates the bole of a tree into three sections based on the generalized relative geometric form. This approach is based on the assumption that in mathematical evaluation a segmented model will best approximate the stem form. A segmented polynomial profile model can be created by grafting three sub-models at two joint points referred to as inflection points (Sharma & Burkhart 2003). Similar to Max & Burkhart (1976) model, Cao (2009) segmented profile model also consists of three sub-models that are joined at two points; each sub-model in Cao (2009) has a form of the modified Goulding & Murray (1976) model, which can be inverted to predict height and volume. Clark et al. (1991) proposed a form-class segmented profile model by incorporating the form quotient of the Girard form class to produce a more robust model. Their model assesses stem structure partially based on the relative measurements involved with the quotient. They compared their model with Max & Burkhart (1976) segmented polynomial equation and found that the form-class segmented model provides more accurate volume estimations. The segmented-profile models including form-class segmented profile models were commonly employed in fitting the stem of hardwood species. For example, Jiang et al. (2005) employed the segmented-profile model by Max & Burkhart (1976) and Clark et al. (1991) in fitting the stem of yellow-poplar; Tian et al. (2017) fitted the invasive tallow stem using both Max & Burkhart (1976) and Cao (2009) as well as Clark et al. (1991). In addition, Alkan & Ozçelik (2020) tested nine taper equations for stem profile and volume of Taurus fir in Taurus Mountains of Turkey and found that among segmented taper models, the Fang et al. (2000) model was the most accurate taper equation for estimating diameter, predicting merchantable, and total volume. Moreover, this model does not require numerical integration like the variable-form equations for volume estimations. Parker (1997) suggested the use of a non-segmented, third degree polynomial for small data sets with limited diameter and height ranges, and developed a computer program for calculating both non-segmented and single or multiple segmented tree profile equations. His equations were based on non-destructive measurements of upper stem diameters of standing trees taken with an optical dendrometer. Overall, segmented polynomial functions composed of a series of sub-models representing various sections of the stem perform better than simple or variable form taper functions and are widely used (Ounekham 2009).
The main objective of this study is to develop and compare the stem profile models for three oak species planted in the Arkansas Delta, located in the lower Mississippi River alluvial valley (LMRAV). Hardwood forests in the LMRAV are a typical example of bottomland hardwood forests in the southeastern US in terms of their ecological features, conversion to croplands, and restoration efforts (MacDonald et al. 1979, Allen 1997). The three oak species are cherry bark oak (CBO - Quercus pagoda Raf), Nuttall oak (NUT - Q. texana Buckley), and Shumard oak (SHU - Q. shumardii Buckl). Different profile models for each species will be estimated and compared. Additionally, wood density will also be measured and compared among the three oak species, indicative of the physical property and carbon content of the oak wood. As natural and primary hardwood forests become increasingly scarce and are protected for conservation purposes, a dominant amount of hardwood timber has come from second growth and plantations. Plantations will become increasingly important in timber supply in the near future. Additionally, stem profile equations are fundamental to quantifying tree biomass which is closely related to the provisions of non-timber ecosystem services. Hence, the results of this study will inform decision making in oak plantation design and management to restore the bottomland hardwood forests for both timber and non-timber objectives while enriching the existing literature in oak stem profile, volume, and wood density.
MethodologyStudy site
The data for this study were collected from an oak restoration plantation established at the University of Arkansas Division of Agriculture Pine Tree Research Station (AgPTRS) in the Arkansas Delta (35° 07′ 48″ N, 90° 55′ 45″ W). The 12.15 ha hardwood plantation was on a prior flood prone and low productive bean field that was enrolled in the Conservation Reserve Program. The oak trees involved in this study were planted in February 2004 after the site was ripped to a depth of 0.91 m in November 2003. Seedlings were hand planted at a 3.05 × 3.05 m spacing in 12 plots with 4 replications. Each of the three red oak species (cherry bark oak, CBO - Quercus pagoda Raf; Nuttall oak, NUT - Quercus texana Buckley; and Shumard oak, SHU - Quercus shumardii Buckl) were randomly assigned to 68.58 × 68.58 m plots. (Fig. 1). No further management was applied to the plots in the 18-year interval between planting and data collection for this study.
Data collection
Destructive sampling procedures are a standard practice in tree profile data collection because stem diameters including diameter-outside-bark (DOB) and diameter-inside-bark (DIB) can easily be measured on felled trees. In addition, stem profile data measured along a felled tree provides reliable data, representing one of the most accurate methods for volume quantification. We employed destructive sampling for accurate and reliable profile data measurements of diameter and height value pairs along the entire length of the tree stem.
A total of 60 oak trees were randomly selected (20 for each species of CBO, NUT, and SHU) equally from all each of the four replications in the study (5 trees from each replication). Prior to felling, total height of the standing sample trees was measured with a clinometer; meanwhile, diameter at breast height (DBH) were also measured for all the 60 selected trees. They were then felled at approximately 0.15 m from the ground, and the bole was cleared of branches to facilitate measurements. Merchantable stem length was recorded for each tree by running a tape measurement along the length of the main apical branch until a 10.16 cm DOB was reached. The stem was divided into 1.83 m sections in the field and disks with 2.5-5.0 cm in thickness were then extracted from the midpoint of each section; DOB and DIB were measured at the lower and upper ends of each section. These collected disks were marked with species, tree number as well as disk number. They were sealed in pre-weighed and marked plastic bags to preserve their green weight. In addition, both DOB and DIB at a selected height (~5.27 m) was also recorded. As a result, DOB and DIB were recorded along the bole at relative heights of 0.15 m, 1.37 m, 1.83 m, 3.66 m, 5.27 m, 5.49 m, 7.32 m etc. until the point of a 10.16 cm diameter outside bark was reached. At each recorded height along the felled stem, both DOB and DIB were measured twice in two different directions for accuracy.
The tree disks were returned to the laboratory and weighed in the pre-weighed bags to determine green weight. The disks were then oven dried at 80 °C and reweighed to determine moisture content. A wedge of wood was then separated from the oven dry disk and the weight of that wedge was recorded. The volume of the wedge was then determined by the displacement method in water at 3 °C after coating the wedge in paraffin wax. The wood density was computed using the dry weight of the wedge divided by the volume of the same wedge, denoted by grams per cubic centimeter (g cm^{-3}).
Profile models
A total of seven profile models were fitted in this study including the segmented polynomial models of Max & Burkhart (1976), Cao (2009), Clark et al. (1991) and optimized Cao (2009); and the form-class models of second- and third-order polynomial equations. To be specific, those seven profile models were: (MB) Max & Burkhart (1976); (CAO) Cao (2009); (CSS) Clark et al. (1991), optimized Cao (2009); (OCAO_F) modified MB to pass through DBH and FC at FC height, optimized Cao (2009); (OCAO_R) modified MB to pass through DBH and FC at relative height; (POLY3) 3^{rd}-order polynomial joined at breast height to 5th order polynomial; and (POLY2) 2^{nd}-order polynomial joined at breast height to 5^{th}-order polynomial. Specific mathematic expressions of these models are shown in the equations below. The notations associated with the seven taper equations are listed in Tab. 1.
The widely used MB profile function is expressed in eqn. 1 and b_{1}, b_{2}, b_{3}, b_{4}, α, and β are parameters to be estimated, α and β represent the joint points as showed in the constraint conditions function. This segmented polynomial model describes a tree stem with three sections (eqn. 1):
The profile model of CAO including the optimized OCAO_F and OCAO_R (Cao 2009) is shown in eqn. 3-5 with parameter b_{1} being modified to parameter b_{1}^{*} for OCAO_R (eqn. 3, eqn. 4, eqn. 5):
The segmented model of CSS describes the stem with four sections composed of the butt section (from stump to breast height of 1.37 m), lower stem (between 1.37 m and 5.27 m), middle stem (from 5.27 m to 40-70% of total height), and upper stem (from 40% to 70% of total height to the tip of the tree - Jiang et al. 2005, Ozçelik et al. 2011). The mathematic expression is presented in eqn. 9, where b_{1}, b_{2}, b_{3} are the parameters for the butt section; b_{4} is the parameter for the lower stem; b_{5}, b_{6} are the coefficients for height above 5.27 m; Z_{1} = 1-(h/TH) = 1 - S; Z_{2} = 1-(4.5/TH) = 1 - S_{3}, Z_{3} = (h-17.3)/(H-17.3).
where b_{1}, b_{2}, b_{3}, b_{4}, b_{5} are the parameters to be estimated.
Statistical analysis
To test for differences in DBH and total height of CBO, NUT, and SHU at a 0.05 significance level, analysis of variance (ANOVA) was employed.
The measurements of DOB, DIB, and height pairs provided the inputs for TProfile (Matney 1996) to estimate the parameters of each profile equation described above. All profile equations fitted in this study were calculated using Tprofile which can fit commonly used profile functions for a variety of hardwood species throughout the Southeast US (Matney 1996, Gordon et al. 2015). One potential problem in fitting the stem profile model is autocorrelation (VanderSchaaf & Burkhart 2007, Yang et al. 2009); to address this concern, Durbin-Watson test for both DOB and DIB data of three oak species were performed. The test statistic values were all around 1.4/1.5 (p-value ~ 0.09), indicating that autocorrelation was not a major concern in the data.
Those fitted equations were evaluated and compared using mean absolute error (MAE), root mean square error (RMSE), and the index of fit (FI). Specifically, the MAE (eqn. 13) was computed along the bole as the average absolute difference between observed and predicted diameter values. The RMSE (eqn. 14) was calculated by the square root of the average of the absolute value of the squared biases whereas FI (eqn. 15) was calculated as one minus the sum of squared errors divided by the total sum of squares.
where Y_{i} denotes the observations/measurements, Ŷ_{i} represents the predictions, bar{Y} is the mean of Y_{i}, and n is the number of observations in the dataset.
Furthermore, we plotted the residuals (computed by observed values minus predicted values) of DOB and DIB against relative height to diagnose how well each model predicted the stem profile data. Together with the computed RMSE and FI values, we determined the best fit profile equation for each oak species.
ResultsDescriptive statistics for three oak species
The descriptive statistics of DBH, total height, and wood density for all three oak species of CBO, NUT, and SHU are summarized in Tab. 2. Average DBH for CBO, NUT, and SHU were 18.44 cm, 16.51 cm, and 16.61 cm, respectively. Total height, on average, was 17.75 m, 15.56 m, and 13.58 m respectively for CBO, NUT, and SHU. The results of ANOVA tests indicated that there was no distinct difference in DBH among CBO, NUT, and SHU; however, a significant difference was found in total height among them. To be specific, CBO was significantly taller than NUT which was taller than SHU, though they were planted in a similar environment at the same time. NUT and SHU have a similar wood density, which is significantly lower than that of CBO.
Profile equation parameters and comparisons
Tprofile estimated parameters for all seven profile models for the three oak species. Tab. S1 (Supplementary material) summarizes parameter estimations of all fitted profile models for both DOB and DIB of each oak species and the number of observations for each model. According to the MAE (lower), RMSE (lower), and FI (higher - Tab. 3) along with graphical evaluation of residual plots (Fig. S1a-S1c in Supplementary material), the CSS profile model was the best fitted for both CBO and NUT; however, the third-order polynomial model (POLY3) was the best fitted for SHU, especially in terms of DIB. Specifically, the FI, MAE, and RMSE for both DOB and DIB of CBO were FI_{OB/IB }= 0.987, MAE_{OB/IB} = 0.025, and RMSE_{OB/IB} = 0.094 with the CSS model; whereas for NUT with the CSS model, the FI, MAE, and RMSE were estimated at FI_{OB} = 0.973, MAE_{OB} = 0.033, and RMSE_{OB} = 0.145 for DOB and at FI_{IB} = 0.967, MAE_{IB} = 0.033, and RMSE_{IB} = 0.150 for DIB, respectively. By contrast, the FI, MAE, and RMSE for SHU with the POLY3 model were estimated at FI_{OB} = 0.976, MAE_{OB} = 0.041, and RMSE_{OB} = 0.140 for DOB, and FI_{IB} = 0.968, MAE_{IB} = 0.064, and RMSE_{IB} = 0.152 for DIB. To be noted, the MAE of CSS for both DOB and DIB profile of CBO species was slightly higher than the MB model (MAE_{MB_OB} = 0.023, MAE_{MB_IB} = 0.015) while it was slightly higher than the CAO model (MAE_{CAO_OB} = 0.028) for DOB profile of NUT. In addition, there was no difference in FI and RMSE for SHU between the POLY3 and CSS models for DOB. Yet, there was a slight difference in FI and RMSE for the DIB profile of SHU between the two models with FI_{IB} = 0.966 and RMSE_{IB} = 0.155 for the CSS model. Tab. S1 (Supplementary material) summarizes the parameters estimated for all seven models for both DOB and DIB profiles together with the total number of observations utilized in the model fitting; Tab. 3summarized all the statistics of model evaluation.
In addition, model performance was further evaluated using average residual plots for all seven models of three oak species (see Fig. S1-S3 in Supplementary material). In these plots, the y-axis represents the average residues and x-axis represents the relative height, h/TH. By comparing seven average residual curves for each species, we found that the best fitted model for both CBO and NUT was the CSS model while POLY3 is the best for SHU.
Finally, graphical profile comparisons were also made to visualize the difference between the average observed and predicted stem forms with 0.1 relative height interval (Fig. 2, Fig. 3, Fig. 4). We constructed the stem profile curves for both DOB and DIB for each species using its best fitted profile equations and the measurement data. Those curves also demonstrated that there was a good fit between the observed and predicted stem diameter values along the stem, which in turn confirmed that those profile equations were able to accurately predict the stem profile.
Discussion
As explained by Avery & Burkhart (2002), numerous factors could affect the performance of profile equations, including species, region, number of sample trees, data collection method among others. Those factors were all carefully considered and taken into account in the development of oak profile equations in this study. As reported by Saarinen et al. (2019), the sample size of destructive stem analysis could affect parametrizing the stem profile equations. A total sample of 60 oak trees, with 20 trees for each species, is considered to be an appropriate sample size for this study given that all oak trees were felled on the same even-aged restoration plantation site. However, since the data comes from a single site, the effect of geographical variability is not addressed. Variables including total height, DBH, various heights along the stem, and DOB and DIB at different heights were measured using the destructive sampling method. This procedure provided reliable data for fitting the oak profile models. Seven different profile models were estimated and compared for each oak species. Multiple indicators including RMSE, index of fit, and residual plots were used to diagnose the fitted models, representing a statistically sound approach for model evaluation.
Our results show that the best fitted model was CSS for both CBO and NUT, whereas it was POLY3 for SHU. Although no difference was found between CSS and POLY3 for fitting the DOB profile of SHU, for the DIB profile of SHU the FI of the CSS model was slightly lower than that of POLY3, while the RMSE of the CSS model was slightly higher than that of POLY3. As reported by Muhairwe et al. (1994), species could be an influencing factor for stem form variation. The MB model did not exhibit better performance than CSS. This may be partly because the MB model is constrained by a lack of representation of upper stem diameter, which is evidenced by an insufficient number of observations in the upper crown area. Moreover, the most important advantage of Clark et al. (1991) is using extra upper diameter measurement at 5.27 m which also contributed to the model performance. Leites & Robinson (2004) reported that crown dimension is another potential factor that contributes to the differences of stem profile. Therefore, future work can incorporate crown length and crown ratio into the profile models. In addition, utilization of modern regression techniques such as the mixed-effects modeling approach can also improve the profile model performance. For example, Gómez-García et al. (2013) developed a mixed-effect variable exponent taper equation for birch trees in northwestern Spain which produced the best fitting statistics. The profile equations that express mathematical relationships between tree height and diameter are essential for quantifying stem volume as well as total tree biomass and carbon content. Based on the profile equations, stem volume can then be calculated (Von Gadow et al. 2001). With properly estimated allometric relationships between the volumes of stem and other parts of the tree, total tree biomass and carbon content can also be derived. Besides, the similarities and differences in DBH and total heights among the three species can aid in selecting a proper species mix to achieve different management objectives. For example, Cherry bark oak is a better species choice than Nuttall or Shumard oak if the management objective is for timber, biomass, or carbon because cherry bark oak grows taller and has a higher wood density than Nuttall or Shumard oak, although they have similar DBH at the same age (Stein et al. 2003). However, because of their different heights, the mix of the three species can create more diverse stand structure and enhance biodiversity. Thus, our results provide a base-station for guiding the restoration and management of bottomland hardwood plantation forests to achieve both timber and non-timber objectives.
Conclusions
We estimated stem profile models for three red oak species (cherry bark, Nuttall, and Shumard oaks) with data collected from 60 felled trees randomly selected from a plantation in the Arkansas Delta in the southern United States. We found the best fitted model was Clark et al. (1991) profile equation for cheery bark and Nuttall oaks and the third-order polynomial model for Shumard oak. Additionally, there is no statistically significant difference in DBH among the three oak species while cherry bark oak is the tallest and has the highest wood density of the three species. The fitted profile equations in this study can be implemented into inventory software such as TCruise (Matney 1996) to improve the stem volume estimates of young oak trees in Arkansas Mississippi Delta. Given that an increasing amount of hardwood timber is supplied from second growth and plantations, our findings will guide decision making in red oaks restoration and management and their wood product marketing. Moreover, our estimated stem profile models and wood density can help quantify total tree biomass and associated carbon content. This, however, will entail the development of allometric relationships between the volumes or biomass quantities of stem and other parts of trees such as branches, leaves, and roots. Future studies on this site will determine if the stem profiles change over time as this plantation matures.
A few caveats of this study should be noted. First, while the results of fitted profile models provide a guideline for estimating the stem volume of oak species, we do not include comparisons with oak volume models given that the focus of this study was to compare stem profiles among different oak species grown on the same site instead of volume prediction. Second, the data were collected from an even-aged plantation and could not symbolize to other natural/planted mixed species and uneven-aged bottomland hardwood forests in this region, further study is needed to assess the viability of widespread application.
Acknowledgments
We are thankful for the support provided by the College of Forestry, Agriculture & Natural Resources at the University of Arkansas at Monticello and the Arkansas Center for Forest Business in completing this study. We also thank Ouname Mhotsha for helping data collection, and Ruxin Tao for his suggestion in making the study site map. Additionally, we are grateful for all field supports including technical and equipment support provided by the Division of Agriculture Pine Tree Research Station.
ReferencesAlkan O, Ozçelik RStem taper equations for diameter and volume predictions of Abies cilicica Carr. in the Taurus Mountains, Turkey. Journal of Mountain Science 17 (12): 3054-3069.2020Allen JAReforestation of bottomland hardwoods and the issue of woody species diversity. Restoration Ecology 5 (2): 125-134.1997Avery TE, Burkhart HEForest Measurements. McGraw-Hill series in forest resources. McGraw-Hill, Boston, MS, USA, pp. 456.2002Brooks JR, Jiang L, Ozcelik JCompatible stem volume and taper equations for Brutian pine, Cedar of Lebanon, and Cilicica fir in Turkey. Forest Ecology and Management 256: 147-151.2008Bruce D, Curtis RO, Vancoevering CDevelopment of a system of taper and volume tables for red alder. Forest Science 14 (3): 339-350.1968Cao QV, Burkhart HE, Max TAEvaluation of two methods for cubic-volume prediction of loblolly pine to any merchantable limit. Forest Science 26: 71-80.1980Cao QVCalibrating a segmented taper equation with two diameter measurements. South Journal of Application Forest 33 (2): 58-61.2009Clark A, Souter RA, Schlaegel BEStem profile equations for southern tree species. Res. Pap. SE-282, USDA Forest Service, Southeastern Forest Experiment Station, Asheville, NC, USA, pp. 113.1991EPABottomland hardwoods. Environmental Protection Agency, Washington, DC, USA, pp. 10.2016Fang Z, Borders BE, Bailey RLCompatible volume taper models for loblolly and slash pine based on system with segmented-stem form factors. Forest Science 46: 1-12.2000Gordon AD, Lundgren C, Hay EDevelopment of a composite taper equation to predict over-and under-bark diameter and volume of Eucalyptus saligna in New Zealand. New Zealand Journal of Forest Science 25 (3): 318-327. -1996Gordon HA, Connor KF, Haywood JDProceedings of the 17th biennial Southern Silvicultural Research Conference. e-Gen. Tech. Rep. SRS-203, USDA Forest Service, Southern Research Station, Asheville, NC, USA, pp. 551.2015Goulding CJ, Murray JCPolynomial taper equations that are compatible with tree volume equations. New Zealand Journal of Forest Science 5 (3): 313-322.1976Gómez-García E, Crecente-Campo F, Diéguez-Aranda USelection of mixed-effects parameters in a variable-exponent taper equation for birch trees in northwestern Spain. Annals of Forest Science 70: 707-715.2013Hilt DETaper-based system for estimating stem volume of upland oaks. Res. Pap. NE-458, USDA Forest Service, Beltsville, MD, USA, pp. 12.1980Husch B, Miller CI, Beers TWForest mensuration (2nd edn). Ronald Press, New York, USA, pp. 410.1972Jiang L, Brooks JR, Wang JCompatible taper and volume equations for yellow-poplar in West Virginia. Forest Ecology Management 213: 399-409.2005Kellison RC, Young MJThe bottomland hardwood forest of the southern United States. Forest Ecology and Management 90 (2-3): 101-115.1997Kozak A, Munro DO, Smith JHGTaper functions and their application in forest inventory. Forest Chronicle 45: 278-283.1969Kozak AMy last words on taper equations. Forestry Chronicle 80 (4): 507-515.2004Lee WK, Seo JH, Son YMModeling stem profiles for Pinus densiflora in Korea. Forest Ecology and Management 172 (1): 69-77.2003Leites LP, Robinson APImproving taper equations of loblolly pine with crown dimensions in a mixed-effects modeling framework. Forest Science 50 (2): 204-212.2004Lockhart BR, Guldin JM, Foti TTree species composition and structure in an old bottomland hardwood forest in South-Central Arkansas. Castanea 75 (3): 315-329.2010MacDonald PO, Frayer WE, Clauser JKDocumentation, chronology, and future projections of bottomland hardwood habitat loss in the Lower Mississippi alluvial Plain (vols. 1 and 2). US Fish and Wildlife Service, Vicksburg, MS, USA, pp. 124.1979Matney TG, Parker RCProfile equations for several hardwood species to variable top diameter limits. Southern Journal of Applied Forestry 16 (2): 75-78.1992Matney TGTProfile. World Wide Heuristic Solutions, Starkville, MS, USA, pp. 24.1996Max TA, Burkhart HESegmented polynomial regression applied to taper equations. Forest Science 22 (3): 283.1976Methol RJComparisons of approaches to modelling tree taper, stand structure and stand dynamics in forest plantations. Ph.D. thesis. University of Canterbury, Christchurch, New Zealand, pp. 298.2001Muhairwe CK, Lemay VM, Kozak AEffects of adding tree, stand, and site variables to Kozak’s variable-exponent taper equation. Canadian Journal of Forest Research 24 (2): 252-259.1994Newnham RA variable-form taper function. Report PI-X-83, Petawawa National Forestry Institute, Canadian Forest Service, Canada, pp. 83.1988Ormerod DWA simple bole model. Forestry Chronicle 49: 136-138.1973Ounekham KDeveloping volume and taper equations for Styrax tonkinensis in Laos. Thesis, University of Canterbury, Christchurch, New Zealand, pp. 77.2009Ozçelik R, Brooks JR, Jiang LModeling stem profile of Lebanon cedar, Brutian pine, and Cilicica fir in Southern Turkey using nonlinear mixed-effects models. European Journal of Forest Research 130: 613-621.2011Parker RCNondestructive sampling applications of the Tele-Relaskop in forest inventory. Southern Journal of Applied Forestry 21 (2): 75-83.1997Rauscher HM, Young MJ, Webb CD, Robison DJTesting the accuracy of growth and yield models for southern hardwood forests. Southern Journal of Applied Forestry 24 (3): 176-185.2000Rupsys P, Petrauskas EDevelopment of q-exponential models for tree height, volume and stem profile. Physical Science International Journal 5 (15): 2369-2378.2010Saarinen N, Kankare V, Pyörälä J, Yrttimaa T, Liang X, Wulder MA, Holopainen M, Hyyppä J, Vastaranta MAssessing the effects of sample size on parametrizing a taper curve equation and the resultant stem-volume estimates. Forests 10 (10): 848.2019Saikku MDown by the riverside: the disappearing bottomland hardwood forest of southeastern North America. Environment and History 2 (1): 77-95.1996SAS Institute Inc.SAS/ETS User’s Guide, Version 9.0. SAS Institute Inc., Cary, NC, USA, pp. 125.2002Sharma M, Burkhart HESelecting a level of conditioning for the segmented polynomial taper equation. Forest Science 49: 324-330.2003Stein J, Binion D, Acciavatti RField guide to native oak species in Eastern North America. USDA Forest Service, Forest Health Technology Enterprise Team, Morgantown, WV, USA, pp. 54.2003Tian N, Fan Z, Matney FG, Schultz EBGrowth and stem profiles of invasive Triadica sebifera in the Mississippi coast of the United States. Forest Science 63 (6): 569-576.2017VanderSchaaf CL, Burkhart HEComparison of methods to estimate Reineke’s maximum size-density relationship species. Forest Science 53: 435-442.2007Von Gadow K, Real P, Alvarez-Gonzalez JGModelizacion del crecimiento y la evolucion de los bosques [Modeling of the growth and evolution of forests]. IUFRO World Series, Vienna, Austria, pp. 242. [in Spanish]2001Westfall JA, Scott CTTaper models for commercial tree species in the northeastern United States. Forest Science 56 (6): 515-528.2010Yang Y, Huang S, Trincado G, Meng SXNonlinear mixed effects modeling of variable-exponent taper equations for lodgepole pine in Alberta, Canada. European Journal of Forest Research 128: 415-429.2009Zakrzewski WT, MacFarlane DWRegional stem profile model for crossborder comparisons of harvested red pine (Pinus resinosa Ait.) in Ontario and Michigan. Forest Science 52: 468-475.2006
Study site and the layout of the measurement plots with all three oak species.
DOB and DIB profiles for cherry bark oak (CBO) using the best fitted model of CSS (Clark et al. 1991).
DOB and DIB profiles for Nuttall oak (NUT) using the best fitted model of CSS (Clark et al. 1991).
DOB and DIB profiles for Shumard oak (SHU) using the best fitted model of POLY3 (3^{rd}-order polynomial joined at breast height to 5^{th}-order polynomial).
Notations associated with stem profile models.
Notations
Definition
TH
The total height of a tree
h_{D}
Breast height (1.37 m)
D
Diameter at breast height (DBH)
d
Diameter at specific height of h, including diameter-outside-bark (DOB) and diameter-inside-bark (DIB)
h
Height above the ground to the measurement point
d^{*}
Calibrated diameter
F
diameter at 5.27 m, including DOB and DIB
Descriptive statistics of DBH and total tree height for CBO, SHU, and NUT. Different letters in each column indicate statistically different values (p<0.05). (CBO): cherry bark oak (Quercus pagoda Raf); (NUT): Nuttall oak (Quercus texana Buckley); (SHU): Shumard oak (Quercus shumardii Buckl); (n): sample size; (SD) standard deviation.
Oak species (n)
Statistics
DBH(cm)
Total Height(m)
Wood density(g cm^{-3})
CBO(20)
Average
18.44^{a}
17.75^{a}
0.73^{a}
Min
13.72
15.24
0.65
Max
23.11
19.81
0.86
SD
2.67
1.09
0.05
NUT(20)
Average
16.51^{a}
15.56^{b}
0.68^{b}
Min
11.68
12.19
0.49
Max
20.60
18.29
0.83
SD
2.59
1.48
0.07
SHU(20)
Average
16.61^{a}
13.58^{c}
0.68^{b}
Min
11.05
11.28
0.59
Max
26.67
16.46
0.77
SD
4.39
1.34
0.04
Model evaluation statistics of MAE, EMSE, and FI of seven profile models for oak species of CBO, NUT, and SHU. (CBO): cherry bark oak (Quercus pagoda Raf); (NUT): Nuttall oak (Quercus texana Buckley); (SHU): Shumard oak (Quercus shumardii Buckl); (MB): Max & Burkhart (1976); (CAO); Cao (2009); (CSS) Clark et al. (1991); (OCAO_F): optimized Cao (2009) modified MB to pass through DBH and FC at FC height; (OCAO_R): optimized Cao (2009) modified MB to pass through DBH and FC at relative height; (POLY3): 3^{rd}-order polynomial joined at breast height to 5^{th}-order polynomial; (POLY2): 2^{nd}-order polynomial joined at breast height to 5^{th}-order polynomial; (DIB): diameter-inside-bark; (DOB): diameter-inside-bark.
Model
Oak species Statistics
CBO
NUT
SHU
DOB
DIB
DOB
DIB
DOB
DIB
MB
RMSE
0.117
0.109
0.234
0.229
0.183
0.188
MAE
0.023
0.015
0.033
0.036
0.089
0.089
Index of Fit
0.979
0.979
0.930
0.925
0.958
0.951
CAO
RMSE
0.135
0.122
0.224
0.244
0.196
0.028
MAE
0.071
0.056
0.028
0.117
0.140
0.254
Index of Fit
0.973
0.975
0.934
0.914
0.953
0.891
CSS
RMSE
0.094
0.094
0.145
0.150
0.140
0.155
MAE
0.025
0.025
0.033
0.033
0.089
0.091
Index of Fit
0.987
0.987
0.973
0.967
0.976
0.966
OCAO_F
RMSE
0.109
0.104
0.201
0.236
0.157
0.221
MAE
0.051
0.043
0.074
0.178
0.079
0.157
Index of Fit
0.982
0.981
0.947
0.920
0.970
0.932
OCAO_R
RMSE
0.099
0.086
0.173
0.191
0.140
0.208
MAE
0.036
0.028
0.048
0.071
0.081
0.160
Index of Fit
0.985
0.986
0.958
0.945
0.974
0.934
POLY3
RMSE
0.102
0.104
0.191
0.229
0.140
0.152
MAE
0.041
0.030
0.036
0.074
0.041
0.064
Index of Fit
0.984
0.984
0.953
0.940
0.976
0.968
POLY2
RMSE
0.104
0.109
0.188
0.206
0.140
0.152
MAE
0.043
0.036
0.036
0.079
0.046
0.064
Index of Fit
0.984
0.982
0.953
0.939
0.976
0.967
Tab. S1 - Parameter estimates for all seven profile models for oak species of CBO, NUT, and SHU.
Fig. S1 - Residue plots of DOB and DIB stem profile of all seven models for species of CBO.
Fig. S2 - Residue plots of DOB and DIB stem profile of all seven models for species of NUT.
Fig. S3 - Residue plots of DOB and DIB stem profile of all seven models for species of SHU.