Allometric equations and wood density parameters for estimating aboveground and woody debris biomass in Cajander larch (Larix cajanderi) forests of northeast Siberia
- Faculty of Science, Vrije Universiteit Amsterdam, Amsterdam, 1081 HV, the Netherlands
Correspondence: Clement Jean Frédéric Delcourt (firstname.lastname@example.org)
Boreal forests are particularly vulnerable to climate warming, which increases the occurrence of natural disturbances, such as fires and insect outbreaks. It is therefore essential to better understand climate-induced changes in boreal vegetation dynamics. This requires accurate estimates of variations in biomass across regions and time. This remains challenging in the extensive larch forests of northeast Siberia because of the paucity of allometric equations and physical properties of woody debris needed for quantifying aboveground biomass pools from field surveys. Our study is the first to present values of mean squared diameter (MSD) and specific gravity that can be used to calculate fine dead and downed woody debris loads in Cajander larch (Larix cajanderi) forests using the line-intersect sampling approach. These values were derived from field measurements collected in 25 forest stands in the Republic of Sakha, Russia, and compared with values reported for other prevalent boreal tree species. We developed allometric equations relating diameter at breast height (DBH, at 1.3 m) to stem wood, stem bark, branches, foliage, and aboveground biomass based on measurements of 63 trees retrieved from previous studies. Differences between our allometric models and existing equations were assessed in predicting larch aboveground biomass in 53 forest stands sampled in the Republic of Sakha. We found that using fine woody debris (FWD) parameters from other boreal tree species and allometric equations developed in other regions may result in significantly lower biomass estimates in larch-dominated forests of northeast Siberia. The FWD parameters and allometric equations presented in our paper can be used to refine estimates of aboveground biomass in Cajander larch forests in northeast Siberia.
The boreal forest represents the largest remaining forest on Earth and covers large areas of the high northern latitudes (FAO, 2020). The boreal forest stores nearly a quarter of the global terrestrial carbon pool (Scharlemann et al., 2014) and as such plays a major role in the global carbon cycle. In addition to its role of a long-term carbon sink (Goodale et al., 2002; Gurney et al., 2002; Pan et al., 2011), boreal forests influence energy and water fluxes through biogeophysical processes associated with surface albedo and evapotranspiration (Bonan et al., 1992; Thomas and Rowntree, 1992). The accelerated warming observed in the boreal regions (Serreze et al., 2009; Cohen et al., 2014) is predicted to increase the extent, frequency, and severity of natural disturbances such as fires (Kasischke and Turetsky, 2006; Flannigan et al., 2009; De Groot et al., 2013) and insect outbreaks (Bale et al., 2002; Gustafson et al., 2010), thereby invoking aboveground biomass losses. These changes in land cover may have significant impacts on multiple forcing agents, including from greenhouse gases and aerosols, and through alterations in surface albedo (Bonan et al., 1995; Randerson et al., 2006). Quantifying variations in forest biomass and productivity across regions, ecosystems, and time is therefore essential for assessing climate-induced changes in boreal vegetation dynamics and feedback mechanisms.
Estimating forest carbon stocks is possible from various techniques from local field inventories to spaceborne assessment at regional to global scales (Picard et al., 2012). Remote sensing is an effective tool for forest structure assessments as it provides repeated information over large areas of interest in a cost-effective manner. Indirect measurements of biomass stocks are based on relationships between remotely sensed proxies and detailed field measurements of biomass. In addition to calibrating and validating remotely sensed measurements and constraining biogeochemical models, field measurements provide a level of detail of the vertical and horizontal structure of biomass that is unmatched by other methods. For example, field-based surveys are required to assess the spatial distribution of dead and downed woody debris (DWD) on the forest floor (Woodall et al., 2008).
DWD is a major component in forest ecosystem functioning. It plays a fundamental role in carbon storage (Woodall et al., 2013), influences other nutrient cycles (Finér et al., 2003; Graham and Cromack, 1982), facilitates tree regeneration (Weaver et al., 2009), and serves as habitat for several plant and animal species (Harmon et al., 1986; Freedman et al., 1996). DWD is also an important controlling factor of forest fire behavior (Rothermel, 1972). DWD fuel load affects surface fire intensity (Byram, 1959; Alexander, 1982) and can trigger the transition between surface and crown fire regimes (Van Wagner, 1977). Fuel size is also critical in fire behavior prediction as the surface area-to-volume ratio determines fuel moisture and combustion efficiency (Byram, 1959). Ignition and consumption are generally larger for smaller diameter pieces that exhibit high surface exposure. Quantification of fine and medium-sized dead and downed woody material is therefore essential for assessing fire effects and forest carbon dynamics.
The line-intersect method is one of the most common field-based approaches to sample woody debris in forest stands (Warren and Olsen, 1964; Van Wagner, 1968, 1982; Brown, 1971; Brown and Roussopoulos, 1974; Brown et al., 1982). This technique allows estimation of the total biomass by measuring the diameter or cross-sectional area of all dead and downed woody pieces at their intersection with a transect line (Warren and Olsen, 1964). As this sampling approach can be tedious and time-consuming for smaller pieces, fine woody debris (FWD, <7 cm in diameter) is generally tallied by diameter size class using a go/no-go sizing gauge, and load for any size class can be retrieved from the number of intercepts over the transect line (Van Wagner, 1982). One of the advantages of this method is that it significantly reduces sampling time in the field, but locally derived species-specific values for specific gravity, tilt angle, and diameter for each size class are needed. Such data have been published for most of the North American boreal forest tree species (e.g., Brown, 1974; Sackett, 1980; Delisle and Woodard, 1988; Nalder et al., 1999). However, to the best of our knowledge, there is no available data for larch species in boreal forests of northeast Siberia.
Destructive harvest is the most accurate approach for estimating tree biomass, yet this method requires intensive sampling efforts, and complete harvesting of a forest stand is time-consuming and not suitable for long-term monitoring. To overcome these limitations, one standard technique is to apply allometric equations to convert easily measured variables, such as tree diameter and height, to biomass for each individual tree in a given area. Allometry examines size-correlated variations in rates of plant growth (Huxley and Teissier, 1936; Niklas, 1994) and can be applied to the different parts of a tree (e.g., stem, foliage, and branches). By sampling individual trees over a size range of interest, allometric equations can be developed using regression models for estimating the relationship between tree component biomass and one or more size parameters. Numerous allometric equations have been published for North American boreal tree species (Smith and Brand, 1983; Penner et al., 1997; Ter-Mikaelian and Korzukhin, 1997). In contrast, allometric relationships are more limited in the larch forests of northeast Siberia, which differ from other boreal forests because they consist of deciduous needleleaf trees from the genus Larix Mill. Characterized by its adaptability to grow on permafrost terrain, Siberian larch species (Larix sibirica Ledeb., Larix gmelinii (Rupr) Rupr., Larix cajanderi Mayr.) have evolved physiological and morphological features to endure the long and extremely low air and soil temperatures as well as vigorous cryogenic processes in permafrost terrain. These forests account approximately for 20 % of the world's boreal forests (Osawa and Zyryanova, 2010); store vast amount of carbon, mainly in organic soils (Matsuura and Hirobe, 2010); and contribute to the stability of the permafrost in a warming climate (Abaimov, 2010; Stuenzi et al., 2021). Despite their critical importance to the global carbon cycle, larch forest ecosystems, especially the remote Cajander larch (Larix cajanderi Mayr.) forests of northeastern Siberia, are understudied, and little is known about their response to climate change.
Most of L. cajanderi habitat is located eastwards of the Verkhoyansk Mountain Range where very large differences between summer and winter temperatures prevail (Abaimov, 2010; Chevychelov and Bosikov, 2010). L. cajanderi can grow on various soil types, from poorly developed stony soils of mountainous terrains to rich alluvial soils, and ranges in size from dwarf and shrub-like forms to 25 m tall trees. It is the least-heat-demanding larch species and is well adapted to the underlying continuous permafrost as it can resist freezing temperatures and develop adventitious roots as the active layer deepens (Abaimov, 2010). Few studies have reported allometric equations for L. cajanderi in northeast Siberia (Kajimoto et al., 2006; Alexander et al., 2012), but local allometric equations have been developed based on limited tree harvests near the forest–tundra ecotone (Chersky, 68.74∘ N, 161.40∘ E; 100 m a.s.l.) and the altitudinal treeline (about 100 km west of Oymyakon, 63.45∘ N, 142.77∘ E; 1160 m a.s.l.). As conifer biomass accumulation and allocation are controlled by climate, soil conditions, and stand characteristics (Gower et al., 1995), there may be limitations to applying site-specific allometric equations over large areas away from the measurement sites.
The objective of this study was to expand the biomass allometry of L. cajanderi for northeast Siberia using raw measurements from a comprehensive dataset of biomass measurements in Eurasia (Schepaschenko et al., 2017a, b, c). We developed allometric equations relating diameter at breast height (DBH) to stem wood, stem bark, foliage, branches, and total aboveground biomass. In addition, we sampled fine woody debris of L. cajanderi using the line-intersect method in Siberian larch-dominated forest stands (Delcourt and Veraverbeke, 2022). We derived diameter and specific gravity values per woody debris size classes. The allometric equations and wood density parameters presented in this work will be of use to researchers who want to quantify aboveground and woody debris biomass in Cajander larch forests of northeast Siberia. These data will facilitate improved quantification and understanding of the dynamics of these two major boreal forest carbon pools in northeast Siberia.
2.1 Fine woody debris sampling
The line-intersect method is a widely used approach to quantify FWD lying on the ground in a forest stand (Warren and Olsen, 1964; Van Wagner, 1968; Brown, 1971). It requires measuring the diameter of each piece of wood at its intersection with a transect line which can be considered as a strip of infinitesimal width containing a series of cross-sectional areas (Van Wagner, 1982). The sum of cross-sectional areas divided by the length of the transect line can be converted to volume by multiplying both numerator and denominator by width. FWD load, or weight per unit ground area, is then obtained from Eq. (1) by multiplying the volume by the specific gravity of wood as follows (Van Wagner, 1982):
where W is the FWD load, is a probability factor that allows us to sum the cross-sectional areas as circles, d is the piece diameter, is the factor required to convert d2 into circular area, L is the length of the transect line, and G is the specific gravity in units of weight per unit volume. Equation (1) assumes that woody pieces are horizontal and does not account for ground slope. To minimize the bias related to tilted pieces that are less likely to be intercepted by the transect line, W can be multiplied by a correction factor equal to the secant of the piece tilt angle relative to horizontal (Brown and Roussopoulos, 1974). Similarly, a correction factor can be calculated from the ground slope angle as follows (Brown, 1974):
where s is the slope correction factor and α is the ground slope angle (degrees). Consequently,
where h is the angle between the piece and the horizontal plane (degrees). Measuring diameter on each intersected piece along a transect line can be tedious and time-consuming, especially where small pieces are abundant. In practice, FWD is tallied by diameter size class using a go/no-go sizing gauge, and the number of intercepts over the transect line is reported for each class (Brown, 1974; McRae et al., 1979). Therefore, the term ∑d2 in Eq. (3) is replaced by , where ni is the number of intercepts over the transect line within the diameter size class i, and Di is the representative class diameter (Van Wagner, 1982). The quadratic mean diameter (QMD) is generally used as the appropriate class diameter so that load for any species can be calculated as follows (Van Wagner, 1982; Nalder et al., 1999):
where Wi is the FWD load (t ha−1), Gi is the specific gravity (g cm−3), hi is the piece tilt angle (degrees) of the diameter size class i, and QMDi is the quadratic mean diameter (cm) of the size class i given by
To determine wood density and size parameters for L. cajanderi, we sampled natural forest stands within the boreal forests of the Republic of Sakha (Russia) in the summer of 2019 (Fig. 1 and Table 1). Our study sites were located on both sides of the Lena River within the Central Yakutian Lowland (0–200 m a.s.l.) and the near-Lena plateau (200–500 m a.s.l.) (Chevychelov and Bosikov, 2010). Larix cajanderi largely prevails in these areas of continuous permafrost with inclusion of Pinus sylvestris communities (Isaev et al., 2010). We selected larch-dominated stands within a gradient of vegetation structure, stand age, and landscape position. At each site, we established a 30 m×30 m quadrant, within which measurements were made along a 2 m×30 m belt transect in the north–south direction intersecting the quadrant's centroid. We recorded latitude, longitude, elevation, slope, aspect, and a general site description. To characterize stand structure and composition, we inventoried every tree within the belt transect by recording species and DBH, and we calculated stem density and basal area. We estimated stand age from basal tree disks or increment cores collected on five trees of the dominant cohort.
FWD was sampled for L. cajanderi using the line-intersect approach. In each site, we inventoried every piece (e.g., twigs, limbs, branches) that intersects the 30 m transect line (west side of our belt transect) in each of the following five roundwood diameter size classes: 0–0.49, 0.50–0.99, 1–2.99, 3–4.99, and 5–6.99 cm (McRae et al., 1979), hereafter referred to as classes I to V. For the first five pieces in each size class, we measured the distance along the transect (m) and their diameters (cm) at their points of intersection with the transect line (Fig. 2) using a caliper to 1 mm. These pieces were systematically collected and returned to the laboratory for specific gravity determination using the water displacement method in accordance with the ASTM International D2395–14 Method B-II (ASTM International, 2014). First, we measured the oven-dry mass of each sample by drying to constant mass (i.e., less than 0.2 % mass change over a 5 h period of drying) at 105 ∘C (Roussopoulos and Johnson, 1973). After drying, the water displacement method requires us to seal the surface of woody pieces to prevent bias in volume determination resulting from water absorption during immersion. Each oven-dried sample was then covered with a thin impermeable layer by immersion in hot liquid paraffin (solidification point 57–60 ∘C, 0.90 g cm−3 at 20 ∘C), and the mass of the coated piece was measured again before volume determination. As suggested by Nalder et al. (1999), we used a paraffin temperature of 180 ∘C to minimize the error associated with the use of a finite-thickness paraffin layer. Using a pointed, slender rod, each sample was submerged in a container of water placed on a balance. The mass of water after immersion in grams corresponds to the volume of the immersed sample in cubic centimeters. Finally, specific gravity of each sample was determined as follows (ASTM International, 2014):
where G is the specific gravity (g cm−3), K is a constant equal to 1 when the mass is in grams and the volume is in cubic centimeters, m0 is the oven-dry mass (g), and V0 is the oven-dry volume (cm3) given by
where mw,disp is the mass of water displaced (g), mcoated is the mass of the oven-dried sample after immersion in the paraffin (g), and ρparaffin is the density of the paraffin wax (g cm−3).
For each site and size class, we calculated the mean specific gravity and mean squared diameter (MSD) using a similar approach to Sackett (1980) and Nalder et al. (1999). Defined as the sum of the squared diameters divided by the number of samples, MSD is a direct measure of cross-sectional area. It is therefore more adequate than QMD to interpret differences in dead and downed woody debris biomass between species and regions (Nalder et al., 1999). All statistical analyses were performed using R statistical software version 4.0.3 (R Core Development Team, 2021). We compared specific gravity and MSD as a function of diameter size class using a one-way analysis of variance (ANOVA) followed by a Tukey–Kramer post hoc analysis for multiple comparisons. Each MSD value was divided by the square of the arithmetic class center (ACC) of the corresponding size class (i.e., 0.25, 0.75, 2, 4, and 6 for classes I to V) to test for the effect of size classes relative to class midpoint. In cases where the assumptions of ANOVA were not met, we used the non-parametric Kruskal–Wallis test and the Wilcoxon rank sum test for pairwise comparisons with Benjamini–Hochberg corrections for multiple testing as implemented in the R “stats” package.
To simplify FWD loads calculations, size-class-specific values of specific gravity, MSD, and tilt angle can be combined into a single factor M (g cm−1) as follows (Nalder et al., 1999):
While tilt bias can be significant in fresh logging slash where smaller pieces are attached to larger ones (e.g., correction factor as high as 1.38) (Brown and Roussopoulos, 1974), a FWD inventory made in harvested and natural Canadian boreal forest sites indicated that it was generally not large (Nalder et al., 1999). To facilitate comparisons with other species and regions, we used a tilt correction factor of 1.13 as suggested by Brown (1974). The single factor M defined in Eq. (8) represents the load per intercept per meter of transect and can then be used to calculate FWD loads for each diameter size class i using the formula
We used Eq. (9) to calculate FWD loads at 47 of 53 forest stands where Larix cajanderi fine woody debris was observed (Table 1) using M factors derived from this study. For each site, we derived the percentage difference between these estimates and FWD loads computed using M factors from other species and boreal regions.
2.2 Biomass data
We retrieved raw biomass data from a comprehensive dataset of destructive sampling measurements compiled by Schepaschenko et al. (2017b, c). This dataset summarizes measurements from about 1200 experiments conducted between 1930 and 2014 in Eurasia. This dataset provides valuable data that were previously difficult to access for the non-Russian scientific community due to limited publications in English language literature. For each sample plot, the dataset consists of forest stand parameters (e.g., tree species composition, stand age, stem density) as well as tree level variables including DBH, stem height, age, and biomass. Tree biomass is reported in kilograms of dry matter for four components, including stem, stem bark, branches, and foliage, as well as for the total aboveground pool. All the data contained in this dataset were collected using similar sampling strategies. We selected raw harvest data from L. cajanderi stands sampled as part of two experiments conducted in the Magadan Oblast and the Ust-Yansky district (Republic of Sakha), Russia (Schepaschenko et al., 2017b, c) (Fig. 1). Table 2 summarizes the characteristics of larch-dominated forest stands sampled in both sites.
2.3 Development of allometric equations
We developed site-specific allometric equations relating tree DBH to component biomass and total aboveground biomass by fitting power functions of the form
where Y is the dry biomass (kg), DBH is the tree diameter at breast height (cm), and a and b are the regression coefficients. Power relationships are commonly used when describing plant allometry and can be approximated following different regression techniques. We first applied a double logarithmic transformation of both biomass and diameter variables as in Eq. (11) and fitted the linear regression using the ordinary least squares method.
The logarithmic transformation introduces a downward bias in the prediction of plant biomass (Baskerville, 1972). A correction factor (CF) was therefore calculated using Eq. (12) for each regression and multiplied to the value of a to correct for this systematic bias (Sprugel, 1983):
where is the residual standard deviation. The log–log form depicted in Eq. (11) has been widely used to model allometric relationships as it corrects for heteroscedasticity and facilitates statistical comparison between several equations compared to curvilinear regressions (Gower et al., 1999). However, the performance of the double logarithmic transformation model may be reduced when predicting biomass over a wide range of tree diameters (Bond-Lamberty et al., 2002).
Equation (10) was then fitted using the weighted nonlinear regression model as implemented by the function “nls” in the R “stats” package. Originally developed in forestry in the 1960s (Cunia, 1964; Wharton and Cunia, 1987), the weighted least squares method can be used to treat models which have non-random residuals (Brown et al., 1989; Parresol, 1999; Moore, 2010). This method involves assigning each sample a positive weight wj as follows (after Picard et al., 2012):
where εj is the residual error and σj is the residual standard deviation for each observation j. In biomass studies, the heteroscedasticity of the residuals is often approximated by a power relationship between the residual variance and the diameter of trees (Picard et al., 2012) such that
where k>0 and c≥0. Consequently,
As shown in Eq. (16), the weighting is defined by the value of the exponent c that must be fixed prior to fitting the nonlinear regression. For each site and biomass component, the value of c was determined by (1) approximation of the conditional variance of biomass and (2) visual assessment of the plots of the weighted residuals against the fitted values (Table B1; Figs. B1 and B2). In the first approach, we followed Picard et al. (2012) by dividing DBH values in K classes centered on DBHk, where . We selected K=5 in this study and visually checked that the power model was appropriate for modeling the residual variance. A linear regression was fitted between the standard deviation of biomass and the median DBH of each class k using a double logarithmic transformation, and the value of the exponent c was approximated as the slope of this regression. We report regression coefficients and root mean square errors (RMSEs) for each linear and nonlinear regression equation. Standard errors of the regression coefficients are also provided for uncertainty propagation purposes. To determine whether there were differences in allometric equations among sites, we fitted one model where regression parameters varied from site to site (reduced model) and another model where the parameters were identical for both sites (full model). We tested the null hypothesis that a and b parameters did not differ between the two sites using the extra sum-of-squares F test (Motulsky and Christopoulos, 2004). Significant test results would suggest site-specific variations in the allometric relationships between DBH and component biomass. Our allometric models were compared with existing biomass equations for L. cajanderi published by Kajimoto et al. (2006) and Alexander et al. (2012). We assessed differences between models in predicting larch aboveground biomass in 53 forest stands located in the study area described in Sect. 2.1 (Table 1). In each stand, we inventoried every tree within the belt transect and measured DBH (at the height of 1.3 m from the base; Delcourt and Veraverbeke, 2022). Significant differences among models were assessed using the non-parametric Kruskal–Wallis test and the Wilcoxon rank sum test for pairwise comparisons (0.05 significance level).
3.1 Fine woody debris
We sampled a total of 25 forest stands, and we measured diameter on 223 pieces and specific gravity on 95 pieces of dead and downed wood of L. cajanderi. Summary values of MSD and specific gravity are reported in Table 3. Thirteen sites were intermediate aged (58.1 ± 6.8 years old) and 11 were mature (123.5 ± 24.5 years old). Basal areas and tree densities ranged from 1.9 to 15.7 m2 ha−1 and 1170 to 69 830 stems ha−1, respectively. The average dominance of L. cajanderi was 64 %, and larch trees co-occurred with either Scots pine (Pinus sylvestris) or silver birch (Betula pendula) or both in 21 of the 25 sampled sites. Specific gravity ranged from 0.625 to 0.672 g cm−3 (Table 3) and differed as a function of diameter size classes (ANOVA, , p=0.026), with Class II being significantly higher than Class III (Tukey, p<0.05). Site-level MSD/ACC2 values were normally distributed (Shapiro–Wilk, p=0.08), but the assumption of homogeneity of variance was violated (Levene, p<0.05). The log-transformed data did not satisfy the assumptions of normality and homoscedasticity (Shapiro–Wilk, p<0.05; Levene, p<0.05). MSD varied significantly across size classes (Kruskal–Wallis, χ2=21.8, df=4, p<0.001), with Class II being significantly different from Class III (Wilcoxon, p<0.05) and Class IV (Wilcoxon, p<0.001).
Specific gravity was the highest in Class II and then decreased with larger pieces until Class V, where it had intermediate values (Table 3). This pattern in the variation of G across diameter size classes is similar to that of Larix laricina (tamarack) in the boreal forests of western Canada (Nalder et al., 1999), although our values for L. cajanderi are higher. Discrepancy in specific gravity between species can be attributed to differences in live tree stem wood specific gravity but also to decomposition rates and growth ring patterns that are both influenced by climate. In particular, the tighter the growth rings, the greater the proportion of cell walls in the wood cells. It is worth noting that our approach to measure specific gravity has some limitations associated with (1) the removal of rotten pieces for which we were not able to determine the volume, (2) the inevitable uptake of moisture between oven-drying and volume determination, and (3) the use of paraffin wax to seal the surface of woody pieces before their immersion in water. The latter may lead to an overestimation of the oven-dry volume resulting from the filling of small cracks with paraffin (Nalder et al., 1999). Furthermore, the number of pieces sampled for specific gravity and MSD determination is relatively small for larger size classes.
Despite these limitations, this study is the first to report measurements of specific gravity and mean squared diameter for Larix cajanderi in the boreal forests of northeast Siberia. Nalder et al. (1999) showed that these two variables can vary substantially between species and regions, emphasizing the need for locally derived species-specific data. When calculated for size classes II, III, IV, and V from values reported in Table 3, the single factor M varies considerably from values derived by Nalder et al. (1999) for common boreal tree species including Picea mariana (black spruce), Picea glauca (white spruce), and Larix laricina in the Canadian Northwest Territories and Saskatchewan (Table 4). For 19 of the 20 comparisons, M was larger for Larix cajanderi than for the other species, with a mean difference of 23 % (range: −2 %–46 %) (Table 4). Fine woody debris biomass averaged 2.19 t ha−1 in the 47 forest stands sampled in the Republic of Sakha when calculated from size-class-specific M values derived in similar Cajander larch forest ecosystems (Table A1). Figure 3 reveals that using M factors derived for other boreal tree species and regions resulted in lower FWD loads between on average 9 % (Picea mariana, Northwest Territories) and 28 % (Larix laricina, Saskatchewan). This suggests that using specific gravity and MSD values from other species and/or regions would significantly underestimate dead and downed wood biomass in larch-dominated boreal forests of northeast Siberia.
3.2 Allometric relationships
In line with previous studies in northeast Siberia (Kajimoto et al., 2006; Alexander et al., 2012), DBH-based site-specific allometry is a reliable approach for estimating aboveground tree biomass in mature stands of L. cajanderi. All equations fitted using linear regressions applied to log-transformed data were statistically significant (p<0.001) (Table C1). At Magadan and Ust-Yansky sites, coefficients of determination (R2) for individual biomass components ranged from 0.84 to 0.96 and 0.93 to 0.99, respectively. Weighted nonlinear regression models resulted in better fits in all cases (Table 5 and Fig. C1). DBH was a significant predictor (p<0.05) of tree components biomass at Magadan (RMSE=2.60–46.29 kg) and Ust-Yansky (RMSE=0.51–5.57 kg) sites (Table 5). Higher RMSE observed at Magadan sites are explained by the larger range of DBHs. Most biomass equations for boreal tree species were derived from linear regressions on log-transformed data (Ter-Mikaelian and Korzukhin, 1997; Bond-Lamberty et al., 2002). However, the slopes and intercepts of the relationships between tree diameter and biomass may differ with stand age (Lieffers and Campbell, 1984; Campbell et al., 1985), which limits the efficiency of the logarithmic transformation for describing plant allometry over a wide range of tree sizes (Bond-Lamberty et al., 2002; Kajimoto et al., 2006). Similar to a previous work on shrub allometry in boreal forests (Berner et al., 2015), our results show that weighted nonlinear regression models can be used to predict plant biomass. Contrary to previous works on Cajander larch allometry, the set of equations presented in this study includes bark as a separate biomass component. Quantifying bark biomass in Siberian larch forests is important, for example, when estimating fuel consumption from wildfires, as bark and wood are characterized by different combustion processes (Lestander et al., 2012).
Our regression models performed poorer in predicting foliage and branches biomass, particularly in the Magadan site (R2<0.90), as shown by the higher errors of the regression coefficients (Table 5). This is in agreement with other studies on shrubs (Buech and Rugg, 1995; Berner et al., 2015) and trees (Gower et al., 1997; Bond-Lamberty et al., 2002) that found weaker allometric models for components related to annual production such as foliage and new growth. It is also clear from Fig. 4b and c and Fig. C2b and c that there is more discrepancy between allometric models for these biomass components, which is likely to be caused by local variations in light, water, and nutrient availability (Bond-Lamberty et al., 2002). While some studies have used crown-base diameter and sapwood area to improve biomass estimation for foliage and branches (Waring et al., 1980; O'Hara, 1988; Comeau and Kimmins, 1989; Osawa, 1990), Kajimoto et al. (2006) found that regression models based on these variables did not always result in better fits in L. gmelinii and L. cajanderi stands. Including tree height in addition to DBH as independent variables is another alternative often used for boreal tree species (Schmitt and Grigal, 1981; Ker, 1984; Harding and Grigal, 1985), including Siberian larch species (Sawamoto et al., 2003). However, it may be more difficult to measure tree height accurately in the field and more time-consuming to apply such an allometric model for biomass estimation, especially in dense forest stands. We tested tree height (H) as a predictor of tree components biomass (Table C2) and found that the improvement in model fits may not be large enough to counterbalance the additional practical difficulties associated with making in situ tree height measurements.
Figure 4d shows differences between our newly developed allometric models and existing aboveground biomass equations developed for Larix cajanderi trees in northeast Siberia (Kajimoto et al., 2006; Alexander et al., 2012). To illustrate these differences, we calculated aboveground biomass at 53 forest stands near Yakutsk (Fig. 1 and Table 1) using each allometric model. Total aboveground biomass estimated by applying site-specific allometry from the Magadan Oblast averaged at 6.58 ± 3.94 kg m−2 (range: 0–17.9 kg m−2) and at 7.21 ± 3.42 kg m−2 (range: 0–15.1 kg m−2) when using the equation developed for the Ust-Yansky district (Fig. 5). Estimates from allometric relationships developed for the Chersky and Oymyakon areas resulted in significant lower biomass estimates (Wilcoxon, p<0.05), ranging from 3.64 to 4.40 kg m−2 (Fig. 5). A similar trend was found when predicting stem biomass (Figs. 4a and C2a). Aboveground biomass is assumed to vary along a latitudinal gradient within larch forests of northeast Siberia (Usoltsev et al., 2002), ranging from 0.1 kg m−2 in northern regions to 18 kg m−2 in more productive southern stands (Usoltsev, 2001; Kajimoto et al., 2010). Using the allometric equations developed in this study, mean aboveground biomass estimates across our larch forest stands were 33 %–48 % higher than predictions from existing biomass equations. Our estimates were in good agreement with allometry-derived values reported for similar larch forest ecosystems located in nearby areas (Fig. D1). Siewert et al. (2015) reported a mean aboveground tree biomass of 7.2 kg m−2 in 12 L. cajanderi dominated forest stands in the Spasskaya Pad/Neleger study area (62.23∘ N, 129.62∘ E). Our study area is in the contact zone between Larix cajanderi and another closely related Siberian larch species, Larix gmelinii (Rupr.) Rupr. (Abaimov, 2010). Schulze et al. (1995) and Sawamoto et al. (2003) used allometric relationships to compute aboveground biomass in L. gmelinii forest stands. Their estimates ranged from 4.4 kg m−2 (49 years old) to 12.0 kg m−2 (125 years old) and from 2.0 kg m−2 (25 years old) to 10.5 kg m−2 (170 years old). The biomass estimates at our study sites are of the same order as those of previous studies in similar mature larch forest ecosystems. Further efforts are needed to verify and validate the utility of our newly developed allometric equations. Part of this effort could focus on cross-comparing biomass estimates derived from allometric equations with remote sensing measurements from optical and light detection and ranging (lidar) sensors (Miesner et al., 2022).
Previous studies have shown that caution is needed when applying regional biomass equations in other regions (Bond-Lamberty et al., 2002; Jenkins et al., 2003; Berner et al., 2015). Differences in the relationship between tree diameter and biomass may relate to differences in local site conditions, including stand age, tree density, growing season length, drainage, and nutrient availability (Schulze et al., 1995; Kajimoto et al., 2010). By influencing these factors, permafrost and fire disturbances may shift biomass allocation of Siberian larch trees (Kajimoto et al., 1999; Alexander et al., 2012). Stand structure and dynamics in the Chersky and Oymyakon areas differ with our study sites in the Magadan Oblast and the Ust-Yansky district (Table 2). Allometric equations provided in Alexander et al. (2012) were derived from trees harvested near Chersky in one young stand (15 years old) and two intermediate-aged stands (60 years old), as well as from a mature stand (155 years old) described in Kajimoto et al. (2006). The mature stand (140 years old) sampled near Oymyakon was located near the altitudinal treeline of L. cajanderi on the middle part of a north-facing slope (10∘) and consisted of relatively small trees (mean DBH=5.93 cm, mean H=4.31 m) (Kajimoto et al., 2006). All sites described in this study are underlain by continuous permafrost; however, differences in active layer dynamics may contribute to differences in allocation patterns. Soil water and nitrogen (N) availability increase as the active layer deepens, which may allow trees to store more aboveground carbon. In contrast, a shift in proportion of carbon investment from the growth of stem and/or branches to roots development has been reported in northeastern Siberian larch forests growing on permafrost soils with thin active layers (Kajimoto et al., 1999, 2006). This root-oriented biomass allocation pattern is primarily explained as a response of larch trees to N-poor soils (Kajimoto et al., 2010; Osawa et al., 2010). The particular soil conditions that prevail on north-facing slopes of mountainous terrain, such as lower soil temperature and shallower active layers, and the limited direct solar radiation may reduce the rate of photosynthesis and respiration of larch trees (Koike et al., 2010). These limited resources of nutrients and light may explain why larch trees at Oymyakon sites (OM) tended to allocate less aboveground biomass per unit DBH (Fig. 4). Our study shows that using our newly developed allometric equations may result in higher biomass estimates compared to values derived from existing DBH-based biomass equations. This may have important implications for understanding changes in boreal vegetation dynamics, carbon, and energy budgets in a warming climate. Our results suggest that site-specific allometric equations should be applied near the sampling sites from which they were developed. Given the vastness of the area covered by Cajander larch and the comparatively limited number of available allometric equations, a practical challenge arises when one would want to measure tree biomass in areas that are farther away from the four sampling locations described in this study (i.e., the Ust-Yansky district, Chersky and Oymyakon areas in the Republic of Sakha, and the Magadan Oblast). In this case, we recommend using site-specific equations when site parameters (e.g., stand age, DBH, and tree height) closely match with one of the sampling sites. When this is not the case, we recommend that one calculates tree biomass from the five different equations (Fig. 4). This will result in a mean estimate and an uncertainty range. We acknowledge that the number of samples resulting in our newly developed allometric equations is still relatively small, and these equations can further be refined in the future when more data becomes available that cover a wider range of climate conditions, tree ages, and stem sizes. Efforts in this direction are already underway (Miesner et al., 2022).
Quantifying regional and temporal variations in biomass of boreal forests is essential for assessing the effects of climate warming on these ecosystems. This study contributes much needed data and equations for estimating aboveground biomass in larch-dominated forests of northeast Siberia. We collected measurements of specific gravity and mean squared diameter from fine dead and downed woody debris in the Republic of Sakha (Russia). These values, the first available for Larix cajanderi, are required to calculate woody debris biomass using the line-intersect method. We showed that using data from other dominant boreal tree species may result in large biases.
Allometric equations relating DBH to aboveground tree biomass, as well as individual tree compartments (i.e., stem wood, stem bark, branches, foliage), were developed for two sites in the Ust-Yansky district and the Magadan Oblast (Russia). In agreement with previous studies, our analysis revealed that DBH-based allometry was a simple and reliable approach for estimating L. cajanderi biomass in Siberian forests. However, we found significant differences between our newly developed allometric relationships and existing equations, which are likely due to differences in stand characteristics and environmental conditions at sampled sites. Our results suggest that site-specific allometric equations should be applied with caution in other regions. While our newly developed allometric equations have further extended the geographical coverage of available allometric equations within the range of the L. cajanderi forests in northeast Siberia, our study also demonstrated the critical need for additional field surveys to reduce uncertainties in biomass pools estimates. The data and equations presented in our study contribute to the quantification of aboveground boreal biomass pools of the larch forests of northeast Siberia.
Codes required to derive the allometric equations and woody debris parameters provided in this study and to reproduce related figures and tables are available from https://github.com/delcourtc/larix-allometry.git (last access: 14 September 2022; https://doi.org/10.5281/zenodo.7079233, Delcourt, 2022).
Tree biomass and size measurements at Magadan and Ust-Yansky sites are freely available from https://doi.org/10.1594/PANGAEA.871491 (Schepaschenko et al., 2017c). Fine woody debris data and tree size measurements collected at Batamay and Yert sites are available from https://doi.org/10.5281/zenodo.7049450 (Delcourt and Veraverbeke, 2022).
CJFD designed the research with input from SV. CJFD organized the field campaign with input from SV and the individuals mentioned in the Acknowledgements. CJFD performed the analysis, created figures, and wrote the manuscript with input from SV. SV acquired funding for this research and supervised the work.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
We would like to thank Trofim C. Maximov and Roman E. Petrov for logistical and field support. We wish to thank Brian Izbicki, Brendan M. Rogers, Rebecca C. Scholten, Tatiana A. Shestakova, and Dave van Wees for contributing to data collection. We are grateful to Hans Cornelissen and Richard van Logtestijn for discussions on woody debris and tree core sampling. We are thankful to Ute Sass-Klaassen and Linar Akhmetzyanov for their input on dendrochronology. We thank Dmitry Schepaschenko for compiling a comprehensive dataset of biomass measurements for Eurasia and making these data available.
This research has been supported by the Aard- en Levenswetenschappen, Nederlandse Organisatie voor Wetenschappelijk Onderzoek (grant no. 016.Vidi.189.070).
This paper was edited by David Medvigy and reviewed by three anonymous referees.
Abaimov, A. P.: Geographical Distribution and Genetics of Siberian Larch Species, in: Permafrost Ecosystems: Siberian Larch Forests (Ecological Studies), edited by: Osawa, A., Zyryanova, O., Matsuura, Y., Kajimoto, T., and Wein, R., Springer, Dordrecht, the Netherlands, 41–58, https://doi.org/10.1007/978-1-4020-9693-8, 2010.
Alexander, H. D., Mack, M. C., Goetz, S., Loranty, M. M., Beck, P. S. A., Earl, K., Zimov, S., Davydov, S., and Thompson, C. C.: Carbon Accumulation Patterns During Post-Fire Succession in Cajander Larch (Larix cajanderi) Forests of Siberia, Ecosystems, 15, 1065–1082, https://doi.org/10.1007/s10021-012-9567-6, 2012.
Alexander, M. E.: Calculating and interpreting forest fire intensities, Can. J. Botany, 60, 349–357, https://doi.org/10.1139/b82-048, 1982.
ASTM International: ASTM D2395-14: Standard test methods for density and specific gravity (relative density) of wood and wood-based materials, ASTM International, West Conshohocken, Pennsylvania, United States, https://doi.org/10.1520/D2395-14, 2014.
Bale, J. S., Masters, G. J., Hodkinson, I. D., Awmack, C., Bezemer, T. M., Brown, V. K., Butterfield, J., Buse, A., Coulson, J. C., Farrar, J., Good, J. E. G., Harrington, R., Hartley, S., Jones, T. H., Lindroth, R. L., Press, M. C., Symrnioudis, I., Watt, A. D., and Whittaker, J. B.: Herbivory in global climate change research: direct effects of rising temperature on insect herbivores, Glob. Change Biol., 8, 1–16, https://doi.org/10.1046/j.1365-2486.2002.00451.x, 2002.
Bartalev, S. A., Belward, A. S., Erchov, D. V., and Isaev, A. S.: A new SPOT4-VEGETATION derived land cover map of Northern Eurasia, Int. J. Remote Sens., 24, 1977–1982, https://doi.org/10.1080/0143116031000066297, 2003.
Baskerville, G. L.: Use of Logarithmic Regression in the Estimation of Plant Biomass, Can. J. Forest Res., 2, 49–53, https://doi.org/10.1139/x72-009, 1972.
Berner, L. T., Alexander, H. D., Loranty, M. M., Ganzlin, P., Mack, M. C., Davydov, S. P., and Goetz, S. J.: Biomass allometry for alder, dwarf birch, and willow in boreal forest and tundra ecosystems of far northeastern Siberia and north-central Alaska, Forest Ecol. Manag., 337, 110–118, https://doi.org/10.1016/j.foreco.2014.10.027, 2015.
Bonan, G. B., Pollard, D., and Thompson, S. L.: Effects of boreal forest vegetation on global climate, Nature, 359, 716–718, https://doi.org/10.1038/359716a0, 1992.
Bonan, G. B., Chapin, F. S., and Thompson, S. L.: Boreal forest and tundra ecosystems as components of the climate system, Climatic Change, 29, 145–167, https://doi.org/10.1007/BF01094014, 1995.
Bond-Lamberty, B., Wang, C., and Gower, S. T.: Aboveground and belowground biomass and sapwood area allometric equations for six boreal tree species of northern Manitoba, Can. J. Forest Res., 32, 1441–1450, https://doi.org/10.1139/x02-063, 2002.
Brown, J. K.: A Planar Intersect Method for Sampling Fuel Volume and Surface Area, Forest Sci., 17, 96–102, 1971.
Brown, J. K.: Handbook for inventorying downed woody material, U. S. Dept. of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station, Ogden, Utah, Gen. Tech. Rep. INT-16, 24 pp., 1974.
Brown, J. K. and Roussopoulos, P. J.: Eliminating biases in the planar intersect method for estimating volumes of small fuels, Forest Sci., 20, 350–356, 1974.
Brown, J. K., Oberheu, R. D., and Johnston, C. M.: Handbook for inventorying surface fuels as biomass in the Interior West, U. S. Dept. of Agriculture, Forest Service, Intermountain Forest and Range Experiment Station, Ogden, Utah, Gen. Tech. Rep. INT-129, 48 pp., https://doi.org/10.2737/INT-GTR-129, 1982.
Brown, S., Gillepsie, A. J. R., and Lugo, A. E.: Biomass estimation methods for tropical forests with applications to forest inventory data, Forest Sci., 35, 881–902, 1989.
Buech, R. R. and Rugg, D. J.: Biomass Relations for Components of Five Minnesota Shrubs, U. S. Dept. of Agriculture, Forest Service, North Central Forest Experiment Station, St. Paul, Minnesota, Res. Pap. NC-325, https://doi.org/10.2737/NC-RP-325, 1995.
Byram, G. M.: Combustion of Forest Fuels, in: Forest fire: control and use, edited by: Davis, K. P., McGraw-Hill, New York, 61–89, ISBN 978-0-07-015484-1, 1959.
Campbell, J. S., Lieffers, V. J., and Pielou, E. C.: Regression equations for estimating single tree biomass of trembling aspen: Assessing their applicability to more than one population, Forest Ecol. Manag., 11, 283–295, https://doi.org/10.1016/0378-1127(85)90106-9, 1985.
Chevychelov, A. P. and Bosikov, N. P.: Natural Conditions, in: The Far North: Plant Biodiversity and Ecology of Yakutia (Plant and Vegetation), edited by: Troeva, E. I., Isaev, A. P., Cherosov, M. M., and Karpov, N. S., Springer, Dordrecht, the Netherlands, 1–23, https://doi.org/10.1007/978-90-481-3774-9, 2010.
Cohen, J., Screen, J. A., Furtado, J. C., Barlow, M., Whittleston, D., Coumou, D., Francis, J., Dethloff, K., Entekhabi, D., Overland, J., and Jones, J.: Recent Arctic amplification and extreme mid-latitude weather, Nat. Geosci., 7, 627–637, https://doi.org/10.1038/ngeo2234, 2014.
Comeau, P. G. and Kimmins, J. P.: Above- and below-ground biomass and production of lodgepole pine on sites with differing soil moisture regimes, Can. J. Forest Res., 19, 447–454, https://doi.org/10.1139/x89-070, 1989.
Cunia, T.: Weighted Least Squares Method and Construction of Volume Tables, Forest Sci., 10, 180–191, 1964.
De Groot, W. J., Cantin, A. S., Flannigan, M. D., Soja, A. J., Gowman, L. M., and Newbery, A.: A comparison of Canadian and Russian boreal forest fire regimes, Forest Ecol. Manag., 294, 23–34, https://doi.org/10.1016/j.foreco.2012.07.033, 2013.
Delcourt, C.: delcourtc/larix-allometry: First release of larix-allometry code (v1.0.0), Zenodo [code], https://doi.org/10.5281/zenodo.7079233, 2022.
Delcourt, C. J. F. and Veraverbeke, S.: Field measurements for estimating aboveground and woody debris biomass in Cajander larch forests of northeast Siberia, Zenodo [data set], https://doi.org/10.5281/zenodo.7049450, 2022.
Delisle, G. P. and Woodard, P. M.: Constants for calculating fuel loads in Alberta, P. M. Canadian Forestry Service, Northern Forestry Centre, Edmonton, Alberta, Forest Management Note 45, https://d1ied5g1xfgpx8.cloudfront.net/pdfs/11688.pdf (last access: 7 September 2022), 1988.
FAO: Global Forest Resources Assessment 2020: Main report, Rome, Italy, https://doi.org/10.4060/ca9825en, 165 pp., 2020.
Finér, L., Mannerkoski, H., Piirainen, S., and Starr, M.: Carbon and nitrogen pools in an old-growth, Norway spruce mixed forest in eastern Finland and changes associated with clear-cutting, Forest Ecol. Manag., 174, 51–63, https://doi.org/10.1016/S0378-1127(02)00019-1, 2003.
Flannigan, M. D., Krawchuk, M. A., De Groot, W. J., Wotton, B. M., and Gowman, L. M.: Implications of changing climate for global wildland fire, Int. J. Wildland Fire, 18, 483–507, https://doi.org/10.1071/WF08187, 2009.
Freedman, B., Zelazny, V., Beaudette, D., Fleming, T., Johnson, G., Flemming, S., Gerrow, J. S., Forbes, G., and Woodley, S.: Biodiversity implications of changes in the quantity of dead organic matter in managed forests, Environ. Rev., 4, 238–265, https://doi.org/10.1139/a96-013, 1996.
Goodale, C. L., Apps, M. J., Birdsey, R. A., Field, C. B., Heath, L. S., Houghton, R. A., Jenkins, J. C., Kohlmaier, G. H., Kurz, W., Liu, S., Nabuurs, G.-J., Nilsson, S., and Shvidenko, A. Z.: Forest Carbon Sinks in the Northern Hemisphere, Ecol. Appl., 12, 891–899, https://doi.org/10.2307/3060997, 2002.
Gower, S. T., Isebrands, J. G., and Sheriff, D. W.: Carbon Allocation and Accumulation in Conifers, in: Resource Physiology of Conifers: Acquisition, Allocation, and Utilization, edited by: Smith, W. and Hinckley, T., Academic Press, 217–254, https://doi.org/10.1016/C2009-0-02454-4, ISBN 978-0-08-092591-2, 1995.
Gower, S. T., Vogel, J. G., Norman, J. M., Kucharik, C. J., Steele, S. J., and Stow, T. K.: Carbon distribution and aboveground net primary production in aspen, jack pine, and black spruce stands in Saskatchewan and Manitoba, Canada, J. Geophys. Res.-Atmos., 102, 29029–29041, https://doi.org/10.1029/97JD02317, 1997.
Gower, S. T., Kucharik, C. J., and Norman, J. M.: Direct and Indirect Estimation of Leaf Area Index, fAPAR, and Net Primary Production of Terrestrial Ecosystems, Remote Sens. Environ., 70, 29–51, https://doi.org/10.1016/S0034-4257(99)00056-5, 1999.
Graham, R. L. and Cromack, K.: Mass, nutrient content, and decay rate of dead boles in rain forests of Olympic National Park, Can. J. Forest Res., 12, 511–521, https://doi.org/10.1139/x82-080, 1982.
Gurney, K. R., Law, R. M., Denning, A. S., Rayner, P. J., Baker, D., Bousquet, P., Bruhwiler, L., Chen, Y. H., Ciais, P., Fan, S., Fung, I. Y., Gloor, M., Heimann, M., Higuchi, K., John, J., Maki, T., Maksyutov, S., Masarie, K., Peylin, P., Prather, M., Pak, B. C., Randerson, J., Sarmiento, J., Taguchi, S., Takahashi, T., and Yuen, C. W.: Towards robust regional estimates of CO2 sources and sinks using atmospheric transport models, Nature, 415, 626–630, https://doi.org/10.1038/415626a, 2002.
Gustafson, E. J., Shvidenko, A. Z., Sturtevant, B. R., and Scheller, R. M.: Predicting global change effects on forest biomass and composition in south-central Siberia, Ecol. Appl., 20, 700–715, https://doi.org/10.1890/08-1693.1, 2010.
Harding, R. B. and Grigal, D. F.: Individual tree biomass estimation equations for plantation-grown white spruce in northern Minnesota, Can. J. Forest Res., 15, 738–739, https://doi.org/10.1139/x85-122, 1985.
Harmon, M. E., Franklin, J. F., Swanson, F. J., Sollins, P., Gregory, S. V., Lattin, J. D., Anderson, N. H., Cline, S. P., Aumen, N. G., Sedell, J. R., Lienkaemper, G. W., Cromack, K., and Cummins, K. W.: Ecology of Coarse Woody Debris in Temperate Ecosystems, Adv. Ecol. Res., 15, 133–302, https://doi.org/10.1016/S0065-2504(08)60121-X, 1986.
Huxley, J. S. and Teissier, G.: Terminology of Relative Growth, Nature, 137, 780–781, https://doi.org/10.1038/137780b0, 1936.
Isaev, A. P., Protopopov, A. V., Protopopova, V. V., Egorova, A. A., Timofeyev, P. A., Nikolaev, A. N., Shurduk, I. F., Lytkina, L. P., Ermakov, N. B., Nikitina, N. V., Efimova, A. P., Zakharova, V. I., Cherosov, M. M., Nikolin, E. G., Sosina, N. K., Troeva, E. I., Gogoleva, P. A., Kuznetsova, L. V., Pestryakov, B. N., Mironova, S. I., and Sleptsova, N. P.: Vegetation of Yakutia: Elements of Ecology and Plant Sociology, in: The Far North: Plant Biodiversity and Ecology of Yakutia (Plant and Vegetation), edited by: Troeva, E. I., Isaev, A. P., Cherosov, M. M., and Karpov, N. S., Springer, Dordrecht, the Netherlands, 143–260, https://doi.org/10.1007/978-90-481-3774-9, 2010.
Jenkins, J. C., Chojnacky, D. C., Heath, L. S., and Birdsey, R. A.: National-scale Biomass Estimators for United States Tree Species, Forest Sci., 49, 12–35, 2003.
Kajimoto, T., Matsuura, Y., Sofronov, M. A., Volokitina, A. V., Mori, S., Osawa, A., and Abaimov, A. P.: Above- and belowground biomass and net primary productivity of a Larix gmelinii stand near Tura, central Siberia, Tree Physiol., 19, 815–822, https://doi.org/10.1093/treephys/19.12.815, 1999.
Kajimoto, T., Matsuura, Y., Osawa, A., Abaimov, A. P., Zyryanova, O. A., Isaev, A. P., Yefremov, D. P., Mori, S., and Koike, T.: Size-mass allometry and biomass allocation of two larch species growing on the continuous permafrost region in Siberia, Forest Ecol. Manag., 222, 314–325, https://doi.org/10.1016/j.foreco.2005.10.031, 2006.
Kajimoto, T., Osawa, A., Usoltsev, V. A., and Abaimov, A. P.: Biomass and Productivity of Siberian Larch Forest Ecosystems, in: Permafrost Ecosystems: Siberian Larch Forests (Ecological Studies), edited by: Osawa, A., Zyryanova, O., Matsuura, Y., Kajimoto, T., and Wein, R., Springer, Dordrecht, the Netherlands, 99–122, https://doi.org/10.1007/978-1-4020-9693-8, 2010.
Kasischke, E. S. and Turetsky, M. R.: Recent changes in the fire regime across the North American boreal region – Spatial and temporal patterns of burning across Canada and Alaska, Geophys. Res. Lett., 33, L09703, https://doi.org/10.1029/2006GL025677, 2006.
Ker, M. F.: Biomass equations for seven major Maritimes species, Canadian Forestry Service, Maritimes Forest Research Centre, Information Report M-X-148, https://d1ied5g1xfgpx8.cloudfront.net/pdfs/7771.pdf (last access: 7 September 2022), 1984.
Koike, T., Mori, S., Zyryanova, O. A., Kajimoto, T., Matsuura, Y., and Abaimov, A. P.: Photosynthetic characteristics of trees and shrubs growing on the north- and south-facing slopes in Central Siberia, in: Permafrost Ecosystems: Siberian Larch Forests (Ecological Studies), edited by: Osawa, A., Zyryanova, O., Matsuura, Y., Kajimoto, T., and Wein, R., Springer, Dordrecht, the Netherlands, 273–287, https://doi.org/10.1007/978-1-4020-9693-8, 2010.
Lestander, T. A., Lundström, A., and Finell, M.: Assessment of biomass functions for calculating bark proportions and ash contents of refined biomass fuels derived from major boreal tree species, Can. J. Forest Res., 42, 59–66, https://doi.org/10.1139/x11-144, 2012.
Lieffers, V. J. and Campbell, J. S.: Biomass and growth of Populus tremuloides in northeastern Alberta: estimates using hierarchy in tree size, Can. J. Forest Res., 14, 610–616, https://doi.org/10.1139/x84-111, 1984.
Matsuura, Y. and Hirobe, M.: Soil Carbon and Nitrogen, and Characteristics of Soil Active Layer in Siberian Permafrost Region, in: Permafrost Ecosystems: Siberian Larch Forests (Ecological Studies), edited by: Osawa, A., Zyryanova, O., Matsuura, Y., Kajimoto, T., and Wein, R., Springer, Dordrecht, the Netherlands, 149–163, https://doi.org/10.1007/978-1-4020-9693-8, 2010.
McRae, D. J., Alexander, M. E., and Stocks, B. J.: Measurement and description of fuels and fire behavior on prescribed burns: A handbook, Canadian Forestry Service, Great Lakes Forest Research Centre, Sault Ste. Marie, Ontario, Information Report O-X-287, 44 pp., https://d1ied5g1xfgpx8.cloudfront.net/pdfs/8963.pdf (last access: 7 September 2022), 1979.
Miesner, T., Herzschuh, U., Pestryakova, L. A., Wieczorek, M., Zakharov, E. S., Kolmogorov, A. I., Davydova, P. V., and Kruse, S.: Forest structure and individual tree inventories of north-eastern Siberia along climatic gradients, Earth Syst. Sci. Data Discuss. [preprint], https://doi.org/10.5194/essd-2022-152, in review, 2022.
Moore, J. R.: Allometric equations to predict the total above-ground biomass of radiata pine trees, Ann. Forest Sci., 67, 806, https://doi.org/10.1051/forest/2010042, 2010.
Motulsky, H. and Christopoulos, A.: Fitting models to biological data using linear and nonlinear regression: A practical guide to curve fitting, Oxford University Press, 352 pp., ISBN 9780198038344, 2004.
Nalder, I. A., Wein, R. W., Alexander, M. E., and de Groot, W. J.: Physical properties of dead and downed round-wood fuels in the Boreal forests of western and Northern Canada, Int. J. Wildland Fire, 9, 85–99, https://doi.org/10.1071/WF00008, 1999.
Niklas, K. J.: Plant Allometry: The Scaling of Form and Process, University of Chicago Press, 412 pp., ISBN 9780226580807, 1994.
O'Hara, K. L.: Stand structure and growing space efficiency following thinning in an even-aged Douglas-fir stand, Can. J. Forest Res., 18, 859–866, https://doi.org/10.1139/x88-132, 1988.
Osawa, A.: Reconstructed development of stem production and foliage mass and its vertical distribution in Japanese larch, Tree Physiol., 7, 189–200, https://doi.org/10.1093/treephys/7.1-2-3-4.189, 1990.
Osawa, A. and Zyryanova, O. A.: Introduction, in: Permafrost Ecosystems: Siberian Larch Forests (Ecological Studies), edited by: Osawa, A., Zyryanova, O., Matsuura, Y., Kajimoto, T., and Wein, R., Springer, Dordrecht, the Netherlands, 3–15, https://doi.org/10.1007/978-1-4020-9693-8, 2010.
Osawa, A., Matsuura, Y., and Kajimoto, T.: Characteristics of permafrost forests in Siberia and potential responses to warming climate, in: Permafrost Ecosystems: Siberian Larch Forests (Ecological Studies), edited by: Osawa, A., Zyryanova, O., Matsuura, Y., Kajimoto, T., and Wein, R., Springer, Dordrecht, the Netherlands, 459–481, https://doi.org/10.1007/978-1-4020-9693-8, 2010.
Pan, Y., Birdsey, R. A., Fang, J., Houghton, R., Kauppi, P. E., Kurz, W. A., Phillips, O. L., Shvidenko, A., Lewis, S. L., Canadell, J. G., Ciais, P., Jackson, R. B., Pacala, S. W., McGuire, A. D., Piao, S., Rautiainen, A., Sitch, S., and Hayes, D.: A large and persistent carbon sink in the world's forests, Science, 333, 988–993, https://doi.org/10.1126/science.1201609, 2011.
Parresol, B. R.: Assessing tree and stand biomass: a review with examples and critical comparisons, Forest Sci., 45, 573–593, 1999.
Penner, M., Power, K., Muhairwe, C., Tellier, R., and Wang, Y.: Canada's Forest Biomass Resources: Deriving Estimates from Canada's Forest Inventory, Canadian Forest Service, Pacific Forestry Centre, Victoria, British Columbia, Information Report BC-X-370, 33 pp., ISBN 0-662-25716-2, 1997.
Picard, N., Saint-André, L., and Henry, M.: Manual for building tree volume and biomass allometric equations: from field measurements to prediction, Food and Agricultural Organization of the United Nations, Rome, Italy, and Centre de Coopération Internationale en Recherche Agronomique pour le Développement, Montpellier, France, 215 pp., E-ISBN 9789251073476, 2012.
R Core Development Team: R: A language and environment for statistical computing, R Foundation for Statistical Computing, Vienna, Austria, https://www.R-project.org/ (last access: 14 March 2022), 2021.
Randerson, J. T., Liu, H., Flanner, M. G., Chambers, S. D., Jin, Y., Hess, P. G., Pfister, G., Mack, M. C., Treseder, K. K., Welp, L. R., Chapin, F. S., Harden, J. W., Goulden, M. L., Lyons, E., Neff, J. C., Schuur, E. A. G., and Zender, C. S.: The impact of boreal forest fire on climate warming, Science, 314, 1130–1132, https://doi.org/10.1126/science.1132075, 2006.
Rothermel, R. C.: A mathematical model for predicting fire spread in wildland fuels, U. S. Dept. of Agriculture, Intermountain Forest and Range Experiment Station, Ogden, Utah, Res. Pap. INT-115, 40 pp., 1972.
Roussopoulos, P. J. and Johnson, V. J.: Estimating slash fuel loading for several Lake States tree species, U. S. Dept. of Agriculture, Forest Service, North Central Forest Experiment Station, St. Paul, Minnesota, Res. Pap. NC-88, 1973.
Sackett, S. S.: Woody fuel particle size and specific gravity of southwestern tree species, U. S. Dept. of Agriculture, Forest Service, Rocky Mountain Forest and Range Experiment Station, Fort Collins, Colorado, Research Note RM-389, 1980.
Sawamoto, T., Hatano, R., Shibuya, M., Takahashi, K., Isaev, A. P., Desyatkin, R. V., and Maximov, T. C.: Changes in net ecosystem production associated with forest fire in taiga ecosystems, near Yakutsk, Russia, Soil Sci. Plant Nutr., 49, 493–501, https://doi.org/10.1080/00380768.2003.10410038, 2003.
Scharlemann, J. P. W., Tanner, E. V. J., Hiederer, R., and Kapos, V.: Global soil carbon: understanding and managing the largest terrestrial carbon pool, Carbon Manag., 5, 81–91, https://doi.org/10.4155/cmt.13.77, 2014.
Schepaschenko, D., Shvidenko, A., Usoltsev, V., Lakyda, P., Luo, Y., Vasylyshyn, R., Lakyda, I., Myklush, Y., See, L., McCallum, I., Fritz, S., Kraxner, F., and Obersteiner, M.: A dataset of forest biomass structure for Eurasia, Sci. Data, 4, 170070, https://doi.org/10.1038/sdata.2017.70, 2017a.
Schepaschenko, D., Shvidenko, A., Usoltsev, V., Lakyda, P., Luo, Y., Vasylyshyn, R., Lakyda, I., Myklush, Y., See, L., Mc-Callum, I., Fritz, S., Kraxner, F., and Obersteiner, M.: Biomass plot data base, PANGAEA [data set], https://doi.org/10.1594/PANGAEA.871465, 2017b.
Schepaschenko, D., Shvidenko, A., Usoltsev, V., Lakyda, P., Luo, Y., Vasylyshyn, R., Lakyda, I., Myklush, Y., See, L., Mc-Callum, I., Fritz, S., Kraxner, F., and Obersteiner, M.: Biomass tree data base, PANGAEA [data set], https://doi.org/10.1594/PANGAEA.871491, 2017c.
Schmitt, M. D. C. and Grigal, D. F.: Generalized biomass estimation equations for Betula papyrifera Marsh., Can. J. Forest Res., 11, 837–840, https://doi.org/10.1139/x81-122, 1981.
Schulze, E.-D., Schulze, W., Koch, H., Arneth, A., Bauer, G., Kelliher, F. M., Hollinger, D. Y., Vygodskaya, N. N., Kusnetsova, W. A., Sogatchev, A., Ziegler, W., Kobak, K. I., and Issajev, A.: Aboveground biomass and nitrogen nutrition in a chronosequence of pristine Dahurian Larix stands in eastern Siberia, Can. J. Forest Res., 25, 943–960, https://doi.org/10.1139/x95-103, 1995.
Serreze, M. C., Barrett, A. P., Stroeve, J. C., Kindig, D. N., and Holland, M. M.: The emergence of surface-based Arctic amplification, The Cryosphere, 3, 11–19, https://doi.org/10.5194/tc-3-11-2009, 2009.
Siewert, M. B., Hanisch, J., Weiss, N., Kuhry, P., Maximov, T. C., and Hugelius, G.: Comparing carbon storage of Siberian tundra and taiga permafrost ecosystems at very high spatial resolution, J. Geophys. Res.-Biogeo., 120, 1973–1994, https://doi.org/10.1002/2015JG002999, 2015.
Smith, B. W. and Brand, G. J.: Allometric biomass equations for 98 species of herbs, shrubs, and small trees, U. S. Dept. of Agriculture, Forest Service, North Central Forest Experiment Station, St. Paul, Minnesota, Research Note NC-299, https://doi.org/10.2737/NC-RN-299, 1983.
Sprugel, D. G.: Correcting for Bias in Log-Transformed Allometric Equations, Ecology, 64, 209–210, https://doi.org/10.2307/1937343, 1983.
Stuenzi, S. M., Boike, J., Gädeke, A., Herzschuh, U., Kruse, S., Pestryakova, L. A., Westermann, S., and Langer, M.: Sensitivity of ecosystem-protected permafrost under changing boreal forest structures, Environ. Res. Lett., 16, 084045, https://doi.org/10.1088/1748-9326/ac153d, 2021.
Ter-Mikaelian, M. T. and Korzukhin, M. D.: Biomass equations for sixty-five North American tree species, Forest Ecol. Manag., 97, 1–24, https://doi.org/10.1016/S0378-1127(97)00019-4, 1997.
Thomas, G. and Rowntree, P. R.: The Boreal Forests and Climate, Q. J. Roy. Meteor. Soc., 118, 469–497, https://doi.org/10.1002/qj.49711850505, 1992.
Usoltsev, V. A.: Forest Biomass of Northern Eurasia: Database and Geography, Russian Academy of Sciences, Ural Branch, Botanical Garden, Yekaterinburg, ISBN 5769112018, 2001.
Usoltsev, V. A., Koltunova, A. I., Kajimoto, T., Osawa, A., and Koike, T.: Geographical Gradients of Annual Biomass Production from Larch Forests in Northern Eurasia, Eurasian J. For. Res., 5, 55–62, 2002.
Van Wagner, C. E.: The Line Intersect Method in Forest Fuel Sampling, Forest Sci., 14, 20–26, 1968.
Van Wagner, C. E.: Conditions for the start and spread of crown fire, Can. J. Forest Res., 7, 23–34, https://doi.org/10.1139/x77-004, 1977.
Van Wagner, C. E.: Practical aspects of the line intersect method, Canadian Forestry Service, Maritimes Forest Research Centre, Fredericton, New Brunswick, Information Report PI-X-12E, 11 pp., ISBN 0-662-11816-2, 1982.
Walker, D. A., Raynolds, M. K., Daniëls, F. J. A., Einarsson, E., Elvebakk, A., Gould, W. A., Katenin, A. E., Kholod, S. S., Markon, C. J., Melnikov, E. S., Moskalenko, N. G., Talbot, S. S., Yurtsev, B. A., and The other members of the CAVM Team: The Circumpolar Arctic vegetation map, J. Veg. Sci., 16, 267–282, https://doi.org/10.1111/j.1654-1103.2005.tb02365.x, 2005.
Waring, R. H., Thies, W. G., and Muscato, D.: Stem Growth per Unit of Leaf Area: A Measure of Tree Vigor, Forest Sci., 26, 112–117, 1980.
Warren, W. G. and Olsen, P. F.: A Line Intersect Technique for Assessing Logging Waste, Forest Sci., 10, 267–276, 1964.
Weaver, J. K., Kenefic, L. S., Seymour, R. S., and Brissette, J. C.: Decaying wood and tree regeneration in the Acadian Forest of Maine, USA, Forest Ecol. Manag., 257, 1623–1628, https://doi.org/10.1016/j.foreco.2009.01.023, 2009.
Wharton, E. H. and Cunia, T.: Estimating tree biomass regressions and their error, proceedings of the workshop on tree biomass regression functions and their contribution to the error, U. S. Dept. of Agriculture, Forest Service, Northeastern Forest Experiment Station, Broomall, Pennsylvania, Gen. Tech. Rep. NE-117, 303 pp., https://doi.org/10.2737/NE-GTR-117, 1987.
Woodall, C. W., Heath, L. S., and Smith, J. E.: National inventories of down and dead woody material forest carbon stocks in the United States: Challenges and opportunities, Forest Ecol. Manag., 256, 221–228, https://doi.org/10.1016/j.foreco.2008.04.003, 2008.
Woodall, C. W., Walters, B. F., Oswalt, S. N., Domke, G. M., Toney, C., and Gray, A. N.: Biomass and carbon attributes of downed woody materials in forests of the United States, Forest Ecol. Manag., 305, 48–59, https://doi.org/10.1016/j.foreco.2013.05.030, 2013.