Next Article in Journal
Coupled Relationship between Soil Physicochemical Properties and Plant Diversity in the Process of Vegetation Restoration
Previous Article in Journal
Evasive Planning for the Management of Eucalyptus Rust Austropuccinia psidii for Espírito Santo State, Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Communication

Stand Age Class to Size Class Crosswalk by Forest Type Group in Minnesota, USA

Department of Forest Resources, University of Minnesota, St. Paul, MN 55108, USA
*
Author to whom correspondence should be addressed.
Forests 2022, 13(5), 647; https://doi.org/10.3390/f13050647
Submission received: 18 March 2022 / Revised: 15 April 2022 / Accepted: 19 April 2022 / Published: 22 April 2022
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)

Abstract

:
Many forestry models require the input of stand size class information, a variable with multiple definitions across forest inventories. This work describes a crosswalk between stand age class and stand size class for several forest type groups in Minnesota to facilitate consistency and availability of the latter information. Refinements to the crosswalk include ratio adjustments that redistribute the number of plots (or hectares) from one estimated size class to another, based on known crosswalk error rates. The results showed that 61.9% of all plots were correctly classified, and 95.5% were within one size class. Correct classifications for individual forest type/size class combinations ranged from 16.7–79.2%. Applying the crosswalk to a validation dataset produced percent errors from −46.4% to 49.4%, but with the ratio adjustment, errors dropped to −20.1% to 14.4%. These results suggest that the crosswalk and ratio adjustments provide a coarse, yet reasonable, substitute to using more complex stand size classification methodologies, particularly when forecasting future stand conditions.

1. Introduction

Stand size class plays a crucial role in many forest and ecological contexts, including when scheduling harvest, prescribing silvicultural treatments, and evaluating wildlife habitat. Descriptive and predictive models of forests may include stand size class information as a response or explanatory variable (e.g., [1,2]). A size class represents the predominant tree size or stage of tree development within a forest stand or landscape (e.g., seedling/sapling class). The U.S. Forest Service, Forest Inventory and Analysis (FIA) program measures forest resources across the entire country and its territories. This program defines three size classes based on the stocking majority of large diameter trees (≥27.94 cm for hardwoods; ≥22.86 cm for softwoods), medium diameter trees (≥12.7 cm and less than large trees), and small diameter trees (<12.7 cm), with diameters generally measured at diameter-breast-height (dbh; 1.37 m above ground) [3]. In addition, most state natural resource agencies (e.g., Minnesota Department of Natural Resources (MNDNR)) record size class during their inventories and use this variable in management decisions [4].
For other operational forest inventories in Minnesota, many do not use size classification methodologies as complex as those used by FIA or MNDNR, which require detailed information and algorithms (see [4,5]). This complexity hinders the addition of such detail in existing and future inventories for other ownerships in the state (e.g., county lands, private landowners, tribal lands). Perhaps more critically, applied modeling and forecasting efforts benefit from straightforward models that increment through time (e.g., [1]). Linking the stand age with other forest attributes (e.g., stand size class) facilitates temporal stand and habitat projections by enabling iterative updates to these variables of interest. Further, the use of stand age alone avoids the complexities of incrementing additional predictor variables before estimating the response. Therefore, this research sought to (1) create an age class to size class crosswalk that allowed for rapid stand size class imputation using stand age class information; (2) provide ratio adjustments that re-assign hectares from one estimated size class to another, based on known crosswalk error rates; and (3) validate the crosswalk and the ratio adjustments using additional observations and quantify the error reduction (or addition) from using the ratio adjustments.

2. Materials and Methods

Per plot, the FIA program uses four 0.0169 ha subplots to measure stand and larger tree (≥12.7 cm dbh) attributes and four 0.0013 ha microplots to measure smaller tree (<12.7 cm dbh) attributes. If a plot encompasses a major shift in certain attributes (e.g., ownership, forest type, stand size class), the plot is subdivided into conditions based on the differences. Annually, FIA measures 20% of the plots in Minnesota, finishing a complete measurement every five years. Training data for this study came from FIA measurement periods from 1999–2013 and comprised 9263 plots [6]. This dataset includes three consecutive statewide inventories (1999–2003, 2004–2008, and 2009–2013). Further, the data was pooled to increase sample sizes [7]. Only sampled plots with one FIA condition and adequate stocking (≥10%) were included to minimize variability and to have an assigned forest type, respectively. For a full description of the FIA database, see [3,8].
Subsequently, an age class to size class map was explored for several forest types by identifying the age class associated with a change in size class. However, small sample sizes prevented determining a reliable map for several forest types. Thus, we focused on four broad forest type groupings to provide sufficient sample sizes: upland conifer, lowland conifer, northern hardwoods, and aspen-birch. These groupings represent distinct habitats and find common reference in Minnesota (e.g., [9]). We then summarized plot counts by forest type group (hereafter referred to as forest types) and five-year age classes and determined the age class where the majority of plots transitioned from one size class to another for each forest type. The selection of the transition points was straightforward, except for a few cases that required discretion to ensure a reasonable choice (e.g., lowland conifer poletimber to sawtimber). Finally, the crosswalk was applied to the training data to compute the overall error rates by forest type and the error rates for individual forest type/size class combinations.
With only three size classes, the crosswalk will assign a plot to either the correct size class or one of two incorrect alternatives. However, the distributions of plots by size class may overlap considerably, leading to high percentages of incorrect assignments. Therefore, the proportions of correct and incorrect assignments were computed within each level of estimated size class. These values comprise a ratio adjustment to the map that redistributes hectares (or the number of plots) from one estimated size class to another, based on the observed crosswalk error rates (Equation (1)). As such, we describe it as the adjusted crosswalk.
h i * = A i h i
where h i * = three element vector of adjusted hectares (or number of plots) by size class, A i = 3 × 3 matrix of ratio adjustments for redistributing estimated size classes, h i = three element vector of unadjusted hectares (or number of plots) by size class, and i = an index for forest type. Note that this adjusted crosswalk only applies to applications using aggregate hectares by forest type and size class. When assigning a size class to an individual plot or stand, the unadjusted crosswalk is appropriate and will be subject to the uncorrected error rates.
The validation of the original and adjusted crosswalk involved using them to estimate size classes for the latest complete FIA inventory (2014–2018) (see [10] for a similar approach). The same filters were applied to the validation data as to the training data, resulting in 2913 plots. Both crosswalk estimates were compared to observed classes via percent errors (see Equation (2)). In addition, percent differences (Equation (3)) between estimates using the original and adjusted crosswalk enabled the quantification of the gains (or losses) from using the adjustment.
e = ( O b s E s t ) / O b s
d = | e c w | | e a c w |
where e = percent error, Obs = observed plot counts (or total hectares) in a forest type/size class combination, Est = estimated plot counts (similar to Obs), d = percent difference, ecw = percent error using the original crosswalk, and eacw = percent error using the ratio adjustment. All analyses were conducted using the R statistical program [11].

3. Results and Discussion

Figure 1 shows the distribution of FIA plots by age class and size class for each forest type group. Where the size class distributions intersect represents the age thresholds between size classes, and Table 1 gives these age classes associated with each size class. Regardless of the forest type, the pattern remains constant: as stands age, they increase in size class. However, for most forest types, the size class distributions overlap considerably (particularly the poletimber and sawtimber classes), suggesting a moderate to weak relationship between age class and size class.
Reapplying the unadjusted crosswalk from Table 1 to the training data showed that 61.9% of all plots were assigned the correct size class, and 95.5% were within one size class (Table 2). However, correct classifications varied within the forest type group. The aspen-birch group had the highest accuracy and lowland conifer the lowest (Table 2). This result corresponds directly to the extent of size class distribution overlap within a forest type (Figure 1). In addition, classification errors showed no strong trend toward larger or smaller size classes. Early attempts to increase the proportion of correct category assignments (including alternative approaches such as regression analysis and mixed-effects modeling) showed marginal improvement while reducing utility.
Table 3 breaks down the error by individual forest type/size class combinations. Within a forest type group, the percentages represent the proportion of observed size classes within each estimated size class, and the diagonal elements give the percentage of plots correctly classified. This table shows that the individual error rates had considerably more variation than the overall error rates in Table 2, ranging from 16.7% to 79.2% correct assignments per forest type/size class combination. Correctly assigning the sawtimber size class to lowland conifer stands proved particularly challenging. These stands develop slowly and may often only produce small to moderate sized trees, thus limiting the number of sawtimber stands for defining a clear size class threshold (as evidenced by the flat sawtimber curve in Figure 1). Northern hardwood poletimber stands were also difficult to classify, as the bulk of these stands were older than the age class intersection between the poletimber and sawtimber distributions. Historically, these forests received low priority management, leading to the majority of stands aging into larger size classes. Upland conifer and aspen-birch provided relatively clean breaks between size class distributions, except that most sawtimber aspen-birch stands were younger than the poletimber/sawtimber intersection. The often intensive management and short rotations of aspen-birch limits the number of sawtimber aspen stands and prevents a clearer threshold.
In order to compensate for the large classification errors, ratio adjustments were computed that would redistribute a portion of plots (or hectares) from one estimated size class to the other size classes, based on misclassification rates (i.e., the adjusted crosswalk). Table 4 shows the ratio adjustments for each forest type/size class combination, with the percentages representing the proportions of an estimated size class within the corresponding observed size classes. For example, for those hectares of upland conifer with an estimated sawtimber size class, 62.3% of the hectares will remain in the sawtimber size class, 19.3% will be assigned the poletimber size class, and 18.4% will be given the seedling/sapling size class. The same interpretation holds for the other forest type/estimated size class combinations.
These ratios suggest that for each level of the estimated size class, the majority of plots were correctly classified (except the poletimber class for lowland conifers). Still, the large variability across ratios (e.g., 0.401 to 0.859 for hectares that remain in the estimated size class) results from the considerable overlap observed in Figure 1 and re-enforces the need for using the adjustments when applicable. Note that although similar, Table 3 and Table 4 provide distinct error rates, with those in the former computed relative to all plots within a true size class and those in the latter computed relative to all plots within a predicted size class.
Table 5 provides the validation statistics when applying the crosswalk with and without the ratio adjustments to the 2014–2018 FIA data. The original crosswalk produced overall percent errors ranging from −46.4% to 49.4% and behaved the poorest for the conifers. Although northern hardwoods had mixed results, the crosswalk did provide more reasonable estimates for aspen-birch. Adding the ratio adjustment (using Equation (1)) improved results considerably, with overall percent errors ranging from −20.1% to 14.4%. Table 5 also provides the percent difference between using the original crosswalk and using the ratio adjustment. The overall differences ranged from −8.6% to 44.4%, with three forest types (upland conifer, lowland conifer, and northern hardwoods) showing a mostly substantial improvement in accuracy after applying the ratios. For aspen-birch, the positive effect was less pronounced, but still with an overall improvement. Three combinations (sawtimber upland conifer, seedling/sapling northern hardwoods, and seedling/sapling aspen-birch) showed slight decreases in accuracy when using the ratio adjustments, likely due to sampling variation. Still, the majority of large, positive differences suggest that using the ratio adjustments should improve results over the unadjusted crosswalk.

4. Conclusions

Ultimately, the age class to size class crosswalk represents a straightforward method for imputing stand size class information. The unadjusted classification errors ranged from −46.4% to 49.4%, but including the ratio adjustments lowered the errors to −20.1% to 14.4%, indicating that the adjusted crosswalk should be used whenever possible. Natural resource practitioners should recognize the coarseness of the map and its adjustment. Still, the crosswalk provides a reasonable substitute for complex stand size classification methodologies, particularly when needing to efficiently forecast or simulate future forest stand conditions.

Author Contributions

Conceptualization, J.M.Z. and A.R.E.; methodology, J.M.Z. and A.R.E.; software, J.M.Z.; validation, J.M.Z.; formal analysis, J.M.Z.; investigation, J.M.Z.; resources, J.M.Z. and A.R.E.; data curation, J.M.Z. and T.S.G.; writing—original draft preparation, J.M.Z.; writing—review and editing, J.M.Z., A.R.E. and T.S.G.; visualization, J.M.Z.; supervision, J.M.Z. and A.R.E.; project administration, J.M.Z.; funding acquisition, A.R.E. and J.M.Z. All authors have read and agreed to the published version of the manuscript.

Funding

Funding for this work was provided by the Interagency Information Cooperative, the Department of Forest Resources, University of Minnesota, and the Minnesota Agricultural Experiment Station and National Institute of Food and Agriculture, U.S. Department of Agriculture, McIntire Stennis capacity grant #1020087. Any opinions, findings, conclusions, or recommendations expressed in this publication are those of the author(s) and do not necessarily reflect the view of the U.S. Department of Agriculture.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data used in this study is publicly available from https://apps.fs.usda.gov/fia/datamart/ (accessed on 26 August 2021).

Acknowledgments

The authors would like to thank the U.S. Forest Service, Forest Inventory and Analysis program for their assistance in the use of the data.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Jaakko Pöyry Consulting, Inc. Forest Wildlife: A Technical Paper for a Generic Environmental Impact Statement on Timber Harvesting and Forest Management in Minnesota; Jaakko Pöyry Consulting, Inc.: Tarrytown, NY, USA, 1992. [Google Scholar]
  2. Frelich, L.E.; Ek, A.R.; Zobel, J.M.; Page, K. Forest Wildlife Habitat Description and Data for Minnesota Species; Staff Paper Series No. 219; University of Minnesota, Department of Forest Resources: St. Paul, MN, USA, 2012. [Google Scholar]
  3. Burrill, E.A.; DiTommaso, A.M.; Turner, J.A.; Pugh, S.A.; Christensen, G.; Perry, C.J.; Conkling, B.L. The Forest Inventory and Analysis Database: Database Description and User Guide for Phase 2, v9.0.1.; USDA Forest Service, Rocky Mountain Research Station: Fort Collins, CO, USA, 2021.
  4. Minnesota Department of Natural Resources (MNDNR). Cooperative Stand Assessment (CSA) Users’ Manual; Minnesota Department of Natural Resources, Division of Forestry: St. Paul, MN, USA; Available online: https://iic.umn.edu/sites/iic.umn.edu/files/files/media/csa_manual_1.pdf (accessed on 4 April 2022).
  5. Arner, S.L.; Woudenberg, S.; Waters, S.; Vissage, J.; MacLean, C.; Thompson, M.; Hansen, M. National Algorithms for Determining Stocking Class, Stand Size Class, and Forest Type for Forest Inventory and Analysis Plots. Unpublished Report, USDA Forest Service, Forest Inventory and Analysis, 2001 (Modified 2003). Available online: https://www.fia.fs.fed.us/library/sampling/docs/supplement4_121704.pdf (accessed on 4 April 2022).
  6. United States Department of Agriculture (USDA). Forest Inventory and Analysis Database; USDA Forest Service, Northern Research Station: St. Paul, MN, USA, 2021. Available online: https://apps.fs.usda.gov/fia/datamart/CSV/datamart_csv.html (accessed on 26 August 2021).
  7. Zobel, J.M.; Ek, A.R.; Burk, T.E. Comparison of Forest Inventory and Analysis surveys, basal area models, and fitting methods for the aspen forest type in Minnesota, For. Ecol. Manag. 2011, 262, 188–194. [Google Scholar] [CrossRef]
  8. United States Department of Agriculture (USDA). Forest Inventory and Analysis National Core Field Guide, Vol. 1: Field Data Collection Procedures for Phase 2 Plots, v9.1.; Internal report; USDA Forest Service, Forest Inventory and Analysis: Washington, DC, USA, 2021.
  9. Zobel, J.M.; Ek, A.R. The Wildlife Habitat Indicator for Native Genera and Species (WHINGS): Methodology and Application; Staff Paper Series No. 231; University of Minnesota, Department of Forest Resources: St. Paul, MN, USA, 2014. [Google Scholar]
  10. Zobel, J.M.; Ek, A.R.; O’Hara, T.J. Quantifying the opportunity cost of extended rotation forestry with cohort yield metrics in Minnesota. For. Sci. 2015, 61, 1050–1057. [Google Scholar] [CrossRef]
  11. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021. [Google Scholar]
Figure 1. The distribution of size classes (seedling/sapling (green), poletimber (red), sawtimber (blue)) by age class for four forest type groups (upland conifer (UC), lowland conifer (LC), northern hardwoods (NH), and aspen-birch (AB)) in Minnesota. Data came from Forest Inventory and Analysis (FIA) records spanning 1999–2013. Due to considerable variation in the sample size, the y-axis units vary between graphs, but the x-axis units are identical. All curves were fit using a smoothing spline (with a smoothing parameter equal to 0.6).
Figure 1. The distribution of size classes (seedling/sapling (green), poletimber (red), sawtimber (blue)) by age class for four forest type groups (upland conifer (UC), lowland conifer (LC), northern hardwoods (NH), and aspen-birch (AB)) in Minnesota. Data came from Forest Inventory and Analysis (FIA) records spanning 1999–2013. Due to considerable variation in the sample size, the y-axis units vary between graphs, but the x-axis units are identical. All curves were fit using a smoothing spline (with a smoothing parameter equal to 0.6).
Forests 13 00647 g001
Table 1. Stand age class (years) to stand size class crosswalk. The table lists the age ranges associated with each size class for four forest type groups in Minnesota based on Forest Inventory and Analysis (FIA) records spanning 1999–2013.
Table 1. Stand age class (years) to stand size class crosswalk. The table lists the age ranges associated with each size class for four forest type groups in Minnesota based on Forest Inventory and Analysis (FIA) records spanning 1999–2013.
Forest TypeSize Class
Seedling/SaplingPoletimberSawtimber
Upland Conifer0–2021–55>55
Lowland Conifer0–8081–140>140
Northern Hardwoods0–3536–60>60
Aspen-Birch0–3031–65>65
Table 2. Age class to size class crosswalk percent accuracy and error for each forest type group. Also included are the error rates across all groups and the sample sizes. Errors are reported as proportions of the observed total number of plots in the associated category. Note that the FIA size classes were coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Thus, negative and positive numbers represent estimates smaller and larger than the actual size class, respectively. Data came from FIA plot records spanning 1999–2013.
Table 2. Age class to size class crosswalk percent accuracy and error for each forest type group. Also included are the error rates across all groups and the sample sizes. Errors are reported as proportions of the observed total number of plots in the associated category. Note that the FIA size classes were coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Thus, negative and positive numbers represent estimates smaller and larger than the actual size class, respectively. Data came from FIA plot records spanning 1999–2013.
Forest TypeError (Observed-Fitted)Sample Size
−2−1012
Upland Conifer0.0060.1340.6070.1800.072777
Lowland Conifer0.0390.2150.5520.1820.0122308
Northern Hardwoods0.0160.1270.5660.2580.0332362
Aspen-Birch0.0100.1280.6950.1450.0223816
All0.0180.1500.6190.1860.0279263
Table 3. Age class to size class crosswalk accuracy and error rates for each forest type/size class combination, along with plot counts for each actual/estimated size class combination. Percentages reported as proportions of the observed total number of plots in a size class (i.e., conditional proportions based on rows). Size classes were coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Data came from FIA plot records spanning 1999–2013.
Table 3. Age class to size class crosswalk accuracy and error rates for each forest type/size class combination, along with plot counts for each actual/estimated size class combination. Percentages reported as proportions of the observed total number of plots in a size class (i.e., conditional proportions based on rows). Size classes were coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Data came from FIA plot records spanning 1999–2013.
Forest TypeObserved Size ClassEstimated Size Class
123
Upland Conifer10.674(190)0.309(87)0.018(5)
20.223(59)0.713(189)0.064(17)
30.243(56)0.352(81)0.404(93)
Lowland Conifer10.167(41)0.472(116)0.362(89)
20.023(17)0.466(347)0.511(380)
30.021(28)0.306(403)0.673(887)
Northern Hardwoods10.771(890)0.196(226)0.033(38)
20.593(538)0.325(295)0.083(75)
30.260(78)0.237(71)0.503(151)
Aspen-Birch10.525(403)0.426(327)0.049(38)
20.211(319)0.683(1034)0.106(161)
30.055(84)0.153(235)0.792(1215)
Table 4. Ratio adjustments for each forest type/estimated size class combination. Values in the table are reported as proportions of the estimated total number of plots in a size class (i.e., conditional proportions based on columns). Each proportion represents how plots or hectares within an estimated size class should either remain in the estimated size class (diagonal elements within each forest type group) or be redistributed to the other classes based on known classification error rates (off-diagonal elements). Size classes were coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Data came from FIA plot records spanning 1999–2013.
Table 4. Ratio adjustments for each forest type/estimated size class combination. Values in the table are reported as proportions of the estimated total number of plots in a size class (i.e., conditional proportions based on columns). Each proportion represents how plots or hectares within an estimated size class should either remain in the estimated size class (diagonal elements within each forest type group) or be redistributed to the other classes based on known classification error rates (off-diagonal elements). Size classes were coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Data came from FIA plot records spanning 1999–2013.
Forest TypeObserved Size ClassEstimated Size Class
123
Upland Conifer10.6230.2440.043
20.1930.5290.148
30.1840.2270.809
Lowland Conifer10.4770.1340.066
20.1980.4010.280
30.3260.4650.654
Northern Hardwoods10.5910.3820.144
20.3570.4980.284
30.0520.1200.572
Aspen-Birch10.5000.2050.027
20.3960.6480.114
30.1040.1470.859
Table 5. Validation statistics (percent error and percent difference) and sample sizes for the age class to size class crosswalk by forest type across all FIA plot records spanning 2014–2018. Sample sizes are those for observed size classes, which are coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Percent errors are given for using the crosswalk alone and with the ratio adjustments. Percent differences quantify the relative change in percent accuracy from including the adjustments, with positive values indicating an improvement.
Table 5. Validation statistics (percent error and percent difference) and sample sizes for the age class to size class crosswalk by forest type across all FIA plot records spanning 2014–2018. Sample sizes are those for observed size classes, which are coded as 1—sawtimber, 2—poletimber, and 3—seedling/sapling. Percent errors are given for using the crosswalk alone and with the ratio adjustments. Percent differences quantify the relative change in percent accuracy from including the adjustments, with positive values indicating an improvement.
Forest TypeObserved Size Class% Error No Adj.% Error Ratio Adj.% Diff.Sample Size
Upland Conifer10.0550.123−0.068109
2−0.464−0.0720.39284
30.493−0.1090.38367
Lowland Conifer10.494−0.0510.44487
2−0.394−0.1010.293231
30.1030.0600.044464
Northern Hardwoods1−0.2740.0950.179419
20.475−0.1060.369255
3−0.075−0.161−0.08680
Aspen-Birch1−0.318−0.2010.116173
20.1620.1440.018469
3−0.044−0.069−0.025475
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zobel, J.M.; Ek, A.R.; Gifford, T.S. Stand Age Class to Size Class Crosswalk by Forest Type Group in Minnesota, USA. Forests 2022, 13, 647. https://doi.org/10.3390/f13050647

AMA Style

Zobel JM, Ek AR, Gifford TS. Stand Age Class to Size Class Crosswalk by Forest Type Group in Minnesota, USA. Forests. 2022; 13(5):647. https://doi.org/10.3390/f13050647

Chicago/Turabian Style

Zobel, John M., Alan R. Ek, and Tyler S. Gifford. 2022. "Stand Age Class to Size Class Crosswalk by Forest Type Group in Minnesota, USA" Forests 13, no. 5: 647. https://doi.org/10.3390/f13050647

APA Style

Zobel, J. M., Ek, A. R., & Gifford, T. S. (2022). Stand Age Class to Size Class Crosswalk by Forest Type Group in Minnesota, USA. Forests, 13(5), 647. https://doi.org/10.3390/f13050647

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop