Allometric Equations for Shrub and Short-Stature Tree Aboveground Biomass within Boreal Ecosystems of Northwestern Canada

Aboveground biomass (AGB) of short-stature shrubs and trees contain a substantial part of the total carbon pool within boreal ecosystems. These ecosystems, however, are changing rapidly due to climate-mediated atmospheric changes, with overall observed decline in woody plant AGB in boreal northwestern Canada. Allometric equations provide a means to quantify woody plant AGB and are useful to understand aboveground carbon stocks as well as changes through time in unmanaged boreal ecosystems. In this paper, we provide allometric equations, regression coefficients, and error statistics to quantify total AGB of shrubs and short-stature trees. We provide species- and genus-specific as well as multispecies allometric models for shrub and tree species commonly found in northwestern boreal forest and peatland ecosystems. We found that the three-dimensional field variable (volume) provided the most accurate prediction of shrub multispecies AGB (R2 = 0.79, p < 0.001), as opposed to the commonly used one-dimensional variable (basal diameter) measured on the longest and thickest stem (R2 = 0.23, p < 0.001). Short-stature tree AGB was most accurately predicted by stem diameter measured at 0.3 m along the stem length (R2 = 0.99, p < 0.001) rather than stem length (R2 = 0.29, p < 0.001). Via the two-dimensional variable cross-sectional area, small-stature shrub AGB was combined with small-stature tree AGB within one single allometric model (R2 = 0.78, p < 0.001). The AGB models provided in this paper will improve our understanding of shrub and tree AGB within rapidly changing boreal environments.


Introduction
Ecosystems in northwestern Canada are changing rapidly due to a warming climate, drier conditions, extended growing season, and climate-mediated increases in frequency and intensity of disturbances, such as wildfire, permafrost thaw, insect and pathogen outbreaks, and anthropogenic natural resource extraction (e.g., [1][2][3]). One of the significant outcomes of climate-mediated change in these environments is the increased abundance of short-stature vegetation, such as shrubs [4,5] and low productive and juvenile trees, in particular where wildfire disturbance sets back ecosystems to an early successional stage post fire [4] or in the rapidly changing transition zones between elevated forests and adjacent peatlands due to permafrost thaw [6]. Increased shrub cover influences important ecosystem functions, such as energy balance and hydrology [6] at local to regional scales and greenhouse gas/carbon-climate cycle feedbacks at national to global scales [6,7], which could ultimately exacerbate these changes [8,9]. Shrubs and short-stature trees contain a substantial part of the total aboveground carbon of unmanaged boreal forest and peatland ecosystems, although specific numbers remain

Study Area
Shrub and tree species examined in this study were destructively sampled across the southern margin of the sporadic to discontinuous permafrost zone of the Taiga Plains and Taiga Shield ecozones ( Figure 1). The mean annual air temperature varies from −2.5 (Fort Simpson, Taiga Plains) to between −3 and −4 (Yellowknife, Taiga Shield) • C, while cumulative annual precipitation is 390 (Fort Simpson) and 360 (Yellowknife) mm, respectively [19,20].
The geology is dominated by lacustrine plains overlain by peatlands with fine-to coarse-textured lacustrine and till in the mid-boreal Taiga Plains, which transitions towards the high boreal Taiga Shield into nearly level to rolling and hilly bedrock [19,20]. In the Taiga Plains, approximately 25%-50% of the land surface area is covered by peatlands, consisting of peat plateaus underlain by permafrost, bogs, and fens [21]. While shrubs dominantly grow in the transition zones between peat Forests 2020, 11, 1207 3 of 16 plateaus and bogs or fens, short-stature trees grow within upland forests or on elevated peat plateaus. In the high-boreal Taiga Shield, peatlands cover an area of less than 5% and occupy hollows within and between bedrock exposures [20]. The dominant vegetation in the mid-boreal Taiga Plains consists of black spruce (Picea mariana) and mixed hardwood and softwood forests containing similar species as found in the Taiga Shield ecozone. Taiga Shield ecosystems are dominated by black spruce and jack pine (Pinus banksiana) with abundance of paper birch (Betula papyrifera) and trembling aspen (Populus tremuloides) [19,20]. The geology is dominated by lacustrine plains overlain by peatlands with fine-to coarsetextured lacustrine and till in the mid-boreal Taiga Plains, which transitions towards the high boreal Taiga Shield into nearly level to rolling and hilly bedrock [19,20]. In the Taiga Plains, approximately 25%-50% of the land surface area is covered by peatlands, consisting of peat plateaus underlain by permafrost, bogs, and fens [21]. While shrubs dominantly grow in the transition zones between peat plateaus and bogs or fens, short-stature trees grow within upland forests or on elevated peat plateaus. In the high-boreal Taiga Shield, peatlands cover an area of less than 5% and occupy hollows within and between bedrock exposures [20]. The dominant vegetation in the mid-boreal Taiga Plains consists of black spruce (Picea mariana) and mixed hardwood and softwood forests containing similar species as found in the Taiga Shield ecozone. Taiga Shield ecosystems are dominated by black spruce and jack pine (Pinus banksiana) with abundance of paper birch (Betula papyrifera) and trembling aspen (Populus tremuloides) [19,20].

Shrub Measurements, Destructive Sampling, and Processing
To determine AGB of small-stature shrubs and trees growing in the mid-boreal Taiga Plains and high-boreal Taiga Shield ecoregions, we derived allometric models for five common shrub genera and species (Alnus spp., Betula spp., Dasiphora fruticosa, Salix spp., and Shepherdia canadensis) and four common tree genera and species (Betula papyrifera, Picea spp., Populus balsamifera, and Populus tremuloides). Plant individuals were measured and destructively sampled within 65 peatlands and upland forest ecosystems distributed across the two ecoregions in late July/early August 2018 and 2019. Field sample locations were situated in late successional sites and in sites disturbed by wildland fire within the last 50 years ( Figure 1, Table 1) in order to represent the high variability of boreal ecosystem disturbance by wildfire and permafrost thaw in our allometric models.

Shrub Measurements, Destructive Sampling, and Processing
To determine AGB of small-stature shrubs and trees growing in the mid-boreal Taiga Plains and high-boreal Taiga Shield ecoregions, we derived allometric models for five common shrub genera and species (Alnus spp., Betula spp., Dasiphora fruticosa, Salix spp., and Shepherdia canadensis) and four common tree genera and species (Betula papyrifera, Picea spp., Populus balsamifera, and Populus tremuloides). Plant individuals were measured and destructively sampled within 65 peatlands and upland forest ecosystems distributed across the two ecoregions in late July/early August 2018 and 2019. Field sample locations were situated in late successional sites and in sites disturbed by wildland fire within the last 50 years ( Figure 1, Table 1) in order to represent the high variability of boreal ecosystem disturbance by wildfire and permafrost thaw in our allometric models. Harvested plants were located within <10 m of field transects (Figure 2a), which were set up randomly within each site. Between one to three transects were installed per site. Field transects were 25 m in length, starting in upland forests and traversing into peatlands and crossing the forest-peatland transition zone perpendicularly. This setup was chosen in order to capture the shrub abundance within the transition zones between peat plateaus and bogs/fens, where the largest and fastest ecosystem changes have been observed [6]. Shrubs were located along the transect with distance from the start of the transect. Shrubs were selected per transect within standing height ranges of 0.5 up to 3.5 m. A shrub individual was selected when it was alive and mostly free of leaf or stem disturbances. Up to five shrub individuals were sampled per transect. Harvested plants were located within <10 m of field transects (Figure 2a), which were set up randomly within each site. Between one to three transects were installed per site. Field transects were 25 m in length, starting in upland forests and traversing into peatlands and crossing the forestpeatland transition zone perpendicularly. This setup was chosen in order to capture the shrub abundance within the transition zones between peat plateaus and bogs/fens, where the largest and fastest ecosystem changes have been observed [6]. Shrubs were located along the transect with distance from the start of the transect. Shrubs were selected per transect within standing height ranges of 0.5 up to 3.5 m. A shrub individual was selected when it was alive and mostly free of leaf or stem disturbances. Up to five shrub individuals were sampled per transect. Measured 1D variables for each individual shrub sample included maximum height (herein "max height" (m)), average maximum height (m), number of individual stems, and basal diameter of each stem (cm). Heights were measured using a tape measure, and diameters were measured with a caliper. Previous studies have examined the relationship between the AGB of a single stem and the 1D independent variable measurement of the same stem [11,16]. However, boreal shrubs are commonly multi-stemmed. In order to analyze how accurately 1D field-measured variables from a single stem can predict the AGB of the entire plant, we tested the basal diameter of the stem with the widest diameter (herein "max basal diameter" (cm)) and the length of the longest stem (herein "max stem length" (m)) within allometric relationships. To test a 2D variable, we transformed each basal diameter to cross-sectional area (cm 2 ) and summed this to total cross-sectional area per shrub individual. For the 3D variable, we measured the extent of the uppermost foliage layer perpendicular to the transect (herein "width" (m)) and parallel to the transect (herein "line-intercept cover" (m)) using a tape measure. The 3D shrub volume (m 3 ) (1) was then calculated as follows: (1) Measured 1D variables for each individual shrub sample included maximum height (herein "max height" (m)), average maximum height (m), number of individual stems, and basal diameter of each stem (cm). Heights were measured using a tape measure, and diameters were measured with a caliper. Previous studies have examined the relationship between the AGB of a single stem and the 1D independent variable measurement of the same stem [11,16]. However, boreal shrubs are commonly multi-stemmed. In order to analyze how accurately 1D field-measured variables from a single stem can predict the AGB of the entire plant, we tested the basal diameter of the stem with the widest diameter (herein "max basal diameter" (cm)) and the length of the longest stem (herein "max stem length" (m)) within allometric relationships. To test a 2D variable, we transformed each basal diameter to cross-sectional area (cm 2 ) and summed this to total cross-sectional area per shrub individual. For the 3D variable, we measured the extent of the uppermost foliage layer perpendicular to the transect (herein "width" (m)) and parallel to the transect (herein "line-intercept cover" (m)) using a tape measure. The 3D shrub volume (m 3 ) (1) was then calculated as follows: Following measurements in situ, shrubs were clipped directly above the soil surface ( Figure 2c) and stored in paper bags for further processing. Dead stems were not harvested. In the laboratory, harvested samples were air dried for up to four months, separated into stem and leaf parts, and oven dried at 60 • C for 48 h [22]. Twigs and fruits were included as leaf parts. The total AGB was determined as dry weight (g) by weighing each shrub part (woody and leafy) and summing the dry weight of all parts per shrub individual.

Tree Measurements, Sampling, and Processing
Small-stature tree AGB was collected in late July 2019 along 20 random transects adjacent to the permanent sample plots set up by the Canadian Forest Service, located near Fort Liard, Northwest Territories (NWT), Canada. Live trees were chosen from the understory and open areas across different height ranges determined in intervals of 0.5 m up to ≤4.5 m. In all cases, samples were mostly free of foliage disturbance/mortality and stem blemishes. In situ 1D tree measurements included stem length (m), and stem diameter (cm) measured at 0.03, 0.15, 0.30, and 1.3 m along the stem length starting from the average ground surface surrounding the tree. Stem diameters were transformed to cross-sectional area (cm 2 ) per tree individual to provide a 2D variable analogous to that for shrubs and thus offer the potential for a joint shrub and juvenile/low productive tree allometric equation. 3D volume was not measured for trees because tree AGB can best be predicted with diameter, stem length, or both variables combined (e.g., [14]). Following measurements in situ, trees were cut as close to the ground surface as possible and packed into large paper bags to be transported back to the University of Lethbridge. In the laboratory, trees were separated into stem, branch, and leaf components after air drying of up to four months. Branches were cut off directly at the stem. Twigs and fruits were included as leaf parts, while bark was included as part of the stem. Dead branches were not included in the analysis. Oven drying and biomass derivation was completed using the methods described above for shrubs.

Derivation of Aboveground Biomass Allometric Equations
In situ structural measurements of harvested shrubs and trees were used to determine the most accurate 1D, 2D, or 3D independent variables to predict AGB. Three different forms of single variable regression analysis, which are most commonly used in biomass allometry [11,16,17,[23][24][25] were tested for each 1D, 2D, and 3D independent variable per shrub and tree genus and species and for the pooled data. These were used to determine the most descriptive regression model of AGB for genus/species, multispecies shrubs, and multispecies trees as well as in general for all trees and shrubs combined. The first allometric biomass model (2) uses linear regression of the log-transformed dependent (y) and independent in situ (x) variables (herein "linear logarithmic regression (LLR)"): where α and β are the regression coefficients. The back-transformation to an arithmetic scale was achieved using Equation (3): However, the back-transformation resulted in a skewed distribution ofŷ [23]. Baskerville [23] reported a general underestimation of 10%-20% when back-transforming the logarithmic regression estimates of AGB without correcting for skewness. To set this into context, we compared the results of LLR with the results of linear logarithmic regression with correction (herein "LLRC") (4-6): where ε represents a multiplicative correction factor of the back-transformation with MSE as the mean square error of the regression [23,25]. This correction removes the bias in Equation (3), which occurs following the back-transformation from a normal distribution of ln(ŷ) for a given ln(x). LLR results are presented only to set LLRC model results into context, and usage is not recommended without the provided correction factor.
To avoid the problem of skewness, a majority of the research on biomass allometry have used the untransformed nonlinear relationship between dependent and independent variables to predict tree [17] and shrub [11,16,17] AGB. These models use iterative nonlinear least squares regression via a power function (herein "nonlinear least squares regression (NLS)") with an additive error term ε: Forests 2020, 11, 1207 6 of 16 The NLS function for biomass prediction is available in statistical software packages (e.g., "nls" function in the "stats" package in R [26]). The default nonlinear equation that is implemented in R assumes a homogeneous variance of regression residuals [25]. However, trees and shrubs can be inherently heteroscedastic, such that the assumption of homogeneous residual variance across the range of the independent variable could lead to biased model predictions [25]. Although weights can be specified in NLS functions, these should be applied when both arithmetic and logarithmic variances do not show uniformity [23]. Our AGB data showed uniform variances on arithmetic scales for most species and on logarithmic scales for all species. To understand the performance of using a non-weighted nonlinear model that does not address potential heteroscedasticity within the data, we compared NLS (7) with the results of the logarithmic-based models LLR (2-3) and LLRC (4-6).
The modeled biomass results of these three allometric models were evaluated using root mean square error (RMSE), coefficient of determination (R 2 ), and regression residual analysis. Residual analysis was performed using visual inspection of the relationships between dependent and independent variables as well as the total percentage error (%) derived via (8): Significance of the differences between the genus/species-specific AGB model means and the multispecies model means were evaluated with the t-test for equal variances and Welch's t-test for unequal variances. Equality of variances was tested with the F-test.

Biomass Allometric Models
In total, we used 1D, 2D, and 3D input variables and three different forms of regression (LLR, LLRC, and NLS) to predict shrub AGB. Regression coefficients and standard errors were calculated for (a) total shrub AGB per genus/species; (b) all shrubs combined per ecoregion; and (c) all shrubs within both ecoregions combined. Here, (c) represents a multispecies shrub allometric equation that can be applied to all shrubs across the southern half of the study region (NWT, Canada; Figure 1). For short-stature trees, we examined 1D and 2D predictor variables (3D tree volume predictors were not measured) and the three different regression equations described above (LLR, LLRC, and NLS). Similar to shrubs, we have provided the allometric models that most accurately predict (a) total AGB for each tree species; (b) all tree species combined per ecoregion; and (c) for all trees within the two ecoregions combined (herein "multispecies trees"). In addition to individual "shrub" and "tree" models, we have provided a 2D input variable allometric model for total AGB prediction for small-stature shrubs and trees combined (herein "general shrubs and trees"). With the "general shrub and tree" model, we explored the utility of a single 1D (stem length) or 2D (cross-sectional area) variable for combined shrub and tree AGB prediction to understand whether these variables scale uniformly between shrubs and trees. Such combined approaches may have utility for rapid assessment in the field or using less invasive observation techniques (e.g., unmanned airborne vehicles or laser scanning) where vegetation species and type may be indeterminate.

Results and Discussion
The ranges and averages (± standard deviation) of the predictor variables for 205 shrubs and 106 trees are provided in Supplementary Material 1, Tables S1 and S2, respectively. Regression coefficients, correction factors, and standard errors are presented in Supplementary Material 2, Tables S3-S5, and can be input into Equations (2)- (7). For all allometric models, regression coefficients were positive, indicating increasing biomass with increasing predictor variable for total AGB, as was found by Lambert et al. [14].
The applicability of multispecies models to predict total AGB for single genus and species has been shown for shrubs by He et al. [16]. Similarly, the differences between our modeled means of total AGB via the multispecies equations and the genus/species-specific modeled means were not Forests 2020, 11, 1207 7 of 16 significant (p > 0.05). Furthermore, the difference between the modeled means of shrub and tree biomass in previously burned sites and the unburned sites was not significant (p > 0.05). Species-specific coefficients for input into the allometric Equations (2)-(7) are provided in Supplementary Material 2, Table S3, using volume as the predictor and Supplementary Material 2, Table S4, using cross-sectional area as the predictor.

Comparison of 1D, 2D, and 3D Variables for Shrub Total AGB Prediction
For the genus/species-specific 1D-based models using max stem length or max basal diameter as input into each of the three allometric models (LLR, LLRC, and NLS), Betula spp., Dasiphora fruticosa, and Salix spp. had lower R 2 (0.005 ≤ R 2 ≤ 0.325) compared to Shepherdia canadensis (0.433 ≤ R 2 ≤ 0.809, Table S6). In addition, Dasiphora fruticosa was the only species where stem length was not significantly related (p > 0.05) to the dependent variable of measured total AGB. For the multispecies shrub models (pooled for all shrub genera and species), the use of 1D predictor variables yielded the lowest model fits (RMSE) ranging from 262 (NLS) to 318 (LLRC) g for max stem length and 252 (NLS) to 388 (LRC) g for max basal diameter ( Table 2). These results are in contrast with previous allometric models for boreal [11,16,17] or subtropical [27] shrubs, where the 1D variable basal diameter of the longest stem had provided the most accurate prediction of total AGB. However, 1D field variables, although related to the dependent variable (p < 0.001) (with the exception of max stem length of Dasiphora fruticosa), did not explain total AGB variability when considering all stems of the entire plant (0.228 ≤R 2 ≤ 0.335, Table 2, Figure 3a,b). This was true for each genus/species model as well as the multispecies equation, with the exception of Sheperdia canadensis. For this species, max basal diameter was a similarly good predictor variable (R 2 = 0.809) to cross-sectional area (R 2 = 0.738) and volume (R 2 = 0.765, Table S6). The performance of total AGB models increased for all other genera and species as well as for all genera and species combined using the 3D predictor variable of volume, with R 2 ranging between 0.684 (Betula spp. and Dasiphora fruticosa) and 0.882 (Alnus spp.). The RMSEs for the multispecies models ranged from 141 (NLS) to 144 (LLR and LLRC) g with R 2 of~0.790 using any of the three models (LLR, LLRC, and NLS; Table 2). Figure 3a-d shows the relationship between the three models for 1D, 2D, and 3D variables for the multispecies shrub AGB.   Of the three model forms tested, the 3D volume predictor produced the most consistent shrub AGB results (Figure 3, Table 3), suggesting the choice of model form is less critical when using a 3D predictor. As expected, the exponent was greater for 1D models and close to unity for the 3D models. This suggests that a simple linear model would result in similar model fits compared to the LLRC or NLS models when using a 3D predictor. However, after testing this (results not shown), the  Of the three model forms tested, the 3D volume predictor produced the most consistent shrub AGB results (Figure 3, Table 3), suggesting the choice of model form is less critical when using a 3D predictor. As expected, the exponent α was greater for 1D models and close to unity for the 3D models. This suggests that a simple linear model would result in similar model fits compared to the LLRC or NLS models when using a 3D predictor. However, after testing this (results not shown), the model results achieved slightly less goodness of fit (R 2 = 0.769, RMSE 148.6 g) compared to LLRC or NLS (Table 2). Decreasing exponents (holding all else equal) can be explained by the nature of allometric scaling between 1D, 2D, and 3D measurements of a plant to its mass (a 3D attribute) via a power function. Assuming no change in the multiplier (β) (which might be considered analogous to a density attribute), scaling from a 1D measurement to a 3D property requires a higher exponent (α) compared to scaling from a 2D or 3D measurement (Table 3). Therefore, predictions that are extrapolated from lower to higher dimensions contain more inherent model-based uncertainty than predictions requiring no dimensional extrapolation. However, field volume observations consisted of three single measurements and therefore might contain a high overall measurement uncertainty compared to a single 1D measurement. The exact quantity of model vs. field measurement error propagation is unknown, but the net outcome of the tests performed shows that 3D volume produced the highest AGB model accuracies, followed by 2D and then 1D models. Better model performance and linearity via the 3D predictor variable of volume may also be explained by the structural variability of multi-stemmed shrubs. For example, shrub stems can grow comparably long while being simultaneously thinner rather than being shorter but thicker, so the total dry weight of the shrub with the longer stems may be lower than the dry weight of a shrub that has shorter but thicker stems. This variability is represented in the scatterplots of Figure 3a,b and illustrates that shrub structural variability cannot be sufficiently explained using a measurement from one single perspective alone. The structural heterogeneity of shrubs is a function of site conditions, such as nutrient, water, and light availability. To capture structural heterogeneity, a measurement is needed that describes the shrub structure from three different perspectives. For example, a taller shrub with a single stem will be narrower in width and cover than a shrub with many stems extending in multiple directions. If we assume that the shrub with many stems has a larger width, then it is also likely that the shrub with many stems will have more biomass (dry weight). We demonstrated that neither max stem length nor max basal diameter could be used to predict the dry weight of multi-stemmed shrubs. Volume, however, captured the shrub extent and directional growth and therefore predicted total AGB with less total model uncertainty. The exception of better model performance using max basal diameter for Shepherdia canadensis can be explained by the comparably low number of stems for each harvested individual (<12 stems per plant) and the observed uniform growth of this species in the areas sampled.
A second alternative to volume is the measurement of basal diameters for all stems per plant, converted to cross-sectional area and summed. This is because (a) stem count is represented and (b) shrubs with larger stem counts usually have greater extents (width and line-intercept cover) and dry weight compared to shrubs with lower stem counts. To improve the accuracy of AGB predictions for northern boreal shrubs, measurements to determine volume (max height, line-intercept cover, and width) are recommended. These can be measured rapidly in the field. 1D measures may take slightly less time for shrub individuals that have developed a low number of stems. However, 1D measurements result in both over-and underestimation of AGB depending on the regression form used, especially for shrubs with a high number of stems. 2D cross-sectional area provides the second-best predictions of shrub AGB, although it is slightly more time intensive to measure the basal diameter of each stem per shrub individual. Here, the time required increases with number of stems.

Comparison of Regression Models for Shrub Total AGB Prediction
With regard to model comparisons (LLR, LLRC, and NLS), NLS produced the best model fits for each shrub genus and species (Supplementary Material 3, Table S6) as well as for the multispecies data (Table 2). Here, we found significant differences using 1D vs. 3D models. RMSE varied between 251.68 (NLS, using max basal diameter) and 317.95 (LLRC, using max stem length) g and were greatly reduced with volume as the input variable (RMSE between 141.04 (NLS) and 144.12 (LLRC) g). For the multispecies data, NLS and LLRC overestimated total AGB for all 1D, 2D, and 3D models, while LLR continuously underestimated total AGB ( Table 2). NLS produced the best model fit, independent of the variable used, while the 3D models of all regression forms resulted in the best AGB predictions (RMSE = 141.04 g, R 2 = 0.790 (NLS), RMSE 144.08 g, R 2 = 0.788 (LLR), RMSE = 144.12 g, R 2 = 0.788 (LLRC); Table 2). However, although NLS produces slightly better model fits, nonlinear models require an even variance of errors across the domain of the predictor variable in order to perform valid comparisons of model uncertainties and regression coefficients amongst datasets (e.g., [23]). Residual analysis of our models showed that errors were free of heteroscedasticity. However, when models are transferred to different areas and data, we recommend using the LLRC-based models. This is because biomass data can contain natural heteroscedastic variation. Heteroscedasticity needs to be accounted for in the model development to ensure that model results do not contain bias [25]. For example, Mascaro et al. [25] reported a bias of~100% overestimation when extending predicted small tree (diameter at breast height (DBH) range 2-12 cm) aboveground biomass to stand level biomass using NLS. For model transfer purposes, we have provided regression coefficients and error statistics not only for our best models based on NLS but also for our LLRC models (Supplementary Material 2, Table S3).

Comparison of 1D and 2D Variables for Tree Total AGB Prediction
For the short-stature tree AGB equations, diameter measured at 0.3 m stem length resulted in better model fits compared to stem length for each genus and species (Supplementary Material 3, Table S7). For the multispecies data, total AGB prediction using stem length (1D variable) provided the lowest goodness of fit (highest RMSE and lowest R 2 , Table 4, Figure 4), while inclusion of stem diameter at 0.3 m, alternatively cross-sectional area at 0.3 m, improved predictions by 87% using LLRC or NLS (Table 4).
For multispecies tree total AGB prediction, the exponent α was greater for the 1D model based on stem length and close to unity for the 1D model diameter at 0.3 m, alternatively 2D model cross-sectional area at 0.3 m ( Table 5). The 2D cross-sectional area model produced equivalent results compared to the 1D diameter model but allowed to effectively reconcile tree with shrub AGB predictions. The combined prediction of shrub and tree AGB has not been developed for this study region yet, but would allow the prediction of AGB for boreal short-stature shrubs and trees (<4.5 m) as well as tall-stature trees (>4.5 m) in just a few steps (see Section 3.5). When transferring these allometric models to different areas within the region, we recommend measuring the diameter at 0.3 m stem length.  For multispecies tree total AGB prediction, the exponent was greater for the 1D model based on stem length and close to unity for the 1D model diameter at 0.3 m, alternatively 2D model crosssectional area at 0.3 m ( Table 5). The 2D cross-sectional area model produced equivalent results compared to the 1D diameter model but allowed to effectively reconcile tree with shrub AGB predictions. The combined prediction of shrub and tree AGB has not been developed for this study region yet, but would allow the prediction of AGB for boreal short-stature shrubs and trees (<4.5 m) as well as tall-stature trees (>4.5 m) in just a few steps (see Section 3.5). When transferring these allometric models to different areas within the region, we recommend measuring the diameter at 0.3 m stem length.

Comparison of Regression Models for Tree Total AGB Prediction
For the prediction of tree total AGB per genus/species and all genera/species combined, NLS achieved the lowest RMSE, highest R 2 , and lowest total percentage error compared to measured biomass (<−5%-4%), while the dependent and independent variables were significantly related (p < 0.001, Table 4 and Table S7). LLR predictions resulted in the highest RMSE, similar R 2 , and highest total percentage error relative to the measured biomass compared to LLRC and NLS for each genus/species and for all data combined. The single exception was for predicting total AGB for Picea spp. based on stem length (Table S7). Using LLR, the prediction based on stem length resulted in an underestimation of total AGB of −51% and an underestimation of −20% when using diameter or cross-sectional area, respectively, as input variable for the multispecies models. LLRC had comparably lower RMSE, similar R 2 , and underestimated total multispecies AGB by −6% (stem length) to <−10% (diameter at 0.3 m, cross-sectional area at 0.3 m). In order to address potential heteroscedasticity effects, we have provided the regression coefficients for both NLS and LLRC models with cross-sectional area (measured at 0.3 m stem length) as predictor variable (Supplementary Material 2, Table S4).

Comparison of Regression Models for General Shrub and Tree Total AGB Prediction
The prediction of total AGB for shrubs and trees combined resulted in similar predictive capability to multispecies shrub and multispecies short-stature trees (R 2 ≥ 0.770, p < 0.001). This was determined using the 2D independent variable of cross-sectional area, measured at the base for shrubs and at 0.3 m stem length for trees (Table 6, Figure 5). Relating modeled total AGB to measured total AGB showed no evident bias in the 2D-based prediction of combined shrub and tree AGB in comparison to 3D multispecies shrub and 2D multispecies tree AGB models. This is depicted in Figure 6, which shows modeled AGB in relation to measured AGB of the 2D general shrub and tree AGB model (Figure 6e,f) in comparison to the 3D multispecies shrub (Figure 6a,b) and 2D multispecies tree (Figure 6b,c) AGB models. Similar to the multispecies shrub models, NLS achieved the lowest RMSE and highest R 2 (RMSE = 94.80 g, R 2 = 0.776). The RMSE of model LLR increased by 53% (RMSE = 202.17 g, R 2 = 0.770) and by 54% for the LLRC model (RMSE = 206.37 g, R 2 = 0.770). Compared to Ali et al. [27], who derived best model fits for combined shrub and tree AGB prediction using diameter of the longest stem and total plant height combined, our model results show that AGB of boreal plants can be predicted with a simpler one-variable model using cross-sectional area. For our AGB models based on cross-sectional area, the exponent α was greater for the 1D model based on stem length and close to unity for the 2D model cross-sectional area at 0.3 m (Table 7). However, stem length of shrubs and trees was also weakly related to measured total AGB (p < 0.001, Table 6) and thus represents an alternative to cross-sectional area, which has potential for use in rapid field measurement or non-invasive observation situations (e.g., remote sensing via airborne lidar) where it may be acceptable to trade accuracy at the individual sample-level for greater overall population representation. For model transfer purposes, we recommend the use of LLRC regression coefficients and correction factor in order to address heteroscedasticity.       Figure 6. Measured total AGB related to modeled total AGB using LLRC and NLS for the 3D-based equations for multispecies shrubs (a,b) and the 2D-based equations for multispecies trees (c,d) as well as general shrubs and trees (e,f).

Conclusions
Our analysis has shown that AGB of shrubs can be modeled with higher accuracies when using a 3D field variable, such as volume. Small-stature tree AGB can be most accurately predicted with the stem diameter measured at 0.3 m stem length. In addition, we found that shrub AGB can be reconciled with small-stature tree AGB when using total stem cross-sectional area as the predictor variable. Based on the two best models, we have provided regression coefficients and error statistics for the modeling of short-stature shrub and tree AGB for the region of sporadic to discontinuous permafrost in NWT, Canada. For model uncertainty propagation to total AGB predictions, we have provided the standard error of each model coefficient. We have provided species-and genus-specific as well as multispecies allometric models for shrubs and trees commonly found in boreal forest and peatland ecosystems. These equations are necessary for improving understanding and quantification of biomass change and the potential implications for carbon pools in northern environments, which are highly susceptible to climate change.
Supplementary Materials: The following are available online at http://www.mdpi.com/1999-4907/11/11/1207/s1, Table S1: Descriptive statistic by shrub species and genus by plant compartment (range of values in parentheses, followed by average ± standard deviation); Table S2: Descriptive statistic by tree species by plant compartment (range of values in parentheses, followed by average ± standard deviation); Table S3: Volume-based regression coefficient estimates with error statistics to be input into Equations (4)-(7) as appropriate to derive shrub AGB; Table S4: Volume-based regression coefficient estimates with error statistic to be input into Equations (4)-(7) as appropriate to calculate tree AGB; Table S5: Volume-based regression coefficients with error statistic to be input into Equations (4)-(7) as appropriate to calculate general shrub and tree AGB; Table S6: Model performances for shrub total ABG prediction per genus/species (all p values < 0.001); Table S7: Model performances for tree total ABG prediction per genus/species (all p values < 0.001).